Patent Publication Number: US-2023162081-A1

Title: Fast preparation of highly-entangled quantum states

Description:
RELATED APPLICATIONS 
     This application claims the benefit of U.S. Provisional Application No. 63/282,873, filed on Nov. 24, 2021, and U.S. Provisional Application No. 63/296,715, filed on Jan. 5, 2022. The entire teachings of the above applications are incorporated herein by reference. 
    
    
     BACKGROUND 
     Highly entangled quantum states are useful in such application as quantum computing ad quantum metrology. Producing such quantum states remain a challenge. 
     SUMMARY OF THE INVENTION 
     In an example embodiment, the present invention is a method of preparing an entangled state of a plurality of qubits. The method comprises: preparing an array of qubits, each qubit being in a ground qubit state and having an excited qubit state, the array comprising: a plurality of first qubits and a plurality of second qubits, each first qubit disposed at a vertex of a first lattice, each second qubit disposed at a vertex of a second lattice, wherein each first qubit has at least one nearest neighbor second qubit, wherein for any two first qubits, a spacing between said first qubits and their respective nearest neighbor second qubits is the same. The method further comprises: -causing a first single-site rotation of the plurality of the first qubits and of the plurality of the second qubits, thereby causing each first and each second qubit to transition into a superposition of the ground qubit state and the excited qubit state; -evolving the array of qubits for a pre-determined time τ, thereby producing a cluster state of the array of qubits; and -measuring the state of each first qubit, wherein the measuring comprises reversing the single-site rotation of the plurality of the first qubits, thereby preparing an entangled state of the plurality of second qubits. 
     In another example embodiment, the present invention is a method of preparing an entangled state of a plurality of qubits. The method comprises preparing an array of qubits, each qubit having a first state and a second state, the array comprising: a plurality of first qubits and a plurality of second qubits, wherein each first qubit has at least one neighboring second qubit, wherein each first qubit is controllably in communication with its respective neighboring second qubits. The method further comprises: -causing each second qubit to transition into a superposition of the first state and the second state, thereby preparing a product state; -applying a controlled-Z gate to each of the first qubits and its nearest neighbor second qubit; and -measuring the state of each first qubit, thereby preparing an entangled state of the plurality of second qubits. 
     In another example embodiment, the present invention is a method of preparing an entangled state of a plurality of qubits. The method comprises preparing an array of qubits, each qubit having a ground qubit state and an excited qubit state, the array comprising: a plurality of first qubits and a plurality of second qubits, each first qubit disposed at a vertex of a first lattice, each second qubit disposed at a vertex of a second lattice, wherein each first qubit has at least one nearest neighbor second qubit, wherein for any two first qubits, a spacing between said first qubits and their respective nearest neighbor second qubits is the same, wherein the plurality of the first qubits is in an “input” state, and each of the plurality of second qubits is in the ground state. The method further comprises: -causing a single-site rotation of the plurality of the second qubits, thereby causing each second qubit to transition into a superposition of the ground qubit state and the excited qubit state; -evolving the array of qubits for a pre-determined time τ, thereby producing an entangled state of the plurality of the first qubits and the plurality of the second qubits; and -measuring the state of each first qubit, wherein the measuring comprises applying a second single-site rotation of the plurality of the first qubits, thereby preparing an output state of the plurality of the second qubits. 
     In another example embodiment, the present invention is a device, comprising an array of qubits, each qubit having a ground qubit state and having an excited qubit state, the array comprising: a plurality of first qubits and a plurality of second qubits, each first qubit disposed at a vertex of a first lattice, each second qubit disposed at a vertex of a second lattice, wherein each first qubit has at least one nearest neighbor second qubit, wherein for any two first qubits, a spacing between said first qubits and their respective nearest neighbor second qubits is the same. At least some of the plurality of second qubits is in a GHZ state, and the fidelity of the Greenberger-Horne-Zeilinger (GHZ) state per site is at least 0.99. 
    
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
       The patent or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawing(s) will be provided by the Office upon request and payment of the necessary fee. 
         FIG.  1 A ,  FIG.  1 B ,  FIG.  1 C , and  FIG.  1 D  collectively illustrate a protocol for fast preparation of GHZ state by measuring 1D cluster state. 
         FIG.  2 A  and  FIG.  2 B , collectively illustrate a transition from toric and color code to non-Abelian topological order. 
         FIG.  3 A  and  FIG.  3 B , collectively, illustrate  FIG.  3    a protocol for realizing a fractal spin liquid with fracton order. 
         FIG.  4 A  and  FIG.  4 B , collectively, illustrate Z 3  and S 3  non-Abelian topological order. 
         FIG.  5 A  and  FIG.  5 B  each depicts a stabilizer useful in implementing a fraction order. 
         FIG.  6    is a plot of a distribution of the correlation length ξ over a distance n= ξ  for t=0.99×t SPT . The black line is the predicted normal distribution. 
         FIG.  7    depicts two sublattices of a square lattice. 
         FIG.  8    depicts four sublattices of a square lattice. 
         FIG.  9    illustrates a cluster state useful for obtaining a toric code on the diamond lattice having two sublattices. 
         FIG.  10 A  and  FIG.  10 B , collectively, show the connectivity of the cluster state used to prepare the X-cube model. The A sites at the center of the cube (red) are connected to the B sites situated at the twelve edges surrounding the cube (blue). Together, they form the FCC lattice. 
         FIG.  11    depicts a Lieb lattice. 
         FIG.  12    shows a stabilizer used for preparation of S 3  topological order. 
         FIG.  13    shows a stabilizer of the ground state of the    3  toric code. 
         FIG.  14    depicts sublattices C (red) and D (purple) which form a shifted Lieb lattice (dotted lines) are added. Each site hosts a qubit. 
         FIG.  15   ,  FIG.  16   , and  FIG.  17    collectively show stabilizers used for preparation of S 3  topological order. 
         FIG.  18   ,  FIG.  19   ,  FIG.  20   ,  FIG.  21   , and  FIG.  22    collectively depicts projectors used for preparation of S 3  topological order. 
         FIG.  23    depicts a pattern of C operators related to the quantum double model D(S3). 
         FIG.  24    and  FIG.  25    collectively depict projectors related to the quantum double model D(S3). 
         FIG.  26    and  FIG.  27    are the vertex term and the plaquette terms, respectively, related to the quantum double model D(S3). 
         FIG.  28    depicts the relation between the Matrix-Product State representation of the cluster state and the Matrix-Product Operator representation of the Kramers-Wannier duality 
         FIG.  29    depicts the protocol to perform the Kramers-Wannier duality using finite-time unitary evolution and measurements, up to single-site rotations. The foregoing will be apparent from the following more particular description of example embodiments of the invention, as illustrated in the accompanying drawings in which like reference characters refer to the same parts throughout the different views. The drawings are not necessarily to scale, emphasis instead being placed upon illustrating embodiments of the present invention. 
     
    
    
     DETAILED DESCRIPTION OF THE INVENTION 
     As used herein, a “lattice” can be a unidimensional, a two-dimensional, or a three-dimensional lattice. 
     As used herein, a “qubit” refers to any physical implementation of a quantum mechanical object that can be in any one of at least two states or any superposition thereof. In certain embodiments, a qubit includes a qutrit, i.e. a quantum mechanical object that can be in any one of three states or any superposition thereof. Physical implementations of a qubit include neutral atoms, ions, superconducting electronic circuit, etc. 
     As used herein, the phrase “nearest neighbor(s)” of a given site in one sublattice is defined as the collection of sites in the other sublattice with the shortest distance to said site. 
     A description of example embodiments of the invention is provided below. 
     Long-range entangled quantum states—like cat states and topological order—are key for quantum metrology and information purposes, but they cannot be prepared by any scalable unitary process. Intriguingly, using measurements as an additional ingredient could circumvent such no-go theorems. However, efficient schemes are known for only a limited class of long-range entangled states, and their implementation on existing quantum devices via a sequence of gates and measurements is hampered by high overheads. Described herein is a method and a device that resolve these problems, implementing scalable preparation of a broad range of long-range entangled states with the use of existing experimental platforms. The two-step process described herein can be implemented in Rydberg atom arrays, only requiring time-evolution under the intrinsic atomic interactions, followed by measuring a single sub-lattice (by using, e.g., two atom species). Remarkably, this protocol can prepare the 1D Greenberger-Horne-Zeilinger (GHZ) ‘cat’ state and 2D toric code with fidelity per site exceeding 0.9999, and a 3D fracton state with fidelity ≳0.998. In light of recent experiments showcasing 3D Rydberg atom arrays, this paves the way to the first experimental realization of fracton order. In addition to the examples based on efficiently preparing and measuring cluster states, also proposed herein is a multi-step procedure to create S 3  and D 4  non-Abelian topological order in Rydberg atom arrays and other quantum devices -offering a route towards universal topological quantum computation. 
     The defining property of long-range entangled states is that they cannot be prepared efficiently by unitary processes, that is, prepared within a fixed time window independent of system size. Remarkably, measurements can circumvent this issue in quantum devices. In particular, certain long-range entangled states can in principle be prepared in a fixed amount of time by first acting with simple quantum gates (typically controlled-Z or Ising)—preparing the so-called cluster state. Subsequently measuring a subset of quantum bits leaves the remainder in a long-range entangled state. 
     Referring now to  FIGS.  1 A through  1 D , which illustrate a protocol for fast preparation of a GHZ state by measuring 1D cluster state. Panel (a) shows that the protocol starts with a product state in the X-basis which we let time-evolve under an Ising interaction. After a time t=t SPT , this produces the cluster state. Measuring the red sites in the X-basis leads to a cat state. Using the information of the measurement outcomes, one can flip the appropriate spins to obtain the GHZ state on the blue sites. Panel (b) If no time-evolution is performed with exactly t=t SPT , the post-measurement state is only an approximate cat state. Nevertheless, the resulting correlation length (within which there is long-range order) is very large even for moderate time deviations. Panel (c) describes a tensor-network simulation of the protocol for Rydberg atoms interacting with a 1/r 6  van der Waals interaction. This simulation was done for Rb 70S 1/2  with two different lattice spacings, which give different results due to incorporating the fact that the initial product state preparation pulses are not instantaneous. Panel (d) describes contributions from longer-range van der Waals interactions can be systematically suppressed by interspersing time-evolution with X-pulses (orange) on particular sublattices. This particular example cancels out couplings at distance r=2a, although this was not used to achieve the results in (c). 
     Some well-known examples are the cluster state on a 1D chain or the 2D Lieb lattice, where single-site measurements transform these into the GHZ ‘cat’ state or the toric code, respectively. However, even these known cases have resisted experimental implementation. Indeed, the infrastructure necessary to apply such quantum gates can require a high overhead, and setting aside a fraction of the qubits to be measured is costly in emerging quantum computers—with recent realizations of the toric code phase on quantum devices opting for a unitary protocol. Moreover, going beyond such known examples, an important open question is what other exotic state with potentially valuable properties could be prepared within a measurement-based protocol. 
     In the present disclosure, both these issues are being addressed. In terms of novel examples, 3D cluster states are being presented where single-site measurements produce the elusive fracton order—an intriguing family of 3D states which have resisted experimental detection and hold promise as optimal quantum hard drives. Moreover, going beyond the cluster state framework, there is disclosed a method of obtaining a non-Abelian topological order via multiple alternations between gates and measurements, establishing the first proposal for efficiently obtaining such an exotic state in any platform. The second issue addressed is that of implementability: all of the aforementioned examples—old and new—can be realized using existing Rydberg atom technology. In particular, the entangling operations are permitted to naturally emerge from the time-evolution under the intrinsic interactions between Rydberg atoms, rather than constructed piecemeal via the gates of a quantum computer, thereby trading in the overhead of the latter for scalability. 
     Rydberg atom arrays offer an ideal platform satisfying the two requirements of Ising interactions and high-resolution measurements. These are tunable lattices of atoms trapped by optical tweezers, each serving as an effective qubit consisting of its ground state (n=0) and a particular Rydberg excited state (n=1). When two atoms are simultaneously excited, they experience a van der Waals interaction V(r)˜1/r 6 , giving rise to a Hamiltonian of the form: 
     
       
         
           
             
               
                 
                   
                     H 
                     = 
                     
                       
                         
                           ∑ 
                           i 
                         
                         
                           
                             
                               Ω 
                               i 
                             
                             2 
                           
                           ⁢ 
                           
                             X 
                             i 
                           
                         
                       
                       - 
                       
                         
                           ∑ 
                           i 
                         
                         
                           
                             δ 
                             i 
                           
                           ⁢ 
                           
                             n 
                             i 
                           
                         
                       
                       + 
                       
                         
                           1 
                           2 
                         
                         ⁢ 
                         
                           
                             ∑ 
                             
                               i 
                               ≠ 
                               j 
                             
                           
                           
                             
                               V 
                               ij 
                             
                             ⁢ 
                             
                               n 
                               i 
                             
                             ⁢ 
                             
                               n 
                               j 
                             
                           
                         
                       
                     
                   
                   , 
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   1 
                 
               
             
           
         
       
     
     where X, Y, Z denote the Pauli matrices. Heref Ω i  and δ i  are tuneable, and by shifting the latter, we obtain an explicit Ising interaction 
     
       
         
           
             
               
                 
                   1 
                   2 
                 
                 ⁢ 
                 
                   
                     ∑ 
                     
                       i 
                       ≠ 
                       j 
                     
                   
                   
                     
                       J 
                       ij 
                     
                     ⁢ 
                     
                       Z 
                       i 
                     
                     ⁢ 
                     
                       Z 
                       j 
                     
                     ⁢ 
                         
                     with 
                     ⁢ 
                         
                     
                       J 
                       ij 
                     
                   
                 
               
               = 
               
                 
                   1 
                   4 
                 
                 ⁢ 
                 
                   V 
                   ij 
                 
               
             
             , 
           
         
       
     
     which will play the role of our ‘entangling gates’. Recent years have seen a rapid growth in available system sizes, exploring the many-body quantum physics of Equation 1 in 1D and 2D arrays of many hundreds of atoms, and even demonstrating 3D arrays. Confirming the platform&#39;s promise for realizing exotic states was recently established the onset of toric code topological order. 
     In the simplest instances, the state preparation procedure is a two-step protocol of time-evolution and measurement. At the end of the first step, high-quality cluster states are obtained despite the longer-range van der Waals corrections. This is due to a judicious choice of lattices, as well as reinterpreting the power of cluster states as being a property of their symmetry-protected topological (SPT) phase. While cluster states are already of interest (serving as resources for measurement-based quantum computation), to produce long-range entanglement an ability to measure only a subsystem of its atoms is needed. Recent breakthrough experiments have established that this is possible in dual-species arrays, where one species is transparent to measurements of the other. This way, the 1D GHZ, the 2D tonic code and 3D fracton states are produced. 
     Non-Abelian D 4  or S 3  topological order can be obtained by using a four-step protocol where one time-evolves and measures twice, creating the color code or Z 3  toric code as an intermediate state. This can be conceptually interpreted as sequentially gauging certain symmetries, thereby expanding the possible states of matter which can be efficiently prepared by finite-time evolution and measurement. 
     The 1D Cluster and Ghz States 
     To describe how to efficiently prepare a Schrödinger cat state, first consider the idealized setting of nearest-neighbor Ising interactions, after which the effect of imperfections can be discussed. Starting with a product state in the Pauli-X basis, a chain of N qubits is time-evolved: 
       |ψ( t )) e   −itJΣ     n     Z     n     Z     n+1   |−) ⊗N .   Equation 2
 
     After a time 
     
       
         
           
             
               
                 t 
                 SPT 
               
               := 
               
                 π 
                 
                   4 
                   ⁢ 
                   
                     
                       ❘ 
                       &#34;\[LeftBracketingBar]&#34; 
                     
                     J 
                     
                       ❘ 
                       &#34;\[RightBracketingBar]&#34; 
                     
                   
                 
               
             
             , 
           
         
       
     
     the resulting wave function forms the so-called 1D cluster state. Its defining characteristic is: 
         Z   n−1   X   n   Z   n+1 |ψ( t   SPT ) =|ψ( t   SPT )).   Equation 3
 
     This implies that the entangled state has the remarkable property that measuring, for example, every even site in the X-basis leaves the remaining qubits in a Schrödinger cat state. Indeed, replacing X 2n →±1 by the measurement outcome, the post-measurement state satisfies 
         Z   2n−1   Z   2n+1 |ψ out   =±|ψ out   
 
     In addition, one can show that the state is symmetric, i.e.,  Z 2n−1   =0. The only state satisfying these properties, is the (disordered) GHZ state, a Schrödinger cat state where neighboring spins on the odd sublattice are perfectly aligned or anti-aligned depending on the measurement outcome of the intervening site, as illustrated in  FIG.  1 A . Furthermore, since this information is available, one can always apply single-site spin flips to bring it to the ideal GHZ state: 
     
       
         
           
             
               
                 
                   
                     
                       
                         
                           ❘ 
                           &#34;\[LeftBracketingBar]&#34; 
                         
                         GHZ 
                       
                       〉 
                     
                     = 
                     
                       
                         
                           
                             ❘ 
                             &#34;\[LeftBracketingBar]&#34; 
                           
                           
                             
                               
                                 ↑ 
                                 ↑ 
                                 … 
                                 ↑ 
                               
                               〉 
                             
                             + 
                             
                               
                                 ❘ 
                                 &#34;\[LeftBracketingBar]&#34; 
                               
                               
                                 ↓ 
                                 ↓ 
                                 … 
                                 ↓ 
                               
                             
                           
                         
                         〉 
                       
                       
                         2 
                       
                     
                   
                   , 
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   4 
                 
               
             
           
         
       
     
     although it is a Schrödinger cat state independent of this procedure. 
     To make this procedure practical, understanding is needed of its sensitivity to deviations from the above idealization, such as evolving with additional couplings or for the incorrect time. The ability to create a cat state via measurement is a stable property of the phase of matter that the cluster state belongs to. More precisely, it is known that in addition to the global Ising symmetry Π n  X n , the cluster state is symmetric under spin-flips on even and odd sites, separately; deformations of the cluster state that preserve this symmetry define its so-called SPT phase. 
     As a first indication of the importance of this symmetry, it is explored what happens when the symmetry is slightly broken by time-evolving for the wrong duration. If t is not a multiple of t SPT , the state does not enjoy the aforementioned sublattice symmetry. Indeed, a straightforward computation shows that (ignoring boundary effects and taking|1−t/t SPT |&lt;1/2): 
     
       
         
           
             
               
                 
                   
                     〈 
                     
                       ψ 
                       ⁡ 
                       ( 
                       t 
                       ) 
                     
                     
                       ❘ 
                       &#34;\[LeftBracketingBar]&#34; 
                     
                     
                       
                         ∏ 
                         n 
                       
                       
                         
                           X 
                           
                             2 
                             ⁢ 
                             n 
                           
                         
                         
                           ❘ 
                           &#34;\[RightBracketingBar]&#34; 
                         
                         
                           ψ 
                           ⁡ 
                           ( 
                           t 
                           ) 
                         
                         〉 
                       
                     
                   
                   = 
                   
                     
                       
                         sin 
                         ⁡ 
                         ( 
                         
                           2 
                           ⁢ 
                           J 
                           ⁢ 
                           t 
                         
                         ) 
                       
                       N 
                     
                     ∼ 
                     
                       
                         e 
                         
                           
                             - 
                             N 
                           
                           / 
                           
                             
                               ξ 
                               s 
                             
                             ( 
                             t 
                             ) 
                           
                         
                       
                       . 
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   5 
                 
               
             
           
         
       
     
     This defines a length scale ξ s (t) within which the state looks symmetric, with ξ s (t SPT )=∞. In line with the idea that the symmetry is the determining factor, it is surmised that after measuring every other site, an imperfect cat state is obtained with correlation length Indeed, it is confirmed by a full solution of the model that the resulting state has long-range order for  Z 2m−1 Z 2n−1    within a correlation length closely matching ξ s , beyond which it decays to zero. This is plotted in  FIG.  1 B , where we see that δ grows rapidly as we approach t→t SPT . E.g., a 0.1% deviation results in long-range order over nearly half a million atoms; a 1% error still gives δ/(2a)≈4000 (the original lattice having spacing a). Hence, even in the presence of timing imperfections, this procedure can produce large cat states in a finite amount of time. 
     A chain of Rydberg atoms gives a good approximation to Equation 9 with 
     
       
         
           
             
               J 
               = 
               
                 
                   V 
                   ⁡ 
                   ( 
                   a 
                   ) 
                 
                 4 
               
             
             , 
           
         
       
     
     with additional longer-range corrections. In case we use a single species of Rydberg atoms, the first correction occurs at distance r=2a with a van der Waals coupling V(2a)=V(a)/2 6 , giving an additional evolution by 
     
       
         
           
             
               e 
               
                 
                   - 
                   
                     1 
                     
                       2 
                       6 
                     
                   
                 
                 ⁢ 
                 itJ 
                 ⁢ 
                 
                   
                     ∑ 
                     n 
                   
                   
                     
                       z 
                       n 
                     
                     ⁢ 
                     
                       z 
                       
                         n 
                         + 
                         2 
                       
                     
                   
                 
               
             
             . 
           
         
       
     
     However, one can show that since this Ising coupling preserves the spin-flip symmetry on even and odd sublattices, measuring the even sites still leads to a cat state with an infinite correlation length. It is thus only the third-nearest-neighbor correction which will lead to an imperfect cat state. Similar to Equation 5, we find that this leads to a correlation length 
     
       
         
           
             
               
                 
                   
                     ξ 
                     
                       2 
                       ⁢ 
                       a 
                     
                   
                   = 
                   
                     
                       1 
                       
                         2 
                         ⁢ 
                         
                           
                             ❘ 
                             &#34;\[LeftBracketingBar]&#34; 
                           
                           
                             ln 
                             ⁢ 
                             
                               cos 
                               ⁡ 
                               ( 
                               
                                 π 
                                 / 
                                 
                                   ( 
                                   
                                     2 
                                     ⨯ 
                                     
                                       3 
                                       6 
                                     
                                   
                                   ) 
                                 
                               
                               ) 
                             
                           
                           
                             ❘ 
                             &#34;\[RightBracketingBar]&#34; 
                           
                         
                       
                     
                     ≈ 
                     
                       
                         2 
                         ⨯ 
                         
                           10 
                           5 
                         
                       
                       . 
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   6 
                 
               
             
           
         
       
     
     A near-perfect cat state is thus obtained. 
     In addition to knowing the size of the cat state that was prepare, it can be asked how close to the fixed-point GHZ limit it is. A useful fact is that, to a very high degree, the GHZ-fidelity of the final state can be equated with the cluster-state-fidelity of the pre-measurement state. To characterize the latter, the fact can be used that the expectation value of the stabilizer Z n−1 X n Z n+1  gives a lower bound on the fidelity. At first sight, this expectation value is affected by the second-nearest neighbor coupling at r=2a. However, in this 1D case, it can be shown that its effect is to simply tilt the Ising order parameter away from Z n , which can be easily corrected by a single-site rotation at the end of the time-evolution. In conclusion, even for a single species the dominant correction is at r=3a, leading to a cluster stabilizer 
     
       
         
           
             
               
                 
                   
                     〈 
                     
                       
                         Z 
                         
                           n 
                           - 
                           1 
                         
                       
                       ⁢ 
                       
                         X 
                         n 
                       
                       ⁢ 
                       
                         Z 
                         
                           n 
                           + 
                           1 
                         
                       
                     
                     〉 
                   
                   ≈ 
                   
                     
                       cos 
                       2 
                     
                     ( 
                     
                       π 
                       
                         2 
                         ⨯ 
                         
                           3 
                           6 
                         
                       
                     
                     ) 
                   
                   ≈ 
                   
                     
                       0 
                       . 
                       9 
                     
                     ⁢ 
                     9 
                     ⁢ 
                     9 
                     ⁢ 
                     9 
                     ⁢ 
                     9 
                     ⁢ 
                     
                       5 
                       . 
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   7 
                 
               
             
           
         
       
     
     This suggests a fidelity per site of 
     
       
         
           
             
               
                 
                   
                     
                       
                         
                           ❘ 
                           &#34;\[LeftBracketingBar]&#34; 
                         
                         
                           〈 
                           
                             cluster 
                             | 
                             
                               ψ 
                               ⁡ 
                               ( 
                               
                                 t 
                                 
                                   S 
                                   ⁢ 
                                   P 
                                   ⁢ 
                                   T 
                                 
                               
                               ) 
                             
                           
                           〉 
                         
                         
                           ❘ 
                           &#34;\[RightBracketingBar]&#34; 
                         
                       
                       2 
                     
                     N 
                   
                   ⁢ 
                   
                     
                       
                         
                           &gt; 
                         
                       
                     
                     
                       
                         
                           ~ 
                         
                       
                     
                   
                   ⁢ 
                   
                     
                       1 
                       + 
                       
                         〈 
                         
                           
                             Z 
                             
                               n 
                               - 
                               1 
                             
                           
                           ⁢ 
                           
                             X 
                             n 
                           
                           ⁢ 
                           
                             Z 
                             
                               n 
                               + 
                               1 
                             
                           
                         
                         〉 
                       
                     
                     n 
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   8 
                 
               
             
           
         
       
       
         
           
             
               
                 
                   ≈ 
                   
                     0.9 
                     9 
                     ⁢ 
                     9 
                     ⁢ 
                     9 
                     ⁢ 
                     9 
                     ⁢ 
                     7 
                     ⁢ 
                     5 
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   9 
                 
               
             
           
         
       
     
     This is such a high fidelity. One might worry that at this point other imperfections—which may usually reasonably be ignored -become dominant. One such imperfection is the fact that the initial product state is prepared by a pulse, which has a finite time on the order of tens of nanoseconds. A full tensor network simulation was performed of the proposed state preparation for the case of Rb 70S 1/2 , incorporating this imperfection. As we see in  FIG.  1 C , very high fidelities can be reached. The fact that the resulting fidelity is one order of magnitude worse than Equation 9 is due to the finite pulse time. Indeed, it was confirmed that artificially reducing the pulse time recovers Equation 9. 
     Toric and Color Code 
     Having established cat states in 1D, we now turn to 2D topological order. 
     Referring now to  FIG.  2 A  and  FIG.  2 B . Panel illustrates (a) placement of Rydberg atoms on the vertices (red) and bonds (blue) of the honeycomb lattice. Initializing into the product state|−   ⊗N  and time-evolving for t=π/V(a), the cluster state is obtained whose two types of stabilizers are depicted. Measuring the Rydberg atoms on the red (vertex) sites in the Y-basis produces the toric code state on the blue (bond) sites up to known single-site spin flips. Corrections to topological order are given by longer-range van der Waals interactions connecting the red (vertex) and blue (bond) sublattices; the leading correction has a pre-factor 1/√{square root over (7)} 6≈ 0.003. Panel (b) illustrates an example of putting atoms on the centers (red, A) and vertices (blue, B) of the hexagons. In this case, the time-evolution and measurement of the red (center, A) sites produces the color code. Subsequently loading atoms on the links (purple, C) can produce D 4  non-Abelian topological order after another round of time-evolution and measuring the blue sites. 
     Consider qubits on the vertices (A sublattice) and bonds (B sublattice) of the honeycomb lattice as shown in  FIG.  2 A . Like in the 1D case, start with a product state|−   ⊗N  and time-evolve evolve with a nearest-neighbor Ising interaction for a time t SPT  to obtain a 2D cluster state which is invariant under the two types of stabilizers shown in  FIG.  2 A . In particular, for a given vertex v, we have Y v     v |ψ(t SPT ) =|(t SPT )  where Y v  is the Pauli-Y on that vertex, and    v =Π bϵv Z b  is a product over the three neighboring bond qubits. Hence, upon measuring the A sublattice in the Y-basis, we obtain a state where    v =±1 depending on the measurement outcome. Moreover, one can also show that for every hexagonal plaquette we have    p =Π bϵp  X b =1. These are exactly the two stabilizers of the toric code. The state has    2  topological order for either    v =±1. If one wishes to use this state for quantum computation purposes, one can simply store these measurement outcomes and use them in post-processing. Alternatively, one can always bring the state to    v =±1 by applying string operators formed by single-site X-rotations. 
     If implement using Rydberg atoms, the time-evolution gets contributions from the longer-range van der Waals interactions. However, like in 1D, the ability to produce topological order via measurement is a robust property of the entire SPT phase of this 2D state. The practical consequence is that the topological phase is stable to all longer-range Ising couplings which do not couple the A and B sublattice. The dominant correction is thus at a distance r=√{square root over (7)}a, shown in  FIG.  2 A  and  FIG.  2 B . This means that the post-measurement state exhibits topological order only over a finite number of atoms on the B sublattice, which we compute (similar to Equation 5 and Equation 6) to be 
     
       
         
           
             
               
                 
                   
                     
                       3 
                       2 
                     
                     ⨯ 
                     
                       1 
                       6 
                     
                     ⨯ 
                     
                       1 
                       
                         
                           ❘ 
                           &#34;\[LeftBracketingBar]&#34; 
                         
                         
                           ln 
                           ⁢ 
                           
                             cos 
                             ⁡ 
                             ( 
                             
                               π 
                               / 
                               
                                 ( 
                                 
                                   2 
                                   × 
                                   
                                     
                                       7 
                                     
                                     6 
                                   
                                 
                                 ) 
                               
                             
                             ) 
                           
                         
                         
                           ❘ 
                           &#34;\[RightBracketingBar]&#34; 
                         
                       
                     
                   
                   ≈ 
                   
                     
                       2 
                       ⨯ 
                       
                         10 
                         4 
                       
                     
                     . 
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                      
                   10 
                 
               
             
           
         
       
     
     Here 3/2 counts the number of B sites per A site, and 6 represents the number of couplings at distance r=√{square root over (7)}a. To exemplify the importance of choosing an optimal lattice: for the square and triangular lattices we find    2  topological order for roughly 3000 and 300 atoms, respectively. 
     In addition to achieving the topological phase in a macroscopic system, a more refined question is how close the described protocol is to the toric code fixed-point limit. Unlike the above discussion, this depends on the strength of the AA and BB couplings. For concreteness, consider dual-species arrays, where these intra-species couplings can be made negligible (alternatively, they can be suppressed using local addressing). The expectation value of the stabilizers is then (see also  FIGS.  2 A , B): 
     
       
         
           
             
               
                 ❘ 
                 &#34;\[LeftBracketingBar]&#34; 
               
               
                 〈 
                 
                   𝒜 
                   v 
                 
                 〉 
               
               
                 ❘ 
                 &#34;\[RightBracketingBar]&#34; 
               
             
             ≈ 
             
               
                 
                   cos 
                   6 
                 
                 ( 
                 
                   π 
                   
                     2 
                     ⁢ 
                     
                       
                         7 
                       
                       6 
                     
                   
                 
                 ) 
               
               ⁢ 
                  
               
                 
                   cos 
                   6 
                 
                 ( 
                 
                   π 
                   
                     2 
                     ⁢ 
                     
                       
                         
                           1 
                           ⁢ 
                           3 
                         
                       
                       6 
                     
                   
                 
                 ) 
               
               ⁢ 
                  
               … 
             
             ≈ 
             0.99993 
           
         
       
       
         
           
             
               〈 
               
                 ℬ 
                 p 
               
               〉 
             
             ≈ 
             
               
                 
                   cos 
                   
                     1 
                     ⁢ 
                     2 
                   
                 
                 ( 
                 
                   π 
                   
                     2 
                     ⁢ 
                     
                       
                         7 
                       
                       6 
                     
                   
                 
                 ) 
               
               ⁢ 
                  
               
                 
                   cos 
                   
                     1 
                     ⁢ 
                     2 
                   
                 
                 ( 
                 
                   π 
                   
                     2 
                     ⁢ 
                     
                       
                         
                           1 
                           ⁢ 
                           3 
                         
                       
                       6 
                     
                   
                 
                 ) 
               
               ⁢ 
                  
               … 
             
             ≈ 
             
               0.99987 
               . 
             
           
         
       
     
     This implies a tonic code fidelity  ≳0.9999, making the state useful for quantum information purposes. Indeed, it is well-known that terminating the lattice with an alternation of rough and smooth boundaries gives a surface code. Applying local X- and Z-pulses in string patterns creates and moves e- and m-anyons, acting as logical operators on the topological qubit. Active error correction requires repeatedly measuring the    v  stabilizers, which is achieved by reloading the A sublattice and repeating the preparation procedure. Similarly, the    p  stabilizers can be measured by performing a X→Z rotation and loading atoms at the center of the plaquettes (for which a square lattice is most natural). 
     This procedure can be generalized to other topological orders. For instance, if the A sublattice lives at the centers of plaquettes of the honeycomb lattice, and B are the vertices (see  FIG.  2 B ), then repeating the above procedure produces the color code on the B sublattice, whose stabilizers are Π vϵp  X v  and Π vϵp  Z v  for every plaquette. For this state—which has significantly more quantum computing power than the toric code—a fidelity per site is  ≳0.997. The 3D toric code can also be prepared -which exhibits more thermal stability—with  ≳0.9998 on the diamond lattice. 
     Fracton Order 
     Referring now to  FIG.  3 A  and  FIG.  3 B . Panel (a) depicts Rydberg atoms placed at the red (“2”) and blue (“1”) sites on the hexagonal prism lattice. After time-evolving under the Rydberg interaction and measuring the red sites, the ground state of the so-called Sierpinski prism model is obtained. Panel (b) depicts quasiparticles that have restricted mobility in the hexagonal layers due to being created at the corners of fractal (Sierpinski) operators. 
     To achieve fracton order—a very exotic version of topological order since it cannot be described by a conventional quantum field theory—we consider qubits on the 3D hexagonal prism lattice shown in  FIG.  3 A , with the red and blue dots defining the A and B sublattices, respectively. As above, time-evolving a product state under a nearest-neighbor Ising interaction prepares the cluster state on this lattice. Upon measuring the red sublattice, we obtain a state which can be equated with the fractal spin liquid introduced by Yoshida. This exhibits fracton order. In particular, while particles can move in the z-direction by the application of a string operator (as in the tonic code), in the hexagonal planes they live at the endpoints of a fractal operator (in particular, the Sierpinski triangle; see  FIG.  3 B ), restricting their mobility. 
     Unlike the previously discussed cases of the GHZ state or topological order, cluster SPT phases giving rise to fracton order are not protected by a global symmetry on the A sublattice, but instead by a subsystem or fractal one. A practical consequence is that AA van der Waals couplings will affect the stability of the fractal spin liquid. First consider the dual species case where these contributions can be neglected. We see in  FIGS.  3 A and  3 B  that the first correction, coupling A to B, occurs at r=2a. From this, the fracton phase is found to persists over approximately 400 atoms. However, by interspersing the time-evolution with four sublattice pulses, this number increases to 1300, which can be systematically improved upon using local addressing. The corresponding expectation value of the stabilizers is 0.997, giving a fidelity  ≳0.998. This establishes a high-fidelity implementation of a fracton state. Moreover, the aforementioned procedure has the added benefit of cancelling nearest-neighbor AA couplings, meaning this could also be successfully realized in single-species arrays if one can measure the A sublattice. 
     Hereinbelow, we discuss a different lattice where measurements lead to an approximate realization of the paradigmatic X-cube model which can be realized with stabilizer expectation values 0.94 in the absence of local addressing. 
     Non-Abelian Defects and Anyons 
     Referring now to  FIG.  4 A  and  FIG.  4 B . Panel (a) depicts time-evolving qutrits (formed by pairs of Rydberg-blockaded atoms) on the vertices and bonds of this 3D configuration and measuring the A sublattice creates the    3  toric code with a non-Abelian defect in the corner. This effectively realizes a square lattice with a disclination; an e-anyon traveling around it transforms into its conjugate ē. Panel (b) illustrates that, starting with the    3  toric code on the B sublattice, Rydberg atoms are loaded on the C and D sublattices. An appropriate sequence of time-evolution and pulses, followed by measuring the C sublattice, produces non-Abelian S 3  topological order. 
     The long-range entangled states we have obtained thus far are all Abelian, i.e., braiding their excitations can at most lead to phase factors. In this last section, it is shown how non-Abelian topological order can arise by using multiple time-evolution and measurement steps. 
     The first type of non-Abelian topological order discuss herein is associated to the group D 4  (i.e., the symmetry group of the square). This can be obtained by combining the toric and color code constructions discussed above. Consider the honeycomb lattice in  FIG.  2 B , with the A sublattice (red) at the center of plaquettes, the B sublattice (blue) at the vertices, and the C sublattice (purple) on the bonds. A strikingly simple representation of a state on the C sublattice with non-Abelian D 4  topological order is found through the following sequence of controlled-Z and single-site rotations: 
     
       
         
           
             
               
                 
                   〈 
                   + 
                   
                     
                       
                         ❘ 
                         &#34;\[LeftBracketingBar]&#34; 
                       
                         
                     
                     AB 
                   
                   ⁢ 
                   
                     
                       ∏ 
                       
                         〈 
                         
                           b 
                           , 
                           c 
                         
                         〉 
                       
                     
                     
                       
                         CZ 
                         
                           b 
                           , 
                           c 
                         
                       
                       ⁢ 
                       
                         
                           ∏ 
                           
                             b 
                             ∈ 
                             B 
                           
                         
                         
                           
                             e 
                             
                               
                                 - 
                                 i 
                               
                               ⁢ 
                               
                                 π 
                                 8 
                               
                               ⁢ 
                               
                                 Y 
                                 b 
                               
                             
                           
                           ⁢ 
                           
                             
                               ∏ 
                               
                                 〈 
                                 
                                   a 
                                   , 
                                   b 
                                 
                                 〉 
                               
                             
                             
                               
                                 CZ 
                                 
                                   a 
                                   , 
                                   b 
                                 
                               
                               ⁢ 
                               
                                 
                                   ❘ 
                                   &#34;\[LeftBracketingBar]&#34; 
                                 
                                 
                                   
                                     + 
                                     
                                       〉 
                                       ABC 
                                     
                                   
                                   , 
                                 
                               
                             
                           
                         
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   11 
                 
               
             
           
         
       
     
     where  b, c  denotes nearest neighbors connecting the B and C sublattices. While Equation 11 is written in terms of post-selection for simplicity, this is not required since measurement outcomes can be corrected. To create this state in Rydberg atom arrays, load the A and B sublattice into a product state, with controlled-Z being naturally generated by time-evolution. As discussed above, measuring the A sublattice produces the color code, whose measurement-induced Abelian anyons can be removed with single-site rotations. Now apply the Y-rotation in Equation 11 and load the C sublattice such that time-evolution induces the second controlled-Z. Finally, measuring the B sublattice can again produce Abelian anyons, removable by single-site X-rotations. 
     Interestingly, if one had measured the A sublattice at the end (or if one had not removed the anyons in the intermediate color code), a density of non-Abelian anyons could have been produced in the resulting D 4  topological order, which are much more difficult to efficiently remove. In this regard, the multi-step process is important, even when implementing Equation 11 in a quantum computer rather than through time-evolution in a Rydberg atom array. 
     The intermediate state will now be the    3  toric code, rather than the color code. For this, qutrits instead of qubits are used. One straightforward option is to consider pairs of qubits where one uses only three states, which we label|0 , |1 , |2 (i.e., never initializing or evolving into the fourth state). In the context of Rydberg atoms, this is naturally achieved by having two atoms in close proximity such that the Rydberg blockade forbids them from being simultaneously excited (although using multiple energy levels in a single atom offers an alternative route to qutrit encoding). Similar to the qubit case, time-evolution can generate the controlled-Z gate for qutrits, defined as 
     
       
         
           
             
               
                 
                   
                     
                       
                         
                           
                             
                               
                                 
                                   
                                     
                                       CZ 
                                       ⁢ 
                                       
                                         
                                           ❘ 
                                           &#34;\[LeftBracketingBar]&#34; 
                                         
                                         i 
                                       
                                     
                                     〉 
                                   
                                   ⊗ 
                                 
                                 
                                   ❘ 
                                   &#34;\[RightBracketingBar]&#34; 
                                 
                               
                               ⁢ 
                               j 
                             
                             〉 
                           
                           = 
                           
                             
                               ω 
                               
                                 i 
                                 ⁢ 
                                 j 
                               
                             
                             ⁢ 
                             
                               
                                 ❘ 
                                 &#34;\[LeftBracketingBar]&#34; 
                               
                               i 
                             
                           
                         
                         〉 
                       
                       ⊗ 
                     
                     
                       ❘ 
                       &#34;\[RightBracketingBar]&#34; 
                     
                   
                   ⁢ 
                   j 
                 
                 〉 
               
               ⁢ 
                   
               with 
               ⁢ 
               
                   
                    
               
               ⁢ 
               ω 
             
             = 
             
               
                 e 
                 
                   
                     2 
                     ⁢ 
                     π 
                     ⁢ 
                     i 
                   
                   3 
                 
               
               . 
             
           
         
       
     
     Hence, placing our qutrits on the lattice in  FIG.  2 A  can generate the    3  toric code on the honeycomb lattice. However, it will be convenient to instead consider the Lieb lattice (i.e., vertices (A) and bonds (B) of the square lattice), thereby producing this topological state on the B sublattice. 
     As an aside, although this topological order is Abelian, it can be used to trap non-Abelian defects. For instance, consider the 3D configuration in  FIG.  4 A , which can be interpreted as a 2D Lieb lattice with a disclination. Repeating the above preparation protocol produces the    3  toric code with such a disclination defect. If an e- or m-anyon moves around the corner in  FIG.  4 A , it will transform into the conjugate anyon (ē or  m ), leading to a quantum dimension d=3 for the defect. 
     Finally, to transform this state into S 3  topological order, consider the C and D sublattices defined in  FIG.  4 B , which are loaded with effective qubits. Through time-evolution and sublattice rotations, a controlled-charge-conjugation can first be implemented for the nearest neighbors connecting the B and C sublattices. Subsequently, controlled-Z can be implemented for the nearest neighbors connecting the C and D sublattices. Remarkably, measuring the C qubits produces the ground state of the S 3  quantum double model on the remaining qutrits on B and qubits on. 
     At a conceptual level, this construction can be interpreted as effectively gauging the    2  charge-conjugation symmetry of the    3  toric code, in line with S 3 =   3       2 . 
     Outlook 
     In summary, a variety of highly sought after long-range entangled states can be efficiently created from simple ingredients: initialization of qubits, unitary time-evolution under Ising interactions, and measurement of a subset of qubits, which can be realized on existing experimental platforms. An implementation in Rydberg atom arrays was described, by identifying lattice structures that leverage their intrinsic interactions and demonstrate high fidelity preparation of 1D GHZ and 2D tonic and color code states, as well as 3D fracton phases. Several possible generalizations of the disclosed protocol have been discussed. For instance, using the Rydberg blockade one can realize long-range entangled qutrits; using 3D configurations one can obtain 2D lattices with exotic defects; and by appealing to multiple measurement steps, one can even realize states hosting non-Abelian anyons. 
     Given that the protocol described herein, it is a worthwhile endeavour to further explore how the resulting states could be used. The 1D GHZ state is well-known to be useful in quantum metrology due to its nonclassical sensitivity to external fields. Creation of surface codes for the toric and color codes has been described herein, even allowing for active error correction. Similarly, the 3D fracton states described herein could be similarly used to perform quantum computation, although its surface code properties have not yet been extensively explored. 
     A particularly exciting prospect is the realization of non-Abelian topological order in engineered quantum systems. Here it is demonstrated how both S 3  and D 4  topological order can be realized in Rydberg atom arrays. The novel schemes proposed herein could be used to realize D 4  topological order in a chip of superconducting qubits using the simple geometry of  FIG.  2 B . Note that Equation 11 requires a circuit which is only ten layers deep, which is achievable with current coherence times. 
     Methods 
     Written out in spin notation, the Hamiltonian of interest is: 
     
       
         
           
             
               
                 
                   H 
                   = 
                   
                     
                       
                         ∑ 
                         i 
                       
                       
                         
                           
                             Ω 
                             i 
                           
                           2 
                         
                         ⁢ 
                         
                           X 
                           i 
                         
                       
                     
                     - 
                     
                       
                         ∑ 
                         i 
                       
                       
                         
                           
                             h 
                             i 
                           
                           2 
                         
                         ⁢ 
                         
                           Z 
                           i 
                         
                       
                     
                     + 
                     
                       
                         1 
                         8 
                       
                       ⁢ 
                       
                         
                           ∑ 
                           
                             i 
                             ≠ 
                             j 
                           
                         
                         
                           
                             V 
                             ij 
                           
                           ⁢ 
                           
                             Z 
                             i 
                           
                           ⁢ 
                           
                             Z 
                             j 
                           
                         
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   12 
                 
               
             
           
         
       
     
     with 
     
       
         
           
             
               h 
               i 
             
             = 
             
               
                 δ 
                 i 
               
               + 
               
                 
                   1 
                   2 
                 
                 ⁢ 
                 
                   
                     ∑ 
                     
                       j 
                       ≠ 
                       i 
                     
                   
                   
                     V 
                     
                       i 
                       ⁢ 
                       j 
                     
                   
                 
               
             
           
         
       
     
     (compared to Equation 1). The extra shift of the longitudinal field—explicitly giving an Ising interaction—is also practically significant, since it is important for the stability of the resulting state. 
     Dual Species and Förster Reasonances 
     One interesting option is to use different Alkali atoms for the two sublattices of the set-ups described herein. In addition to making it possible to measure a single sublattice, it can be used to suppress the intra-species interaction. In particular, it is known that by using (near-) Förster resonances, one can choose the Rydberg levels of each atom such that the intra-species interaction is negligible. For instance, the inter-species interaction between the 48S 1/2  state in Rubidium and the 51S 1/2  state in Cesium is roughly 100 times larger than the intra-species coupling, which could potentially be further optimized by exploring higher Rydberg states or other atoms. 
     The tunability of the dual species interaction derives from finding Rydberg states such that there are nearly-degenerate states (e.g., for a given choice of S-states in Rubidium and Cesium, the energy might be approximately unchanged if the Rubidium (Cesium) atom in an S-state transitions to a lower (higher) energy P-state). It is desirable to place the atoms far apart enough such that this transition is considerably off-resonant; beyond the so-called ‘crossover distance’ R c , the interaction is effectively van der Waals (for distances shorter than R c , the interaction decays as ˜1/r 3  and additional spin-flop terms would appear in the Hamiltonian, which should be avoided). For instance, for the aforementioned example of the Rb 48S 1/2  and Cs 51S 1/2  states, R c ≈7.5 μm. For a spacing a=10 μm, the inter-species interaction is still at the MHz scale (with the intra-species coupling being negligible). The proposed preparation time is 
     
       
         
           
             π 
             
               V 
               ⁡ 
               ( 
               a 
               ) 
             
           
         
       
     
     which is thus on me order of microseconds—fast enough to neglect decoherence effects which arise at ˜100 μs. The preparation time can be reduced by considering higher principal quantum numbers. 
     Lower Bound on Stabilizer Fidelity Pre-Measurement 
     It has been shown that the quantum state|ψ  on N qubits has a fidelity  ≥1−Nε/2 for a stabilizer state corresponding to stabilizer operators K if  ψ|K|ψ =1−ε. If Nε&lt;&lt;1, this implies that the fidelity density satisfies the lower bound  ≥1−ε/2. Since the fidelity density will converge a finite value as N→∞(at least for a state with a finite correlation length), it is expected that the condition Nε&lt;&lt;1 does not need to be explicitly enforced (although this is of course the region where the bound is the most useful). Indeed, in one spatial dimension one can straightforwardly prove using matrix product state techniques that the bound  ≥ε/2 is always satisfied for  ψ|K|ψ =1−ε, even in the thermodynamic limit. 
     Lower Bound on Stabilizer Fidelity Post-Measurement 
     The above gives a lower bound on the fidelity for the cluster state, based on the expectation value of the stabilizer. Of ultimate interest is the fidelity with respect to the long-range entangled state for the post-measurement state. Fortunately, these two can be equated for the following reason. Firstly, one can show that  X n   =0, making both measurement outcomes equally likely. Secondly, the correlations between X-measurements on different sites is completely negligible (e.g., in the 1D case we find that  X n X n+2   &lt;10 −10  even upon including all longer-range couplings in the single-species set-up). Hence, to a very good degree, the post-measurement state is given by: 
     
       
         
           
             
               
                 
                   
                     
                       
                         | 
                         
                           ψ 
                           
                             o 
                             ⁢ 
                             u 
                             ⁢ 
                             t 
                           
                         
                       
                       〉 
                     
                     = 
                     
                       
                         ∏ 
                         
                           a 
                           ∈ 
                           A 
                         
                       
                       
                         
                           ( 
                           
                             
                               
                                 X 
                                 a 
                               
                               + 
                               
                                 s 
                                 a 
                               
                             
                             
                               2 
                             
                           
                           ) 
                         
                            
                         | 
                         ψ 
                         〉 
                       
                     
                   
                   , 
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   13 
                 
               
             
           
         
       
     
     where s a =±1 denotes the measurement outcome. From the above explicit formula, it is straightforward to confirm our claim about the equality of fidelities. 
     Numerical Simulation of Cluster State Preparation 
     Consider  87 Rb in the Rydberg level 70S 1/2 . The intrinsic van der Waals interaction is H int =UΣ i&lt;j n i n j /|r i −r j | 6  with U≈5THz μm 6 . If a is the lattice spacing, then the nearest-neighbor interaction is V(a)=U/a 6 . For a=10 μm (a=12 μm) we have V(a)≈5 MHz (V(a)≈1.7 MHz). All atoms are initialized in their ground state. An X-pulse is applied for a time t=t pulse  with strength 
     
       
         
           
             Ω 
             = 
             
               
                 π 
                 
                   2 
                   ⁢ 
                   
                     t 
                     pulse 
                   
                 
               
               . 
             
           
         
       
     
     Taking t pulse =20 ns, corresponding to Ω≈78 MHz, a Z-pulse is then applied for the same time with 
     
       
         
           
             h 
             = 
             
               - 
               
                 
                   π 
                   
                     2 
                     ⁢ 
                     
                       t 
                       pulse 
                     
                   
                 
                 . 
               
             
           
         
       
     
     (Note that one can avoid the need for this separate pulse: since Z commutes with the Rydberg interaction, one can perform both at the same time.) After these pulses, the state|−   ⊗N  is realized (up to corrections due to the interactions which are present during the short pulse; the simulations include these interactions). Next, time-evolution is performed under the natural Ising interaction (Ω=h=0) for a time 
     
       
         
           
             
               t 
               SPT 
             
             = 
             
               
                 π 
                 
                   V 
                   ⁡ 
                   ( 
                   a 
                   ) 
                 
               
               - 
               
                 
                   5 
                   2 
                 
                 ⁢ 
                 
                   
                     t 
                     pulse 
                   
                   . 
                 
               
             
           
         
       
     
     After this time-evolution, a final X-pulse is applied of strength 
     
       
         
           
             Ω 
             = 
             
               - 
               
                 
                   
                     π 
                     
                       2 
                       ⁢ 
                       
                         t 
                         pulse 
                       
                     
                   
                   ⨯ 
                   
                     1 
                     
                       2 
                       6 
                     
                   
                 
                 . 
               
             
           
         
       
     
     The simulations are performed using the density matrix renormalization group method for infinite systems, thereby directly working in the thermodynamic limit. It was found that small bond dimensions,  =20−40, were already sufficient to guarantee convergence. (Note: if we the pulse time can be reduced to e.g., t pulse =5 ns, then the fidelity goes up another order of magnitude: we find 0.999997 (for a=12 μm), in close agreement with the prediction in Equation 9). 
     Dual Stabilizer 
     When measuring X b  on a given site of the B sublattice, then any stabilizer that does not commute with it is no longer a stabilizer of the post-measurement state. However, sometimes a product of such stabilizers commutes with all X-measurements, thereby providing a constraint on the post-measurement state. E.g., in the 1D case, the product of X 2n Z 2n+1 X 2n+2  (each individual factor not commuting with the X-measurements) gives the constraint that the resulting post-measurement state has the    2  symmetry Π n  Z 2n+1 . In the 2D toric code, a local constraint is obtained instead: the product of the cluster term on four bonds around a plaquette gives the condition    p =+1. 
     The Effect of Symmetry Preserving Couplings 
     For the case of the GHZ and toric code states, it is shown that AA and BB Ising couplings do not affect the resulting phase of matter. The argument for BB is simple to see: these gates commute with the measurement and hence dress the post-measurement stat with a finite-depth circuit (which cannot destroy the long-range entangled states). The AA couplings are more interesting and subtle: in fact, they also push through to just being effective finite-depth circuits on the post-measurement states, which one can argue by using the stabilizer property of the cluster state, which allows one to replace it by X (or a product thereof) on the B sublattice. This is derived in explicit detail for the 1D case in the Supplemental Material. In conclusion, only the longer-range AB couplings affect the quality of the resulting phase of matter. 
     Improvements Using Local Addressing 
     The protocol described herein can be further improved using spatially-dependent fields, called local addressing. This tool has already been used in large Rydberg arrays. 
     To illustrate how this can suppress or even cancel unwanted van der Waals interactions, consider again the 1D case. Using spatially-dependent X-pulses, qubits can be flipped on a sublattice of arbitrary choice. Let    k  correspond to flipping every fourth site starting from k=1,2,3,4. If U(t) denotes the time-evolution, then consider the following alternation of time-evolution and sublattice pulses sketched in  FIG.  1 D : 
           2     3   U ( t   SPT /2)   1     3   U ( t   SPT /2)   1     2   U ( t   SPT ).  Equation 14
 
     The additional procedure does not affect the net time-evolution for any sites separated by an odd multiple of lattice spacings; in particular, the nearest-neighbor evolution still prepares the cluster state. However, the second-nearest neighbor coupling drops out completely, making r=3a the dominant correction. (The fourth-nearest neighbor experiences an effective doubled time-evolution.) Similarly, local addressing can suppress unwanted contributions for the fracton model. 
     This cancellation procedure can also be useful after the state preparation. Indeed, the remaining B sublattice will continue to time-evolve under the long-range Ising interaction. By construction, this coupling is small, especially in the dual-species set-up. Nevertheless, such effects can accumulate. One can then apply local pulses as above to effectively change the sign of a problematic coupling, such that further time-evolution will undo its effects. In the 1D case this is not necessary, since the GHZ state is an eigenstate of the van der Waals interactions. 
     Fracton Order 
     Upon measuring the red sublattice of the cluster state on the 3D lattice in  FIG.  3 A , we obtain the stabilizers shown in  FIG.  5 A  and  FIG.  5 B . 
     These can be equated with the stabilizers of the fractal spin liquid introduced in B. Yoshida, Exotic topological order in fractal spin liquids, Phys. Rev. B 88, 125122 (2013). 
     The following expectation value can be obtained for the stabilizers in the dual-species set-up: 
     
       
         
           
             
               
                 
                   
                     〈 
                     
                       𝒜 
                       v 
                     
                     〉 
                   
                   ≈ 
                   
                     
                       
                         cos 
                         
                           1 
                           ⁢ 
                           5 
                         
                       
                       ( 
                       
                         π 
                         
                           2 
                           · 
                           
                             2 
                             6 
                           
                         
                       
                       ) 
                     
                     ⁢ 
                     
                       
                         cos 
                         6 
                       
                       ( 
                       
                         π 
                         
                           2 
                           · 
                           
                             
                               5 
                             
                             6 
                           
                         
                       
                       ) 
                     
                     ⁢ 
                        
                     … 
                   
                   ≈ 
                   
                     0 
                     .995 
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   15 
                 
               
             
           
         
       
       
         
           
             
               
                 
                   
                     〈 
                     
                       ℬ 
                       p 
                     
                     〉 
                   
                   ≈ 
                   
                     
                       
                         cos 
                         
                           4 
                           ⁢ 
                           2 
                         
                       
                       ( 
                       
                         π 
                         
                           2 
                           · 
                           
                             2 
                             6 
                           
                         
                       
                       ) 
                     
                     ⁢ 
                     
                       
                         cos 
                         
                           1 
                           ⁢ 
                           2 
                         
                       
                       ( 
                       
                         π 
                         
                           2 
                           · 
                           
                             
                               5 
                             
                             6 
                           
                         
                       
                       ) 
                     
                     ⁢ 
                        
                     … 
                   
                   ≈ 
                   
                     
                       0 
                       . 
                       9 
                     
                     ⁢ 
                     8 
                     ⁢ 
                     6 
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   16 
                 
               
             
           
         
       
     
     This result can be improved by using local addressing: the out-of-plane contributions at r=2a can be cancelled by noticing that the red sublattice is tripartite (forming a triangular lattice in each layer). Correspondingly, split A up into A 1 ∪ A 2 ∪ A 3 , and let    A     i    denote flipping qubits on the A 1  sublattice via an X-pulse. Consider 
           A     3     U ( t )   A     3       A     2     U ( t )   A     2       A     1     U ( t )   A1   U ( t )  Equation 17
 
     Setting t=t SPT /2 will create the cluster state, and effectively remove the aforementioned couplings. The following improved results have been obtained: 
     
       
         
           
             
               
                 
                   
                     〈 
                     
                       𝒜 
                       v 
                     
                     〉 
                   
                   ≈ 
                   
                     
                       
                         cos 
                         3 
                       
                       ( 
                       
                         π 
                         
                           2 
                           · 
                           
                             2 
                             6 
                           
                         
                       
                       ) 
                     
                     ⁢ 
                     
                       
                         cos 
                         6 
                       
                       ( 
                       
                         π 
                         
                           2 
                           · 
                           
                             
                               5 
                             
                             6 
                           
                         
                       
                       ) 
                     
                     ⁢ 
                        
                     … 
                   
                   ≈ 
                   
                     0 
                     .9985 
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   18 
                 
               
             
           
         
       
       
         
           
             
               
                 
                   
                     〈 
                     
                       ℬ 
                       p 
                     
                     〉 
                   
                   ≈ 
                   
                     
                       
                         cos 
                         6 
                       
                       ( 
                       
                         π 
                         
                           2 
                           · 
                           
                             2 
                             6 
                           
                         
                       
                       ) 
                     
                     ⁢ 
                     
                       
                         cos 
                         
                           1 
                           ⁢ 
                           2 
                         
                       
                       ( 
                       
                         π 
                         
                           2 
                           · 
                           
                             
                               5 
                             
                             6 
                           
                         
                       
                       ) 
                     
                     ⁢ 
                        
                     … 
                   
                   ≈ 
                   
                     
                       0 
                       . 
                       9 
                     
                     ⁢ 
                     9 
                     ⁢ 
                     7 
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   19 
                 
               
             
           
         
       
     
     D 4  Topological Order 
     To prove that Equation 11 exhibits D 4  non-Abelian topological order, first consider 
     
       
         
           
             
               
                 
                   〈 
                   
                     + 
                     
                       A 
                     
                   
                   
                     
                       ❘ 
                       &#34;\[LeftBracketingBar]&#34; 
                     
                          
                     
                       
                         ∏ 
                         
                           〈 
                           
                             a 
                             , 
                             b 
                           
                           〉 
                         
                       
                       
                         C 
                         ⁢ 
                            
                         
                           Z 
                           
                             a 
                             , 
                             b 
                           
                         
                       
                     
                     
                       ❘ 
                       &#34;\[RightBracketingBar]&#34; 
                     
                   
                   + 
                   
                     
                       〉 
                       
                         A 
                         , 
                         B 
                       
                     
                     . 
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   20 
                 
               
             
           
         
       
     
     As discussed above, this is the color code on the B sublattice. Next, use the fact that gauging the    2  Hadamard symmetry of the color code produces D 4  topological order. Usually, gauging is a nonlocal transformation which would require a linear-depth circuit. However, it is known that an arbitrary initial state can be gauged by a finite-depth circuit (made of controlled-Z gates) and single-site projections or measurements, and this general fact can be used to produce D 4  topological order. In particular, first perform a Y basis rotation to transform the Hadamard symmetry into Π bϵb  X b  after which entangle the B degrees of freedom with a product state on the C sublattice through controlled-Z: 
     
       
         
           
             
               
                 
                   
                     
                       
                         
                           
                             
                               
                                 ∏ 
                                 
                                   
                                     ( 
                                     
                                       b 
                                       , 
                                       c 
                                     
                                   
                                   〉 
                                 
                               
                               
                                 C 
                                 ⁢ 
                                    
                                 
                                   Z 
                                   
                                     b 
                                     , 
                                     c 
                                   
                                 
                                 ⁢ 
                                 
                                   
                                     ∏ 
                                     
                                       b 
                                       ∈ 
                                       B 
                                     
                                   
                                   
                                     
                                       e 
                                       
                                         
                                           - 
                                           i 
                                         
                                         ⁢ 
                                         
                                           π 
                                           8 
                                         
                                         ⁢ 
                                         
                                           Y 
                                           b 
                                         
                                       
                                     
                                     〈 
                                   
                                 
                               
                             
                             + 
                           
                           
                             ❘ 
                             &#34;\[RightBracketingBar]&#34; 
                           
                         
                         
                           A 
                         
                       
                       ⁢ 
                       
                         
                           ∏ 
                           
                             
                               ( 
                               
                                 a 
                                 , 
                                 b 
                               
                             
                             〉 
                           
                         
                         
                           C 
                           ⁢ 
                           
                             Z 
                             
                               a 
                               , 
                               b 
                             
                           
                         
                       
                     
                     
                       ❘ 
                       &#34;\[RightBracketingBar]&#34; 
                     
                   
                   + 
                   
                     
                       〉 
                       
                         A 
                         , 
                         B 
                       
                     
                     . 
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   21 
                 
               
             
           
         
       
     
     Now projecting the B sites into  +1 will effectively implement the gauging of the original color code, thereby producing D 4  topological order. 
         3  Toric Code, its Defects, and S 3  Topological Order 
     To create the    3  toric code, place the qutrits on the vertices and bonds of the square lattice, forming the A and B sublattices respectively. One can show that combining the natural Rydberg interaction with several sublattice pulses gives the following effective nearest-neighbor interaction: 
         H   a,b   =U ( n   1,a   n   1,b   +n   2,a   n   2,b   −n   1,a   n   2,b   −n   2,a   n   1,b ),  Equation 22
 
     where n i,λ  denotes whether atom i on site λ is excited. Hence, the time-evolution 
     
       
         
           
             e 
             
               
                 - 
                 
                   
                     4 
                     ⁢ 
                     π 
                     ⁢ 
                     i 
                   
                   
                     3 
                     ⁢ 
                     U 
                   
                 
               
               ⁢ 
               
                 H 
                 
                   a 
                   , 
                   b 
                 
               
             
           
         
       
     
     implements the controlled-Z gate defined in the main text, producing the topological state. Note that the first correction connecting different sublattices has a prefactor 
     
       
         
           
             
               
                 
                   ( 
                   
                     
                       r 
                       
                         A 
                         ⁢ 
                         B 
                       
                     
                     
                       r 
                       
                         AB 
                         ′ 
                       
                     
                   
                   ) 
                 
                 6 
               
               = 
               
                 
                   1 
                   
                     
                       5 
                     
                     6 
                   
                 
                 ≈ 
                 
                   
                     0 
                     . 
                     0 
                   
                   ⁢ 
                   0 
                   ⁢ 
                   8 
                 
               
             
             , 
           
         
       
     
     justifying a nearest-neighbor approximation. 
     If the lattice being used has a disclination, it can trap a non-Abelian defect. Indeed, the usual    3  cluster state is prepared by having a pattern of CZ and (CZ) †  being unitarily equivalent to using only CZ if we act with charge-conjugation on, for example, every other A site along with the B sites immediately above and to the right. However, a disclination destroys the bipartiteness of A, leading to a charge-conjugation defect. This kind of defect has a quantum dimension d=3. 
     Described in the sections below are the details of how this    3  toric code can be transformed into non-Abelian topological order. The key steps are summarized below. First, starting with the    3  toric code on the B sublattice in  FIG.  4 B , load the C sublattice with Rydberg atoms serving as qubits. Their natural interaction is: 
     
       
         
           
             
               
                 
                   
                     H 
                     
                       B 
                       ⁢ 
                       C 
                     
                   
                   = 
                   
                     
                       U 
                       
                         B 
                         ⁢ 
                         C 
                       
                     
                     ⁢ 
                     
                       
                         ∑ 
                         
                           〈 
                           
                             b 
                             , 
                             c 
                           
                           〉 
                         
                       
                       
                         
                           ( 
                           
                             
                               
                                 n 
                                 
                                   1 
                                   , 
                                   b 
                                 
                               
                               ⁢ 
                               
                                 n 
                                 c 
                               
                             
                             + 
                             
                               
                                 n 
                                 
                                   2 
                                   , 
                                   b 
                                 
                               
                               ⁢ 
                               
                                 n 
                                 c 
                               
                             
                           
                           ) 
                         
                         . 
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   23 
                 
               
             
           
         
       
     
     Combining time-evolution for a time 
     
       
         
           
             t 
             = 
             
               π 
               
                 U 
                 
                   B 
                   ⁢ 
                   C 
                 
               
             
           
         
       
     
     with an appropriate pulse effectively implements controlled-charge-conjugation for nearest-neighbors. Then load the D sublattice with qubits, and let C and D interact under 
     
       
         
           
             
               
                 
                   
                     H 
                     CD 
                   
                   = 
                   
                     
                       U 
                       CD 
                     
                     ⁢ 
                     
                       
                         ∑ 
                         
                           
                             ( 
                             
                               c 
                               , 
                               d 
                             
                           
                           〉 
                         
                       
                       
                         
                           n 
                           c 
                         
                         ⁢ 
                         
                           n 
                           d 
                         
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   24 
                 
               
             
           
         
       
     
     for time 
     
       
         
           
             
               t 
               = 
               
                 π 
                 
                   U 
                   
                     C 
                     ⁢ 
                     D 
                   
                 
               
             
             , 
           
         
       
     
     combined with a pulse that cancels out coupling to B sublattice, implementing a controlled-Z. Finally, measuring the C qubits produces the ground state of the S 3  quantum double model on the remaining qutrits on B and qubits on D. 
     Supplemental Materials 
     Preparation of Cat State with Timing Imperfections 
     Consider the state|ψ(t) =e −iHt |−) ⊗N  obtained by evolving the product state with the Ising Hamiltonian H=JΣZ n Z n+1 . (For the special time 
     
       
         
           
             
               
                 t 
                 
                   S 
                   ⁢ 
                   P 
                   ⁢ 
                   T 
                 
               
               = 
               
                 π 
                 
                   4 
                   ⁢ 
                   
                     
                       ❘ 
                       &#34;\[LeftBracketingBar]&#34; 
                     
                     J 
                     
                       ❘ 
                       &#34;\[RightBracketingBar]&#34; 
                     
                   
                 
               
             
             , 
           
         
       
     
     this coincides with the cluster state.) 
     Theorem. Suppose one measures the even sites of |ψ(t)  in the Pauli-X basis, producing the state|ψ out (t)  for a given measurement outcome. Then the remaining odd sites have a correlation function 
         ψ out ( t )| Z   2m−1   Z   2(m+n)−1 |ψ out ( t ))˜ e   −n/ξ   Equation 25
 
     Define 
     
       
         
           
             
               α 
               = 
               
                 
                   
                     
                       
                         cos 
                         2 
                       
                       ( 
                       
                         2 
                         ⁢ 
                         t 
                       
                       ) 
                     
                     . 
                         
                     If 
                   
                   ⁢ 
                       
                   n 
                 
                 ≫ 
                 
                   
                     1 
                     - 
                     α 
                   
                   
                     1 
                     + 
                     α 
                   
                 
               
             
             , 
           
         
       
     
     then the correlation length ξ is normally distributed (corresponding to the randomness of measurement), ξ˜N( ξ ,σ ξ ), with the following mean and standard deviation: 
     
       
         
           
             
               
                 
                   
                     ξ 
                     _ 
                   
                   = 
                   
                     
                       
                         2 
                         
                           
                             ( 
                             
                               1 
                               + 
                               α 
                             
                             ) 
                           
                           ⁢ 
                              
                           ln 
                           ⁢ 
                              
                           
                             ( 
                             
                               
                                 1 
                                 + 
                                 α 
                               
                               
                                 1 
                                 - 
                                 α 
                               
                             
                             ) 
                           
                         
                       
                       ⁢ 
                           
                       and 
                       ⁢ 
                            
                       
                         σ 
                         ξ 
                       
                     
                     = 
                     
                       
                         
                           ξ 
                           ¯ 
                         
                         
                           n 
                         
                       
                       ⁢ 
                       
                         
                           
                             
                               1 
                               - 
                               α 
                             
                             
                               1 
                               + 
                               α 
                             
                           
                         
                         . 
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   26 
                 
               
             
           
         
       
     
     Before proving this result, let us make a few observations:
         1. For the special case t=0 (i.e., α=1), one recovers  ξ =0 for the product state, whereas the cluster state at t=t SPT  (i.e., α=0) gives rise to a true cat state with  ξ =∞.       

     2. Expanding around t=t SPT +δt (i.e., around α=0), one find 
     
       
         
           
             ξ 
             = 
             
               
                 1 
                 α 
               
               - 
               1 
               + 
               
                 
                   O 
                   ⁡ 
                   ( 
                   α 
                   ) 
                 
                 . 
               
             
           
         
       
     
     The correlation length thus blows up as 
     
       
         
           
             ξ 
             ∼ 
             
               
                 
                   4 
                   ⁢ 
                   
                     t 
                     
                       S 
                       ⁢ 
                       P 
                       ⁢ 
                       T 
                     
                     2 
                   
                 
                 
                   π 
                   2 
                 
               
               ⁢ 
               
                 1 
                 
                   
                     ( 
                     
                       δ 
                       ⁢ 
                       t 
                     
                     ) 
                   
                   2 
                 
               
             
           
         
       
     
     upon approaching t→t SPT .
         3. For all t, we see that σ ξ →0 as n→∞. However, for practical purposes it can be meaningful to set n= ξ  to get a sense of how the correlation length will vary over physically relevant length scales. An example is shown in  FIG.  6    for t=0.99t SPT , where one finds that the deviations around  ξ ≈4000 are small.   4. The plot in  FIG.  6    corresponds to  ξ  with a thickness of one standard deviation where n= ξ  has been chosen.       

     Proof Firstly, since the time-evolution is implemented using a purely-diagonal nearest-neighbor Ising Hamiltonian, the X-correlations are uncorrelated beyond two sites. More precisely,  ψ(t) X 2n |ψ(t) =   ψ(t)|X 2n |ψ(t)  for any set of operators that does not include all even operators. Hence, if one measures X on all even sites (leaving out one for convenience), then the state with outcome X 2n =s 2n  can be written as 
     
       
         
           
             
               
                 
                   
                     
                       
                         
                           ❘ 
                           &#34;\[LeftBracketingBar]&#34; 
                         
                         
                           
                             ψ 
                             
                               o 
                               ⁢ 
                               u 
                               ⁢ 
                               t 
                             
                           
                           ⁢ 
                           
                             ( 
                             t 
                             ) 
                           
                         
                       
                       〉 
                     
                     = 
                     
                       
                         ∏ 
                         n 
                       
                       
                         
                           
                             
                               X 
                               
                                 2 
                                 ⁢ 
                                 n 
                               
                             
                             + 
                             
                               s 
                               
                                 2 
                                 ⁢ 
                                 n 
                               
                             
                           
                           
                             
                               2 
                               - 
                               
                                 2 
                                 ⁢ 
                                 α 
                                 ⁢ 
                                 
                                   s 
                                   
                                     2 
                                     ⁢ 
                                     n 
                                   
                                 
                               
                             
                           
                         
                         ⁢ 
                         
                           
                             ❘ 
                             &#34;\[LeftBracketingBar]&#34; 
                           
                           
                             ψ 
                             ⁡ 
                             ( 
                             t 
                             ) 
                           
                         
                       
                     
                   
                   〉 
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   27 
                 
               
             
           
         
       
       
         
           
             
               where 
               ⁢ 
                   
               α 
             
             := 
             
               
                 
                   ❘ 
                   &#34;\[LeftBracketingBar]&#34; 
                 
                 
                   〈 
                   X 
                   〉 
                 
                 
                   ❘ 
                   &#34;\[RightBracketingBar]&#34; 
                 
               
               = 
               
                 
                   cos 
                   2 
                 
                 ⁢ 
                    
                 
                   
                     ( 
                     
                       2 
                       ⁢ 
                       t 
                     
                     ) 
                   
                   . 
                 
               
             
           
         
       
     
     From this, one can derive that 
     
       
         
           
             
               
                 
                   
                     
                       
                         
                           
                             ❘ 
                             &#34;\[LeftBracketingBar]&#34; 
                           
                           
                             〈 
                             
                               ψ 
                               
                                 o 
                                 ⁢ 
                                 u 
                                 ⁢ 
                                 t 
                               
                             
                           
                           
                             ❘ 
                             &#34;\[RightBracketingBar]&#34; 
                           
                         
                         ⁢ 
                         
                           Z 
                           
                             
                               2 
                               ⁢ 
                               m 
                             
                             - 
                             1 
                           
                         
                         ⁢ 
                         
                           Z 
                           
                             
                               2 
                               ⁢ 
                               n 
                             
                             + 
                             1 
                           
                         
                         ⁢ 
                         
                           
                             ❘ 
                             &#34;\[LeftBracketingBar]&#34; 
                           
                           
                             ψ 
                             
                               o 
                               ⁢ 
                               u 
                               ⁢ 
                               t 
                             
                           
                         
                       
                       〉 
                     
                     
                       ❘ 
                       &#34;\[RightBracketingBar]&#34; 
                     
                   
                   = 
                   
                     
                       ∏ 
                       
                         k 
                         = 
                         m 
                       
                       n 
                     
                     
                       
                         
                           ( 
                           
                             
                               1 
                               - 
                               α 
                             
                             
                               1 
                               + 
                               α 
                             
                           
                           ) 
                         
                         
                           
                             1 
                             - 
                             
                               s 
                               
                                 2 
                                 ⁢ 
                                 k 
                               
                             
                           
                           2 
                         
                       
                       . 
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   28 
                 
               
             
           
         
       
     
     For instance, this follows from the factorization property  Z 2m−1 Z k Z k Z 2n+1   = Z 2m−1 Z k     Z k Z 2n+1    and a direct calculation for the simplest case n=m. 
     Hence, it can be seen that the correlation length on the odd sites is given by 
     
       
         
           
             
               
                 
                   ξ 
                   = 
                   
                     
                       1 
                       
                         ln 
                         ⁢ 
                            
                         
                           ( 
                           
                             
                               1 
                               + 
                               α 
                             
                             
                               1 
                               - 
                               α 
                             
                           
                           ) 
                         
                       
                     
                     × 
                     
                       1 
                       ρ 
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   29 
                 
               
             
           
         
       
     
     where ρ is the fraction of measured sites between the two Z-operators which have an outcome s=−1. If one considers all possible measurement outcomes, then ρ×n is described by a binomial distribution B (n,p) where n is the number of even sites between two Z-operators, and 
     
       
         
           
             p 
             = 
             
               
                 1 
                 + 
                 α 
               
               2 
             
           
         
       
     
     is the probability of having outcome s=−1. Hence, for large n, p is described by the normal distribution N (ρ, σ) with standard deviation 
     
       
         
           
             σ 
             = 
             
               
                 
                   
                     p 
                     ⁡ 
                     ( 
                     
                       1 
                       - 
                       p 
                     
                     ) 
                   
                   n 
                 
               
               . 
             
           
         
       
     
     Hence, if 
     
       
         
           
             
               σ 
               ≪ 
               
                 p 
                 ⁡ 
                 ( 
                 
                   
                     
                       i 
                       . 
                       e 
                       . 
                           
                       n 
                     
                     ≫ 
                     
                       
                         1 
                         p 
                       
                       - 
                       1 
                     
                   
                   = 
                   
                     
                       1 
                       - 
                       α 
                     
                     
                       1 
                       + 
                       α 
                     
                   
                 
                 ) 
               
             
             , 
           
         
       
     
     then 1/ρ is normally distributed as N(1/p, σ/p 2 ). 
     Details for Preparing Cluster and GHZ State in ID 
     Consider 1D chain with N atoms with spacing a: 
     
       
         
           
             
               
                 
                   
                     
                       H 
                       ⁢ 
                       
                         ( 
                         
                           Ω 
                           , 
                           h 
                         
                         ) 
                       
                     
                     = 
                     
                       
                         
                           Ω 
                           2 
                         
                         ⁢ 
                         
                           
                             ∑ 
                             
                               n 
                               = 
                               1 
                             
                             N 
                           
                           
                             X 
                             n 
                           
                         
                       
                       - 
                       
                         
                           h 
                           2 
                         
                         ⁢ 
                         
                           
                             ∑ 
                             
                               n 
                               = 
                               1 
                             
                             N 
                           
                           
                             Z 
                             n 
                           
                         
                       
                       + 
                       
                         
                           
                             V 
                             ⁡ 
                             ( 
                             a 
                             ) 
                           
                           4 
                         
                         ⁢ 
                         
                           
                             ∑ 
                             
                               n 
                               = 
                               1 
                             
                             N 
                           
                           
                             
                               ∑ 
                               
                                 k 
                                 ≥ 
                                 1 
                               
                             
                             
                               
                                 v 
                                 k 
                               
                               ⁢ 
                               
                                 Z 
                                 n 
                               
                               ⁢ 
                               
                                 Z 
                                 
                                   n 
                                   + 
                                   k 
                                 
                               
                             
                           
                         
                       
                     
                   
                   ⁢ 
                       
                     
                   
                     
                       with 
                       ⁢ 
                         
                       h 
                     
                     = 
                     
                       δ 
                       + 
                       
                         V 
                         ⁢ 
                         
                           ( 
                           a 
                           ) 
                         
                         ⁢ 
                         
                           
                             ∑ 
                             
                               k 
                               ≥ 
                               1 
                             
                           
                           
                             v 
                             k 
                           
                         
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   30 
                 
               
             
           
         
       
     
     In the nearest-neighbor model described herein, v n&gt;1 =0; in the single-species model v n =1/n 6 ; in the dual species (where even and odd sites have different species), v 2n−1 =1/(2n−1) 6  and v 2n =0. However, the formulas will be valid for general v k . 
     Of interest is the time-evolution with the interaction term over a time 
     
       
         
           
             
               t 
               
                 S 
                 ⁢ 
                 P 
                 ⁢ 
                 T 
               
             
             = 
             
               
                 π 
                 
                   V 
                   ⁡ 
                   ( 
                   a 
                   ) 
                 
               
               . 
             
           
         
       
     
     In particular, consider the following unitary operator: 
     
       
         
           
             
               
                 
                   
                     U 
                     
                       S 
                       ⁢ 
                       P 
                       ⁢ 
                       T 
                     
                   
                   := 
                   
                     
                       e 
                       
                         
                           - 
                           i 
                         
                         ⁢ 
                         
                           t 
                           SPT 
                         
                         ⁢ 
                         
                           H 
                           ⁡ 
                           ( 
                           
                             0 
                             , 
                             0 
                           
                           ) 
                         
                       
                     
                     = 
                     
                       e 
                       
                         
                           - 
                           
                             
                               i 
                               ⁢ 
                               π 
                             
                             4 
                           
                         
                         ⁢ 
                         
                           
                               
                             ∑ 
                           
                           n 
                         
                         ⁢ 
                             
                         
                           
                             ∑ 
                               
                           
                           
                             k 
                             ≥ 
                             1 
                           
                         
                         ⁢ 
                         
                           v 
                           k 
                         
                         ⁢ 
                         
                           Z 
                           n 
                         
                         ⁢ 
                         
                           Z 
                           
                             n 
                             + 
                             k 
                           
                         
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   31 
                 
               
             
           
         
       
     
     This is referred to as the SPT-pulse. Indeed, in the nearest-neighbor model v k =δ k,1 , this would create the ideal cluster SPT phase, a state which is denoted by|cluster) and which is the ground state of the cluster model, H cluster =−Σ n  Z n−1 X n Z n+1 . 
     Realizing the Cluster State 
     Starting from the ground state of H(+∞, 0) and then applying the SPT-pulse, the state being creating is: 
     
       
         
           
             
               
                 
                   
                     
                       
                         
                           
                             
                               
                                 ❘ 
                                 &#34;\[LeftBracketingBar]&#34; 
                               
                               ψ 
                             
                             〉 
                           
                           = 
                           
                             
                               U 
                               SPT 
                             
                             ⁢ 
                             
                               
                                 ❘ 
                                 &#34;\[LeftBracketingBar]&#34; 
                               
                               - 
                             
                           
                         
                         〉 
                       
                       
                         ⊗ 
                         N 
                       
                     
                     = 
                     
                       
                         ∏ 
                         n 
                       
                       
                         
                           ∏ 
                           
                             k 
                             ≥ 
                             2 
                           
                         
                         
                           
                             [ 
                             
                               
                                 cos 
                                 ⁡ 
                                 ( 
                                 
                                   
                                     π 
                                     ⁢ 
                                     
                                       v 
                                       k 
                                     
                                   
                                   4 
                                 
                                 ) 
                               
                               - 
                               
                                 i 
                                 ⁢ 
                                 
                                   sin 
                                   ⁡ 
                                   ( 
                                   
                                     
                                       π 
                                       ⁢ 
                                       
                                         v 
                                         k 
                                       
                                     
                                     4 
                                   
                                   ) 
                                 
                                 ⁢ 
                                 
                                   Z 
                                   n 
                                 
                                 ⁢ 
                                 
                                   Z 
                                   
                                     n 
                                     + 
                                     k 
                                   
                                 
                               
                             
                             ] 
                           
                           ⁢ 
                           
                             
                               ❘ 
                               &#34;\[LeftBracketingBar]&#34; 
                             
                             cluster 
                           
                         
                       
                     
                   
                   〉 
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   32 
                 
               
             
           
         
       
     
     From this, various exact properties of the resulting state can be derived: 
     
       
         
           
             
               
                 
                                   
                   
                     
                       ℱ 
                       SPT 
                     
                     = 
                     
                       
                         
                           
                             
                               ❘ 
                               &#34;\[LeftBracketingBar]&#34; 
                             
                             
                               〈 
                               
                                 cluster 
                                 ❘ 
                                 ψ 
                               
                               〉 
                             
                             
                               ❘ 
                               &#34;\[RightBracketingBar]&#34; 
                             
                           
                           2 
                         
                         N 
                       
                       ≥ 
                       
                         
                           ∏ 
                           
                             k 
                             = 
                             2 
                           
                           ∞ 
                         
                           
                         
                           
                             cos 
                             2 
                           
                           ( 
                           
                             
                               π 
                               ⁢ 
                               
                                 v 
                                 k 
                               
                             
                             4 
                           
                           ) 
                         
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   33 
                 
               
             
             
               
                 
                                   
                   
                     
                       〈 
                       
                         ψ 
                         ⁢ 
                         
                           
                             ❘ 
                             &#34;\[LeftBracketingBar]&#34; 
                           
                           
                             
                               Z 
                               
                                 m 
                                 - 
                                 1 
                               
                             
                             ⁢ 
                             
                               X 
                               m 
                             
                             ⁢ 
                             
                               Z 
                               
                                 m 
                                 + 
                                 1 
                               
                             
                           
                           
                             ❘ 
                             &#34;\[RightBracketingBar]&#34; 
                           
                         
                         ⁢ 
                         ψ 
                       
                       〉 
                     
                     = 
                     
                       
                         ∏ 
                         
                           k 
                           = 
                           2 
                         
                         ∞ 
                       
                         
                       
                         
                           cos 
                           2 
                         
                         ⁢ 
                         
                           ( 
                           
                             
                               π 
                               ⁢ 
                               
                                 v 
                                 k 
                               
                             
                             2 
                           
                           ) 
                         
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   34 
                 
               
             
             
               
                 
                                   
                   
                     
                       〈 
                       
                         ψ 
                         ⁢ 
                         
                           
                             ❘ 
                             &#34;\[LeftBracketingBar]&#34; 
                           
                           
                             X 
                             n 
                           
                           
                             ❘ 
                             &#34;\[RightBracketingBar]&#34; 
                           
                         
                         ⁢ 
                         ψ 
                       
                       〉 
                     
                     = 
                     
                       
                         〈 
                         
                           ψ 
                           ⁢ 
                           
                             
                               ❘ 
                               &#34;\[LeftBracketingBar]&#34; 
                             
                             
                               Y 
                               n 
                             
                             
                               ❘ 
                               &#34;\[RightBracketingBar]&#34; 
                             
                           
                           ⁢ 
                           ψ 
                         
                         〉 
                       
                       = 
                       
                         
                           〈 
                           
                             ψ 
                             ⁢ 
                             
                               
                                 ❘ 
                                 &#34;\[LeftBracketingBar]&#34; 
                               
                               
                                 Z 
                                 n 
                               
                               
                                 ❘ 
                                 &#34;\[RightBracketingBar]&#34; 
                               
                             
                             ⁢ 
                             ψ 
                           
                           〉 
                         
                         = 
                         0 
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   35 
                 
               
             
             
               
                 
                   
                     〈 
                     
                       ψ 
                       ⁢ 
                       
                         
                           ❘ 
                           &#34;\[LeftBracketingBar]&#34; 
                         
                         
                           
                             X 
                             n 
                           
                           ⁢ 
                           
                             X 
                             
                               n 
                               + 
                               2 
                             
                           
                         
                         
                           ❘ 
                           &#34;\[RightBracketingBar]&#34; 
                         
                       
                       ⁢ 
                       ψ 
                     
                     〉 
                   
                   = 
                   
                     
                       
                         
                           c 
                           2 
                           2 
                         
                         ⁢ 
                         
                           t 
                           3 
                           2 
                         
                       
                       2 
                     
                     ⁢ 
                     
                       
                         ∑ 
                         
                           
                             σ 
                             = 
                             0 
                           
                           , 
                           1 
                         
                       
                       
                         
                           
                             ( 
                             
                               - 
                               1 
                             
                             ) 
                           
                           σ 
                         
                         ⁢ 
                         
                           
                             ∏ 
                             
                               k 
                               = 
                               3 
                             
                             ∞ 
                           
                             
                           
                             
                               [ 
                               
                                 
                                   
                                     c 
                                     k 
                                     4 
                                   
                                   ( 
                                   
                                     1 
                                     + 
                                     
                                       
                                         
                                           ( 
                                           
                                             - 
                                             1 
                                           
                                           ) 
                                         
                                         σ 
                                       
                                       ⁢ 
                                       
                                         t 
                                         
                                           k 
                                           - 
                                           1 
                                         
                                       
                                       ⁢ 
                                       
                                         t 
                                         
                                           k 
                                           + 
                                           1 
                                         
                                       
                                     
                                   
                                   ) 
                                 
                                 2 
                               
                               ] 
                             
                             . 
                           
                         
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   36 
                 
               
             
           
         
       
     
     The inequality in the first formula is true for v k  which decay fast enough with k (such as in the Rydberg case v k =1/k 6 ). The other three formulas are true for all v k . In the latter formula, we the shorthand is introduced c k : =cos(πv k /2) and t k : =tan(πv k /2). 
     These formulas are straightforwardly evaluated. First consider the single-species model v k =1/k 6 . In this case, a fidelity per site of    SPT ≳0.99985, i.e., which is very close to the ideal cluster state. E.g., even for N=100 atoms, the fidelity| cluster|ψ | 2 ≈0.985. Relatedly, the stabilizer defining the cluster state also has a very large value:  ψ|Z m−1 X n Z m+1 |ψ ≈0.9994. Moreover, 
         ψ| X   n |ψ =0and | ψ| X   n   X   n+2 |ψ |&lt;10 −10 .  Equation 37
 
     This means that measuring in the X-basis on, say, the even sites essentially corresponds to a completely uncorrelated flat distribution (to an extremely good approximation). 
     Creating the GHZ State Through Measuring Half the System 
     Turning now to measuring X n  on every odd site. Due to Equation 37, the resulting state is 
     
       
         
           
             
               
                 
                   
                     
                       
                         
                           
                             ❘ 
                             &#34;\[LeftBracketingBar]&#34; 
                           
                           
                             ψ 
                             proj 
                           
                         
                         〉 
                       
                       ≈ 
                       
                         
                           ∏ 
                           n 
                         
                         
                           
                             
                               
                                 s 
                                 
                                   
                                     2 
                                     ⁢ 
                                     n 
                                   
                                   - 
                                   1 
                                 
                               
                               + 
                               
                                 X 
                                 
                                   
                                     2 
                                     ⁢ 
                                     n 
                                   
                                   - 
                                   1 
                                 
                               
                             
                             
                               2 
                             
                           
                           ⁢ 
                           
                             
                               ❘ 
                               &#34;\[LeftBracketingBar]&#34; 
                             
                             ψ 
                           
                         
                       
                     
                     〉 
                   
                   , 
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   38 
                 
               
             
           
         
       
     
     where s n  is the measurement outcome. Note that after judiciously chosen single-site flips, one can write 
     
       
         
           
             
               
                 
                   
                     
                       
                         
                           
                             ❘ 
                             &#34;\[LeftBracketingBar]&#34; 
                           
                           
                             ψ 
                             proj 
                           
                         
                         〉 
                       
                       ≈ 
                       
                         
                           ∏ 
                           n 
                         
                         
                           
                             
                               1 
                               + 
                               
                                 X 
                                 
                                   
                                     2 
                                     ⁢ 
                                     n 
                                   
                                   - 
                                   1 
                                 
                               
                             
                             
                               2 
                             
                           
                           ⁢ 
                           
                             
                               ❘ 
                               &#34;\[LeftBracketingBar]&#34; 
                             
                             ψ 
                           
                         
                       
                     
                     〉 
                   
                   . 
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   39 
                 
               
             
           
         
       
     
     From this, it follows that | GHZ|ψ proj   | 2 ≥| SPT|ψ | 2 . I.e., a good fidelity for the cluster state implies a good fidelity for the GHZ state after preparation. E.g., for the above cluster fidelity of ≈0.985 for 100 atoms, we find a lower bound| GHZ|ψ proj   | 2 ≥0.97. However, it turns out that this underestimates the true result. In addition, this can be improved by noting that the correction due to the second-nearest neighbor can be removed by a single-site rotation. 
     Modified Stabilizer 
     Consider the effect of v 2 , leading to the state 
     
       
         
           
             
               
                 
                   
                     
                       
                         
                           
                             ❘ 
                             &#34;\[LeftBracketingBar]&#34; 
                           
                           ψ 
                         
                         〉 
                       
                       = 
                       
                         
                           exp 
                           ⁡ 
                           ( 
                           
                             
                               - 
                               
                                 
                                   i 
                                   ⁢ 
                                   π 
                                   ⁢ 
                                   
                                     v 
                                     2 
                                   
                                 
                                 4 
                               
                             
                             ⁢ 
                             Σ 
                             ⁢ 
                             
                               Z 
                               n 
                             
                             ⁢ 
                             
                               Z 
                               
                                 n 
                                 + 
                                 2 
                               
                             
                           
                           ) 
                         
                         ⁢ 
                         
                           
                             ❘ 
                             &#34;\[LeftBracketingBar]&#34; 
                           
                           cluster 
                         
                       
                     
                     〉 
                   
                   . 
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   40 
                 
               
             
           
         
       
     
     This is still a    2 ×   2  SPT, so measuring, for example, the even sites should give a cat state on the odd sites. To see this, rewrite 
     
       
         
           
             
               
                 
                   
                     
                       
                         
                           ❘ 
                           &#34;\[LeftBracketingBar]&#34; 
                         
                         ψ 
                       
                       〉 
                     
                     = 
                     
                       
                         ∏ 
                         n 
                       
                       
                         
                           exp 
                           ⁡ 
                           ( 
                           
                             
                               - 
                               
                                 
                                   i 
                                   ⁢ 
                                   π 
                                   ⁢ 
                                   
                                     v 
                                     2 
                                   
                                 
                                 4 
                               
                             
                             ⁢ 
                             
                               Z 
                               
                                 
                                   2 
                                   ⁢ 
                                   n 
                                 
                                 - 
                                 1 
                               
                             
                             ⁢ 
                             
                               Z 
                               
                                 
                                   2 
                                   ⁢ 
                                   n 
                                 
                                 + 
                                 1 
                               
                             
                           
                           ) 
                         
                         ⁢ 
                         
                           
                             ∏ 
                             n 
                           
                           
                             
                               exp 
                               ⁡ 
                               ( 
                               
                                 
                                   - 
                                   
                                     
                                       i 
                                       ⁢ 
                                       π 
                                       ⁢ 
                                       
                                         v 
                                         2 
                                       
                                     
                                     4 
                                   
                                 
                                 ⁢ 
                                 
                                   X 
                                   
                                     
                                       2 
                                       ⁢ 
                                       n 
                                     
                                     + 
                                     1 
                                   
                                 
                               
                               ) 
                             
                             ⁢ 
                             
                               
                                 ❘ 
                                 &#34;\[LeftBracketingBar]&#34; 
                               
                               cluster 
                             
                           
                         
                       
                     
                   
                   〉 
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   41 
                 
               
             
           
         
       
     
     Note there are no operators acting on the even sites. Hence, if one measures on the even sites, the projectors that project the state onto our measurement outcome commute with the above operators. The result is thus 
     
       
         
           
             
               
                 
                   
                     
                       
                         
                           
                             
                               
                                 ❘ 
                                 &#34;\[LeftBracketingBar]&#34; 
                               
                               
                                 ψ 
                                 out 
                               
                             
                             〉 
                           
                           = 
                           
                             
                               ∏ 
                               n 
                             
                             
                               
                                 exp 
                                 ⁡ 
                                 ( 
                                 
                                   
                                     - 
                                     
                                       
                                         i 
                                         ⁢ 
                                         π 
                                         ⁢ 
                                         
                                           v 
                                           2 
                                         
                                       
                                       4 
                                     
                                   
                                   ⁢ 
                                   
                                     Z 
                                     
                                       
                                         2 
                                         ⁢ 
                                         n 
                                       
                                       - 
                                       1 
                                     
                                   
                                   ⁢ 
                                   
                                     Z 
                                     
                                       
                                         2 
                                         ⁢ 
                                         n 
                                       
                                       + 
                                       1 
                                     
                                   
                                 
                                 ) 
                               
                               ⁢ 
                               
                                 
                                   ∏ 
                                   n 
                                 
                                 
                                   
                                     exp 
                                     ⁡ 
                                     ( 
                                     
                                       
                                         - 
                                         
                                           
                                             i 
                                             ⁢ 
                                             π 
                                             ⁢ 
                                             
                                               v 
                                               2 
                                             
                                           
                                           4 
                                         
                                       
                                       ⁢ 
                                       
                                         X 
                                         
                                           
                                             2 
                                             ⁢ 
                                             n 
                                           
                                           + 
                                           1 
                                         
                                       
                                     
                                     ) 
                                   
                                   ⁢ 
                                   
                                     
                                       ❘ 
                                       &#34;\[LeftBracketingBar]&#34; 
                                     
                                     cat 
                                   
                                 
                               
                             
                           
                         
                         〉 
                       
                       = 
                       
                         U 
                         ⁢ 
                         
                           
                             ❘ 
                             &#34;\[LeftBracketingBar]&#34; 
                           
                           cat 
                         
                       
                     
                     〉 
                   
                   , 
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   42 
                 
               
             
           
         
       
     
     where the |cat  state is in the diagonal basis (its precise configuration depending on the measurement outcome). This is a demonstration of the fact that Ising couplings on a given sublattice still give rise to a true cat state. 
     Furthermore, note that this cat state is an eigenstate of the two-point function of 
     
       
         
           
             
               
                 
                   
                     
                       UZ 
                       
                         
                           2 
                           ⁢ 
                           m 
                         
                         - 
                         1 
                       
                     
                     ⁢ 
                     
                       U 
                       † 
                     
                   
                   = 
                   
                     
                       Z 
                       
                         
                           2 
                           ⁢ 
                           m 
                         
                         - 
                         1 
                       
                     
                     ⁢ 
                     
                       
                         ∏ 
                         n 
                       
                       
                         
                           exp 
                           ⁡ 
                           ( 
                           
                             
                               - 
                               
                                 
                                   i 
                                   ⁢ 
                                   π 
                                   ⁢ 
                                   
                                     v 
                                     2 
                                   
                                 
                                 4 
                               
                             
                             ⁢ 
                             
                               Z 
                               
                                 
                                   2 
                                   ⁢ 
                                   n 
                                 
                                 - 
                                 1 
                               
                             
                             ⁢ 
                             
                               Z 
                               
                                 
                                   2 
                                   ⁢ 
                                   n 
                                 
                                 + 
                                 1 
                               
                             
                           
                           ) 
                         
                         ⁢ 
                             
                         
                           exp 
                           ⁡ 
                           ( 
                           
                             
                               
                                 i 
                                 ⁢ 
                                 π 
                                 ⁢ 
                                 
                                   v 
                                   2 
                                 
                               
                               2 
                             
                             ⁢ 
                             
                               X 
                               
                                 
                                   2 
                                   ⁢ 
                                   m 
                                 
                                 - 
                                 1 
                               
                             
                           
                           ) 
                         
                         ⁢ 
                         
                           
                             ∏ 
                             n 
                           
                           
                             exp 
                             ⁡ 
                             ( 
                             
                               
                                 
                                   i 
                                   ⁢ 
                                   π 
                                   ⁢ 
                                   
                                     v 
                                     2 
                                   
                                 
                                 4 
                               
                               ⁢ 
                               
                                 Z 
                                 
                                   
                                     2 
                                     ⁢ 
                                     n 
                                   
                                   - 
                                   1 
                                 
                               
                               ⁢ 
                               
                                 Z 
                                 
                                   
                                     2 
                                     ⁢ 
                                     n 
                                   
                                   + 
                                   1 
                                 
                               
                             
                             ) 
                           
                         
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   43 
                 
               
             
           
         
       
       
         
           
             
               
                 
                   = 
                   
                     
                       Z 
                       
                         
                           2 
                           ⁢ 
                           m 
                         
                         - 
                         1 
                       
                     
                     ( 
                     
                       
                         cos 
                         ⁡ 
                         ( 
                         
                           π 
                           ⁢ 
                           
                             v 
                             2 
                           
                           / 
                           2 
                         
                         ) 
                       
                       + 
                       
                         i 
                         ⁢ 
                         
                           sin 
                           ⁡ 
                           ( 
                           
                             π 
                             ⁢ 
                             
                               v 
                               2 
                             
                             / 
                             2 
                           
                           ) 
                         
                         ⁢ 
                         
                           X 
                           
                             
                               2 
                               ⁢ 
                               m 
                             
                             - 
                             1 
                           
                         
                         ⁢ 
                         
                           exp 
                           ⁡ 
                           ( 
                           
                             
                               
                                 i 
                                 ⁢ 
                                 π 
                                 ⁢ 
                                 
                                   v 
                                   2 
                                 
                               
                               2 
                             
                             ⁢ 
                             
                               Z 
                               
                                 
                                   2 
                                   ⁢ 
                                   m 
                                 
                                 - 
                                 3 
                               
                             
                             ⁢ 
                             
                               Z 
                               
                                 
                                   2 
                                   ⁢ 
                                   m 
                                 
                                 - 
                                 1 
                               
                             
                           
                           ) 
                         
                         ⁢ 
                             
                         
                           exp 
                           ⁡ 
                           ( 
                           
                             
                               
                                 i 
                                 ⁢ 
                                 π 
                                 ⁢ 
                                 
                                   v 
                                   2 
                                 
                               
                               2 
                             
                             ⁢ 
                             
                               Z 
                               
                                 
                                   2 
                                   ⁢ 
                                   m 
                                 
                                 - 
                                 1 
                               
                             
                             ⁢ 
                             
                               Z 
                               
                                 
                                   2 
                                   ⁢ 
                                   m 
                                 
                                 + 
                                 1 
                               
                             
                           
                           ) 
                         
                       
                     
                     ) 
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   44 
                 
               
             
           
         
       
       
         
           
             
               
                 
                                   
                   
                     = 
                     
                       
                         cZ 
                         
                           
                             2 
                             ⁢ 
                             m 
                           
                           - 
                           1 
                         
                       
                       - 
                       
                         
                           
                             sY 
                             
                               
                                 2 
                                 ⁢ 
                                 m 
                               
                               - 
                               1 
                             
                           
                           ( 
                           
                             c 
                             + 
                             
                               
                                 isZ 
                                 
                                   
                                     2 
                                     ⁢ 
                                     m 
                                   
                                   - 
                                   3 
                                 
                               
                               ⁢ 
                               
                                 Z 
                                 
                                   
                                     2 
                                     ⁢ 
                                     m 
                                   
                                   - 
                                   1 
                                 
                               
                             
                           
                           ) 
                         
                         ⁢ 
                         
                           ( 
                           
                             c 
                             + 
                             
                               
                                 isZ 
                                 
                                   
                                     2 
                                     ⁢ 
                                     m 
                                   
                                   - 
                                   1 
                                 
                               
                               ⁢ 
                               
                                 Z 
                                 
                                   
                                     2 
                                     ⁢ 
                                     m 
                                   
                                   - 
                                   3 
                                 
                               
                             
                           
                           ) 
                         
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   45 
                 
               
             
           
         
       
       
         
           
             
               
                 
                                   
                   
                     = 
                     
                       
                         cZ 
                         
                           
                             2 
                             ⁢ 
                             m 
                           
                           - 
                           1 
                         
                       
                       - 
                       
                         sY 
                         
                           
                             2 
                             ⁢ 
                             m 
                           
                           - 
                           1 
                         
                       
                       + 
                       
                         O 
                         ⁡ 
                         ( 
                         
                           s 
                           2 
                         
                         ) 
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   46 
                 
               
             
           
         
       
     
     where c=cos(πv 2 /2) and s=sin(πv 2 /2). This means that the Ising order parameter has simply been rotated, which can be undone by applying a field. 
     Preparation of Other Long-Range Entangled States 
     Here other states of interest are presented which can be prepared with the disclosed protocol. The first example is discussed in detail. The other examples work analogously. 
     Xu-Moore Model (Subsystem Symmetry Breaking) 
     Consider the square lattice with the two sublattices (A is red, B is blue) shown in  FIG.  7   . 
     The other states of interest are presented which can be prepared with the protocol described herein. The first example is discussed in detail. The other examples work analogously. 
     Start with a product state|+   ⊗(N     A     +N     B     )  and time-evolve it with the nearest-neighbor Ising Hamiltonian H=JΣ (n,m) Z n Z m  for a time 
     
       
         
           
             
               t 
               SPT 
             
             = 
             
               π 
               
                 4 
                 ⁢ 
                 
                   
                     ❘ 
                     &#34;\[LeftBracketingBar]&#34; 
                   
                   J 
                   
                     ❘ 
                     &#34;\[RightBracketingBar]&#34; 
                   
                 
               
             
           
         
       
     
     (for simplicity we now set J=1); the resulting state is a cluster state on the above graph, with stabilizer X v Π v′∈∂v Z v′ =1 (where v′∈∂v denotes that v and v′ are nearest neighbors). Upon measuring the A sublattice in the X basis, the Xu-Moore state is obtained on the B sublattice with Π v∈  Z v =±1, depending on the measurement outcome. Note that this state corresponds to spontaneous breaking of subsystem symmetries independent of measurement outcome. Acting with judiciously chosen ΠX v  can make it into the homogeneous Xu-Moore state with Π v∈  Z v =+1. 
     Despite not having topological order, this is an interesting state with exotic correlations. In particular, there is long-range order for four-point functions of Z if they form a rectangle on the tilted square lattice, whereas one-, two- and three-point functions will vanish. This is thus an analogue of the Schrödinger cat state with infinitely many branches in the wavefunction. 
     Rydberg Implementation: Single-Species 
     If a single species of Rydberg atoms is used (all targeted for the same Rydberg state), then the Rydberg excited states at different sites will interact with a potential V(r)=U/r 6 . Hence, the first correction to the cluster stabilizer is at the diagonal distance r=√{square root over (2)}a, giving a stabilizer 
     
       
         
           
             
               
                 
                   
                     〈 
                     
                       
                         X 
                         v 
                       
                       ⁢ 
                       
                         
                           ∏ 
                           
                             
                               v 
                               ⁢ 
                               ′ 
                             
                             ∈ 
                             
                               ∂ 
                               v 
                             
                           
                         
                         
                           Z 
                           
                             v 
                             ⁢ 
                             ′ 
                           
                         
                       
                     
                     〉 
                   
                   ≈ 
                   
                     
                       cos 
                       4 
                     
                     ( 
                     
                       π 
                       
                         2 
                         ⁢ 
                         
                           
                             2 
                           
                           6 
                         
                       
                     
                     ) 
                   
                   ≈ 
                   
                     0.93 
                     . 
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   47 
                 
               
             
           
         
       
     
     However, this can be drastically improved by using local addressing, where it is presumed that individual pulses can be applied on the four sublattices shown in  FIG.  8    (here A is subdivided into A and D; similarly we split B into B and C). 
     Eventually, the A and D atoms are measured after preparing the cluster state. Spin flips    λ =Π v∈λ X v  (for any of the sublattices A=A, B, C, D) can be effectively generated by time-evolving the X-term by π/2. If U(t) denotes time-evolving with the Ising interaction of the system, consider the following combined evolution: 
           A     C   U (π/8)   A     D   U (π/8)   C     D   U (π/4)  Equation 48
 
     This effectively generates a t=π/4 evolution for the Ising couplings connecting AB, AC, CD, BD, a t=0 evolution for AD and BC couplings, and a t=π/2 evolution for AA, BB, CC, DD couplings. In particular, this entirely eliminates the diagonal corrections (of length √{square root over (2)}). The result is a stabilizer 
     
       
         
           
             
               
                 
                   
                     〈 
                     
                       
                         X 
                         v 
                       
                       ⁢ 
                       
                         
                           ∏ 
                           
                             
                               v 
                               ⁢ 
                               ′ 
                             
                             ∈ 
                             
                               ∂ 
                               v 
                             
                           
                         
                         
                           Z 
                           
                             v 
                             ⁢ 
                             ′ 
                           
                         
                       
                     
                     〉 
                   
                   ≈ 
                   
                     
                       cos 
                       4 
                     
                     ( 
                     
                       2 
                       × 
                       
                         π 
                         
                           2 
                           × 
                           
                             2 
                             6 
                           
                         
                       
                     
                     ) 
                   
                   ≈ 
                   
                     0.9952 
                     . 
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   49 
                 
               
             
           
         
       
     
     The extra factor of two arose from these contributions adding up over all time-evolutions. Including further corrections, a more precise answer is obtained: 
     
       
         
           
             
               
                 
                   
                     
                       〈 
                       
                         
                           X 
                           v 
                         
                         ⁢ 
                         
                           
                             ∏ 
                             
                               
                                 v 
                                 ⁢ 
                                 ′ 
                               
                               ∈ 
                               
                                 ∂ 
                                 v 
                               
                             
                           
                           
                             Z 
                             
                               v 
                               ⁢ 
                               ′ 
                             
                           
                         
                       
                       〉 
                     
                     ≈ 
                     
                       
                         
                           cos 
                           4 
                         
                         ( 
                         
                           2 
                           × 
                           
                             π 
                             
                               2 
                               × 
                               
                                 2 
                                 6 
                               
                             
                           
                         
                         ) 
                       
                       × 
                       
                         
                           cos 
                           8 
                         
                         ( 
                         
                           π 
                           
                             2 
                             ⁢ 
                             
                               
                                 5 
                               
                               6 
                             
                           
                         
                         ) 
                       
                       × 
                       
                         
                           cos 
                           4 
                         
                         ( 
                         
                           2 
                           × 
                           
                             π 
                             
                               2 
                               ⁢ 
                               
                                 
                                   8 
                                 
                                 6 
                               
                             
                           
                         
                         ) 
                       
                     
                     ≈ 
                     0.9945 
                   
                   , 
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   50 
                 
               
             
           
         
       
     
     confirming that the result is barely changed by including further corrections. 
     Rydberg Implementation: Dual-Species 
     Now consider a dual-species implementation, where the same-sublattice (AA and BB) couplings can be made negligible. Then the first correction is at a distance r=√{square root over (5)}a, leading to the cluster stabilizer: 
     
       
         
           
             
               
                 
                   
                     〈 
                     
                       
                         X 
                         v 
                       
                       ⁢ 
                       
                         
                           ∏ 
                           
                             
                               v 
                               ⁢ 
                               ′ 
                             
                             ∈ 
                             
                               ∂ 
                               v 
                             
                           
                         
                         
                           Z 
                           
                             v 
                             ⁢ 
                             ′ 
                           
                         
                       
                     
                     〉 
                   
                   ≈ 
                   
                     
                       
                         cos 
                         8 
                       
                       ( 
                       
                         π 
                         
                           2 
                           ⁢ 
                           
                             
                               5 
                             
                             6 
                           
                         
                       
                       ) 
                     
                     × 
                     
                       
                         cos 
                         4 
                       
                       ( 
                       
                         π 
                         
                           2 
                           × 
                           
                             3 
                             6 
                           
                         
                       
                       ) 
                     
                   
                   ≈ 
                   
                     0.9994 
                     . 
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   51 
                 
               
             
           
         
       
     
     2D Color Code 
     A protocol for obtaining the toric code on the honeycomb lattice is discussed above (which required putting qubits on the sites and bonds of the honeycomb lattice). Similarly, the color code can be created on the honeycomb lattice, which requires putting qubits on the sites of the dice lattice. The resulting stabilizers for the dual-species set-up are: 
     
       
         
           
             
               
                 
                   
                     〈 
                     
                       
                         X 
                         v 
                       
                       ⁢ 
                       
                         
                           ∏ 
                           
                             b 
                             ∈ 
                             v 
                           
                         
                         
                           Z 
                           b 
                         
                       
                     
                     〉 
                   
                   ≈ 
                   
                     
                       
                         cos 
                         6 
                       
                       ( 
                       
                         π 
                         
                           2 
                           ⁢ 
                           
                             Δ2 
                             6 
                           
                         
                       
                       ) 
                     
                     × 
                     
                       
                         cos 
                         6 
                       
                       ( 
                       
                         π 
                         
                           2 
                           × 
                           
                             
                               7 
                             
                             6 
                           
                         
                       
                       ) 
                     
                     × 
                     ⋯ 
                   
                   ≈ 
                   0.998 
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   52 
                 
               
             
             
               
                 
                   
                     〈 
                     
                       B 
                       p 
                     
                     〉 
                   
                   ≈ 
                   
                     
                       
                         cos 
                         18 
                       
                       ( 
                       
                         π 
                         
                           2 
                           × 
                           
                             2 
                             6 
                           
                         
                       
                       ) 
                     
                     × 
                     
                       
                         cos 
                         36 
                       
                       ( 
                       
                         π 
                         
                           2 
                           × 
                           
                             
                               7 
                             
                             6 
                           
                         
                       
                       ) 
                     
                     × 
                     ⋯ 
                   
                   ≈ 
                   0.994 
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   53 
                 
               
             
           
         
       
     
     3D Toric Code on the Diamond Lattice 
     Consider the 3D toric code. Instead of placing it on the usual cubic lattice, excellent results are obtained for the diamond lattice. In particular, referring to  FIG.  9   , consider the A sublattice to be the vertices of the diamond lattice, and the B sublattice the bonds. In the dual-species set-up, we then find: 
     
       
         
           
             
               
                 
                   
                     〈 
                     
                       
                         X 
                         v 
                       
                       ⁢ 
                       
                         
                           ∏ 
                           
                             b 
                             ∈ 
                             tet 
                           
                         
                         
                           Z 
                           b 
                         
                       
                     
                     〉 
                   
                   ≈ 
                   
                     
                       
                         cos 
                         12 
                       
                       ( 
                       
                         π 
                         
                           2 
                           · 
                           
                             
                               
                                 19 
                                 / 
                                 3 
                               
                             
                             6 
                           
                         
                       
                       ) 
                     
                     × 
                     
                       
                         cos 
                         24 
                       
                       ( 
                       
                         π 
                         
                           2 
                           · 
                           
                             
                               
                                 35 
                                 / 
                                 3 
                               
                             
                             6 
                           
                         
                       
                       ) 
                     
                     × 
                     
                       
                         cos 
                         24 
                       
                       ( 
                       
                         π 
                         
                           2 
                           · 
                           
                             
                               
                                 51 
                                 / 
                                 3 
                               
                             
                             6 
                           
                         
                       
                       ) 
                     
                     × 
                     ⋯ 
                   
                   ≈ 
                   0.9998 
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   54 
                 
               
             
             
               
                 
                   
                     〈 
                     
                       B 
                       p 
                     
                     〉 
                   
                   ≈ 
                   
                     
                       
                         cos 
                         18 
                       
                       ( 
                       
                         π 
                         
                           2 
                           · 
                           
                             
                               
                                 19 
                                 / 
                                 3 
                               
                             
                             6 
                           
                         
                       
                       ) 
                     
                     × 
                     
                       
                         cos 
                         36 
                       
                       ( 
                       
                         π 
                         
                           2 
                           · 
                           
                             
                               
                                 35 
                                 / 
                                 3 
                               
                             
                             6 
                           
                         
                       
                       ) 
                     
                     × 
                     
                       
                         cos 
                         36 
                       
                       ( 
                       
                         π 
                         
                           2 
                           · 
                           
                             
                               
                                 51 
                                 / 
                                 3 
                               
                             
                             6 
                           
                         
                       
                       ) 
                     
                     ⁢ 
                     
                       ( 
                       
                         
                           × 
                           ⋯ 
                         
                         ≈ 
                         0.9996 
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   55 
                 
               
             
           
         
       
     
     To obtain the toric code on the diamond lattice, the cluster state is created on the above lattice and then measure the A (red) sublattice. 
     X-Cube Model 
     To realize the X-cube model, the A and B sublattices are taken to be the vertices and the faces centers of the FCC lattice, respectively (see  FIG.  10 A  and  FIG.  10 B ). Defining A c =Π b∈cube Z b , we obtain for the dual-species set-up: 
     
       
         
           
             
               
                 
                   
                     〈 
                     
                       
                         X 
                         v 
                       
                       ⁢ 
                       
                         A 
                         c 
                       
                     
                     〉 
                   
                   ≈ 
                   
                     
                       
                         cos 
                         24 
                       
                       ( 
                       
                         π 
                         
                           2 
                           · 
                           
                             
                               3 
                             
                             6 
                           
                         
                       
                       ) 
                     
                     × 
                     
                       
                         cos 
                         24 
                       
                       ( 
                       
                         π 
                         
                           2 
                           · 
                           
                             
                               5 
                             
                             6 
                           
                         
                       
                       ) 
                     
                   
                   ≈ 
                   0.957 
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   56 
                 
               
             
           
         
       
       
         
           
             
               
                 
                   
                     〈 
                     
                       B 
                       p 
                     
                     〉 
                   
                   ≈ 
                   
                     
                       
                         cos 
                         32 
                       
                       ( 
                       
                         π 
                         
                           2 
                           · 
                           
                             
                               3 
                             
                             6 
                           
                         
                       
                       ) 
                     
                     × 
                     
                       
                         cos 
                         32 
                       
                       ( 
                       
                         π 
                         
                           2 
                           · 
                           
                             
                               5 
                             
                             6 
                           
                         
                       
                       ) 
                     
                   
                   ≈ 
                   0.944 
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   57 
                 
               
             
           
         
       
     
     Details for Preparing S 3  Topological Order 
     A formal outline of the preparation of S 3  topological order is provided below, followed by a discussion of the physical implementation in Rydberg atom arrays. 
     Formal Preparation 
     The preparation consists of two steps. Preparing the    3  toric code by gauging a product state with    3  symmetry, and further gauging the    2  charge conjugation symmetry to obtain the S 3  topological order. 
     The lattice we will use in the first step is the so-called Lieb lattice with the A, B sublattice on the vertices and bonds of the square lattice, respectively (see  FIG.  11   ). We denote 
     
       
         
           
             
               
                 
                   
                     
                       
                         𝒳 
                         = 
                         
                           ( 
                           
                             
                               
                                 0 
                               
                               
                                 0 
                               
                               
                                 1 
                               
                             
                             
                               
                                 1 
                               
                               
                                 0 
                               
                               
                                 0 
                               
                             
                             
                               
                                 0 
                               
                               
                                 1 
                               
                               
                                 0 
                               
                             
                           
                           ) 
                         
                       
                     
                     
                       
                         𝒵 
                         = 
                         
                           ( 
                           
                             
                               
                                 1 
                               
                               
                                 0 
                               
                               
                                 0 
                               
                             
                             
                               
                                 0 
                               
                               
                                 ω 
                               
                               
                                 0 
                               
                             
                             
                               
                                 0 
                               
                               
                                 0 
                               
                               
                                 
                                   ω 
                                   _ 
                                 
                               
                             
                           
                           ) 
                         
                       
                     
                     
                       
                         𝒞 
                         = 
                         
                           ( 
                           
                             
                               
                                 1 
                               
                               
                                 0 
                               
                               
                                 0 
                               
                             
                             
                               
                                 0 
                               
                               
                                 0 
                               
                               
                                 1 
                               
                             
                             
                               
                                 0 
                               
                               
                                 1 
                               
                               
                                 0 
                               
                             
                           
                           ) 
                         
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   58 
                 
               
             
           
         
       
     
     to be the shift, clock, and charge conjugation of qutrits respectively. Here ω=e 2π/3 . 
     A qutrit is initialized in the state 
     
       
         
           
             
               
                 
                   
                     ❘ 
                     &#34;\[LeftBracketingBar]&#34; 
                   
                   + 
                 
                 〉 
               
               = 
               
                 
                   
                     
                       
                         
                           
                             
                               ❘ 
                               &#34;\[LeftBracketingBar]&#34; 
                             
                             0 
                           
                           〉 
                         
                         + 
                         
                           
                             ❘ 
                             &#34;\[LeftBracketingBar]&#34; 
                           
                           1 
                         
                       
                       〉 
                     
                     + 
                     
                       
                         ❘ 
                         &#34;\[LeftBracketingBar]&#34; 
                       
                       2 
                     
                   
                   〉 
                 
                 
                   3 
                 
               
             
             , 
           
         
       
     
     which is an eigenstate of X. Thus the state is uniquely specified by stabilizers   for every site. Next, an evolution is performed that creates the    3  cluster state on this lattice. This cluster state can be viewed as an SPT protected by a    3  global symmetry and a    3  1-form symmetry. The unitary is given by 
     
       
         
           
             
               
                 
                   
                     U 
                     AB 
                   
                   = 
                   
                     
                       ∏ 
                       
                         〈 
                         
                           a 
                           , 
                           b 
                         
                         〉 
                       
                     
                     
                       C 
                       ⁢ 
                       
                         𝒵 
                         ab 
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   59 
                 
               
             
           
         
       
     
     where C =ω n     a     n     b   , is the controlled—Z gate. Conjugating by this unitary, the stabilizers are given by  FIG.  12   . 
     Next a measurement is performed on the A sublattice in the   basis. In the ideal case, when all measurement outcomes are|+ , the ground state of the    3  toric code is obtained, given by stabilizers shown in  FIG.  13   . The other measurement outcomes which are eigenstates of   signify the presence of e and ē anyons of the toric code, which can be paired up and annihilated via local unitaries. 
     The toric code wavefunction respects a global charge conjugation symmetry Π B    . In order to gauge this symmetry, add two additional sublattices C and D as shown in  FIG.  14   . The additional sublattices consist of qubits initialized with stabilizers X. C is used to perform a symmetry transformation so that the charge conjugation symmetry acts as Π c  X. From this, this    2  symmetry can be gauged by coupling the C and D lattices using the    2  cluster state entangler and measuring the C sublattice. 
     The symmetry transformation we will perform will exchange the following symmetries 
     
       
         
           
             
               
                 
                   
                     
                       ∏ 
                       B 
                     
                     X 
                   
                   ↔ 
                   
                     
                       ∏ 
                       B 
                     
                     
                       X 
                       ⁢ 
                       
                         
                           ∏ 
                           C 
                         
                         𝒞 
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   60 
                 
               
             
           
         
       
     
     The latter enriches the    3  toric code by charge conjugating both e and m anyons. The unitary that achieves this transformation is given by 
     
       
         
           
             
               
                 
                   
                     U 
                     BC 
                   
                   = 
                   
                     
                       ∏ 
                       
                         〈 
                         
                           b 
                           , 
                           c 
                         
                         〉 
                       
                     
                     
                       C 
                       ⁢ 
                       
                         𝒞 
                         cb 
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   61 
                 
               
             
           
         
       
     
     where  bc  denotes nearest neighbor sites of b ∈ B and c ∈ C, and C   cb  is the controlled-  gate where c is the control. The stabilizers are now given by  FIG.  15   . 
     We then couple the C and D sublattices with a    2  cluster state. This is achieved via 
     
       
         
           
             
               
                 
                   
                     U 
                     CD 
                   
                   = 
                   
                     
                       ∏ 
                       
                         〈 
                         
                           c 
                           , 
                           d 
                         
                         〉 
                       
                     
                     
                       CZ 
                       cd 
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   62 
                 
               
             
           
         
       
     
     At this stage the stabilizers of the wavefunction are given by  FIG.  16   . 
     Now, we perform a measurement on the C sublattice. It is clear that the following vertex and plaquette terms constructed from the order two stabilizers are stabilizers of the final model, as shown in  FIG.  17   . 
     On the other hand, the order three stabilizers do not commute. Call the first    3  stabilizer A v . We note that 
     
       
         
           
             
               1 
               + 
               
                 A 
                 v 
               
               + 
               
                 A 
                 v 
                 † 
               
             
             3 
           
         
       
     
     is a projector for the wavefunction that commutes with the measurement X on the C sublattice. Therefore, the resulting projector after the measurement can be written down. Denote the number operator of the four qutrits according to the cardinals n N , n E , n W , n S . Since the operator is diagonal, it acts on the basis states as 
     
       
         
           
             
               
                 
                   
                     A 
                     v 
                   
                   = 
                   
                     ω 
                     
                       
                         
                           n 
                           N 
                         
                         ⁢ 
                         
                           Z 
                           1 
                         
                       
                       + 
                       
                         
                           n 
                           E 
                         
                         ⁢ 
                         
                           Z 
                           1 
                         
                       
                       + 
                       
                         
                           n 
                           S 
                         
                         ⁢ 
                         
                           Z 
                           2 
                         
                       
                       + 
                       
                         
                           n 
                           W 
                         
                         ⁢ 
                         
                           Z 
                           3 
                         
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   63 
                 
               
             
           
         
       
       
         
           
             Therefore 
             , 
           
         
       
       
         
           
             
               
                 
                   
                     
                       A 
                       v 
                     
                     + 
                     
                       A 
                       v 
                       † 
                     
                   
                   = 
                   
                     2 
                     ⁢ 
                     
                       cos 
                       [ 
                       
                         
                           
                             2 
                             ⁢ 
                             π 
                           
                           3 
                         
                         ⁢ 
                         
                           
                             Z 
                             1 
                           
                           ( 
                           
                             
                               n 
                               N 
                             
                             + 
                             
                               n 
                               E 
                             
                             + 
                             
                               
                                 n 
                                 S 
                               
                               ⁢ 
                               
                                 Z 
                                 1 
                               
                               ⁢ 
                               
                                 Z 
                                 2 
                               
                             
                             + 
                             
                               
                                 n 
                                 W 
                               
                               ⁢ 
                               
                                 Z 
                                 1 
                               
                               ⁢ 
                               
                                 Z 
                                 3 
                               
                             
                           
                           ) 
                         
                       
                       ] 
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   64 
                 
               
             
           
         
       
     
     where we have pulled out a common factor Z 1 . Since this operator has eigenvalues±1, it can be removed from the cosine. We are left with 
     
       
         
           
             
               
                 
                   
                     
                       A 
                       v 
                     
                     + 
                     
                       A 
                       v 
                       † 
                     
                   
                   = 
                   
                     2 
                     ⁢ 
                     
                       cos 
                       [ 
                       
                         
                           
                             2 
                             ⁢ 
                             π 
                           
                           3 
                         
                         ⁢ 
                         
                           
                             Z 
                             1 
                           
                           ( 
                           
                             
                               n 
                               N 
                             
                             + 
                             
                               n 
                               E 
                             
                             + 
                             
                               
                                 n 
                                 S 
                               
                               ⁢ 
                               
                                 Z 
                                 1 
                               
                               ⁢ 
                               
                                 Z 
                                 2 
                               
                             
                             + 
                             
                               
                                 n 
                                 W 
                               
                               ⁢ 
                               
                                 Z 
                                 1 
                               
                               ⁢ 
                               
                                 Z 
                                 3 
                               
                             
                           
                           ) 
                         
                       
                       ] 
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   65 
                 
               
             
           
         
       
       
         
           
             
               
                 
                   = 
                   
                     
                       ω 
                       
                         
                           n 
                           N 
                         
                         + 
                         
                           n 
                           E 
                         
                         + 
                         
                           
                             n 
                             S 
                           
                           ⁢ 
                           
                             Z 
                             1 
                           
                           ⁢ 
                           
                             Z 
                             2 
                           
                         
                         + 
                         
                           
                             n 
                             W 
                           
                           ⁢ 
                           
                             Z 
                             1 
                           
                           ⁢ 
                           
                             Z 
                             3 
                           
                         
                       
                     
                     + 
                     
                       h 
                       . 
                       c 
                       . 
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   66 
                 
               
             
           
         
       
     
     The result is schematically represented in  FIG.  18   . 
     Therefore, combining with the order two stabilizers, the projector after the measurement is found to be the one schematically represented in  FIG.  19   . 
     Similarly, the projector that results from the order three plaquette term after the measurement is schematically represented in  FIG.  20   . 
     To conclude, the final wavefunction is not described by a stabilizer code, but instead by the following commuting projectors schematically represented in  FIG.  21    and  FIG.  22   . 
     Relation to the Quantum Double Model D(S 3 ) 
     To make a connection to the quantum double model, the two square lattices are shifted so that they coincide. Each edge now has a Hilbert space of dimension six composed of a qutrit and a qubit. It is shown that the projectors can be mapped exactly to those of the quantum double model defined on the square lattice. To make the explicit connection, the following transformations (which mutually commute) are made:
         1. Swap   and   on all qutrits;   2. Swap Z and X on all qubits;   3. Apply charge conjugation   in the pattern shown in  FIG.  23   .       

     After the above transformations, the projectors take the form depicted schematically in  FIG.  24    and  FIG.  25   . 
     These are exactly the projectors of the quantum double model of S 3 . The projectors for the quantum double model are the vertex term shown in  FIG.  26    and the plaquette terms shown in  FIG.  27   , 
     where L g  and R g  are left and right multiplication operators 
         L   g   |h   =| gh     R   g   |h   =| h   g     Equation 67
 
     Represent the group S 3  using the generators r, s that satisfy r 3 =s 2 =srsr=1. In the regular representation, the left and right multiplications of S 3  are 6×6 matrices, which can be explicitly written in terms of a qubit and a qutrit as 
         L   r   =I⊗X L   s   =X⊗C   Equation 68
 
         R   r   =     t⊕     =     −Z   R   s   =X⊗I   Equation 69
 
     Therefore, the product of our two vertex projectors is exactly A. Similarly, one can show that the product of our plaquette projectors is exactly B p . 
     Implementation with Rydberg Atoms 
     To realize a qutrit on the A and B sublattices, each site is endowed with two Rydberg atoms, which will be called atoms 1 and 2. Due to the proximity of the atoms, the state where both atoms are excited is prohibited, giving an effective three-level system. Label the empty state as |0  and the state where atom 1 or 2 is excited as |1  or |2  respectively. On the other hand, the C and D sublattice will be the usual Rydberg two-level system: |0  for the ground state and |1  for the excited state. 
     Assuming local addressing, the following pulses for the qutrits can be performed: 
     
       
         
           
             
               
                 
                   
                     
                       X 
                       1 
                     
                     = 
                     
                       ( 
                       
                         
                           
                             0 
                           
                           
                             1 
                           
                           
                             0 
                           
                         
                         
                           
                             1 
                           
                           
                             0 
                           
                           
                             0 
                           
                         
                         
                           
                             0 
                           
                           
                             0 
                           
                           
                             0 
                           
                         
                       
                       ) 
                     
                   
                   , 
                   
                     
                       Z 
                       1 
                     
                     = 
                     
                       ( 
                       
                         
                           
                             1 
                           
                           
                             0 
                           
                           
                             0 
                           
                         
                         
                           
                             0 
                           
                           
                             
                               - 
                               1 
                             
                           
                           
                             0 
                           
                         
                         
                           
                             0 
                           
                           
                             0 
                           
                           
                             0 
                           
                         
                       
                       ) 
                     
                   
                   , 
                   
 
                   
                     
                       X 
                       2 
                     
                     = 
                     
                       ( 
                       
                         
                           
                             0 
                           
                           
                             0 
                           
                           
                             1 
                           
                         
                         
                           
                             0 
                           
                           
                             0 
                           
                           
                             0 
                           
                         
                         
                           
                             1 
                           
                           
                             0 
                           
                           
                             0 
                           
                         
                       
                       ) 
                     
                   
                   , 
                   
                     
                       Z 
                       2 
                     
                     = 
                     
                       ( 
                       
                         
                           
                             1 
                           
                           
                             0 
                           
                           
                             0 
                           
                         
                         
                           
                             0 
                           
                           
                             0 
                           
                           
                             0 
                           
                         
                         
                           
                             0 
                           
                           
                             0 
                           
                           
                             
                               - 
                               1 
                             
                           
                         
                       
                       ) 
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   70 
                 
               
             
           
         
       
     
     which are the result of projecting the spin rotation operators into the qutrit subspace. These four generators are enough to generate arbitrary SU(3) rotations. 
     The preparation procedure can be summarized as follows:
         1. Load the A and B sublattices and perform the    3  fourier transform;   2. Implement U AB  Equation 59;   3. Perform the inverse    3  fourier transform on the A sublattice and measure the occupancy (   3  toric code is obtained at this step);   4. Perform single site rotations to pair up charges in the    3  toric code;   5. Load the C sublattice and perform Hadamard;   6. Implement U BC  Equation 61;   7. Load the D sublattice and perform Hadamard;   8. Implement U CD  Equation 62;   9. Perform Hadamard on the C sublattice and measure. (S 3  topological order is obtained at this step);   10. (optional) Perform single site rotations to pair up charges in the S 3  topological order.       

     Implementation of U AB , U BC , and U CD  is discussed below. 
     U AB    
     To create the    3  cluster state, generate the unitary U AB  in Equation 59. 
     The innate two-body interaction has dominated by the nearest neighbor sites between sublattices A and B, and is given by 
     
       
         
           
             
               
                 
                   
                     
                       
                         H 
                         AB 
                       
                       ( 
                       
                         U 
                         , 
                         
                           U 
                           ′ 
                         
                       
                       ) 
                     
                     = 
                     
                       
                         ∑ 
                         
                           〈 
                           
                             a 
                             , 
                             b 
                           
                           〉 
                         
                       
                       
                         [ 
                         
                           
                             U 
                             ⁡ 
                             ( 
                             
                               
                                 
                                   n 
                                   
                                     1 
                                     , 
                                     a 
                                   
                                 
                                 ⁢ 
                                 
                                   n 
                                   
                                     1 
                                     , 
                                     b 
                                   
                                 
                               
                               + 
                               
                                 
                                   n 
                                   
                                     2 
                                     , 
                                     a 
                                   
                                 
                                 ⁢ 
                                 
                                   n 
                                   
                                     2 
                                     , 
                                     b 
                                   
                                 
                               
                             
                             ) 
                           
                           + 
                           
                             
                               U 
                               ′ 
                             
                             ( 
                             
                               
                                 
                                   n 
                                   
                                     1 
                                     , 
                                     a 
                                   
                                 
                                 ⁢ 
                                 
                                   n 
                                   
                                     2 
                                     , 
                                     b 
                                   
                                 
                               
                               + 
                               
                                 
                                   n 
                                   
                                     2 
                                     , 
                                     a 
                                   
                                 
                                 ⁢ 
                                 
                                   n 
                                   
                                     1 
                                     , 
                                     b 
                                   
                                 
                               
                             
                             ) 
                           
                         
                         ] 
                       
                     
                   
                   , 
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   71 
                 
               
             
           
         
       
     
     where it is presumed 1↔1 and 2↔2 are equidistant, as are 1↔2 and 2↔1. (One way of achieving this is by having the two atoms separated in the third direction.) Here,  a, b  denote nearest neighbors, and n 1  and n 2  denotes the occupancy of atoms 1 and 2, respectively. Note that by conjugating the number operator n i  with a 
     
       
         
           
             π 
             2 
           
         
       
     
     pulse of X j  for j=1,2 gives 
     
       
         
           
             
               
                 
                   
                     
                       e 
                       
                         
                           
                             π 
                             ⁢ 
                             i 
                           
                           2 
                         
                         ⁢ 
                         
                           X 
                           j 
                         
                       
                     
                     ⁢ 
                     
                       n 
                       j 
                     
                     ⁢ 
                     
                       e 
                       
                         
                           - 
                           
                             
                               π 
                               ⁢ 
                               i 
                             
                             2 
                           
                         
                         ⁢ 
                         
                           X 
                           j 
                         
                       
                     
                   
                   = 
                   
                     1 
                     - 
                     
                       n 
                       1 
                     
                     - 
                     
                       
                         n 
                         2 
                       
                       . 
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   72 
                 
               
             
           
         
       
     
     One can show that by successive applications of this substitution, as well as the freedom of tuning single-site chemical potentials, one can generate an effective evolution where U′=−U. In fact, if U′&gt;U in Equation 71, the effective evolution will have U=−U′&gt;0. If U′&lt;U, one can simply exchange the 1↔2 labels on the A sublattice to reduce back to U′&gt;U. In conclusion, we thus obtain the C  gate U AB  via 
         U   AB   =e   −4πi/(3U)H(U,−U) .  Equation 73
 
     U BC    
     To implement U BC  in Equation 61, a basis transformation on the B sublattice is performed first, such that   is diagonal. A choice is 
     
       
         
           
             
               
                 
                   U 
                   = 
                   
                     
                       e 
                       
                         
                           
                             π 
                             ⁢ 
                             i 
                           
                           4 
                         
                         ⁢ 
                         
                           X 
                           1 
                         
                       
                     
                     ⁢ 
                     
                       e 
                       
                         
                           
                             π 
                             ⁢ 
                             i 
                           
                           4 
                         
                         ⁢ 
                         
                           Z 
                           1 
                         
                       
                     
                     ⁢ 
                     
                       e 
                       
                         
                           
                             π 
                             ⁢ 
                             i 
                           
                           2 
                         
                         ⁢ 
                         
                           X 
                           2 
                         
                       
                     
                     ⁢ 
                     
                       e 
                       
                         
                           - 
                           
                             
                               π 
                               ⁢ 
                               i 
                             
                             2 
                           
                         
                         ⁢ 
                         
                           Z 
                           2 
                         
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   74 
                 
               
             
           
         
       
       
         
           
             
               
                 
                   
                     
                       𝒞 
                       ~ 
                     
                     ≡ 
                     
                       U 
                       ⁢ 
                       𝒞 
                       ⁢ 
                       
                         U 
                         † 
                       
                     
                   
                   = 
                   
                     ( 
                     
                       
                         
                           
                             - 
                             1 
                           
                         
                         
                           0 
                         
                         
                           0 
                         
                       
                       
                         
                           0 
                         
                         
                           1 
                         
                         
                           0 
                         
                       
                       
                         
                           0 
                         
                         
                           0 
                         
                         
                           1 
                         
                       
                     
                     ) 
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   75 
                 
               
             
           
         
       
     
     The leading order interactions between the B and C sublattices are of the form 
     
       
         
           
             
               
                 
                   
                     H 
                     
                       B 
                       ⁢ 
                       C 
                     
                   
                   = 
                   
                     
                       
                         ∑ 
                         
                           〈 
                           bc 
                           〉 
                         
                       
                       
                         
                           n 
                           
                             1 
                             , 
                             b 
                           
                         
                         ⁢ 
                         
                           n 
                           c 
                         
                       
                     
                     + 
                     
                       
                         n 
                         
                           2 
                           , 
                           b 
                         
                       
                       ⁢ 
                       
                         n 
                         c 
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   76 
                 
               
             
           
         
       
     
     Since 
         e   πin     c     (n     1,b     +n     2,b)   =diag(1,1,1,1,−1,−1)= Z   c   ·C     c,b   Equation 77
 
     where C  is the controlled   gate, 
     
       
         
           
             
               
                 
                   
                     e 
                     
                       π 
                       ⁢ 
                       i 
                       ⁢ 
                       
                         H 
                         
                           B 
                           ⁢ 
                           C 
                         
                       
                     
                   
                   = 
                   
                     
                       ∏ 
                       
                         〈 
                         
                           b 
                           , 
                           c 
                         
                         〉 
                       
                     
                     
                       C 
                       ⁢ 
                       
                         
                           𝒞 
                           ˜ 
                         
                         
                           c 
                           , 
                           b 
                         
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   78 
                 
               
             
           
         
       
     
     where we used the property that for a fixed c ∈ C, there are only two nearest neighbors b ∈ B, thus cancelling away Z c . To conclude, 
     
       
         
           
             
               
                 
                   
                     U 
                     
                       B 
                       ⁢ 
                       C 
                     
                   
                   = 
                   
                     
                       
                         ∏ 
                         
                           〈 
                           
                             b 
                             , 
                             c 
                           
                           〉 
                         
                       
                       
                         C 
                         ⁢ 
                         
                           𝒞 
                           
                             c 
                             , 
                             b 
                           
                         
                       
                     
                     = 
                     
                       U 
                       ⁢ 
                       
                         e 
                         
                           π 
                           ⁢ 
                           i 
                           ⁢ 
                           
                             H 
                             
                               B 
                               ⁢ 
                               C 
                             
                           
                         
                       
                       ⁢ 
                       
                         U 
                         † 
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   79 
                 
               
             
           
         
       
     
     U CD    
     U CD  in Equation 61 is the    2  cluster state entangler between the C and D sublattices, so this can be implemented via the ordinary time evolution. However, one has to ensure that there is no net evolution between the pairs BC and BD. This can be done by flipping all C and D qubits halfway through the evolution. That is, 
     
       
         
           
             
               
                 
                   
                     U 
                     
                       C 
                       ⁢ 
                       D 
                     
                   
                   = 
                   
                     
                       ( 
                       
                         
                           ∏ 
                           
                             C 
                             , 
                             D 
                           
                         
                         X 
                       
                       ) 
                     
                     ⁢ 
                     
                       
                         e 
                         
                           
                             
                               π 
                               ⁢ 
                               i 
                             
                             2 
                           
                           ⁢ 
                           H 
                         
                       
                       ( 
                       
                         
                           ∏ 
                           
                             C 
                             , 
                             D 
                           
                         
                         X 
                       
                       ) 
                     
                     ⁢ 
                     
                       e 
                       
                         
                           
                             π 
                             ⁢ 
                             i 
                           
                           2 
                         
                         ⁢ 
                         H 
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   80 
                 
               
             
           
         
       
     
     Detailed Experimental Protocol in 1D 
     Start with a state|ψ(0) , which corresponds to all Rydberg atoms in the ground state. I.e., in the spin notation, |ψ(0) =|↓   ⊗N . First apply an X-pulse to rotate this into a desired product state, after which apply the SPT-pulse. However, due to the interactions, this is not a pure on-site pulse. So to get a quantitatively correct description, one should first investigate the effects of this X-pulse. A useful identity: 
       exp( A+εB )=exp(ε∫ 0   1   e   sA   Be   −sA   ds+O (ε 2 ))×expA  Equation 81
 
     We define the following quantities: 
     
       
         
           
             
               t 
               
                 p 
                 ⁢ 
                 r 
                 ⁢ 
                 o 
                 ⁢ 
                 d 
               
             
             = 
             
               π 
               
                 2 
                 ⁢ 
                 Ω 
               
             
           
         
       
       
         
           
             e 
             
               
                 - 
                 i 
               
               ⁢ 
               
                 H 
                 ⁡ 
                 ( 
                 
                   Ω 
                   , 
                   0 
                 
                 ) 
               
               ⁢ 
               t 
             
           
         
       
       
         
           
             A 
             = 
             
               
                 - 
                 i 
               
               ⁢ 
               
                 π 
                 4 
               
               ⁢ 
               
                 
                   ∑ 
                   n 
                 
                 
                   X 
                   n 
                 
               
             
           
         
       
       
         
           
             B 
             = 
             
               
                 
                   - 
                   i 
                 
                 ⁢ 
                 
                   
                     ∑ 
                     n 
                   
                   
                     
                       Z 
                       n 
                     
                     ⁢ 
                     
                       Z 
                       
                         n 
                         + 
                         1 
                       
                     
                   
                 
               
               + 
               … 
             
           
         
       
       
         
           
             ε 
             = 
             
               
                 
                   
                     V 
                     ⁡ 
                     ( 
                     a 
                     ) 
                   
                   ⁢ 
                   t 
                 
                 4 
               
               = 
               
                 
                   
                     π 
                     ⁢ 
                     
                       V 
                       ⁡ 
                       ( 
                       a 
                       ) 
                     
                   
                   
                     8 
                     ⁢ 
                     Ω 
                   
                 
                 ≪ 
                 1 
               
             
           
         
       
     
     Product State Pulse 
     So at leading order in E: 
     
       
         
           
             
               
                 
                   
                     exp 
                     ⁡ 
                     ( 
                     
                       
                         - 
                         i 
                       
                       ⁢ 
                       
                         H 
                         ⁡ 
                         ( 
                         
                           Ω 
                           , 
                           0 
                         
                         ) 
                       
                       ⁢ 
                       
                         t 
                         
                           p 
                           ⁢ 
                           r 
                           ⁢ 
                           o 
                           ⁢ 
                           d 
                         
                       
                     
                     ) 
                   
                   = 
                   
                     
                       exp 
                       ⁡ 
                       ( 
                       
                         
                           - 
                           i 
                         
                         ⁢ 
                         ε 
                         ⁢ 
                         
                           
                             ∑ 
                             n 
                           
                           
                             
                               ∫ 
                               0 
                               1 
                             
                             
                               
                                 e 
                                 
                                   
                                     - 
                                     i 
                                   
                                   ⁢ 
                                   s 
                                   ⁢ 
                                   
                                     π 
                                     4 
                                   
                                   ⁢ 
                                   
                                     
                                       ∑ 
                                       m 
                                     
                                     
                                       X 
                                       m 
                                     
                                   
                                 
                               
                               ⁢ 
                               
                                 Z 
                                 n 
                               
                               ⁢ 
                               
                                 Z 
                                 
                                   n 
                                   + 
                                   1 
                                 
                               
                               ⁢ 
                               
                                 e 
                                 
                                   i 
                                   ⁢ 
                                   s 
                                   ⁢ 
                                   
                                     π 
                                     4 
                                   
                                   ⁢ 
                                   
                                     
                                       ∑ 
                                       m 
                                     
                                     
                                       X 
                                       m 
                                     
                                   
                                 
                               
                               ⁢ 
                               d 
                               ⁢ 
                               s 
                             
                           
                         
                       
                       ) 
                     
                     × 
                     
                       exp 
                       ⁡ 
                       ( 
                       
                         
                           - 
                           i 
                         
                         ⁢ 
                         
                           π 
                           4 
                         
                         ⁢ 
                         
                           
                             ∑ 
                             n 
                           
                           
                             X 
                             n 
                           
                         
                       
                       ) 
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   82 
                 
               
             
           
         
       
       
         
           
             
               
                 
                   = 
                   
                     
                       exp 
                       ⁡ 
                       ( 
                       
                         
                           - 
                           i 
                         
                         ⁢ 
                         ε 
                         ⁢ 
                         
                           
                             ∑ 
                             n 
                           
                           
                             
                               Z 
                               n 
                             
                             ⁢ 
                             
                               Z 
                               
                                 n 
                                 + 
                                 1 
                               
                             
                             ⁢ 
                             
                               
                                 ∫ 
                                 0 
                                 1 
                               
                               
                                 
                                   e 
                                   
                                     i 
                                     ⁢ 
                                     s 
                                     ⁢ 
                                     
                                       π 
                                       2 
                                     
                                     ⁢ 
                                     
                                       ( 
                                       
                                         
                                           X 
                                           n 
                                         
                                         + 
                                         
                                           X 
                                           
                                             n 
                                             + 
                                             1 
                                           
                                         
                                       
                                       ) 
                                     
                                   
                                 
                                 ⁢ 
                                 d 
                                 ⁢ 
                                 s 
                               
                             
                           
                         
                       
                       ) 
                     
                     × 
                     
                       exp 
                       ⁡ 
                       ( 
                       
                         
                           - 
                           i 
                         
                         ⁢ 
                         
                           π 
                           4 
                         
                         ⁢ 
                         
                           
                             ∑ 
                             n 
                           
                           
                             X 
                             n 
                           
                         
                       
                       ) 
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   83 
                 
               
             
           
         
       
       
         
           
             	 
             
               Note 
               : 
             
           
         
       
       
         
           
             
               
                 
                   
                     
                       ∫ 
                       0 
                       1 
                     
                     
                       
                         e 
                         
                           i 
                           ⁢ 
                           s 
                           ⁢ 
                           
                             π 
                             2 
                           
                           ⁢ 
                           
                             ( 
                             
                               
                                 X 
                                 n 
                               
                               + 
                               
                                 X 
                                 
                                   n 
                                   + 
                                   1 
                                 
                               
                             
                             ) 
                           
                         
                       
                       ⁢ 
                       d 
                       ⁢ 
                       s 
                     
                   
                   = 
                   
                     
                       ∫ 
                       0 
                       1 
                     
                     
                       
                         
                           [ 
                           
                             
                               cos 
                               ⁡ 
                               ( 
                               
                                 
                                   s 
                                   ⁢ 
                                   π 
                                 
                                 2 
                               
                               ) 
                             
                             + 
                             
                               i 
                               ⁢ 
                               
                                 sin 
                                 ⁡ 
                                 ( 
                                 
                                   
                                     s 
                                     ⁢ 
                                     π 
                                   
                                   2 
                                 
                                 ) 
                               
                               ⁢ 
                               
                                 X 
                                 n 
                               
                             
                           
                           ] 
                         
                         [ 
                         
                           
                             cos 
                             ⁡ 
                             ( 
                             
                               
                                 s 
                                 ⁢ 
                                 π 
                               
                               2 
                             
                             ) 
                           
                           + 
                           
                             i 
                             ⁢ 
                             
                               sin 
                               ⁡ 
                               ( 
                               
                                 
                                   s 
                                   ⁢ 
                                   π 
                                 
                                 2 
                               
                               ) 
                             
                             ⁢ 
                             
                               X 
                               
                                 n 
                                 + 
                                 1 
                               
                             
                           
                         
                         ] 
                       
                       ⁢ 
                       ds 
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   84 
                 
               
             
           
         
       
       
         
           
             
               
                 
                   = 
                   
                     
                       ∫ 
                       0 
                       1 
                     
                     
                       
                         [ 
                         ⁠ 
                         
                           
                             
                               cos 
                               2 
                             
                             ( 
                             
                               
                                 s 
                                 ⁢ 
                                 π 
                               
                               2 
                             
                             ) 
                           
                           - 
                           
                             
                               
                                 sin 
                                 2 
                               
                               ( 
                               
                                 
                                   s 
                                   ⁢ 
                                   π 
                                 
                                 2 
                               
                               ) 
                             
                             ⁢ 
                             
                               X 
                               n 
                             
                             ⁢ 
                             
                               X 
                               
                                 n 
                                 + 
                                 1 
                               
                             
                           
                           + 
                           
                             i 
                             ⁢ 
                             
                               sin 
                               ⁡ 
                               ( 
                               
                                 
                                   s 
                                   ⁢ 
                                   π 
                                 
                                 2 
                               
                               ) 
                             
                             ⁢ 
                             
                               cos 
                               ⁡ 
                               ( 
                               
                                 
                                   s 
                                   ⁢ 
                                   π 
                                 
                                 2 
                               
                               ) 
                             
                             ⁢ 
                             
                               ( 
                               
                                 
                                   X 
                                   n 
                                 
                                 + 
                                 
                                   X 
                                   
                                     n 
                                     + 
                                     1 
                                   
                                 
                               
                               ) 
                             
                           
                         
                         ] 
                       
                       ⁢ 
                       ds 
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   85 
                 
               
             
           
         
       
       
         
           
             
               
                 
                                    
                   
                     = 
                     
                       
                         1 
                         2 
                       
                       - 
                       
                         
                           1 
                           2 
                         
                         ⁢ 
                         
                           X 
                           n 
                         
                         ⁢ 
                         
                           X 
                           
                             n 
                             + 
                             1 
                           
                         
                       
                       + 
                       
                         i 
                         ⁢ 
                         
                           1 
                           π 
                         
                         ⁢ 
                         
                           ( 
                           
                             
                               X 
                               n 
                             
                             + 
                             
                               X 
                               
                                 n 
                                 + 
                                 1 
                               
                             
                           
                           ) 
                         
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   86 
                 
               
             
           
         
       
     
     So at leading order in ε: 
     
       
         
           
             
               
                 
                   
                     e 
                     
                       
                         - 
                         i 
                       
                       ⁢ 
                       
                         H 
                         ⁡ 
                         ( 
                         
                           Ω 
                           , 
                           0 
                         
                         ) 
                       
                       ⁢ 
                       
                         t 
                         
                           p 
                           ⁢ 
                           r 
                           ⁢ 
                           o 
                           ⁢ 
                           d 
                         
                       
                     
                   
                   = 
                   
                     
                       exp 
                       ⁡ 
                       ( 
                       
                         
                           - 
                           i 
                         
                         ⁢ 
                         
                           ε 
                           2 
                         
                         ⁢ 
                         
                           
                             ∑ 
                             n 
                           
                           
                             [ 
                             
                               
                                 
                                   Z 
                                   n 
                                 
                                 ⁢ 
                                 
                                   Z 
                                   
                                     n 
                                     + 
                                     1 
                                   
                                 
                               
                               + 
                               
                                 
                                   Y 
                                   n 
                                 
                                 ⁢ 
                                 
                                   Y 
                                   
                                     n 
                                     + 
                                     1 
                                   
                                 
                               
                               - 
                               
                                 
                                   2 
                                   π 
                                 
                                 ⁢ 
                                 
                                   ( 
                                   
                                     
                                       
                                         Y 
                                         n 
                                       
                                       ⁢ 
                                       
                                         Z 
                                         
                                           n 
                                           + 
                                           1 
                                         
                                       
                                     
                                     + 
                                     
                                       
                                         Z 
                                         n 
                                       
                                       ⁢ 
                                       
                                         Y 
                                         
                                           n 
                                           + 
                                           1 
                                         
                                       
                                     
                                   
                                   ) 
                                 
                               
                             
                             ] 
                           
                         
                       
                       ) 
                     
                     × 
                     
                       exp 
                       ⁡ 
                       ( 
                       
                         
                           - 
                           i 
                         
                         ⁢ 
                         
                           π 
                           4 
                         
                         ⁢ 
                         
                           
                             ∑ 
                             n 
                           
                           
                             X 
                             n 
                           
                         
                       
                       ) 
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   87 
                 
               
             
           
         
       
       
         
           
             
               
                 
                   = 
                   
                     
                       ∏ 
                       n 
                     
                     
                       
                         exp 
                         ⁡ 
                         ( 
                         
                           
                             - 
                             i 
                           
                           ⁢ 
                           
                             
                               ε 
                               2 
                             
                             [ 
                             
                               
                                 
                                   Z 
                                   n 
                                 
                                 ⁢ 
                                 
                                   Z 
                                   
                                     n 
                                     + 
                                     1 
                                   
                                 
                               
                               + 
                               
                                 
                                   Y 
                                   n 
                                 
                                 ⁢ 
                                 
                                   Y 
                                   
                                     n 
                                     + 
                                     1 
                                   
                                 
                               
                               - 
                               
                                 
                                   2 
                                   π 
                                 
                                 ⁢ 
                                 
                                   ( 
                                   
                                     
                                       
                                         Y 
                                         n 
                                       
                                       ⁢ 
                                       
                                         Z 
                                         
                                           n 
                                           + 
                                           1 
                                         
                                       
                                     
                                     + 
                                     
                                       
                                         Z 
                                         n 
                                       
                                       ⁢ 
                                       
                                         Y 
                                         
                                           n 
                                           + 
                                           1 
                                         
                                       
                                     
                                   
                                   ) 
                                 
                               
                             
                             ] 
                           
                         
                         ) 
                       
                       × 
                       
                         
                           ∏ 
                           n 
                         
                         
                           exp 
                           ⁡ 
                           ( 
                           
                             
                               - 
                               i 
                             
                             ⁢ 
                             
                               π 
                               4 
                             
                             ⁢ 
                             
                               X 
                               n 
                             
                           
                           ) 
                         
                       
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   88 
                 
               
             
           
         
       
     
     In the last step, Trotter factorization is used (which is accurate to order ε). (Note that the above directly extends to all interaction distances weighted by v k !) The resulting state is 
     
       
         
           
             
               
                 
                   
                     
                       
                         ❘ 
                         &#34;\[LeftBracketingBar]&#34; 
                       
                       
                         ψ 
                         ⁡ 
                         ( 
                         
                           t 
                           
                             p 
                             ⁢ 
                             r 
                             ⁢ 
                             o 
                             ⁢ 
                             d 
                           
                         
                         ) 
                       
                     
                     〉 
                   
                   = 
                   
                     
                       e 
                       
                         
                           - 
                           i 
                         
                         ⁢ 
                         
                           H 
                           ⁡ 
                           ( 
                           
                             Ω 
                             , 
                             0 
                           
                           ) 
                         
                         ⁢ 
                         
                           t 
                           
                             p 
                             ⁢ 
                             r 
                             ⁢ 
                             o 
                             ⁢ 
                             d 
                           
                         
                       
                     
                     ⁢ 
                     
                       
                         ❘ 
                         &#34;\[LeftBracketingBar]&#34; 
                       
                       ↓ 
                     
                   
                 
                 〉 
               
               
                 ⊗ 
                 N 
               
             
             = 
             
               
                 
                   ∏ 
                   n 
                 
                 
                   
                     exp 
                     ⁡ 
                     ( 
                     
                       
                         - 
                         i 
                       
                       ⁢ 
                       
                         
                           ε 
                           2 
                         
                         [ 
                         
                           
                             
                               Z 
                               n 
                             
                             ⁢ 
                             
                               Z 
                               
                                 n 
                                 + 
                                 1 
                               
                             
                           
                           + 
                           
                             
                               Y 
                               n 
                             
                             ⁢ 
                             
                               Y 
                               
                                 n 
                                 + 
                                 1 
                               
                             
                           
                           - 
                           
                             
                               2 
                               π 
                             
                             ⁢ 
                             
                               ( 
                               
                                 
                                   
                                     Y 
                                     n 
                                   
                                   ⁢ 
                                   
                                     Z 
                                     
                                       n 
                                       + 
                                       1 
                                     
                                   
                                 
                                 + 
                                 
                                   
                                     Z 
                                     n 
                                   
                                   ⁢ 
                                   
                                     Y 
                                     
                                       n 
                                       + 
                                       1 
                                     
                                   
                                 
                               
                               ) 
                             
                           
                         
                         ] 
                       
                     
                     ) 
                   
                   ⁢ 
                   
                     
                       ( 
                       
                         
                           
                             
                               
                                 
                                   
                                     
                                       ❘ 
                                       &#34;\[LeftBracketingBar]&#34; 
                                     
                                     ↑ 
                                   
                                   〉 
                                 
                                 n 
                               
                               + 
                               
                                 i 
                                 ⁢ 
                                 
                                   
                                     ❘ 
                                     &#34;\[LeftBracketingBar]&#34; 
                                   
                                   ↓ 
                                 
                               
                             
                             〉 
                           
                           n 
                         
                         
                           2 
                         
                       
                       ) 
                     
                     
                       ⊗ 
                       N 
                     
                   
                 
               
               = 
               
                 
                   exp 
                   ⁡ 
                   ( 
                   
                     
                       - 
                       i 
                     
                     ⁢ 
                     
                       
                         2 
                         ⁢ 
                         ε 
                       
                       
                         V 
                         ⁡ 
                         ( 
                         a 
                         ) 
                       
                     
                     ⁢ 
                     
                       H 
                       ⁡ 
                       ( 
                       
                         0 
                         , 
                         0 
                       
                       ) 
                     
                   
                   ) 
                 
                 × 
                 
                   
                     ∏ 
                     n 
                   
                   
                     
                       exp 
                       ⁡ 
                       ( 
                       
                         
                           - 
                           i 
                         
                         ⁢ 
                         
                           ε 
                           π 
                         
                         ⁢ 
                         
                           ( 
                           
                             
                               
                                 Y 
                                 n 
                               
                               ⁢ 
                               
                                 Z 
                                 
                                   n 
                                   + 
                                   1 
                                 
                               
                             
                             + 
                             
                               
                                 
                                   Z 
                                   n 
                                 
                                 
                                   ︸ 
                                   
                                     
                                       = 
                                       
                                         : 
                                         
                                           
                                             ❘ 
                                             &#34;\[LeftBracketingBar]&#34; 
                                           
                                           φ 
                                         
                                       
                                     
                                     〉 
                                   
                                 
                               
                               ⁢ 
                               
                                 Y 
                                 
                                   n 
                                   + 
                                   1 
                                 
                               
                             
                           
                           ) 
                         
                       
                       ) 
                     
                     ⁢ 
                     
                       
                         ( 
                         
                           
                             
                               
                                 
                                   
                                     
                                       
                                         ❘ 
                                         &#34;\[LeftBracketingBar]&#34; 
                                       
                                       ↑ 
                                     
                                     〉 
                                   
                                   n 
                                 
                                 + 
                                 
                                   i 
                                   ⁢ 
                                   
                                     
                                       ❘ 
                                       &#34;\[LeftBracketingBar]&#34; 
                                     
                                     ↓ 
                                   
                                 
                               
                               〉 
                             
                             n 
                           
                           
                             2 
                           
                         
                         ) 
                       
                       
                         ⊗ 
                         N 
                       
                     
                   
                 
               
             
           
         
       
     
     This is a  =2 state! Note if c=cos(ε/π) and s=sin(ε/π), then 
     
       
         
           
             
               
                 
                   ❘ 
                   &#34;\[LeftBracketingBar]&#34; 
                 
                 ϕ 
               
               〉 
             
             = 
             
               
                 
                   ∏ 
                   n 
                 
                 
                   
                     ( 
                     
                       c 
                       - 
                       
                         i 
                         ⁢ 
                         s 
                         ⁢ 
                         
                           Y 
                           
                             n 
                             - 
                             1 
                           
                         
                         ⁢ 
                         
                           Z 
                           n 
                         
                       
                     
                     ) 
                   
                   ⁢ 
                   
                     ( 
                     
                       c 
                       - 
                       
                         i 
                         ⁢ 
                         s 
                         ⁢ 
                         
                           Z 
                           n 
                         
                         ⁢ 
                         
                           Y 
                           
                             n 
                             + 
                             1 
                           
                         
                       
                     
                     ) 
                   
                   ⁢ 
                   
                     
                       ( 
                       
                         
                           
                             
                               
                                 
                                   
                                     
                                       ❘ 
                                       &#34;\[LeftBracketingBar]&#34; 
                                     
                                     ↑ 
                                   
                                   〉 
                                 
                                 n 
                               
                               + 
                               
                                 i 
                                 ⁢ 
                                 
                                   
                                     ❘ 
                                     &#34;\[LeftBracketingBar]&#34; 
                                   
                                   ↓ 
                                 
                               
                             
                             〉 
                           
                           n 
                         
                         
                           2 
                         
                       
                       ) 
                     
                     
                       ⊗ 
                       N 
                     
                   
                 
               
               = 
               
                 
                   ∏ 
                   n 
                 
                 
                   
                     ( 
                     
                       
                         c 
                         2 
                       
                       - 
                       
                         
                           s 
                           2 
                         
                         ⁢ 
                         
                           Y 
                           
                             n 
                             - 
                             1 
                           
                         
                         ⁢ 
                         
                           Y 
                           
                             n 
                             + 
                             1 
                           
                         
                       
                       - 
                       
                         i 
                         ⁢ 
                         c 
                         ⁢ 
                         s 
                         ⁢ 
                         
                           
                             Z 
                             n 
                           
                           ( 
                           
                             
                               Y 
                               
                                 n 
                                 - 
                                 1 
                               
                             
                             + 
                             
                               Y 
                               
                                 n 
                                 + 
                                 1 
                               
                             
                           
                           ) 
                         
                       
                     
                     ) 
                   
                   ⁢ 
                   
                     
                       ( 
                       
                         
                           
                             
                               
                                 
                                   
                                     
                                       ❘ 
                                       &#34;\[LeftBracketingBar]&#34; 
                                     
                                     ↑ 
                                   
                                   〉 
                                 
                                 n 
                               
                               + 
                               
                                 i 
                                 ⁢ 
                                 
                                   
                                     ❘ 
                                     &#34;\[LeftBracketingBar]&#34; 
                                   
                                   ↓ 
                                 
                               
                             
                             〉 
                           
                           n 
                         
                         
                           2 
                         
                       
                       ) 
                     
                     
                       ⊗ 
                       N 
                     
                   
                 
               
             
           
         
       
     
     From this, we see that: 
     
       
         
           
             
               
                 
                   
                     〈 
                     φ 
                   
                   
                     ❘ 
                     &#34;\[RightBracketingBar]&#34; 
                   
                 
                 ⁢ 
                 
                   Y 
                   m 
                 
                 ⁢ 
                 
                   
                     ❘ 
                     &#34;\[LeftBracketingBar]&#34; 
                   
                   ϕ 
                 
               
               〉 
             
             = 
             
               
                 cos 
                 ⁡ 
                 ( 
                 
                   
                     4 
                     ⁢ 
                     ε 
                   
                   π 
                 
                 ) 
               
               = 
               
                 1 
                 - 
                 
                   
                     8 
                     ⁢ 
                     
                       ε 
                       2 
                     
                   
                   
                     π 
                     2 
                   
                 
                 + 
                 
                   O 
                   ⁡ 
                   ( 
                   
                     ε 
                     4 
                   
                   ) 
                 
               
             
           
         
       
     
     at leading order in ε. Also, 
     
       
         
           
             
               
                 
                   
                     
                       ❘ 
                       &#34;\[LeftBracketingBar]&#34; 
                     
                   
                   
                     〈 
                     
                       
                         + 
                         y 
                       
                       
                         | 
                         φ 
                       
                     
                     〉 
                   
                   
                     ❘ 
                     &#34;\[RightBracketingBar]&#34; 
                   
                 
                 2 
               
               N 
             
             = 
             
               
                 
                   cos 
                   2 
                 
                 ( 
                 
                   
                     2 
                     ⁢ 
                     ε 
                   
                   π 
                 
                 ) 
               
               = 
               
                 1 
                 - 
                 
                   
                     4 
                     ⁢ 
                     
                       ε 
                       2 
                     
                   
                   
                     π 
                     2 
                   
                 
                 + 
                 
                   O 
                   ⁡ 
                   ( 
                   
                     ε 
                     4 
                   
                   ) 
                 
               
             
           
         
       
     
     SPT Pulse 
     Application of the SPT-pulse is now desirable over a time 
     
       
         
           
             
               t 
               = 
               
                 
                   
                     t 
                     
                       S 
                       ⁢ 
                       P 
                       ⁢ 
                       T 
                     
                   
                   - 
                   
                     
                       2 
                       ⁢ 
                       ε 
                     
                     
                       V 
                       ⁡ 
                       ( 
                       a 
                       ) 
                     
                   
                 
                 = 
                 
                   
                     π 
                     - 
                     
                       2 
                       ⁢ 
                       ε 
                     
                   
                   
                     V 
                     ⁡ 
                     ( 
                     a 
                     ) 
                   
                 
               
             
             . 
           
         
       
     
     For the end result, one expects a state whose fidelity (per site) for the cluster state is approximately 
     
       
         
           
             ℱ 
             = 
             
               
                 
                   
                     
                       
                         〈 
                         cluster 
                       
                       
                         ❘ 
                         &#34;\[RightBracketingBar]&#34; 
                       
                     
                     ⁢ 
                     
                       ψ 
                       ⁡ 
                       ( 
                       
                         
                           t 
                           
                             p 
                             ⁢ 
                             r 
                             ⁢ 
                             o 
                             ⁢ 
                             d 
                           
                         
                         + 
                         
                           t 
                           
                             S 
                             ⁢ 
                             P 
                             ⁢ 
                             T 
                           
                         
                         - 
                       
                       ) 
                     
                   
                   〉 
                 
                 N 
               
               ≈ 
               
                 
                   ( 
                   
                     1 
                     - 
                     
                       
                         4 
                         ⁢ 
                         
                           ε 
                           2 
                         
                       
                       
                         π 
                         2 
                       
                     
                   
                   ) 
                 
                 ⁢ 
                 
                   
                     ∏ 
                     
                       k 
                       = 
                       2 
                     
                     ∞ 
                   
                   
                     
                       
                         cos 
                         2 
                       
                       ( 
                       
                         
                           π 
                           ⁢ 
                           
                             v 
                             k 
                           
                         
                         4 
                       
                       ) 
                     
                     . 
                   
                 
               
             
           
         
       
     
     For the Rydberg case v k =1/k 6 , we expect  ≈0.9958 for the choice ε=0.1, and  ≈0.99981 for ε=0.01. Doing the full numerical tensor-network simulation of the state preparation:  ≈0.9956 for ε=0.1 and  ≈0.99981 for ε=0.01. Note, the above fidelity is for the ground state of H=Σ n Z n−1 Y n Z n+1 , not H=−Σ n Z n−1 X n Z n+1 . 
     Pre-Measurement Pulse 
     Before measuring the even sites, rotation of those sites back to the z-basis is needed. I.e., time-evolution is desirable for a time 
     
       
         
           
             
               t 
               
                 m 
                 ⁢ 
                 e 
                 ⁢ 
                 a 
                 ⁢ 
                 s 
               
             
             = 
             
               
                 
                   
                     4 
                     ⁢ 
                     ε 
                   
                   
                     V 
                     ⁡ 
                     ( 
                     a 
                     ) 
                   
                 
                 ⁢ 
                     
                 where 
                 ⁢ 
                     
                 Ω 
               
               = 
               
                 - 
                 
                   
                     π 
                     ⁢ 
                     
                       V 
                       ⁡ 
                       ( 
                       a 
                       ) 
                     
                   
                   
                     8 
                     ⁢ 
                     ε 
                   
                 
               
             
           
         
       
     
     on even sites and Ω=0 on odd sites. 
     The time-evolution operator takes the form: 
     
       
         
           
             
               
                 
                   
                     exp 
                     ⁡ 
                     ( 
                     
                       
                         - 
                         iH 
                       
                       ⁢ 
                       
                         t 
                         
                           m 
                           ⁢ 
                           e 
                           ⁢ 
                           a 
                           ⁢ 
                           s 
                         
                       
                     
                     ) 
                   
                   = 
                     
                   
                     
                       exp 
                       ⁡ 
                       ( 
                       
                         i 
                         ⁢ 
                         
                           π 
                           4 
                         
                         ⁢ 
                         
                           
                             ∑ 
                             n 
                           
                           
                             X 
                             
                               2 
                               ⁢ 
                               n 
                             
                           
                         
                       
                       ) 
                     
                     × 
                     
                       exp 
                       ⁡ 
                       ( 
                       
                         
                           - 
                           i 
                         
                         ⁢ 
                         ε 
                         ⁢ 
                         
                           
                             ∑ 
                             n 
                           
                           
                             
                               ∫ 
                               0 
                               1 
                             
                             
                               
                                 e 
                                 
                                   
                                     - 
                                     i 
                                   
                                   ⁢ 
                                   s 
                                   ⁢ 
                                   
                                     π 
                                     4 
                                   
                                   ⁢ 
                                   
                                     
                                       ∑ 
                                       m 
                                     
                                     
                                       X 
                                       
                                         2 
                                         ⁢ 
                                         m 
                                       
                                     
                                   
                                 
                               
                               ⁢ 
                               
                                 Z 
                                 n 
                               
                               ⁢ 
                               
                                 Z 
                                 
                                   n 
                                   + 
                                   1 
                                 
                               
                               ⁢ 
                               
                                 e 
                                 
                                   i 
                                   ⁢ 
                                   s 
                                   ⁢ 
                                   
                                     π 
                                     4 
                                   
                                   ⁢ 
                                   
                                     
                                       ∑ 
                                       m 
                                     
                                     
                                       X 
                                       
                                         2 
                                         ⁢ 
                                         m 
                                       
                                     
                                   
                                 
                               
                               ⁢ 
                               d 
                               ⁢ 
                               s 
                             
                           
                         
                       
                       ) 
                     
                   
                 
               
             
             
               
                 
                   = 
                     
                   
                     
                       exp 
                       ⁡ 
                       ( 
                       
                         i 
                         ⁢ 
                         
                           π 
                           4 
                         
                         ⁢ 
                         
                           
                             ∑ 
                             n 
                           
                           
                             X 
                             
                               2 
                               ⁢ 
                               n 
                             
                           
                         
                       
                       ) 
                     
                     × 
                     
                       exp 
                       ⁡ 
                       ( 
                       
                         
                           - 
                           i 
                         
                         ⁢ 
                         ε 
                         ⁢ 
                         
                           
                             ∑ 
                             n 
                           
                           
                             
                               
                                 Z 
                                 
                                   2 
                                   ⁢ 
                                   n 
                                 
                               
                               ( 
                               
                                 
                                   Z 
                                   
                                     
                                       2 
                                       ⁢ 
                                       n 
                                     
                                     - 
                                     1 
                                   
                                 
                                 + 
                                 
                                   Z 
                                   
                                     
                                       2 
                                       ⁢ 
                                       n 
                                     
                                     + 
                                     1 
                                   
                                 
                               
                               ) 
                             
                             ⁢ 
                             
                               
                                 ∫ 
                                 0 
                                 1 
                               
                               
                                 
                                   e 
                                   
                                     i 
                                     ⁢ 
                                     s 
                                     ⁢ 
                                     
                                       π 
                                       2 
                                     
                                     ⁢ 
                                     
                                       X 
                                       
                                         2 
                                         ⁢ 
                                         n 
                                       
                                     
                                   
                                 
                                 ⁢ 
                                 d 
                                 ⁢ 
                                 s 
                               
                             
                           
                         
                       
                       ) 
                     
                   
                 
               
             
             
               
                 
                   = 
                     
                   
                     
                       exp 
                       ⁡ 
                       ( 
                       
                         i 
                         ⁢ 
                         
                           π 
                           4 
                         
                         ⁢ 
                         
                           
                             ∑ 
                             n 
                           
                           
                             X 
                             
                               2 
                               ⁢ 
                               n 
                             
                           
                         
                       
                       ) 
                     
                     × 
                     
                       exp 
                       ⁡ 
                       ( 
                       
                         
                           - 
                           i 
                         
                         ⁢ 
                         ε 
                         ⁢ 
                         
                           
                             ∑ 
                             n 
                           
                           
                             
                               
                                 Z 
                                 
                                   2 
                                   ⁢ 
                                   n 
                                 
                               
                               ( 
                               
                                 
                                   Z 
                                   
                                     
                                       2 
                                       ⁢ 
                                       n 
                                     
                                     - 
                                     1 
                                   
                                 
                                 + 
                                 
                                   Z 
                                   
                                     
                                       2 
                                       ⁢ 
                                       n 
                                     
                                     + 
                                     1 
                                   
                                 
                               
                               ) 
                             
                             ⁢ 
                             
                               2 
                               π 
                             
                             ⁢ 
                             
                               ( 
                               
                                 1 
                                 + 
                                 
                                   i 
                                   ⁢ 
                                   
                                     X 
                                     
                                       2 
                                       ⁢ 
                                       n 
                                     
                                   
                                 
                               
                               ) 
                             
                           
                         
                       
                       ) 
                     
                   
                 
               
             
             
               
                 
                   = 
                     
                   
                     
                       exp 
                       ⁡ 
                       ( 
                       
                         i 
                         ⁢ 
                         
                           π 
                           4 
                         
                         ⁢ 
                         
                           
                             ∑ 
                             n 
                           
                           
                             X 
                             
                               2 
                               ⁢ 
                               n 
                             
                           
                         
                       
                       ) 
                     
                     × 
                     
                       exp 
                       ⁡ 
                       ( 
                       
                         i 
                         ⁢ 
                         
                           
                             2 
                             ⁢ 
                             ε 
                           
                           π 
                         
                         ⁢ 
                         
                           
                             ∑ 
                             n 
                           
                           
                             
                               Y 
                               
                                 2 
                                 ⁢ 
                                 n 
                               
                             
                             ( 
                             
                               
                                 Z 
                                 
                                   
                                     2 
                                     ⁢ 
                                     n 
                                   
                                   - 
                                   1 
                                 
                               
                               + 
                               
                                 Z 
                                 
                                   
                                     2 
                                     ⁢ 
                                     n 
                                   
                                   + 
                                   1 
                                 
                               
                             
                             ) 
                           
                         
                       
                       ) 
                     
                     × 
                     
                       exp 
                       ⁡ 
                       ( 
                       
                         
                           - 
                           i 
                         
                         ⁢ 
                         
                           
                             8 
                             ⁢ 
                             ε 
                           
                           
                             π 
                             ⁢ 
                             
                               V 
                               ⁡ 
                               ( 
                               a 
                               ) 
                             
                           
                         
                         ⁢ 
                         
                           H 
                           ⁡ 
                           ( 
                           
                             0 
                             , 
                             0 
                           
                           ) 
                         
                       
                       ) 
                     
                   
                 
               
             
           
         
       
     
     Here all equalities are up to order ε. 
     It can be seen that the product state pulse applies an SPT pulse tor a time 
     
       
         
           
             
               2 
               ⁢ 
               ε 
             
             
               V 
               ⁡ 
               ( 
               a 
               ) 
             
           
         
       
     
     and the pre-measurement pulse applies an SPT pulse for a time 
     
       
         
           
             
               
                 8 
                 ⁢ 
                 ε 
               
               
                 π 
                 ⁢ 
                 
                   V 
                   ⁡ 
                   ( 
                   a 
                   ) 
                 
               
             
             . 
           
         
       
     
     Hence, if the goal is not to prepare the best cluster state, but instead to prepare the best GHZ state (after measuring the even sites), one should only apply the SPT pulse for a total time 
     
       
         
           
             s 
             = 
             
               
                 
                   t 
                   
                     S 
                     ⁢ 
                     P 
                     ⁢ 
                     T 
                   
                 
                 - 
                 
                   
                     2 
                     ⁢ 
                     ε 
                   
                   
                     V 
                     ⁡ 
                     ( 
                     a 
                     ) 
                   
                 
                 - 
                 
                   
                     8 
                     ⁢ 
                     ε 
                   
                   
                     π 
                     ⁢ 
                     
                       V 
                       ⁡ 
                       ( 
                       a 
                       ) 
                     
                   
                 
               
               = 
               
                 
                   
                     π 
                     2 
                   
                   - 
                   
                     
                       ( 
                       
                         
                           2 
                           ⁢ 
                           π 
                         
                         + 
                         8 
                       
                       ) 
                     
                     ⁢ 
                     ε 
                   
                 
                 
                   π 
                   ⁢ 
                   
                     V 
                     ⁡ 
                     ( 
                     a 
                     ) 
                   
                 
               
             
           
         
       
     
     E.g., for ε=0.1, this reduces the total time from the naive 
     
       
         
           
             
               π 
               
                 V 
                 ⁡ 
                 ( 
                 a 
                 ) 
               
             
             ≈ 
             
               3 
               .14 
               / 
               
                 V 
                 ⁡ 
                 ( 
                 a 
                 ) 
               
             
           
         
       
     
     to a time≈2.69/V(a), i.e., a reduction by more than 14%, which is significant. (For ε=0.01, and a reduction by slightly over 1.4%). 
     Doing the full numerical simulation, for ε=0.01, a state was found with (Z 2n−1 Z 2n Z 2n+1 )≈0.9993 and (Z 2n )&lt;10 −5  and (Z 2n Z 2n+2 )&lt;10 −7 . Measuring the even sites, and presuming it can be approximated as an equal weight distributions (as suggested by the aforementioned correlations),  ≈0.99984. 
     Addendum: SPT Model and Properties 
     Parent Hamiltonian for Equation 32 is 
     
       
         
           
             
               
                 
                   H 
                   = 
                     
                   
                     
                       
                         ∑ 
                         n 
                       
                       
                         
                           U 
                           
                             S 
                             ⁢ 
                             P 
                             ⁢ 
                             T 
                           
                         
                         ⁢ 
                         
                           X 
                           n 
                         
                         ⁢ 
                         
                           U 
                           
                             S 
                             ⁢ 
                             P 
                             ⁢ 
                             T 
                           
                           † 
                         
                       
                     
                     = 
                     
                       - 
                       
                         
                           ∑ 
                           n 
                         
                         
                           
                             Z 
                             
                               n 
                               - 
                               1 
                             
                           
                           ⁢ 
                           
                             X 
                             n 
                           
                           ⁢ 
                           
                             Z 
                             
                               n 
                               + 
                               1 
                             
                           
                           ⁢ 
                           
                             e 
                             
                               i 
                               ⁢ 
                               
                                 π 
                                 2 
                               
                               ⁢ 
                               
                                 Z 
                                 n 
                               
                               ⁢ 
                               
                                 
                                   ∑ 
                                   
                                     k 
                                     ≥ 
                                     2 
                                   
                                 
                                 
                                   
                                     v 
                                     k 
                                   
                                   ( 
                                   
                                     
                                       Z 
                                       
                                         n 
                                         - 
                                         k 
                                       
                                     
                                     + 
                                     
                                       Z 
                                       
                                         n 
                                         + 
                                         k 
                                       
                                     
                                   
                                   ) 
                                 
                               
                             
                           
                         
                       
                     
                   
                 
               
             
             
               
                 
                   = 
                     
                   
                     
                       - 
                       𝒵 
                     
                     ⁢ 
                     
                       
                         ∑ 
                         n 
                       
                       
                         
                           Z 
                           
                             n 
                             - 
                             1 
                           
                         
                         ⁢ 
                         
                           X 
                           n 
                         
                         ⁢ 
                         
                           Z 
                           
                             n 
                             + 
                             1 
                           
                         
                         ⁢ 
                         
                           
                             ∏ 
                             
                               k 
                               ≥ 
                               2 
                             
                           
                           
                             ( 
                             
                               1 
                               - 
                               
                                 
                                   t 
                                   k 
                                   2 
                                 
                                 ⁢ 
                                 
                                   Z 
                                   
                                     n 
                                     - 
                                     k 
                                   
                                 
                                 ⁢ 
                                 
                                   Z 
                                   
                                     n 
                                     + 
                                     k 
                                   
                                 
                               
                               + 
                               
                                 i 
                                 ⁢ 
                                 
                                   t 
                                   k 
                                 
                                 ⁢ 
                                 
                                   
                                     Z 
                                     n 
                                   
                                   ( 
                                   
                                     
                                       Z 
                                       
                                         n 
                                         - 
                                         k 
                                       
                                     
                                     + 
                                     
                                       Z 
                                       
                                         n 
                                         + 
                                         k 
                                       
                                     
                                   
                                   ) 
                                 
                               
                             
                             ) 
                           
                         
                       
                     
                   
                 
               
             
             
               
                 
                   = 
                     
                   
                     
                       - 
                       𝒵 
                     
                     ⁢ 
                     
                       
                         ∑ 
                         n 
                       
                       
                         [ 
                         
                           
                             
                               Z 
                               
                                 n 
                                 - 
                                 1 
                               
                             
                             ⁢ 
                             
                               X 
                               n 
                             
                             ⁢ 
                             
                               Z 
                               
                                 n 
                                 + 
                                 1 
                               
                             
                           
                           + 
                           
                             
                               t 
                               2 
                             
                             ⁢ 
                             
                               Z 
                               
                                 n 
                                 - 
                                 2 
                               
                             
                             ⁢ 
                             
                               Z 
                               
                                 n 
                                 - 
                                 1 
                               
                             
                             ⁢ 
                             
                               Y 
                               n 
                             
                             ⁢ 
                             
                               Z 
                               
                                 n 
                                 + 
                                 1 
                               
                             
                           
                           + 
                           
                             
                               t 
                               2 
                             
                             ⁢ 
                             
                               Z 
                               
                                 n 
                                 - 
                                 1 
                               
                             
                             ⁢ 
                             
                               Y 
                               n 
                             
                             ⁢ 
                             
                               Z 
                               
                                 n 
                                 + 
                                 1 
                               
                             
                             ⁢ 
                             
                               Z 
                               
                                 n 
                                 + 
                                 2 
                               
                             
                           
                           - 
                           
                             
                               t 
                               2 
                               2 
                             
                             ⁢ 
                             
                               Z 
                               
                                 n 
                                 - 
                                 2 
                               
                             
                             ⁢ 
                             
                               Z 
                               
                                 n 
                                 - 
                                 1 
                               
                             
                             ⁢ 
                             
                               X 
                               n 
                             
                             ⁢ 
                             
                               Z 
                               
                                 n 
                                 + 
                                 1 
                               
                             
                             ⁢ 
                             
                               Z 
                               
                                 n 
                                 + 
                                 2 
                               
                             
                           
                           + 
                           
                             O 
                             ⁡ 
                             ( 
                             
                               t 
                               3 
                             
                             ) 
                           
                         
                         ] 
                       
                     
                   
                 
               
             
           
         
       
       
         
           
             	 
             
               Z 
               = 
               
                 
                   ∏ 
                   
                     k 
                     = 
                     2 
                   
                   ∞ 
                 
                 
                   
                     cos 
                     2 
                   
                   ( 
                   
                     π 
                     ⁢ 
                     
                       v 
                       k 
                     
                     / 
                     2 
                   
                   ) 
                 
               
             
           
         
       
     
     The first term that does not commute with    2 ×   2  has prefactor t 3 =tan(πv 3 /2); for v k =1/k 6 , we have t 3 ≈0.002. This is thus an SPT to a very good approximation. One can perform MBQC up to this threshold. In particular, the exact  =4 MPS state is up to the approximation where one neglects t 3  and higher orders. 
     Neglecting t 3  and higher orders, an exact SPT corresponding to the stabilizer is 
     
       
         
           
             
               Z 
               
                 n 
                 - 
                 1 
               
             
             ⁢ 
             
               X 
               n 
             
             ⁢ 
             
               Z 
               
                 n 
                 + 
                 1 
               
             
             ⁢ 
             
               exp 
               ⁡ 
               ( 
               
                 i 
                 ⁢ 
                 
                   π 
                   
                     1 
                     ⁢ 
                     2 
                     ⁢ 
                     8 
                   
                 
                 ⁢ 
                 
                   ( 
                   
                     
                       
                         Z 
                         
                           n 
                           - 
                           2 
                         
                       
                       ⁢ 
                       
                         Z 
                         n 
                       
                     
                     + 
                     
                       
                         Z 
                         n 
                       
                       ⁢ 
                       
                         Z 
                         
                           n 
                           + 
                           2 
                         
                       
                     
                   
                   ) 
                 
               
               ) 
             
           
         
       
     
     One can ask what the fidelity is for the experimentally prepared state with respect to this X=4 SPT state. Approximately 0.999999. State: 
     
       
         
           
             
               A 
               ↑ 
             
             = 
             
               
                 
                   ( 
                   
                     
                       
                         ω 
                       
                       
                         0 
                       
                       
                         0 
                       
                       
                         0 
                       
                     
                     
                       
                         
                           i 
                           ⁢ 
                           
                             ω 
                             ¯ 
                           
                         
                       
                       
                         0 
                       
                       
                         0 
                       
                       
                         0 
                       
                     
                     
                       
                         0 
                       
                       
                         
                           
                             - 
                             i 
                           
                           ⁢ 
                           ω 
                         
                       
                       
                         0 
                       
                       
                         0 
                       
                     
                     
                       
                         0 
                       
                       
                         
                           - 
                           
                             ω 
                             ¯ 
                           
                         
                       
                       
                         0 
                       
                       
                         0 
                       
                     
                   
                   ) 
                 
                 ⁢ 
                     
                 and 
                 ⁢ 
                     
                 
                   A 
                   ↓ 
                 
               
               = 
               
                 
                   
                     ( 
                     
                       
                         
                           0 
                         
                         
                           0 
                         
                         
                           
                             ω 
                             ¯ 
                           
                         
                         
                           0 
                         
                       
                       
                         
                           0 
                         
                         
                           0 
                         
                         
                           
                             i 
                             ⁢ 
                             ω 
                           
                         
                         
                           0 
                         
                       
                       
                         
                           0 
                         
                         
                           0 
                         
                         
                           0 
                         
                         
                           
                             
                               - 
                               i 
                             
                             ⁢ 
                             
                               ω 
                               ¯ 
                             
                           
                         
                       
                       
                         
                           0 
                         
                         
                           0 
                         
                         
                           0 
                         
                         
                           
                             - 
                             ω 
                           
                         
                       
                     
                     ) 
                   
                   ⁢ 
                       
                   with 
                   ⁢ 
                       
                   ω 
                 
                 = 
                 
                   
                     exp 
                     ⁡ 
                     ( 
                     
                       
                         2 
                         ⁢ 
                         π 
                         ⁢ 
                         i 
                       
                       
                         5 
                         ⁢ 
                         1 
                         ⁢ 
                         2 
                       
                     
                     ) 
                   
                   ≈ 
                   
                     0.9999 
                     + 
                     
                       
                         0 
                         . 
                         0 
                       
                       ⁢ 
                       1 
                       ⁢ 
                       2 
                       ⁢ 
                       3 
                       ⁢ 
                       i 
                     
                   
                 
               
             
           
         
       
       
         
           
             	 
             
               
                 Then 
                 ⁢ 
                     
                 
                   
                     〈 
                     
                       χ 
                       = 
                       
                         4 
                         ⁢ 
                         SPT 
                         ⁢ 
                         
                           
                             ❘ 
                             &#34;\[LeftBracketingBar]&#34; 
                           
                           ψ 
                         
                       
                     
                     〉 
                   
                   N 
                 
               
               ≈ 
               
                 
                   0 
                   . 
                   9 
                 
                 ⁢ 
                 9 
                 ⁢ 
                 9 
                 ⁢ 
                 9 
                 ⁢ 
                 9 
                 88. 
               
             
           
         
       
     
     Close to 
     
       
         
           
             
               
                 cos 
                 2 
               
               ( 
               
                 π 
                 
                   
                     2 
                     2 
                   
                   × 
                   
                     3 
                     6 
                   
                 
               
               ) 
             
             ≈ 
             
               1 
               - 
               
                 
                   
                     π 
                     2 
                   
                   
                     
                       2 
                       4 
                     
                     × 
                     
                       3 
                       
                         1 
                         ⁢ 
                         2 
                       
                     
                   
                 
                 . 
               
             
           
         
       
     
     Note: stabilizer and string order can also be measured in two-species experimental set-up. 
     What is the length-scale within which the experimentally-prepared state looks like a genuine SPT state? This is set by the failure of the state being symmetric under e.g. Π n X 2n ; it will mean that string order parameter will eventually decay exponentially to zero, with some length scale ξ SPT . This can be obtained by looking at the largest eigenvalue of the mixed transfer matrix. The absolute value of this eigenvalue was found to be roughly 0.9999953, implying ξ SPT &gt;200000 in units of two-site unit cell, i.e., ξ SPT  is bigger than four hundred thousand lattice sites. In fact, in the approximation where we set v k =0 for k≥5, one can exactly calculate the above eigenvalue to be λ=cos(πv 3 /2) (for v 3 ≤1/2), so for v 3 =1/3 6  one obtains 
     
       
         
           
             
               ξ 
               
                 S 
                 ⁢ 
                 P 
                 ⁢ 
                 T 
               
             
             = 
             
               
                 1 
                 
                   
                     ❘ 
                     &#34;\[LeftBracketingBar]&#34; 
                   
                   
                     ln 
                     ⁢ 
                     
                       cos 
                       ⁡ 
                       ( 
                       
                         π 
                         
                           2 
                           × 
                           
                             3 
                             6 
                           
                         
                       
                       ) 
                     
                   
                   
                     ❘ 
                     &#34;\[RightBracketingBar]&#34; 
                   
                 
               
               ≈ 
               
                 
                   
                     2 
                     3 
                   
                   × 
                   
                     3 
                     
                       1 
                       ⁢ 
                       2 
                     
                   
                 
                 
                   π 
                   2 
                 
               
               ≈ 
               
                 4 
                 ⁢ 
                 3 
                 ⁢ 
                 0 
                 ⁢ 
                 7 
                 ⁢ 
                 7 
                 ⁢ 
                 
                   0 
                   . 
                 
               
             
           
         
       
     
     i.e., roughly half a million sites. 
     String correlation function for the  =4 MPS is 
     
       
         
           
             
               〈 
               
                 
                   Z 
                   
                     
                       2 
                       ⁢ 
                       m 
                     
                     - 
                     1 
                   
                 
                 ⁢ 
                 
                   X 
                   
                     2 
                     ⁢ 
                     m 
                   
                 
                 ⁢ 
                 
                   X 
                   
                     
                       2 
                       ⁢ 
                       m 
                     
                     + 
                     2 
                   
                 
                 ⁢ 
                     
                 … 
                 ⁢ 
                     
                 
                   X 
                   
                     
                       2 
                       ⁢ 
                       n 
                     
                     - 
                     2 
                   
                 
                 ⁢ 
                 
                   X 
                   
                     2 
                     ⁢ 
                     n 
                   
                 
                 ⁢ 
                 
                   Z 
                   
                     
                       2 
                       ⁢ 
                       n 
                     
                     + 
                     1 
                   
                 
               
               〉 
             
             = 
             
               
                 
                   cos 
                   2 
                 
                 ( 
                 
                   
                     π 
                     ⁢ 
                     
                       v 
                       2 
                     
                   
                   2 
                 
                 ) 
               
               = 
               
                 
                   
                     cos 
                     2 
                   
                   ( 
                   
                     π 
                     
                       2 
                       7 
                     
                   
                   ) 
                 
                 ≈ 
                 
                   
                     0 
                     . 
                     9 
                   
                   ⁢ 
                   9 
                   ⁢ 
                   9 
                   ⁢ 
                   4 
                 
               
             
           
         
       
     
     for any m≤n. Including all v k , one expects to see the above string correlator for beyond which it decays exponentially to zero. (In contrast, the trivial string order decays to zero even for short distances.) 
     Kramers-Wannier and General Initial States 
     Referring to  FIG.  28   , the KW MPO is obtained by starting with the MPS of the 1D cluster state flipping the legs on the B (blue) sublattice. Generalized KW dualities can be similarly obtained by a cluster state which is a non-trivial SPT protected by the desired symmetries on a bipartite lattice. 
     The protocol disclosed herein applies to arbitrary initial states, the claim being that by applying the Ising time-evolution discussed in the main text and subsequently measuring the A sublattice will leave the remaining B sites in the Kramers-Wannier transformation of the initial state. (The special case of the initial state being a product state gives many of the examples discussed above, and the special case of the initial state being Z3 top order gives us S3 top order as discussed above.) This general procedure is outlined in a circuit picture in  FIG.  29   . The essential reason for why this works is as follows. First, note that after measuring the A sites, the cluster entangler can be interpreted as creating the cluster state on the red sites (top left of circuit picture) and blue sites (bottom right of circuit picture). The remaining claim is now that this tensor network (by interpreting certain legs as input and the other legs as output) effectively encodes the Kramers-Wannier transformation. For clarify, consider the 1D case, although the explanation readily generalizes to cluster states defined on bipartite lattices in all dimensions. 
     Consider a one-dimensional lattice of 2N qubits. Identify two sublattices A and B corresponding to the odd and even sites of the lattice, respectively. The 1D cluster state can be expressed using a Matrix Product State (MPS) as 
     
       
         
           
             
               
                 
                   
                     
                       
                         
                           ❘ 
                           &#34;\[LeftBracketingBar]&#34; 
                         
                         ψ 
                       
                       〉 
                     
                     = 
                     
                       
                         ∑ 
                         
                           { 
                           s 
                           } 
                         
                       
                       
                         
                           Tr 
                           [ 
                           
                             
                               C 
                               
                                 s 
                                 1 
                               
                             
                             ⁢ 
                             
                               C 
                               
                                 s 
                                 2 
                               
                             
                             ⁢ 
                                 
                             … 
                             ⁢ 
                                 
                             
                               C 
                               
                                 s 
                                 
                                   2 
                                   ⁢ 
                                   N 
                                 
                               
                             
                           
                           ] 
                         
                         ⁢ 
                         
                           
                             ❘ 
                             &#34;\[LeftBracketingBar]&#34; 
                           
                           
                             
                               s 
                               1 
                             
                             , 
                             
                               s 
                               2 
                             
                             , 
                             … 
                                 
                             , 
                             
                               s 
                               
                                 2 
                                 ⁢ 
                                 N 
                               
                             
                           
                         
                       
                     
                   
                   〉 
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   89 
                 
               
             
           
         
       
     
     where s n =0,1 are Z-basis states and the tensor C is defined as 
     
       
         
           
             
               
                 
                   C 
                   = 
                   
                     
                       1 
                       
                         2 
                       
                     
                     ⁢ 
                     
                       ( 
                       
                         
                           
                             
                               
                                 〈 
                                 0 
                               
                               
                                 ❘ 
                                 &#34;\[RightBracketingBar]&#34; 
                               
                             
                           
                           
                             
                               
                                 〈 
                                 0 
                               
                               
                                 ❘ 
                                 &#34;\[RightBracketingBar]&#34; 
                               
                             
                           
                         
                         
                           
                             
                               
                                 〈 
                                 1 
                               
                               
                                 ❘ 
                                 &#34;\[RightBracketingBar]&#34; 
                               
                             
                           
                           
                             
                               
                                 - 
                                 
                                   〈 
                                   1 
                                 
                               
                               
                                 ❘ 
                                 &#34;\[RightBracketingBar]&#34; 
                               
                             
                           
                         
                       
                       ) 
                     
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   90 
                 
               
             
           
         
       
     
     To turn this into a Matrix Product Operator (MPO), first double the unit cell to get an MPS with double the physical legs. 
     
       
         
           
             
               
                 
                   
                     C 
                     ⊗ 
                     C 
                   
                   = 
                   
                     ( 
                     
                       
                         
                           
                             
                               〈 
                               
                                 0 
                                 + 
                               
                             
                             
                               ❘ 
                               &#34;\[RightBracketingBar]&#34; 
                             
                           
                         
                         
                           
                             
                               〈 
                               
                                 1 
                                 - 
                               
                             
                             
                               ❘ 
                               &#34;\[RightBracketingBar]&#34; 
                             
                           
                         
                       
                       
                         
                           
                             
                               〈 
                               
                                 0 
                                 - 
                               
                             
                             
                               ❘ 
                               &#34;\[RightBracketingBar]&#34; 
                             
                           
                         
                         
                           
                             
                               〈 
                               
                                 1 
                                 + 
                               
                             
                             
                               ❘ 
                               &#34;\[RightBracketingBar]&#34; 
                             
                           
                         
                       
                     
                     ) 
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   91 
                 
               
             
           
         
       
     
     Flipping the leg of the first entry upwards yields the MPO 
     
       
         
           
             
               
                 
                   σ 
                   = 
                   
                     ( 
                     
                       
                         
                           
                             
                               
                                 
                                   
                                     ❘ 
                                     &#34;\[LeftBracketingBar]&#34; 
                                   
                                   0 
                                 
                                 〉 
                               
                               ⁢ 
                               
                                 〈 
                                 + 
                               
                             
                             
                               ❘ 
                               &#34;\[RightBracketingBar]&#34; 
                             
                           
                         
                         
                           
                             
                               
                                 
                                   
                                     ❘ 
                                     &#34;\[LeftBracketingBar]&#34; 
                                   
                                   1 
                                 
                                 〉 
                               
                               ⁢ 
                               
                                 〈 
                                 - 
                               
                             
                             
                               ❘ 
                               &#34;\[RightBracketingBar]&#34; 
                             
                           
                         
                       
                       
                         
                           
                             
                               
                                 
                                   
                                     ❘ 
                                     &#34;\[LeftBracketingBar]&#34; 
                                   
                                   0 
                                 
                                 〉 
                               
                               ⁢ 
                               
                                 〈 
                                 - 
                               
                             
                             
                               ❘ 
                               &#34;\[RightBracketingBar]&#34; 
                             
                           
                         
                         
                           
                             
                               
                                 
                                   
                                     ❘ 
                                     &#34;\[LeftBracketingBar]&#34; 
                                   
                                   1 
                                 
                                 〉 
                               
                               ⁢ 
                               
                                 〈 
                                 + 
                               
                             
                             
                               ❘ 
                               &#34;\[RightBracketingBar]&#34; 
                             
                           
                         
                       
                     
                     ) 
                   
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   92 
                 
               
             
           
         
       
     
     This is exactly the Kramers-Wannier duality. For example, plugging in the 1+) product state, one gets the MPS for the GHZ state 
     
       
         
           
             
               
                 
                   ( 
                   
                     
                       
                         
                           
                             〈 
                             0 
                           
                           
                             ❘ 
                             &#34;\[RightBracketingBar]&#34; 
                           
                         
                       
                       
                         0 
                       
                     
                     
                       
                         0 
                       
                       
                         
                           
                             〈 
                             1 
                           
                           
                             ❘ 
                             &#34;\[RightBracketingBar]&#34; 
                           
                         
                       
                     
                   
                   ) 
                 
               
               
                 
                   Equation 
                   ⁢ 
                       
                   93 
                 
               
             
           
         
       
     
     The teachings of all patents, published applications and references cited herein are incorporated by reference in their entirety. 
     While this invention has been particularly shown and described with references to example embodiments thereof, it will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the scope of the invention encompassed by the appended claims.