Symmetry protection of measurement-based quantum computation in ground states

The two-dimensional cluster state, a universal resource for measurement-based quantum computation, is also the gapped ground state of a short-ranged Hamiltonian. Here, we examine the effect of perturbations to this Hamiltonian. We prove that, provided the perturbation is sufficiently small and respects a certain symmetry, the perturbed ground state remains a universal resource. We do this by characterising the operation of an adaptive measurement protocol throughout a suitable symmetry-protected quantum phase, relying on generic properties of the phase rather than any analytic control over the ground state.


I. INTRODUCTION
A quantum computer relies on quantum entanglement to achieve computational speedups. In the traditional, circuit-based model for quantum computation, the required entanglement is built up throughout the course of the computation through application of entangling gates coupling two or more qubits at a time. Alternatively, in the model of measurement-based quantum computation (MBQC) [1,2], universal quantum computation is achieved solely through single-particle operations (specifically, single-particle measurements) on a fixed entangled resource state, independent of the quantum algorithm being performed.
Since the initial discovery that the 2-D cluster state is a universal resource for MBQC [1], much effort has been devoted to characterizing other universal resource states. Many of the universal resource states so far identified [1,[3][4][5][6] have been projected entangled pair states (PEPS) [7] of small bond dimension. The tensor network structure of these states facilitates the analysis of measurements, which might otherwise be an intractable problem. Another advantage of such states is that under appropriate conditions [8], they are unique (possibly gapped) ground states of local frustration-free Hamiltonians on spin lattices. This suggests a method of constructing the resource state by cooling an appropriate interacting spin system [9,10].
However, if we wish to adopt this viewpoint of the resource state for MBQC as the ground state of a quantum spin system, it would be too restrictive to confine ourselves to states in which the effect of measurements can be determined analytically from the tensor-network structure. A generic local Hamiltonian, or even an arbitrarily small generic local perturbation to a PEPS parent Hamiltonian, will not have such a property. Therefore, it is desirable to develop an understanding of MBQC in ground states of spin systems that does not rely on analytic control of the ground state. For this reason, there has been an interest in relating MBQC to forms of quantum order which, as parameters of the Hamiltonian are varied, can disappear only at a quantum phase transition [11][12][13].
In this paper, we will use such a connection between MBQC and quantum order to give a precise characterization of the operation of MBQC in the ground states of a large class of perturbations to the 2-D cluster model. This will allow us to give a rigorous proof that such perturbed ground states remain universal resources for MBQC provided that the perturbation is sufficiently small. Our proof relies in part on an extension of the the relationship introduced in [13] between MBQC and symmetry-protected topological (SPT) order [14][15][16], a form of quantum order characterizing quantum systems which cannot be smoothly deformed into a product state while a certain symmetry is enforced. If the perturbation to the 2-D cluster model respects an appropriate symmetry, then the perturbed ground state will still possess non-trivial SPT order, and we will show that this gives us sufficient information about the ground state to characterize the implications of the perturbation for MBQC. Our result therefore holds independently of any analytic solution for the perturbed ground state. Our proof of universality is in the same spirit as [17]. There, it was shown that, whereas measurements on the cluster state simulate quantum circuits, measurements on a noisy cluster state simulate the same circuits, but with added noise. Here, our task is complicated by the highly correlated nature of the "errors" in the resource state that from result a change in the Hamiltonian. Nevertheless, we will show how to exploit the additional structure resulting from SPT order to establish an effective noise model for ground states of appropriate perturbed cluster models. Therefore, universal quantum computation can be achieved (for sufficiently small perturbations, corresponding to sufficiently weak noise in the effective circuit model) by choosing a measurement protocol which simulates a fault-tolerant quantum circuit. The universality is then a consequence of the threshold theorem [18] for fault-tolerant quantum computation with noisy quantum circuits.

A. Summary of results
Our ultimate goal in this paper is to prove the universality for a MBQC of a class of perturbations of the 2-D cluster state. However, in order to reach this goal, arXiv:1207.4805v3 [quant-ph] 11 Nov 2012 most of this paper will be devoted to a further elucidation of the relationship between SPT order and MBQC. For simplicity of presentation, we will first explore this relationship in one-dimensional systems. It has already been shown that in a class of quantum phases characterized by SPT order, the structure implied by SPT order leads to the perfect operation of the identity gate in MBQC [13]. Here, we consider the 1-D cluster model, which lies in the simplest of the SPT phases considered in [13], and characterize the operation of non-trivial (i.e., not the identity) gates in the presence of a perturbation which respects the symmetry protecting this SPT phase. We obtain the following: Theorem 1 (Effective noise model in one dimension). Consider a measurement protocol which in the exact 1-D cluster model would simulate a sequence of gates. In the perturbed resource state, the same measurement protocol simulates the same gate sequence, but with additional noise associated with each non-trivial gate. So long as the non-trivial gates are sufficiently separated from each other by identity gates, this effective noise has no correlations between different time steps, i.e. it is Markovian.
The proof of Theorem 1 will be divided into two stages. First, in Section II we will establish Theorem 1 for ground states which are pure finitely-correlated states (pFCS), a special case of matrix-product states (MPS). For such states, both the manifestations of SPT order [15,16], and the effect of measurements [3] can be understood straightforwardly in terms of the tensor-network structure. The ideas leading to Theorem 1 can thus be understood most directly in this context. Second, in Section III we will prove Theorem 1 for arbitrary ground states within the SPT phase.
The extension of these ideas to the 2-D cluster model will be considered in Section IV. We will construct an appropriate symmetry group, such that the following result is satisfied for symmetry-respecting perturbations.
Theorem 2 (Effective noise model in two dimensions). Consider a measurement protocol which in the exact 2-D cluster model would simulate a sequence of gates. In the perturbed resource state, the same measurement protocol simulates the same gate sequence, but with additional noise associated with each gate. So long as the non-trivial gates are sufficiently separated from each other by identity gates, this effective noise has no correlations between different time steps, or between different gates taking place at the same time step, i.e., it is local and Markovian.
Combined with the existing results on fault tolerance in the circuit model [18], Theorem 2 will imply the main result of this paper: Theorem 3. For sufficiently small symmetry-respecting perturbations, the perturbed ground state remains a universal resource for measurement-based quantum computation.

II. THE EFFECTIVE NOISE MODEL CONSTRUCTION: FINITELY-CORRELATED STATES
In this section, we will prove our effective noise model result, Theorem 1, for a restricted class of ground states of infinite one-dimensional chains. Specifically, we consider pure finitely-correlated states (pFCS) [19,20]. A pFCS can be considered as the thermodynamic limit of the translationally-invariant matrix-product states (MPS) |Ψ N generated by the MPS tensor A, on finite chains of N sites with periodic boundary conditions, e.g.
(here and throughout this paper, we use a graphical notation to represent the contraction of tensors, e.g. see [21,22]). The MPS tensor A must satisfy an additional condition known as injectivity, which is related to the exponential decay of the correlation functions; each of the finite-chain states |Ψ N (for sufficiently large N ) is then the unique gapped ground state of a local frustration-free Hamiltonian [23]. We have several motivations for considering this class of ground states. First, it is widely believed that pFCS capture the essential physics of gapped ground states of infinite one-dimensional translationally-invariant spin chains in general. (Note that, although the theorem regarding the efficient approximation of ground states of finite spin chains by MPS [24,25] could be regarded as supporting this belief, we cannot use this theorem to draw any rigorous conclusions for our purposes here, since it does not hold that the MPS tensor A can be kept fixed independently of the system size for a constant accuracy.) Second, the ideas leading to our effective noise model result find their simplest and most physically meaningful expression in this context. Finally, the proof presented here will play a dual role in our paper, as it can also be applied to arbitrary quantum states, provided that they satisfy a few extra criteria in common with pFCS. Thus, in order to establish the effective noise model result for general ground states, which we do in Section III, it will suffice to provide a separate proof of these criteria.
The outline of this section is as follows. We begin in Secs. II A and II B by reviewing the properties of the 1-D cluster model and the nature of the SPT phase in which it is contained. In Sec. II C, we review the results of [13] regarding the structure shared by pFCS ground states throughout the whole SPT phase. In Sec. II D, we prove a key result: the standard adaptive measurement protocol acting on a ground state in the phase is equivalent to a non-adaptive dual process acting on a 'topologically disentangled' version of the ground state, which we refer to as the dual state. In Sec. II E, we give a characterization of the dual process in the case that the original resource state is the exact cluster state. Finally, in Sec. II F we . . . . . .
The generators of the on-site Z2 × Z2 symmetry of the 1-D cluster model. Each dot is a qubit, and the shaded areas constitute two-qubit sites.
exploit the short range of the correlations in pFCS to construct the effective noise model for any pFCS ground states within the SPT phase, establishing Theorem 1 for the case of pFCS ground states.
A. The 1-D cluster model in the absence of perturbations Here we recall the properties of the 1-D cluster model in the absence of perturbations. The Hamiltonian is where X i denotes the Pauli X operator acting on the ith site, and similarly for Z i . With appropriate boundary conditions, the system has a unique ground state (the cluster state), and an energy gap of 2, independent of the system size. This model has a global Z 2 ×Z 2 symmetry generated by the symmetry operations i even X i and i odd X i . We consider this symmetry to be on-site, which is to say it acts on states as a unitary representation U (g) of the symmetry group G = Z 2 × Z 2 = {1, x, y, z} (with y = xz), such that U (g) acts as U (g) = [u(g)] ⊗N , where N is the number of sites (we group qubits into two-qubit sites in order to ensure this condition is satisfied; see Fig. 1). As we will see in Sec. II B, the cluster model lies in a nontrivial SPT phase with respect to this symmetry, so that the cluster state cannot be smoothly deformed into a product state without breaking the symmetry [26].
The 1-D cluster state can be represented as a pFCS [3]. For our purposes we will take the MPS tensor A C to have the form This is expressed with respect to a particular basis for a two-qubit site, where |± = (|0 ± |1 )/ √ 2. Here, and throughout this paper, we use the notation A[ψ] to refer to the linear operator obtained from the MPS tensor A by interpreting A ψ * as a linear operator (acting on states from the right), where ψ * denotes the tensor obtained by complex conjugation from the rank-1 tensor corresponding to the state |ψ .
The MPS representation for the cluster state plays a crucial role in the correlation space picture [3] for the operation of the cluster state as a quantum computational wire [27]. When a projective measurement is performed on a site, giving the outcome |ψ , this is interpreted as inducing an evolution A[ψ] on a "correlation system". In the case of the cluster state, for any qubit rotation U about the x-or z-axis, one can find a product basis {|α } for a two-qubit site, such that where B α is an outcome-dependent unitary byproduct operator. This byproduct can be accounted for by adjusting the basis for future measurements depending on the outcome of the current one.
B. Symmetry-protected topological order in finitely-correlated states Here, we will review the results of [15,16] on the manifestation of SPT order in pFCS, and demonstrate that the 1-D cluster model indeed lies in a nontrivial SPT phase with respect to the Z 2 × Z 2 symmetry.
Consider some ground state which is invariant under the on-site representation U (g) = [u(g)] ⊗N of some symmetry group G, and which can be represented as a pFCS, as in Eq. (1). The tensor A can be taken to obey a symmetry condition [16,28] where β(g) is a one-dimensional linear representation of the symmetry group G, and W (g) is a projective unitary representation of the symmetry group G. This means that for some function ω, called the factor system of the projective representation, which maps pairs of group elements to phase factors. By blocking sites, we can ensure that β(g) = 1 (however, for simplicity we will assume that β(g) = 1 without blocking). Eq. (5) can then be represented graphically as Observe that W (g) can be multiplied by a g-dependent phase factor without affecting Eq. (7); a set of factor systems related by such a transformation is referred to as a cohomology class. The arguments of [15,16] show that two such pFCS ground states correspond to the same cohomology class if and only if they are in the same symmetry-protected phase. Nontrivial cohomology classes [those not containing the trivial factor system ω(g 1 , g 2 ) = 1] correspond to phases with nontrivial SPT order.
As an example, consider the cluster model, and its Z 2 × Z 2 symmetry. The on-site representation u(g) of the symmetry is generated by and the MPS tensor is given by Eq. (3). It can be shown that the symmetry condition Eq. (7) is satisfied with the projective representation W = V P , where V P is the Pauli representation This projective representation has nontrivial cohomology class, so that the cluster model lies in a non-trivial symmetry-protected phase.
C. Symmetry-respecting perturbations to the cluster state Suppose we now consider a perturbation to the cluster Hamiltonian Eq. (2), such that the perturbed model still respects the Z 2 × Z 2 symmetry and admits a pFCS ground state. Unless the perturbation is large enough to induce a phase transition, the MPS tensor A corresponding to the perturbed pFCS ground state should still satisfy the symmetry condition Eq. (7), for some projective representation W (g) with the same factor system as the Pauli projective representation V [Eq. (10)].
The general form of the MPS satisfying these symmetry conditions was established in [13]. Here we will briefly review the relevant results from [13]. We observe that the Pauli representation satisfies a property which we refer to as maximal non-commutativity: Definition 1. A projective representation W (g) of an abelian group G is called maximally non-commutative if the subgroup Z W (G) ≡ {g ∈ G : W (g) commutes with W (h) for all h ∈ G} (which we can think of as the "projective centre" of G) is the trivial subgroup.
Notice that the subgroup Z W (G) is actually determined by the factor system ω, since W (g)W (h) = ω(g, h)W (gh) = ω(g, h)ω(h, g) −1 W (h)W (g). Furthermore, it is the same for all factor systems within a given cohomology class. Much of the discussion in this paper can be applied to any SPT phase characterized by a finite abelian symmetry group and a maximally noncommutative cohomology class.
An important consequence of maximal noncommutativity of a factor system is [13] Lemma 1. For each maximally non-commutative factor system ω of a finite abelian group G, there exists a unique (up to unitary equivalence) irreducible projective representation with factor system ω. The dimension of this irreducible representation is |G|.
Specifically, the Pauli representation V P of Z 2 × Z 2 is the unique irreducible projective representation corresponding to its factor system. In general, throughout this paper, we will use V (g) to denote the unique irreducible projective representation for the current factor system. A consequence of Lemma 1 is that, for a tensor satisfying the symmetry condition Eq. (7), the bond space decomposes as a tensor product of a |G|-dimensional protected subsystem in which W (g) acts irreducibly as V (g) and a junk subsystem in which W (g) acts trivially, i.e.
Thus the tensor A appearing the MPS representation of ground states in the symmetry-protected phase satisfies the symmetry condition Here we use a thick line ( ) to represent the protected subsystem, and a dashed line ( ) to represent the junk subsystem. The protected subsystem enjoys several nice properties for storing and manipulating logical information in a quantum computation, as we now show.
Suppose we perform a projective measurement on one site in a simultaneous eigenbasis {|i } (which is |++ , |+− , |−+ , |−− for the Z 2 ×Z 2 cluster state symmetry), and obtain the outcome |i . Then the resulting state on the remaining sites is found by replacing the original MPS tensor at the measured site by Now we make use of another consequence of maximal non-commutativity [13]: Lemma 2. Let u(g) be a linear on-site representation of a finite abelian symmetry group G, and let ω be a maximally non-commutative factor system of G. Then for each basis element |i in a simultaneous eigenbasis {|i } of u(g), there exists a group element g i such that for any projective representation V (g) with factor system ω, where χ i (g) is the scalar representation of G such that u(g)|i = χ i (g)|i .
For the example of the cluster state symmetry, we have g ++ = 1, g +− = x, g −+ = z, g −− = y, as can readily be verified directly.
As was shown in [13], Lemma 2 in conjunction with the symmetry condition Eq. (12) implies the decomposition for some tensor A. Another way of writing this result is that where we have defined the tensor Note that, from a quantum circuit perspective, this tensor can also be interpreted as a unitary controlled operation i |i i|⊗V (g i ) coupling a site to an ancilla particle; hence the choice of notation. Conversely, any MPS tensor of the form Eq. (16) for some tensor A will satisfy the symmetry condition Eq. (12). Following [29], we refer to the tensor A as as the degeneracy tensor ; and to the tensor of Eq. (17), which is determined entirely by the symmetry, as the structural tensor. From Eq. (15), we see that, in the correlation space picture, measuring in a simultaneous eigenbasis {|i } leads to an evolution on the protected subsystem of correlation space given by an outcome-dependent unitary V (g i ); this evolution is determined by the symmetry (hence the same throughout the SPT phase), and decoupled from the junk subsystem. Viewing the unitaries V (g i ) as outcomedependent byproducts of the measurements, which can straightforwardly be accommodated in a deterministic evolution using the standard techniques of measurementbased quantum computation, we say that the identity gate operates perfectly throughout the SPT phase. However, the result of measurement in any other basis is not fixed by the symmetry, and in general leads to the protected subsystem being coupled to the junk subsystem, so that the operation of other measurement-based gates will not be a robust property of the symmetry-protected phase.
D. The dual picture for MBQC on a 1D resource state In order to deal with the randomness of measurement outcomes, the measurement protocol for MBQC with the cluster state needs to be adaptive: the outcome of the measurement on one site will affect the measurement basis on other sites arbitrarily far away. In analysing the effect of this protocol when acting on a perturbed resource state, we would like to make an argument based on the locality of the perturbed Hamiltonian, but the non-local adaptivity of the measurement protocol poses a difficulty. Therefore, in this section, we develop an alternate characterization of the effect of the cluster state adaptive measurement protocol acting on a ground state in the symmetry-protected phase. We will show that this protocol is equivalent to a dual process acting on a related state, which we call the dual state. We will show that this dual process simply consists of a sequence of unitary interactions between selected sites (those corresponding to the locations of non-trivial gates) and an ancilla particle, with no adaptivity.
In our discussion of the dual process, we will represent a pFCS ground state on an infinite chain as a formal tensor network This is not, of course, the mathematically rigorous way to describe pFCS, but we find it useful for facilitating understanding. In Appendix A, we will discuss how to formulate similar arguments in the rigorous pFCS framework. Later on (in Section III), we will also be interested in finite chains; in that case, the arguments of this subsection can be applied more directly, given appropriate boundary conditions [specifically, the boundary conditions at the right edge should be as depicted in Eq. (50)].
The usefulness of the exact 1-D cluster state (with MPS tensor A C ) as a quantum computational wire results from the fact that, for each gate U in a certain set, there exists a basis {|α } such that where B α is the outcome-dependent unitary byproduct operator. When we measure one site projectively and obtain the outcome |α , the original MPS tensor A is replaced at the measured site by Eq. (19) in the tensornetwork description of the resultant state.
In the case of the exact cluster state, the effect of the byproduct operator can be accounted for by adjusting the measurement basis for future measurements. This fact turns out to be closely related to the nontrivial SPT order, as we now demonstrate. Our discussion relies on the observation that, in the cluster state, the byproduct operators are Pauli operators. That is to say, it is always the case that B α is a scalar multiple of V (g α ) for some g α ∈ Z 2 × Z 2 . Hence, we can make use of the symmetry condition [which can be derived from Eq. (12)] where b α = u(g α ). Applying this condition repeatedly shows that the byproduct operator can be displaced arbitrarily far to the right. In our formal tensor-network picture for an infinite chain, we consider that this process is continued indefinitely, so that the byproduct operator "disappears out to infinity", and is replaced with b α acting on all sites to the right of the one on which the measurement took place, i.e.
Hence, whenever we obtain the "wrong" outcome for a measurement (i.e. the corresponding byproduct operator B α is not the identity), we can recover the "correct" resultant state by applying the correction b † α to all the remaining sites on the right (equivalently, we can simply adjust the measurement basis for measurements on those sites).
Let us now examine what happens when we perform the same adaptive measurement protocol on a resource state that is not the exact cluster state. Consider a pFCS ground state contained with the same SPT phase as the cluster state, characterised by the Pauli representation of the group Z2 × Z2. We will keep using the same measurement protocol as for the exact cluster state. (Our argument could be generalised to any pFCS ground state contained within any SPT phase characterized by a finite abelian symmetry group G and a maximally noncommutative cohomology class, so long as the the adaptive correction appearing in the measurement protocol takes the same form as for the cluster state, i.e. application of u(g α ) to the sites on the right for some group elements g α ∈ G.) The resource state is then of the form with the MPS tensor A of the form Eq. (16). We now repeat the above argument, in reverse. We make use of the symmetry condition Eq. (12) in the form from which we obtain Therefore, we have shown that the process we actually perform, i.e. applying the measurement-dependent correction to the sites on the right of the one measured, is equivalent to a different process, in which the measurement-dependent correction is applied in the internal bond space of the MPS, as depicted in the righthand side of Eq. (24). In a physical system, of course, we do not have direct access to the internal bonds of a tensor network state, so we could never perform the latter process directly; nevertheless, the two are equivalent. Following the measurement and the adaptive correction, which we think of as being performed internally, as in the right-hand side of Eq. (24), the outcome of the measurement can be "forgotten", i.e. we describe the resultant state of the system as the mixture of the righthand side of Eq. (24) for all possible measurement outcomes. Without affecting the reduced state on the remaining unmeasured sites, for notational convenience we replace this mixture with a coherent superposition, i.e.
where we have defined the tensor which we can also interpret as a unitary coupling coupling a site to an ancilla particle. We now use the k index throughout the paper to distinguish the "G"'s resulting from different measurement operations. (The label k refers to the site at which the measurement is being performed; we include this label to reflect the dependence on the measurement basis {|α } and byproduct operators B α , which will in general be different for each site at which a measurement is performed.) Now, using the expression Eq. (16) for the MPS tensor A, we can write where we have defined the tensor G k (which can also be interpreted as a unitary coupling between a site and the ancilla particle) according to After the adaptive measurement sequence, we can treat the resultant state as having the form shown. This is equivalent (b) to building |Ψ from the dual state | Ψ [the shaded box; see Eq. (29)] by unitary couplings to an ancilla particle.
We are now in a position to define our dual process. Suppose we perform a sequence of such adaptive measurements at successive sites, which at each site is described by the insertion of the tensor Eq. (26), as in Eq. (25). There will be a different coupling G k associated with each site k. As shown in Figure 2, we find that the original adaptive measurement process, applied to the resource state |Ψ , is equivalent to a dual process applied to the dual state | Ψ . The dual state | Ψ is defined to be the state built from the degeneracy tensor A, with the structural tensor discarded: The dual process comprises a series of consecutive unitary interactions G k between individual sites k and an ancilla particle.
There are several reasons why this "dual picture" is a useful way to understand the operation of MBQC in one-dimensional ground states. First, the dual process lacks the long-range measurement adaptivity which is a characteristic of the original adaptive measurement protocol. Second, the perfect operation of the identity gate is automatically built in, because, for sites k at which the adaptive measurement process at the given site is the one that corresponds in the exact cluster state to the identity gate, [i.e. the measurement basis is the simultaneous eigenbasis {|i } of the symmetry, and the byproduct operators are B i = V (g i ), where the g i are the group elements appearing in Eq. (15)], the corresponding coupling is trivial, G k = I.
The final motivation for the dual picture is that the dual state on which it is based has some physical significance in its own right, and retains some key properties of the the original resource state. For example, if the original MPS tensor A generates a pFCS, then so does A, and the respective correlation lengths obey the inequality ξ ≤ ξ (see Appendix A). Additionally, in Appendix B we will show how our dual state can be obtained from the original ground state through a generalization of the unitary that was introduced by Kennedy and Tasaki [30] to transform the SPT Haldane phase [31,32] into a local symmetry-breaking phase; this unitary has recently been described as a "topological disentangler" [33], and in some sense we can think of the dual state as being a topologically disentangled version of the original resource state.

The dual process for initialization and readout in the 1-D cluster state
Above we only discussed measurement sequences corresponding to unitary gates in correlation space. A complete scheme for using a 1-D resource state as a quantum computational wire also includes measurement sequences corresponding to initialization (i.e. discarding the current state in correlation space and replacing it with a fixed state), and readout (i.e. making the state in correlation space available as the physical state of one qubit). We now describe briefly how the measurement protocols used on the 1-D cluster state for these purposes can be accommodated in our framework.
Initialization.-The initialization procedure involves measuring a site in the computational basis {|00 , |01 , |10 , |11 }. In this basis, the MPS tensor A C for the exact cluster state takes the form The randomness of measurement outcomes can therefore be accounted for by applying the appropriate outcomedependent correction operator in correlation space following the measurement: B † 00 = B † 10 = I, B † 01 = B † 11 = X. Since the correction operators are Pauli operators, the above discussion applies without change.
Readout.-The standard readout procedure for the cluster state involves measuring the second qubit of a two-qubit site in the computational basis, then applying an outcome-dependent correction operator to the first qubit, which acts as the output. Provided that we are only interested in the final state of the output qubit, this procedure is equivalent to a coherent correction operator coupling the two qubits in the site (specifically, it is a controlled-Z gate C Z = |0 0| ⊗ I + |1 1| ⊗ Z). Carrying through a similar argument to that given above for unitary gates, Sec. II D, we obtain the same result, but with the interaction G k in the dual process between the site k in question and the ancilla particle given by E. MPS of minimal bond dimension and the dual picture As an example of the general formalism introduced in Sec. II D, here we will examine the form of the couplings G k appearing in the dual process [Eq. (28)], in the particular case that the resource state is an MPS with bond dimension D = |G|, where G is the symmetry group characterizing the symmetry-protected phase. Given that the dimension of the protected subsystem is |G| (by Lemma 1), this is the smallest possible value of D, and corresponds to the absence of a junk subsystem (or, more precisely, a junk subsystem of dimension 1). In particular, the 1D cluster state is of this type. In general, the MPS tensor A for such an MPS must be of the form where the A[i] here are scalars. It follows that the dual of such a state is a product state, with each site in the state |φ = i A[i]|i . (We choose the normalization for the MPS tensor A so that φ|φ = 1.) Therefore, the effect of the dual process acting on the dual state results from a series of independent interactions of the form We recall that, in the correlation space picture of quantum computational wires, a quantum state can serve as a resource for executing a unitary gate U if there exists some basis {|α } such that for some set of unitary byproduct operators B α and scalars β α . We will now show how this property manifests itself in the dual picture, for the class of states considered here. We make use of the representation for the MPS tensor A as It follows that, at a site k measured in the basis {|α }, with the byproduct operators B α , we have where |ϕ = α β α |α . (It can be shown that our choice of normalization ensures that ϕ|ϕ = 1.) Here the first equality follows from Eq. (38) and the definition of G; and the second inequality follows by Eq. (37). Thus, we have shown that in the dual picture the gate U simply acts on the ancilla particle. Next we will do a similar analysis for the initialization and readout procedures specific to the 1-D cluster state.
Initialization.-Recall the discussion of initialization in Sec. II D. We make use of the form of the MPS tensor A C in the computational basis, Eqs. (30)(31)(32)(33), multiplied by the appropriate normalization factor as discussed above. Thus, for the site k at which initialization takes place, we find that where Γ = 1 √ 2 (|00 0| + |01 0| + |10 1| − |11 1|) (thanks to our choice of normalization, we find that Γ is an isometry, i.e. Γ † Γ = I). Therefore, applying the measurement sequence for initialization leads to the ancilla system getting initialized in the state |0 , as we would expect.
Readout.-From the definition of the operator G k in the case of sites k at which readout takes place [Eq. (34) As long as the locations of non-trivial gates (sites k and l in this diagram) are separated by a distance much greater than the correlation length ξ, the reduced state on those locations will be approximately a product state, and the dual process then reduces (b) to a sequence of independent interactions.
we find that (here we separate a site into its two constituent qubits, each denoted by a thick line) where |I = (1/ √ 2) (|00 + |11 ) is the canonical maximally-entangled state. Thus the state of the ancilla qubit indeed gets transferred onto the output qubit.

F. MBQC on a perturbed resource state simulates a noisy quantum circuit
In the previous subsection, we saw how measurements on an MPS of minimal bond dimension correspond to quantum gates. Now we will consider what happens when we perform the same measurement sequences on a perturbed resource state, assuming that the perturbed state remains within the same SPT phase. We will find that measurements on such a perturbed cluster state simulate the same quantum circuit, but with noisy gates. The noise is described by application of a completely positive, trace preserving (CPTP) noise superoperator following each gate.
In Section II E, we were able to treat each gate independently in the case of the unperturbed cluster state because the dual state was a product state, | Ψ = |φ ⊗N . This will no longer be true once we introduce perturbations, but we still want to treat gates independently. Towards this end, we recall that for a site for which the corresponding sequence is that for the identity gate, the associated coupling G k in the dual picture between that site and the ancilla particle is trivial. Therefore, such a site can be traced out from the beginning without affecting the final state of the output qubit. That is to say, we only need to consider the reduced state ρ = Tr trivial sites | Ψ Ψ| on the remaining sites, which are those corresponding to non-identity gates (we refer to these as the non-trivial sites). We are free to choose our measurement protocol to ensure that the distance between any two non-trivial sites is much greater than the correlation length ξ. For pFCS, it is then straightforward to show that ρ is approximately a product state ρ prod = k ρ k over the non-trivial sites, or more precisely where . 1 denotes the trace norm, m is the number of non-trivial sites, R is the minimum distance between any two non-trivial sites, and f (R) is a function related to the transfer channel of the pFCS, behaving asymptotically as where ξ is the correlation length associated with the pFCS. We first consider the case where ρ = ρ prod exactly. Because ρ is then a product state, we find, as in the previous subsection, that the dual process acting on the dual state is again effectively a sequence of independent interactions, this time of the form Fig. 3). Let us suppose that G k results from the measurement sequence corresponding to a unitary gate U k . Then, after tracing out the physical site, Eq. (46) corresponds to an evolution on the ancilla qubit described by the CPTP map As we saw in Sec. II E, in the absence of perturbations to the cluster state, ρ k = |φ φ| and where . ♦ is the diamond norm on superoperators [34].
The cases when G k corresponds to initialization or readout are analogous. Therefore we have shown (in the case ρ = ρ prod ) that the measurement protocol on the perturbed cluster state reproduces the desired quantum circuit, except that each gate (as well as the initialization and readout steps) is accompanied by some associated noise. Furthermore, if the perturbation is sufficiently small, then the reduced states ρ k will be close to |φ φ| (see Appendix D for the proof), so that the noise will be weak, in the sense that E k is close to the identity superoperator in the diamond norm.
In the general case, in which ρ and ρ prod are not equal, but are -close in the trace distance, we just need to observe that the reduced state of the output qubit following the dual process can be obtained from ρ by application of some CPTP superoperator, which we call B. From the contractivity property of the trace distance, it follows that B( ρ) − B( ρ prod ) 1 ≤ ρ − ρ prod 1 ≤ . Therefore, the effective noisy quantum circuit description correctly describes the final state of the output qubit up to an accuracy . Note that, because the bound Eq. (45) depends on the number of non-trivial gates m, it will be necessary to have the separation R scale with m in order to obtain a fixed accuracy , but only logarithmically; indeed, the minimum separation required to achieve an accuracy scales like R min / ξ = O[log(m/ )].

G. Summary of Section II
In Section II, we have presented, within the context of pure finitely-correlated states, the main ideas leading to our effective noise model construction. Our discussion has hinged around the "dual state" which we associated with each ground state carrying the appropriate SPT order. Loosely speaking, we can think of the entanglement in SPT-ordered ground states as comprising "topological" and "non-topological" components intertwined. The topological component is fixed throughout the phase and is responsible for the distinctive characteristics of the SPT phase, such as the degeneracy in the entanglement spectrum [35], the diverging localizable entanglement length [36,37], and the perfect operation of the identity gate. One can think of the dual state as being obtained from the original ground state by a topological disentangler, "separating out" the topological component of the entanglement and leaving only the non-topological component [33].
In this paper, the importance of the dual state is due to the following fact, which we established in Sec. II D: the cluster state adaptive measurement protocol, when applied to an SPT-ordered ground state, couples in a natural way to the topological component of the entanglement, and the effect is thus equivalent to a "dual process" (with a simpler structure) acting on the dual state. This result gives rise to an effective quantum circuit description describing the outcome of the measurement protocol applied to any SPT-ordered ground state [e.g. see Fig. 3(a)]. The action of non-trivial gates is determined by an interaction with a single site in the dual state, and perturbations to the dual state give rise to noisy gates. As long as the locations of non-trivial gates are sufficiently separated, the reduced state on the sites relevant for the gate operation will be a product state, and this corresponds to independent noise acting on each gate in the effective circuit description.

III. THE EFFECTIVE NOISE MODEL CONSTRUCTION: GENERAL GROUND STATES
In this section, we will extend our characterization of the effective noise model to any ground state within the symmetry-protected phase, without reference to finitely correlated states. Instead of starting from scratch, we will build on the results of Section II, as follows. We formulate a condition which we believe (on physical grounds) to be satisfied for any system within the symmetry-protected phase. We will show that this condition leads to a construction for the dual state of any ground state in the SPT phase, independently of the pFCS formalism. Furthermore, given an exact MPS representation for the dual state (which always exists, albeit possibly with a bond dimension exponentially large in the system size), we show that one can construct a corresponding MPS representation for the original ground state, such that the arguments of Section II can be applied without significant change. In order to establish the approximate factorization condition Eq. (45) in the case of general ground states, we will show that the dual state is (like the original resource state) the gapped ground state of a local Hamiltonian, which can be constructed in a straightforward way from the original Hamiltonian. This will allow us to establish the approximate factorization condition Eq. (45) without assuming that the dual state has a pFCS structure.
A. Symmetry-protected topological order and boundary conditions Because we are considering general ground states, we can no longer make direct use of the characterization of SPT order in finitely-correlated states of [15,16]. Instead, we adopt the perspective in which SPT order is related to the fractionalized edge modes associated with open boundary conditions [32]. Our discussion will, out of necessity, be physically motivated rather than mathematically rigorous, but will suggest the formulation of the precise assumptions under which the rigorous results of this paper can be proven.
Consider a 1-D chain with open boundary conditions, with symmetry-respecting interactions such that, in the bulk, there is no symmetry-breaking and a finite energy gap for excitations. It is still possible that the energy gap for edge excitations may be much smaller than the bulk gap (or even zero). We denote by P the subspace comprising the low-lying edge states. If the chain is sufficiently long, we expect that the gap in the bulk should ensure that the edges are non-interacting, so that P decomposes a tensor product of degrees of free-dom associated with the left and right edges respectively, P = P l ⊗ P r , and the restriction H P of the Hamiltonian H onto this subspace is a sum of terms acting on each edge: H P = h l ⊗ I + I ⊗ h r .
For our purposes, we are interested in how the edge states transform under the symmetry. Because P is spanned by a set of energy eigenspaces, and the Hamiltonian commutes with the representation of the symmetry U (g), it follows that P must be an invariant subspace for U (g). We write the operation of U (g) restricted to P as U P (g). We expect that, for sufficiently long chains, the symmetry should act independently on the respective edge states, i.e.
By assumption, U (g), and hence U P (g), is a linear representation of the symmetry group G. It follows that V l (g) and V r (g) are in general projective representations of G, and if V l (g) has factor system ω, then V r (g) must have factor system ω −1 . We are free to trans-  [38,39].
In non-trivial SPT phases, the edge interactions h l and h r (and therefore the overall Hamiltonian H) will always have degenerate ground states, due to the fact that non-trivial projective representations cannot be onedimensional. On the other hand, we expect, at least in the case that the symmetry group G is abelian, that a non-degenerate ground state can be recovered by introducing terminating particles at the left and right edges, transforming projectively under the symmetry with factor systems ω −1 and ω respectively (see Figure 4). This is because, loosely speaking, these terminating particles can couple to the edge modes, with the composite system at each edge transforming under a linear representation (and therefore, in the case of an abelian symmetry group, generically having a non-degenerate ground state). For example, the ground state of a spin chain in the Haldane phase can be made non-degenerate through coupling to spin-1/2 particles at the edges. Conversely, if the terminating particles do not transform with the cohomology classes [ω −1 ] and [ω] respectively, then the degeneracy cannot be removed completely because there is still a non-trivial projective symmetry transformation at each edge. Thus, we have arrived at yet another characterization of SPT order, which we state as a conjecture in the absence of a rigorous proof: In any case, in the remainder of this section, we will consider systems satisfying Condition 1. Specifically, all the results will apply to finite chains with the appropriate edge couplings imposed to ensure a non-degenerate gapped ground state. This will prove convenient for our analysis, but the properties of the system in the bulk should not, of course, depend on the boundary conditions.
Note also that, in the case of a system with the interactions governed by the parent Hamiltonian of a pFCS [generated by an MPS tensor satisfying the symmetry condition Eq. (7) corresponding to the symmetry-protected phase], Condition 1 can easily be established directly. Furthermore, the stability theorem of [40] ensures that Condition 1 remains true for sufficiently small symmetryrespecting perturbations of such models, regardless of the validity of Conjecture 1.
The "topological disentangler" D applied to the original ground state (a) turns it into the dual state (b), leaving the terminating particles maximally entangled.

B. The general construction for the dual state; exact MPS representation of SPT-ordered ground states
Recall that in Section II D, we defined the dual state in the context of pFCS. Here, we will give an analogous construction for the dual state corresponding to a general ground state within a symmetry-protected phase, provided that the phase is characterized by a finite abelian group G and a maximally non-commutative cohomology class [ω]. The construction applies to a finite chain, with the appropriate boundary conditions as discussed in Section III A. This construction will then allow us to express the original ground state as an MPS, with the MPS tensors satisfying an appropriate symmetry condition.
We consider a finite chain coupled to terminating particles, such that the overall system is invariant under the symmetry U (g) = V * (g) ⊗ [u(g)] ⊗N ⊗ V (g). Here we have taken the right terminating particle to transform under V (g), the unique irreducible projective representation with factor system ω; and the left terminating particle under V * (g) [V * (g) is the operator obtained from V (g) by complex conjugation of the matrix elements in some basis; observe that V * (g) is a projective representation of G with factor system ω −1 ].
The natural analogues in the current setting (finite chains, with the specific choice of boundary conditions) of the pFCS ground states which we considered in Section II are states of the form for some end vectors |L and R|, and where the MPS tensor A satisfies the symmetry condition Eq. (12) [which can be shown to ensure the invariance of the state under U (g)]. Given the decomposition Eq. (16) for the MPS tensor A, it follows that the dual state can be obtained from the original ground state by a sequence of unitary interactions between individual sites and the terminating particle on the right (see Fig. 5); we can think of the overall unitary transformation D as a "topological disentangler". Specifically, we have D|Ψ = | Ψ ⊗ |I , where |I = D k=1 |k |k is the canonical maximally-entangled state between the terminating particles.
We will now show that, for a general gapped symmetryrespecting ground state |Ψ [not necessarily in the MPS form Eq. (50)], it remains the case that D|Ψ = | Ψ ⊗ |I for some state | Ψ on the non-terminating sites; this will serve as the definition of the dual state | Ψ for general ground states.
We observe that the original ground state |Ψ must be invariant under the global symmetry operation U (g), i.e.
(Without loss of generality, we can set α(g) = 1 by absorbing it into into the definition of the symmetry 1 .) This implies that the state D|Ψ is invariant under DU (g)D † .
Let us examine what this 'dual' symmetry looks like. We observe that where is the interaction from which D is built; to get to Eq. (53), we made use of Eq. (14). From this, one can show that It it straighforward to show (using the irreducibility of V ) that invariance of a state under the right-hand side of Eq. (56) implies that it must be of the form D|Ψ = | Ψ ⊗ |I for some state | Ψ , as required.
It is now straightforward to construct an appropriate MPS representation for a general ground state. Indeed, let us consider an MPS representation for the dual state | Ψ of the form We choose this representation to be exact; this may require the bond dimension to be very large (scaling exponentially in the system size), but that is of no importance to us. Then we have Now we can apply the inverse transformation D † to obtain This is a representation of |Ψ as an MPS, with each of the shaded regions corresponding to an MPS tensor A of the form Eq. (16), and hence satisfying the symmetry condition corresponding to the symmetry-protected phase. In addition, we should take note of the boundary conditions at the right edge. These boundary conditions ensure that the arguments of Section II D apply without any need to invoke an infinite-system limit.

C. The dual state as the ground state of a local Hamiltonian
In the previous subsection, we have constructed the dual state for any ground state in the symmetryprotected phase. The original ground state is, by assumption, the gapped ground state of a local Hamiltonian. In this subsection we will show that this is also true of the dual state. That is, starting from the original Hamiltonian H, we construct another local Hamiltonian H for which the dual state is the gapped ground state.
We start by proving a useful fact about the unitary transformation D introduced in the previous section: although it is in general non-local, it maps symmetryrespecting local observables (i.e. those supported on a small set of sites of finite size) to local observables. Indeed, let us consider some local observable h; we will show that DhD † is also local. For concreteness, we suppose that h acts on two adjacent sites somewhere in the bulk. Now, observe that By means of Eqs. (62) and (54), it can be verified that if h commutes with the symmetry, i.e.
Since V (g) is an irreducible projective representation, Eq. (64) implies (by Schur's Lemma) that h acts trivially on the terminating particle, i.e.
Now, using Eqs. (61) and (65), we find that DhD † = h, where h acts on the same two sites as h (see Fig. 6). Thus, although the duality transformation D is non-local, we have shown that it maps local symmetry-respecting operators to local operators, as promised. The exception is operators h at the left edge, which act non-trivially on the left terminating particle; in that case, the above argument breaks down, but we can observe directly from the structure of D that DhD † is supported on the union of the support of h and the right terminating particle. For operators h acting non-trivially on the right terminating particle, the argument must be adjusted, but the conclusion that DhD † is supported on the support of h still holds.
We are now in a position to construct the Hamiltonian for which the dual state | Ψ is the gapped ground state. We observe that DHD † has |I ⊗ | Ψ as its gapped ground state; however, it includes terms acting nontrivially on the terminating particles. We define a Hamiltonian acting only on the intermediate sites according to H = I|DHD † |I ≡ F(H); by the locality result proven above, each local interaction term in H corresponds to a local term in H supported on the same set of sites. It can be shown that | Ψ is the unique ground state of H, and that the gap is at least as large as that of DHD † , or equivalently H.

D. The factorization condition for general ground states
Recall that the other condition that needed to be satisfied in order to apply the arguments of Section II for general ground states was that the factorization condition for the reduced density operator ρ on the non-trivial sites in the dual state, should be satisfied when the non-trivial sites are sufficiently separated from each other. Recall that, for the case of pFCS, one can prove the bound with m the number of non-trivial sites, and f (R) a function scaling asymptotically as where ξ is the correlation length in the dual state. We conjecture that Eq. (67) should be a general property of all gapped ground states of local Hamiltonians. However, we have only been able to rigorously prove the weaker bound where f (R) is as before, and d is the dimension of the Hilbert space at each site; see Appendix C for the proof. Note that if we assume only the weaker bound Eq. (68), then the separation between non-trivial sites will need to scale more rapidly with the number of gates m; we find that the minimum separation R min required for an accuracy scales like This still implies that the number of measurements need scale only as a polynomial in the number of non-trivial gates.

E. Nonzero temperature
The formulation of the dual state as the ground state of a dual Hamiltonian extends naturally to nonzero temperature: under the topological disentangler D, the thermal state of the original Hamiltonian H maps to the thermal state of a dual Hamiltonian H. Furthermore, it can be shown that an appropriate adaptive measurement protocol acting on the thermal state of the original Hamiltonian is equivalent to a non-adaptive dual process (of the same form as in the zero-temperature case), acting on the thermal state of the dual Hamiltonian.
However, it does not appear possible to construct a Markovian effective noise model for nonzero temperature using the same techniques as for zero temperature. The reason is that our arguments were based on the assumption that the reduced state ρ k on each of the non-trivial sites in the dual state does not differ greatly from its value in the dual of the unperturbed resource state. This is indeed the case for small local perturbations to the Hamiltonian (as we prove in Appendix D), but it need not be true for nonzero temperature. For example, consider the one-dimensional Ising model, with Hamiltonian (we include the Z 1 term to select out a unique ground state). In this model, it can be shown (e.g. using the transfer matrix method) that the reduced state on a single spin changes discontinuously as soon as the temperature is switched on (this is closely related to the disappearance of the magnetic order in the 1-D Ising model at nonzero temperature). Given the structure of the dual Hamiltonian as discussed in Appendix B, there is good reason to believe that it will exhibit a similar phenomenon.
The difficulty of treating thermal states in our framework should not be surprising, as the dual process has the perfect operation of the identity gate built in, whereas the cluster model is not expected to have a long-range identity gate at nonzero temperature. On the other hand, there exists a measurement protocol for a threedimensional cluster model which retains the perfect operation of the identity gate at sufficiently small nonzero temperatures [41]. Therefore, if the dual process description could be extended to measurement protocols such as this one, then it might be expected that the dual Hamiltonian would possess an ordered phase that persists at nonzero temperature, such that the local reduced state varies continuously with temperature up to the phase transition.

IV. TWO-DIMENSIONAL SYSTEMS AND FAULT TOLERANCE
The equivalence we demonstrated in Sections II and III, between MBQC on perturbed resource states and noisy quantum circuits, opens up the possibility of exploiting the results in the literature on fault-tolerant quantum computation with noisy quantum circuits. Here, we will extend the results of the previous sections to the 2D cluster model, which, unlike the 1D models considered previously, is a universal resource for quantum computation. We will again find that, provided the perturbation to the Hamiltonian respects a certain symmetry, MBQC using the perturbed ground state as a resource is equivalent to a noisy quantum circuit. We will show that the noise in this effective circuit description has no correlations in time (as in the previous section), nor any correlations in space. This reduction to local, Markovian noise will allow us to invoke the threshold theorem to deduce that, provided the perturbation respects the symmetry and is sufficiently small, the perturbed ground state remains a universal resource for MBQC. It should be emphasized that, although we make use of the theory of fault-tolerant quantum computation, our final result cannot be described as a fault tolerance result for MBQC, since it applies only to symmetry-respecting perturbations, and we assume noiseless operation of the measurement protocol.

A. The 'quasi-1D' model
Here, we make a first attempt at generalizing the 1D results to a 2D model which is universal for quantum computation. The ground state of the model we introduce here is not strictly a universal resource for MBQC unless we allow non-single-qubit measurements; however the discussion here will serve as a stepping stone for consideration of the 2D cluster model in Sec. IV B.
In the absence of perturbations, the 2D model we consider involves N uncoupled 1D cluster states arranged in the second dimension, as shown in Fig. 7. The Hamiltonian acting on each chain is simply the 1D cluster Hamiltonian. For generality we also assume the existence of some uncoupled qubits, each with an associated term −X in the Hamiltonian (i.e. the ground state is |+ ). In order to treat this 2D model within the same framework which we have developed for 1D systems, we will consider an entire N v × 2 block (where N v is the extent in the vertical direction) to be a single 'site', as shown in Fig. 7(a); hence we can consider the lattice to comprise a 1D chain of such 'sites'. The unperturbed ground state, which we denote |Ψ Q , then has an MPS representation which is essentially a tensor product of several copies of the 1D cluster state MPS representation, with a correlation system comprising N qubits.  FIG. 8. In the 'quasi-1D' model, two-qubit gates in correlation space cannot be done with single-qubit measurements. However, after applying controlled-Z gates between neighbouring qubits in order to couple two chains, an entangling gate can be performed in correlation space by means of singlequbit measurements.
Each chain contributes a separate Z 2 × Z 2 symmetry, so that the model is invariant under a symmetry group The projective representation of this symmetry in correlation space is the N -qubit generalization of the Pauli representation, namely where V P is the single-qubit Pauli representation of Z 2 × Z 2 , given by Eq. (10). It can be checked that this projective representation is maximally noncommutative, and because it has dimension 2 N = |G|, it must be the unique irreducible projective representation corresponding to its factor system (by Lemma 1 from Sec. II C). Now, it is easy to see that, where S is the set of gates which can be executed in correlation space by a single-site measurement (up to Pauli byproducts) in the 1D cluster state, we can execute any tensor product s 1 ⊗ s 2 ⊗ · · · ⊗ s N , s 1 , · · · , s N ∈ S (72) (up to Pauli byproducts) in correlation space by a singlesite measurement in our 2D model; we just do the corresponding measurements on each chain separately. We can also find a measurement basis for a columnar site which induces entangling gates between two qubits in correlation space; however, this measurement basis clearly cannot correspond to single-qubit measurements, since the two chains would then remain uncoupled. For reasons that will become clear when we consider the relation of the present model to the 2-D cluster state in Sec. IV B, we will only consider entangling gates between nearestneighbour qubits in correlation space, for which we construct the measurement basis in a particular way, as follows.
We define the on-site unitary u, which involves applying controlled-Z gates between neighbouring qubits to turn our original resource state |Ψ Q into another graph state |Ψ Q in which the two chains of interest are coupled, as shown in Figure 8. If A[·] is the MPS tensor for |Ψ Q at the given site, then A [·] = A[u † (·)] is the MPS tensor for |Ψ Q . Using the measurement sequences described in [42], it can be shown that there exists a measurement basis {|α } for a columnar site, corresponding to single-qubit measurements, such that A [α] = B α U , where U is an entangling two-qubit gate, and the B α are outcome-dependent Pauli byproducts. It follows that this two-qubit gate can be performed in correlation space (up to the same Pauli byproducts) by measuring in the basis {u † |α }.
From the above considerations, we see that the model we are discussing can be considered as a generalization of the 1-D cluster state in which N qubits can be propagated in correlation space, acted on by entangling gates between nearest neighbour qubits as well as single-qubit gates. In the presence of symmetry-respecting perturbations to the Hamiltonian, the arguments of Sections II and III can still be applied for any finite N . However, if we want to exploit the locality of the perturbation in the vertical direction as well as the horizontal, we need to make some additional arguments. First, we observe that (by Lemma 1) the protected subsystem of correlation space (which corresponds to the ancilla system appearing in the dual picture of MBQC) will have dimension 2 N , and by identifying the action of the symmetry within the protected subsystem with Eq. (71), we can decompose the protected subsystem into N qubits, one associated with each chain.
Our argument now hinges on two observations. First, the dual Hamiltonian of which the dual state is the gapped ground state, as constructed in Sec. III C, is in fact a sum of interactions that are local on the original two-dimensional lattice. Second, the unitary couplings G k appearing in the dual process, which a priori could couple an entire columnar site to the entire N -qubit ancilla system, in fact acts trivially outside an appropriately localized area (see Fig. 10). These observations both follow from the form of the interaction between a columnar site and the ancilla system. (Recall that D 1 and its inverse appeared in the development of the dual picture in Sec. II D, as well as in the construction of the duality transformation D from which the dual Hamiltonian H was obtained in Sec. III C.) It is easily seen that in the present quasi-1D setup, D 1 simply comprises a number of applications of the corresponding operator D (1) 1 for the one-dimensional cluster chain (see Fig. 9). Now, consider a quantum circuit comprising a sequence of gates, and let Q k be the sets of physical (not ancilla) qubits acted on by the corresponding couplings G k in the dual process. Just as in the one-dimensional case, we expect that if R ≡ min k1,k2 dist(Q k1 , Q k2 ) is much larger than the correlation length ξ for the dual state, then the reduced state Tr ( k Q k ) c | Ψ Ψ| on k Q k should be approximately a product state over the Q k 's. Thus, arguing as in the one-dimensional case (see Sec. II F), we find that  performing the measurement sequence on a perturbed resource state corresponds to a noisy quantum circuit, with the noise described by a noise superoperator E k following each gate. Furthermore, E k acts non-trivially only on the same qubits that were acted on by the corresponding gate in the original noiseless quantum circuit. The strength of the noise, as given by E k − I ♦ , is determined by the deviation (in the trace norm) of the reduced density operator on Q k from its unperturbed value [see Eq. (48) in Sec. II F], which should be small for small perturbations.
Let us now estimate the required scale-up in the size of the resource state. We only consider in detail the case of local quantum circuits (i.e. containing only gates acting between nearest-neighbour qubits). As in the onedimensional case (Sec. III D), according to the rigorous factorization result proved in Appendix C, the minimum separation R min between any of the Q k 's required for an accuracy scales like The required scale-up can be expressed in terms of R min , as follows. First, we must ensure that, at each time step, all non-trivial gates are separated by a distance of at least R min . This leads to a scale-up by a factor of ∼ R min in the number of time steps. Then, the buffering between horizontal locations at which nontrivial gates take place implies another factor of R min scale-up. Hence, the total scale-up factor is s ∼ R 2 min . On the other hand, if the quantum circuit that we want to simulate is not already local, then translating it into a local circuit will introduce additional overhead (still scaling at worst polynomially in the number of qubits in the quantum circuit).
B. The 2D cluster model Now we will return to the model we are actually interested in: the 2D cluster model on a square lattice. Investigations of the effect on this model of perturbations [11,[43][44][45] have demonstrated a variety of results depending on the perturbation. Here, we will focus on perturbations respecting an appropriate symmetry. When this symmetry is enforced, the cluster model lies in a robust SPT phase, within which the identity gate is protected and the effective noise model construction of this paper can be applied.
In order to achieve our goal, we will establish an equivalence between the 2D cluster model and a 'quasi-1D' model of the type considered in the previous section. The basic idea is to define a (local) duality transformation U (not the same as the duality transformation D which we have considered previously) which relates the two models. Specifically, we define where (CZ) ij is the controlled-Z gate acting on qubits i and j, and the product is over an appropriate set L comprising nearest-neighbour pairs of qubits. By an appropriate choice of L, we can ensure that applying U to the 2D cluster Hamiltonian turns it into a model of the type we considered in the previous section. Now, let H be a perturbation to the 2-D cluster Hamiltonian. Then UHU † is a perturbation to the quasi-1D model, and the arguments of the preceding section can be applied provided that the perturbation respects the appropriate symmetry. Furthermore, the result (in terms of statistics of measurement outcomes) of performing the adaptive measurement protocol described in Section IV A on the ground state of UHU † , involving measuring the observablesô, must be the same as the effect of performing the same protocol on the ground state of H, but measuring the observables U †ô U. We will now examine in detail this corresponding measurement protocol for perturbations of the 2-D cluster model.
In the quasi-1D resource state, there is a set of 'redundant' qubits which never need to be measured. It turns out to be convenient to assume, however, that we do measure those qubits, in the z basis, and that we do this before any other measurements. We observe that all the measurements that are performed on the quasi-1D state (after applying controlled-Z gates to couple chains where we want to perform a two-qubit gate) are all single-qubit, and are either in the z basis, i.e. measuring Z, or in the x-y plane, i.e. measuring σ θ = (cos θ)X + (sin θ)Y for some angle θ. Hence, the corresponding observables to measure in the 2D cluster state are either of the form Z j (for some qubit j), or for some qubit j and angle θ j , and where N j is some set of neighbouring redundant qubits. But, since we measured the redundant qubits first, they are all now in eigenstates of Z. Therefore, labelling the measured values of Z on the redundant qubits by {z j }, we see that measuring Uô j U † is equivalent to measuring σ (j) θj k∈Nj z k , which in turn is equivalent to measuring σ (j) θj and reinterpreting the measurement outcomes based on the value of k∈Nj z k . Therefore, we have shown that the measurement protocol on the 2D cluster state can be implemented using only single-qubit measurements and adaptivity. It can be checked that the measurement protocol so constructed is essentially the same as the usual one for the 2D cluster state on a square lattice, which is described, e.g. in [1,42].
Finally, let us discuss the required symmetry. The duality transformation U can be used to relate the (Z 2 × Z 2 ) ×N symmetry which protects the quasi-1D model to a corresponding one in the 2D cluster model. The form of the generators of the latter symmetry is shown in Figure  11(a). Let us remark that we can also make similar arguments in the case that the 1D chains are arranged on the 2D square lattice in an unconventional way, for example diagonally [11,46] as shown in Figure 11(b). The advantage of the diagonal layout is that the symmetry [one of the generators of which is depicted in Figure 11(b)] takes a particularly simple form, due to the fact that every non-chain qubit neighbours an even number of chain qubits, and so the Z's that would normally appear on non-chain qubits all cancel. In particular, this symmetry commutes with an especially simple and physically meaningful perturbation, namely a uniform magnetic field in the x direction, i.e. V = B i X i . (The effect of such a perturbation has been studied numerically in [45,47]; the SPT cluster phase persists up to a first-order phase transition at |B| = 1.)

C. Perturbed ground states are universal resources
Let us summarize the conclusions which we obtain from the considerations in Sections IV A and IV B by stating them as a theorem. In combination with the threshold theorem of fault-tolerant quantum computation, this theorem will allow us to deduce that sufficiently small symmetry-respecting perturbations to the 2-D cluster Hamiltonian retain ground states which can serve as universal resources for MBQC.  11. Two possible layouts for the 1D chains (red lines) on a 2D lattice. Each layout is associated with a measurement protocol for MBQC on the 2D cluster state, and with a symmetry group (a representative generator of which is shown in each case). For any layout, we can construct an effective noise model corresponding to performing the associated measurement protocol on a perturbed cluster state, provided that the perturbation respects the associated symmetry.
We consider perturbations to the exact cluster Hamiltonian H C on the 2-D square lattice, which we can take to be a sum of local commuting terms, with an energy gap to the first excited state of 2. Suppose now we consider a perturbed Hamiltonian H = H C + V , and V is a perturbation of the form where Λ is the set of all lattice sites, and each V u is an interaction term supported on the set B(u, r) of sites within some fixed distance r (more generally, interactions decaying exponentially with distance would not present an obstacle to our arguments). We define the local strength of the perturbation by The cluster Hamiltonian H C belongs to a class of Hamiltonians for which it has been shown [48] that the gap is stable to local perturbations, i.e. there exists a threshold η > 0 (depending only on r), such that the gap of the perturbed Hamiltonian is at least 1, provided that J ≤ η. Let us assume that the perturbation V respects an appropriate symmetry group, constructed according to the procedure described in Sec. IV B (such as the one of the symmetry groups depicted in Fig. 11). Suppose then we want to use the perturbed ground state |Ψ to simulate a local quantum circuit containing N qubits, T time steps, and m gates, with the gates drawn from the gate set S comprising single-qubit rotations, a two-qubit entangling gate (as constructed in Sec. IV A), and the non-unitary RESTART gate (which corresponds to the reinitialization of a qubit). We obtain the result Theorem 4. Provided J ≤ η, we can find an appropriate measurement protocol on the ground state |Ψ such that the final reduced state on the output qubits is -close in the trace norm to the outcome of the quantum circuit, with added noise. In each time step t of the equivalent circuit process, the appropriate gates are applied, followed by a noise process described by a superoperator E t . This superoperator can be written as a tensor product E t = A E t,A , where the product is over 'locations', i.e. sets of qubits coupled by a gate in the time step t (each qubit not coupled by a gate in the time step t also counts as a location, but E t,A = I in that case). Thus, the noise has no correlations in space (other than those due to gates acting between qubits) or time. Furthermore, the noise operator E t,A at each location and time is close to the identity superoperator in the diamond norm: for some constant c (dependent only on r). The number of qubits measured n satisfies Proof. The only aspect that we have not previously discussed is the bound Eq. (79). Following the same argument as in the one-dimensional case (Sec. II F), we find [using the analogue of Eq. (48)] that the deviation E t,A − I ♦ is bounded above by ∆ X ≡ ρ X − ρ X,0 1 , where X is the set of qubits in the 2-D lattice that affect the operation of the gate in question, and ρ X and ρ X,0 are the reduced states on X of the perturbed dual state | Ψ and the unperturbed dual state | Ψ 0 respectively. Physically, it should be clear that ∆ X will be small for small perturbations; in Appendix D we demonstrate that, so long as J ≤ η, the inequality ∆ X ≤ cJ holds for some constant c depending only on r.
Now that we have shown that perturbations in the Hamiltonian correspond to noisy quantum circuits, we can invoke the threshold theorem of fault-tolerant quantum computation [18,49,50]. For our purposes, the most suitable version is Theorem 13 of Ref. [18], which we can state as follows: Theorem 5. Let us assume a noise model as described in Theorem 4. Then there exists a threshold η > 0 and a constant α such that, so long as E t,A − I ♦ ≤ η for all A,t, then the following propeties hold. For any > 0, and any local quantum circuit C made from gates drawn from S (with N qubits, T time steps, and m gates), there exists another local circuit C with gates drawn from S, such that C with noise produces the same result (in terms of the probability distribution for the final readout, and up to an error ) as C without noise. The scale-up factors for the number of qubits, the number of time steps, and the number of gates are all bounded by (const.) × log α (m/ ). Combining Theorems 4 and 5, we obtain: Theorem 6. Consider the perturbed model H = H C + V as described above. Then there exists some threshold η = min{η, η /c} > 0 (depending only on r) with the following property. Provided that J < η , then for any local quantum circuit C (with N qubits, T time steps, and m gates), with gates drawn from S, we can find an appropriate measurement protocol on the perturbed ground state |Ψ such that the result is equivalent (in terms of the probability distribution for the final readout, and up to an error which can be made arbitrarily small) to the outcome of the original quantum circuit. As m → ∞ with held fixed, the number of measured qubits n satisfies This is sufficient to show that the perturbed ground states remain universal resources, which is Theorem 3 as stated in Sec. I A; it is the main result of this paper.

V. CONCLUSION
In this paper, we have developed a framework to characterize the effectiveness of measurement protocols for MBQC with SPT-ordered ground states of quantum spin systems. This has allowed us to prove the universality for MBQC of the ground states of perturbed versions of the 2-D cluster Hamiltonian, provided that the perturbation is sufficiently small and respects an appropriate symmetry.
The type of SPT order that we have presented here is that present in one-dimensional systems, which is related to a nontrivial factor system (also known as a 2cocycle). It is for this reason that, in order to establish universality in two-dimensional systems, we had to treat them as 'quasi-one-dimensional' and assume an extensive symmetry group (Z 2 × Z 2 ) ×N , which grows with the vertical extent of the system. For standard, non-extensive symmetries in two dimensions, SPT orders can be related to 3-cocycles [51,52], but it remains to be seen whether similar connections can be drawn between such two-dimensional SPT order and MBQC.
Finally, we note that if MBQC in ground states of quantum spin systems is to be a robust form of quantum computation, then it must be possible in the presence of arbitrary (not necessarily symmetry-respecting) local perturbations to the Hamiltonian, as well as at nonzero temperature. Non-symmetry-respecting perturbations break the symmetry that is essential to our argument; the difficulty of extending our treatment to nonzero temperature was discussed in Sec. III E. Nor have we considered the effect of non-ideal measurements, or of decoherence of the resource state taking place during the course of the measurement protocol. Therefore, it remains an open question whether fault-tolerant MBQC is possible with such imperfections. g ∈ G. Therefore, since Λ is the unique eigenvector of A † with eigenvalue 1, it must satisfy V g (Λ) = χ(g)Λ for some scalars χ(g). The fact that V g is trace preserving ensures that χ(g) = 1. Hence we find that Λ commutes with V (g) ⊗ I; it follows by Schur's Lemma that Λ factorizes as Λ = Ω ⊗ Λ for some density operator Λ, and where Ω = I/ |G| is the maximally-mixed state on the protected subsystem (recall that the |G| is the dimension of the protected subsystem). It follows that there is a |G|-fold degeneracy in the entanglement spectrum throughout the SPT phase, generalizing the 2-fold degeneracy in the (Z 2 × Z 2 )-protected Haldane phase [35].
We can define the transfer channel corresponding to the dual FCS (generated by A) according to Observe that A unital implies that A is also unital, and that for any operator σ acting on the junk subsystem. Hence, any eigenvalue of A must also be an eigenvalue of A. From this we can see that the dual FCS is also a pure FCS, and Λ is the unique fixed point of A † ; thus, the entanglement spectrum of the dual state is the same as that of the original state, but with the |G|-fold degeneracy removed 2 . The respective correlation lengths obey the inequality ξ ≤ ξ. We now outline how the arguments of Section II D can be put on a rigorous footing within the pFCS formalism. Given everything that we have established so far, it can be shown that the reduced state ρ n of the original FCS on a block of n adjacent sites can be obtained from the corresponding reduced state ρ n of the dual FCS according to the quantum circuit shown in Fig. 12(a). Assuming that we choose n large enough that all measurements take place within this block of n sites, we can then make arguments similar to those of Sec. II D, and we find that the result of the adaptive measurement protocol is equivalent to a sequence of interactions between the dual state and an ancilla particle as shown in Fig.  12(b), with the unitary interactions G k defined as they were previously.

Appendix B: Connection with the Kennedy-Tasaki transformation
The Kennedy-Tasaki (KT) transformation [30] is a non-local unitary transformation which transforms a 2 A similar property was found numerically for the Kennedy-Tasaki transformation in [33]. We will discuss the connection between that transformation and our "dual state" in Appendix B. 12. (a) The quantum circuit generating the reduced state on n sites for the original FCS. The ancilla particle is initially in the maximally-mixed state Ω, and the symbol indicates that it should be traced out at the end. (b) The adaptive measurement protocol acting on the original FCS is equivalent to a dual process acting on the dual state, comprising a sequence of unitary interactions G k coupling the dual state to an ancilla particle.
spin-1 chain in the Z 2 × Z 2 symmetry-protected Haldane phase into a system where the symmetry is spontaneously broken in the bulk. In this section, adapting Ref. [33], we will define a generalized version of the Kennedy-Tasaki transformation, which can be applied to any system in the non-trivial SPT phase with respect to an on-site representation U (g) = [u(g)] ⊗N of the symmetry group G = Z 2 × Z 2 (like the original KT transformation, our generalization is defined for finite chains with open boundary conditions and no terminating particles). We will show that, when the ground state of the original system can be expressed as a pFCS, the ground state of the KT-transformed system is essentially the same as the state which we have referred to throughout this paper as the "dual state". We expect that for general ground states the situation should be qualitatively similar.
Observe that the symmetry group is generated by two commuting elements x and z; hence for any g ∈ G, we can write g = x m(g) z n(g) for some m(g) and n(g) taking values of 0 or 1. We will choose to write the unique non-trivial irreducible projective representation as V (g) = X m(g) Z n(g) , where X and Z are the appropriate Pauli operators. This is a rephasing of the Pauli representation V P defined in Eq. (10); thus the factor system is different to, but in the same cohomology class as, that of V P . (The construction of the dual state does depend on the specific choice of representative factor system for a cohomology class, although in a fairly trivial way; the present choice is the one that will ensure that the Kennedy-Tasaki transformation reproduces the dual state exactly.) The generalized Kennedy-Tasaki transformation D KT is then defined as follows: where D kl is a unitary coupling the two sites k and l according to Notice that all the operators D kl appearing in the product Eq. (B1) commute. In the case that the particles are spin-1's, with the Z 2 ×Z 2 rotation symmetry consisting of π rotations about the x, y, and z particles (the simultaneous eigenbasis of the symmetry for a single site is then {|x , |y , |z }, where |α , α = x, y, z is the 0 eigenstate of the spin component operator S α ), the transformation D KT reduces to the standard Kennedy-Tasaki transformation. Specifically, Eqs. (B1) and (B2) correspond to Eqs. (4) and (5) in Ref. [33]. For open boundary conditions, we expect there to be an approximate (becoming exact in the limit as the length of the chain goes to infinity) four-fold degeneracy, due to the two-fold degeneracy associated with each edge. An appropriate analogue of the SPT pFCS for this choice of boundary conditions is as follows: the low-energy subspace is spanned by states of the form where the MPS tensor A obeys the usual symmetry condition, where |L * and R * | are fixed end vectors, while |L and R| are allowed to vary (within the twodimensional irrep space) in order to generate the fourdimensional low-energy subspace. Now we want to examine what happens to a state of the form Eq. (B3) under D KT . Suppose we set R| = +|, |L = |0 . It is then straightforward to show [using the symmetry condition Eq. (12) and the decomposition Eq. (15)] that applying all the pairwise interactions D 1l involving the first site gives the result l>1 D l1 |Ψ(0, +) = Continuing this process, we find that which is the dual state. To obtain the other states within the low-energy subspace, it is sufficient to observe that, using the symmetry condition on the tensor A, we have and that D KT commutes with u(g) for all g ∈ Z 2 × Z 2 . Hence, the other basis states for the Kennedy-Tasaki transformed system can be obtained from the dual state by application of a symmetry operation [u(g)] ⊗N for some g ∈ G. The dual state thus represents one of the four distinct symmetry-breaking states within the lowenergy subspace of the transformed system. Let us also note that, although they relate to different boundary conditions, the unitary transformation D that we introduced in Sec. III B is equivalent to the KT transformation D KT , in the following sense. Indeed, an important property of D KT is that for any local symmetryrespecting observableô, D KTô D † KT remains local and symmetry-respecting. It turns out that this also holds for D, in the case of observablesô acting in the bulk (the fact that DôD † is local was established in Sec. III C; it can be shown that DôD † still respects the on-site symmetry as well). Thus, like D KT , the transformation D can be applied to yield a local, symmetry-respecting Hamiltonian H in the bulk. It can be shown H is precisely the KT transformed Hamiltonian [53]. Thus, in line with the results for pFCS described above, we expect H to have four degenerate, locally distinguishable symmetrybreaking ground states in the bulk. On the other hand, when the transformation D is applied to the edge interactions (those which couple the ends of the chain to the terminating particles), the result need not respect the symmetry. Therefore, including the edge interactions favours one of the symmetry-breaking ground states over the others and leads to the non-degeneracy of the dual state | Ψ .
Appendix C: The factorization condition for ground states of a local Hamiltonian In this Appendix, we will show how to derive the approximate factorization result Eq. (68) for a nondegenerate gapped ground state |Ψ of a local Hamiltonian. Our main tool is the existing theorem on the exponential decay of correlation functions for such a ground state [54][55][56]. This theorem states that there exists a correlation length ξ and a function f (x), with f (x) = O[exp(−x/ξ)] as x → ∞, such that for any sets of lattice sites X and Y , and observables A X and B Y supported on X and Y respectively, we have where · Ψ denotes the expectation value of an observable with respect to |Ψ , and |X| denotes the number of sites contained in the set X. Suppose that instead of two sets of lattice sites, we have m disjoint sets X 1 , . . . , X m . Let N = m k=1 |X k | be the total number of lattice sites contained within all of the X k 's. We can obtain the following corollary: Lemma 3. For any observables A X1 , . . . , A Xm supported on the respective sets, where R is the smallest distance between any two of the X k 's, i.e. R ≡ min k =l dist(X k , X l ).
Proof. Without loss of generality, we can assume that A X k = 1 for all k. Observing that | A X1 | ≤ A X1 = 1, we have that The two-body result implies that the first term in the right-hand side of Eq. (C3) is bounded by |X 1 |f (R). Continuing in this way, the lemma follows by induction. Now we want to show that the reduced state on the union of the X k 's is close to a product state. To do this, we make use of the following lemma. We remind the reader that here we make use of both the trace norm A 1 = Tr √ A † A and the spectral norm A = max ψ|ψ =1 A|ψ . They are both distinct from the norm induced by the Hilbert-Schmidt inner product.
Lemma 4. Consider a vector space of dimension D. Then there exists a basis {E i |i = 1, ..., D 2 } for the space of linear operators supported on the site, orthonormal with respect to the Hilbert-Schmidt inner product A, B = Tr(A † B), and comprising Hermitian operators E i such that such that E i 1 E i = 1.
Theorem 7. Let |Ψ be the non-degenerate gapped ground state of a local Hamiltonian. Let ρ be the reduced state of |Ψ on m k=1 X k , and let ρ k be the reduced state on X k . Then where ρ prod = k ρ k .