# Patent

## Title
Systems and methods for performing quantum computations

## Abstract

Apparatus and methods for performing quantum computations are disclosed. Such apparatus and methods may include identifying a first quantum state of a lattice having a system of quasi-particles disposed thereon, moving the quasi-particles within the lattice according to at least one predefined rule, identifying a second quantum state of the lattice after the quasi-particles have been moved, and determining a computational result based on the second quantum state of the lattice.

## Background
CROSS-REFERENCE TO RELATED APPLICATIONS
This application is a continuation of U.S. patent application Ser. No. 10/909,005, filed Jul. 30, 2004, now U.S. Pat. No. 7,109,593. The subject matter disclosed and claimed herein is related to the subject matter disclosed and claimed in U.S. patent application Ser. Nos. 10/930,640, 10/931,083, and 10/931,082, all of which were filed on Aug. 31, 2004. The disclosure of each of the above-referenced U.S. patent applications is incorporated herein in its entirety.


FIELD OF THE INVENTION
This invention relates in general to the field of quantum computing. More particularly, this invention relates to topological quantum computing.
BACKGROUND OF THE INVENTION
Since the discovery of the fractional quantum Hall effect in 1982, topological phases of electrons have been a subject of great interest. Many abelian topological phases have been discovered in the context of the quantum Hall regime. More recently, high-temperature superconductivity and other complex materials have provided the impetus for further theoretical studies of and experimental searches for abelian topological phases. The types of microscopic models admitting such phases are now better understood. Much less is known about non-abelian topological phases. They are reputed to be obscure and complicated, and there has been little experimental motivation to consider non-abelian topological phases. However, non-abelian topological states would be an attractive milieu for quantum computation.
It has become increasingly clear that if a new generation of computers could be built to exploit quantum mechanical superpositions, enormous technological implications would follow. In particular, solid state physics, chemistry, and medicine would have a powerful new tool, and cryptography also would be revolutionized.
The standard approach to quantum computation is predicated on the quantum bit (“qubit”) model in which one anticipates computing on a local degree of freedom such as a nuclear spin. In a qubit computer, each bit of information is typically encoded in the state of a single particle, such as an electron or photon. This makes the information vulnerable. If a disturbance in the environment changes the state of the particle, the information is lost forever. This is known as decoherence—the loss of the quantum character of the state (i.e., the tendency of the system to become classical). All schemes for controlling decoherence must reach a very demanding and possibly unrealizable accuracy threshold to function.
Topology has been suggested to stabilize quantum information. A topological quantum computer would encode information not in the conventional zeros and ones, but in the configurations of different braids, which are similar to knots but consist of several different threads intertwined around each other. The computer would physically weave braids in space-time, and then nature would take over, carrying out complex calculations very quickly. By encoding information in braids instead of single particles, a topological quantum computer does not require the strenuous isolation of the qubit model and represents a new approach to the problem of decoherence.
In 1997, there were independent proposals by Kitaev and Freedman that quantum computing might be accomplished if the “physical Hilbert space” V of a sufficiently rich TQFT (topological quantum field theory) could be manufactured and manipulated. Hilbert space describes the degrees of freedom in a system. The mathematical construct V would need to be realized as a new and remarkable state for matter and then manipulated at will.
In 2000, Freedman showed that some extraordinarily complicated local Hamiltonian H can be written down whose ground state is V. But this H is an existence theorem only, and is far too complicated to be the starting point for a physical realization.
In 2002, Freedman showed a Hamiltonian involving four-body interactions and stated that after a suitable perturbation, the ground state manifold of H will be the desired state V. This H is less complex than the previously developed H, but it is still only a mathematical construct. One does not see particles, ions, electrons, or any of the prosaic ingredients of the physical world in this prior art model. A Hamiltonian is an energy operator that describes all the possible physical states (eigenstates) of the system and their energy values (eigenvalues).
Freedman further defined the notion of d-isotopy, and showed that if it can be implemented as a ground state of a reasonable Hamiltonian, then this would lead to V and to topological quantum computation. Isotopy is defined as deformation, and two structures that are isotopic are considered to be the same. As shown in the toruses 1 and 2 of FIGS. 1A and 1B, respectively, for example, X and X′ are isotopic, because one may be gradually deformed into the other. In d-isotopy, small circles can be absorbed as a factor=d. Such closed curves are referred to as multicurves or multiloops. Loop X′′ in FIG. 1C (winding around torus 3) is not d-isotopic to X or X′. Loops that are unimportant (because, e.g., they comprise a contractible circle) are called trivial loops and it is desirable to remove, as well as count them. Whenever a trivial loop is removed, the picture is multiplied by “d”. In other words, if two multiloops are identical except for the presence of a contractible circle, then their function values differ by a factor of d, a fixed positive real number. It has been shown that d=2 cos π/(k+2), where k is a level such as 1, 2, 3, etc. which is a natural parameter of Cherns-Simons theory.
According to Freedman, the parameter d can take on only the “special” values: 1, root2, golden ratio, root3 . . . 2 cosπ/(k+2) (where k is a natural number). At d=1, the space V becomes something already known, if not observed in solid state physics. For d>1, V is new to the subject. Freedman, et al., later showed that d-isotopy is explicable by field theory and that multiloops as domain walls can be alternately interpreted as Wilson loop operators. Thus, d-isotopy is a mathematical structure that can be imposed on the multiloops, and is based on Cherns-Simons theory.
An exotic form of matter is a fractional quantum Hall fluid. It arises when electrons at the flat interface of two semiconductors are subjected to a powerful magnetic field and cooled to temperatures close to absolute zero. The electrons on the flat surface form a disorganized liquid sea of electrons, and if some extra electrons are added, quasi-particles called anyons emerge. Unlike electrons or protons, anyons can have a charge that is a fraction of a whole number.
The fractional quantum Hall fluids at one-third filling (of the first Landau level) are already a rudimentary (abelian) example of the V of a TQFT. To effect quantum computation, it would be desirable to construct states more stable and more easily manipulated than FQHE (fractional quantum Hall effect) fluids.
One property of anyons is that when they are moved around each other, they remember in a physical sense the knottedness of the paths they followed, regardless of the path's complexity. It is desirable to use anyons in a system with complex enough transformations, called non-abelian transformations, to carry out calculations in a topological quantum computation system.
In view of the foregoing, there is a need for systems and methods that overcome the limitations and drawbacks of the prior art.
SUMMARY OF THE INVENTION
Apparatus and methods according to the invention may include identifying a first quantum state of a lattice having a quasi-particle disposed thereon, moving the quasi-particle within the lattice according to at least one predefined rule, identifying a second quantum state of the lattice after the quasi-particle has been moved, and determining a computational result based on the second quantum state of the lattice. The quasi-particle may be a non-abelian anyon, for example.
The quasi-particle may be an excitation of a least energy state of a system of real particles. The least energy state and the quasi-particle may be determined by a Hamiltonian operator that is defined from interactions of real particles. The Hamiltonian operator may induce a process on the multi-loops that induces rules for creating, deforming, and annihilating loops. The real particles may define a first dimer covering of an underlying lattice (perhaps the triangular lattice), which, in combination with a second, fixed, background dimer covering, defines one or more multi-loops. The least energy state may be a superposition of the multi-loops, and the quasi-particle may be a canonical excitation of the superposition.
The lattice may be the triangular lattice whose edge center for a Kagome lattice comprises a plurality of hexagons. Each hexagon may contain exactly one real particle. The Kagome lattice may include a plurality of triangular sub-lattices, wherein edges of the triangular sub-lattices are sites of the Kagome lattice. The lattice may include a plurality of lattice sites, where none of the lattice sites hosts more than one dimer.
The predefined rules may include one or more combinatorial moves. The combinatorial moves may include, without limitation, a bow-tie move, a triangle move, and a rhombus flip. Moving the quasi-particle may include moving the quasi-particle relative to a second quasi-particle to cause a quantum braid to be formed in the 2D+1-dimensional space-time of the lattice. The computational result may be based on the quantum braid, which may provide an indication as to how the quasi-particle was moved relative to the second quasi-particle.
Additional features and advantages of the invention will be made apparent from the following detailed description of illustrative embodiments that proceeds with reference to the accompanying drawings.



BRIEF DESCRIPTION OF THE DRAWINGS
The foregoing summary, as well as the following detailed description of preferred embodiments, is better understood when read in conjunction with the appended drawings. For the purpose of illustrating the invention, there is shown in the drawings exemplary constructions of the invention; however, the invention is not limited to the specific methods and instrumentalities disclosed. In the drawings:
FIGS. 1A, 1B, and 1C are diagrams useful in describing isotopy;
FIG. 2 is a diagram of an exemplary Kagome lattice in accordance with the present invention;
FIG. 3 is a diagram of an exemplary Kagome lattice in accordance with the present invention;
FIG. 4 is a diagram of an exemplary lattice that is useful for describing aspects of the present invention;
FIG. 5 is a diagram of an exemplary lattice useful for describing dimer moves in accordance with the present invention; and
FIG. 6 is a block diagram showing an exemplary computing environment in which aspects of the invention may be implemented.



DETAILED DESCRIPTION OF ILLUSTRATIVE EMBODIMENTS
Apparatus and methods according to the invention include identifying a first quantum state of a lattice having a quasi-particle disposed thereon, moving the quasi-particle within the lattice according to at least one predefined rule, identifying a second quantum state of the lattice after the quasi-particle has been moved, and determining a computational result based on the second quantum state of the lattice.
In accordance with the present invention, “realistic” microscopics can be provided for d-isotopy, wherein “realistic” refers to atoms and electrons, for example, as opposed to loops and curves. More particularly, realistic means physical degrees of freedom, with local interaction (the particles are nearby to each other and know of each other's presence). The interactions are potential energy costs of bringing atoms near each other. Existing physical relationships, such as van der waals forces, may be used, along with other physical characteristics, such as tunneling amplitudes. It is desirable to obtain a physical embodiment of the mathematical construct of d-isotopy. In other words, it is desirable to turn the abstract description of the d-isotopy into known physical processes (e.g., tunneling, repulsion, Coulomb interaction).
To go from abstract d-isotopy to real physics, multiloops (i.e., multicurves) are implemented. These multiloops are desirably implemented as broken curves (i.e., “dimer covers”, which is a standard term of art in physics). The following rules are desirably implemented as well: (1) isotopy rule—the amplitude of a picture (or a portion of picture) does not change as it is bent; and (2) d rule—to account for the value associated with a small loop. The rules are turned into “fluctuations”, which are terms in the Hamiltonian that relate to different states (e.g., flickering between two states).
An exemplary embodiment is directed to an extended Hubbard model with a two-dimensional Kagome lattice and a ring-exchange term, described further below. A dimer cover (such as a one-sixth filled Kagome lattice) along with a topology (such as loops) and rules (such as the d-isotopy rule) are desirably comprised within an exemplary extended Hubbard model. Exemplary particles used in conjunction with an exemplary extended model can be bosons or spinless fermions.
At a filling fraction of one-sixth, the model is analyzed in the lowest non-vanishing order of perturbation theory. For a lattice populated with a certain percentage of particles (e.g., electrons), the particles will naturally dissipate to form the “perfect” arrangement. Thus, if one-sixth of the edges of the triangular lattice are filled with particles (e.g., electrons), a perfect matching of the electrons results (e.g., as the result of the Coulomb repulsion).
An exactly soluble point is determined whose ground state manifold is the extensively degenerate “d-isotopy space”, a precondition for a certain type of non-Abelian topological order. Near the values d=2 cos π/(k+2), this space collapses to a stable topological phase with anyonic excitations closely related to SU(2) Chern-Simons theory at level k.
A class of Hamiltonians produces a ground state manifold V=d-isotopy (in some cases “weak d-isotopy”) from standard physical processes. The Hamiltonian is an “extended Hubbard model” in certain cases with an additional “ring exchange term”. A parameter domain is described along with a lattice on which an exemplary model operates. An exemplary lattice is the Kagome lattice. An exemplary model as described herein can serve as a blueprint for the construction of phases of matter Vd, which is useful for quantum information processing. The present invention opens the topological path to a quantum computer.
The extended Hubbard model formalizes the Hamiltonian (kinetic+potential energy) as a tunneling term plus energy costs which are dependent on individual particle location (“on sight potential”) and pair-wise locations (e.g., Coulomb repulsion). For this exemplary model, a ring exchange term may be added, which models collective rotations of large groups of particles (e.g., groups of four particles). According to exemplary embodiments, the cases in which the particles are fermions or bosons are treated separately. The particles can be electrons, Cooper pairs, neutral atoms, phonons, or more exotic “electron fractions” such as chargeons, for example. It is noted that spin is not taken into account in the exemplary embodiments, and so if the particle is specialized to be an electron, either spin must be frozen with a magnetic field or a hyperfine splitting of levels tolerated. It is further noted that a formal transformation may be performed such that these occupation based models may be regarded as purely spin models.
The Kagome lattice is well known in condensed matter physics to support highly “frustrated” spin models with rather mysterious ground states. A Kagome lattice is formed using the centers of the edges of the triangular lattice, and resembles a plurality of hexagons with triangles connecting them. An exemplary Kagome lattice used herein is shown with distinguished sublattices, as shown in FIG. 2.
The Hamiltonian given by:




  
    
      
        H
        =
        
          
            
              ∑
              i
            
            ⁢
            
              
                μ
                t
              
              ⁢
              
                n
                t
              
            
          
          +
          
            
              U
              0
            
            ⁢
            
              
                ∑
                i
              
              ⁢
              
                n
                t
                2
              
            
          
          +
          
            U
            ⁢
            
              
                ∑
                
                  
                    (
                    
                      i
                      ,
                      j
                    
                    )
                  
                  ∈
                  O
                
              
              ⁢
              
                
                  n
                  i
                
                ⁢
                
                  n
                  j
                
              
            
          
          +
          
            
              ∑
              
                
                  
                    (
                    
                      i
                      ,
                      j
                    
                    )
                  
                  ∈
                  
                    ⊳
                    ⊲
                  
                
                ,
                
                  ∉
                  O
                
              
            
            ⁢
            
              
                V
                ij
              
              ⁢
              
                n
                i
              
              ⁢
              
                n
                j
              
            
          
          -
          
            
              ∑
              
                (
                
                  i
                  ,
                  j
                
                )
              
            
            ⁢
            
              
                t
                ij
              
              ⁡
              
                (
                
                  
                    
                      c
                      i
                      †
                    
                    ⁢
                    
                      c
                      j
                    
                  
                  +
                  
                    
                      c
                      j
                      †
                    
                    ⁢
                    
                      c
                      i
                    
                  
                
                )
              
            
          
          +
          
            Ring
            .
          
        
      
    
    
      
        (
        1
        )
      
    
  




is the occupation number on site i, and μi is the corresponding chemical potential. U0 is the usual onsite Hubbard energy U0 (superfluous for spinless fermions). U is a (positive) Coulomb penalty for having two particles on the same hexagon while Vij represents a penalty for two particles occupying the opposite corners of “bow-ties” (in other words, being next-nearest neighbors on one of the straight lines). Allowing for the possibility of inhomogeneity, not all Vij are assumed equal. Specifically, define vcab=Vij, where a is the color of site (i), b is the color of (j), and c is the color of the site between them. In the lattice of FIG. 2, each of the following is possible distinct, vgbb, vbbb, vgbg, vbrb, and vbrg, where r∈R, g∈G, and b∈B=K\(R∪G). tij is the usual nearest-neighbor tunnelling amplitude which is also assumed to depend only on the color of the environment: tij≡tcab where c now refers to the color of the third site in a triangle. “Ring” is a ring exchange term—an additional kinetic energy term which is added to the Hamiltonian on an ad hoc basis to allow correlated multi-particle hops which “shift” particles along some closed paths.

The onsite Hubbard energy U0 is considered to be the biggest energy in the problem, and it is set to infinity, thereby restricting the attention to the low-energy manifold with sites either unoccupied or singly-occupied. The rest of the energies satisfy the following relations: U>>tij, Vij, μi.
Equations may be derived to the second order in perturbation theory for the ground state manifold of this extended Hubbard model. Such equations are provided below. The solutions describe parameter regimes within the Hubbard model for the existence of the phase called d-isotopy. Technically, the dotted sublattice of the Kagome is desirably altered by a system of defects, or alternatively a “ring exchange term” is desirably introduced to achieve d-isotopy; the bare model yields the ground state manifold “weak d-isotopy”, but this distinction is only technical. Both “weak” and “ordinary” behave similarly at the next step.
It is known from the theory of C*-algebras that for “special” d=2 cosπ/(k+2), d-isotopy has a unique symmetry which if broken will relieve the extensive degeneracy of the ground state manifold to a mere finite degeneracy (which depends on topology and boundary conditions). It is further known that once this symmetry is broken, the result is the topological phase Vd. Vd functions as a universal quantum computer. In accordance with the present invention, under a large domain of perturbations, this symmetry will be broken. Thus, the Hamiltonians in accordance with the present invention are blueprints for d-isotopy and hence, after perturbation, for a universal quantum computational system. Such a system has a degree of topological protection from decoherence, which is an enemy of quantum computation. The strength of this protection will depend on the size of the spectral gap above Vd—a quantity which is difficult to compute but can be bounded from the energy scales of any given implantation. There is no a priori basis (such as exists in the FQHE) for asserting that this spectral gap will be small.
Although the equations first solve for an occupation model with Kagome geometry, the invention extends to any physical implementation of the quantum doubles of SU(2) Chern-Simons theory (and their Kauffman-Turaev-Viro) variants which proceed by breaking the symmetry inherent in d-isotopy at its “special” values. Such implementations may be based on lattice models with spin or occupation degrees of freedom or based on field theory.
The non-abelian topological phases which arise are related to the doubled SU(2) Chern-Simons theories described in the prior art. These phases are characterized by (k+1)2-fold ground state degeneracy on the torus T2 and should be viewed as a natural family containing the topological (deconfined) phase of Z2 gauge theory as its initial element, k=1. For k≧2, the excitations are non-abelian. For k=3 and k≧5, the excitations are computationally universal.
The conditions that a microscopic model should satisfy in order to be in such a topological phase are described. It is useful to think of such a microscopic model as a lattice regularization of a continuum model whose low energy Hilbert space may be described as a quantum loop gas. More precisely, a state is defined as a collection of non-intersecting loops. A Hamiltonian acting on such state can do the following: (i) the loops can be continuously deformed—this “move” is referred to as an isotopy move; (ii) a small loop can be created or annihilated—the combined effect of this move and the isotopy move is referred to as ‘d-isotopy’; and (iii) when exactly k+1 strands come together in some local neighborhood, the Hamiltonian can cut them and reconnect the resulting “loose ends” pairwise so that the newly-formed loops are still non-intersecting.
More specifically, in order for this model to be in a topological phase, the ground state of this Hamiltonian should be a superposition of all such pictures with the additional requirements that (i) if two pictures can be continuously deformed into each other, they enter the ground state superposition with the same weight; (ii) the amplitude of a picture with an additional loop is d times that of a picture without such loop; and (iii) this superposition is annihilated by the application of the Jones-Wenzl (JW) projector that acts locally by reconnecting k+1 strands. It should be noted that, as described herein, the particular form of these projectors is highly constrained and leads to a non-trivial Hilbert space only for special values of d=±2 cosπ/(k+2). A Hamiltonian is constructed which enforces d-isotopy for its ground state manifold (GSM).
An exemplary model is defined on the Kagome lattice shown in FIG. 3, which is similar to that shown in FIG. 2. The sites of the lattice are not completely equivalent, and two sublattices are shown, as represented by R (red) and G (green). In FIG. 3, solid dots and dashed lines represent sites and bonds of the Kagome lattice K with the special sublattices R and G. Solid lines define the surrounding triangular lattice.
A dimer cover refers to every vertex meeting exactly one (i.e., one and only one) edge. Dimer covers arise physically (naturally) by repulsion. For example, electrons repel each other to form dimer covers.
As shown, green is a perfect match (as defined above), and red is a second perfect match. Two perfect matches result in multiloops that alternate in color (green, red, green, red, etc.).
It is possible that a green dimer and a red dimer cover the same edge. A solution is to consider it a very short loop of length two (a red and green—assume one of the red or green is slightly displaced so it is a “flat” loop). The multiloops will have alternated these dimers.
Encoding is desirable for an exemplary quantum computing model. The loops are encoded in dimers, and the rules are encoded in particle (e.g., electronic) interactions, such as repulsion or tunneling, for example. In certain situations, when electrons get too close, the system passes from the ground state to the excited state and back to the ground state. This is a virtual process seen in second order perturbation theory and is used to build fluctuations; it generates a deformation.
The “undoped” system corresponds to the filling fraction one-sixth (i.e., Np≡Σi ni=N/6, where N is the number of sites in the lattice). The lowest-energy band then consists of configurations in which there is exactly one particle per hexagon, hence all U-terms are set to zero. These states are easier to visualize if a triangular lattice T is considered whose sites coincide with the centers of hexagons of K, where K is a surrounding lattice for T. Then a particle on K is represented by a dimer on T connecting the centers of two adjacent hexagons of K.
The condition of one particle per hexagon translates into the requirement that no dimers share a site. In the 1⁄6-filled case, this low-energy manifold coincides with the set of all dimer coverings (perfect matchings) of T. The “red” bonds of T (the ones corresponding to the sites of sublattice R) themselves form one such dimer covering, a so-called “staggered configuration”. This particular covering is special: it contains no “flippable plaquettes”, or rhombi with two opposing sides occupied by dimers. See FIG. 4, which shows a triangular lattice T obtained from K by connecting the centers of adjacent hexagons. The bonds corresponding to the special sublattices R and G are shown in dashed and dotted lines, respectively. Triangles with one red side are shaded.
Therefore, particles live on bonds of the triangular lattice and are represented as dimers. In particular, a particle hop corresponds to a dimer “pivoting” by 60 degrees around one of its endpoints. Vij=Vcab is now a potential energy of two parallel dimers on two opposite sides of a rhombus, with c being the color of its short diagonal.
Desirably, the triangular lattice is not bipartite. On the edges of a bipartite lattice, the models will have an additional, undesired, conserved quantity (integral winding numbers, which are inconsistent with the JW projectors for k>2), so the edge of the triangular lattice gives a simple realization.
Because a single tunneling event in D leads to dimer “collisions” (two dimers sharing an endpoint) with energy penalty U, the lowest order at which the tunneling processes contribute to the effective low-energy Hamiltonian is 2. At this order, the tunneling term leads to two-dimer “plaquette flips” as well as renormalization of bare onsite potentials due to dimers pivoting out of their positions and back.
By fixing R as in FIG. 3, without small rhombi with two opposite sides red, as the preferred background dimerization, the fewest equations are obtained along with ergodicity under a small set of moves. Unlike in the usual case, the background dimerization R is not merely a guide for the eyes, it is physically distinguished: the chemical potentials and tunneling amplitudes are different for bonds of different color.
Exemplary elementary dimer moves that preserve the proper dimer covering condition includen plaquette flips, triangle moves, and bow-tie moves. A plaquette (rhombus) flip is a two-dimer move around a rhombus made of two lattice triangles. Depending on whether a “red” bond forms a side of such a rhombus, its diagonal, or is not found there at all, the plaquettes are referred to, respectively, as type 1 (or 1′), 2, or 3 (see the lattice diagram of FIG. 5). FIG. 5 shows an overlap of a dimer covering of T (shown in thick black line) with the red covering shown in dashed line) corresponding to the special sublattice R. Shaded plaquettes correspond to various dimer moves described herein.
The distinction between plaquettes of type 1 and 1′ is purely directional: diagonal bonds in plaquettes of type 1 are horizontal, and for type 1′ they are not. This distinction is desirable because the Hamiltonian breaks the rotational symmetry of a triangular (or Kagome) lattice. A triangle move is a three-dimer move around a triangle made of four elementary triangles. One such “flippable” triangle is labelled 4 in FIG. 5. A bow-tie move is a four-dimer move around a “bow-tie” made of six elementary triangles. One such “flippable” bow tie is labelled 5 in FIG. 5.
To make each of the above moves possible, the actual dimers and unoccupied bonds desirably alternate around a corresponding shape. For both triangle and bow-tie moves, the cases when the maximal possible number of “red” bonds participate in their making (2 and 4 respectively) are depicted. Note that there are no alternating red/black rings of fewer than 8 lattice bonds (occupied by at most 4 non-colliding dimers). Ring moves only occur when red and black dimers alternate; the triangle labelled 4 in FIG. 5 does not have a ring term associated with it, but the bow-tie labelled 5 does.
The correspondence between the previous smooth discussion and rhombus flips relating dimerizations of J of T is now described. The surface is now a planar domain with, possibly, periodic boundary conditions (e.g., a torus). A collection of loops is generated by R∪J, with the convention that the dimers of R∩J be considered as length 2 loops or bigons). Regarding isotopy, move 2 is an isotopy from R∪J to R∪J′ but by itself, it does almost nothing. It is impossible to build up large moves from type 2 alone. So it is a peculiarity of the rhombus flips that there is no good analog of isotopy alone but instead go directly to d-isotopy. The following relations associated with moves of type 5 and 1 (1′) are imposed:

d3Ψ−Ψ=0,  (2a)

dΨ−Ψ=0  (2b)

because from one to four loops in Equation (2a) is passed, and zero to one loop is passed in Equation (2b).

Having stated the goal, the effective Hamiltonian

 H:D→D

    
    
        is a 2×2 matrix corresponding to a dimer move in the two-dimensional basis of dimer configurations connected by this move. ΔIJ=1 if the dimerizations I, J∈D are connected by an allowed move, and ΔIJ=0 otherwise.
    
    


Therefore, it suffices to specify these 2×2 matrices for the off-diagonal processes. For moves of types (1)-(3), they are given below:




  
    
      
        
          
            
              H
              ~
            
            
              (
              1
              )
            
          
          =
          
            
              (
              
                
                  
                    
                      v
                      gb
                      b
                    
                  
                  
                    
                      
                        -
                        2
                      
                      ⁢
                      
                        t
                        rb
                        b
                      
                      ⁢
                      
                        t
                        gb
                        b
                      
                    
                  
                
                
                  
                    
                      
                        -
                        2
                      
                      ⁢
                      
                        t
                        rb
                        b
                      
                      ⁢
                      
                        t
                        gb
                        b
                      
                    
                  
                  
                    
                      v
                      rb
                      b
                    
                  
                
              
              )
            
            =
            
              (
              
                
                  
                    
                      v
                      gb
                      b
                    
                  
                  
                    
                      
                        -
                        2
                      
                      ⁢
                      
                        c
                        0
                      
                      ⁢
                      
                        ε
                        2
                      
                    
                  
                
                
                  
                    
                      
                        -
                        2
                      
                      ⁢
                      
                        c
                        0
                      
                      ⁢
                      
                        ε
                        2
                      
                    
                  
                  
                    
                      v
                      rb
                      b
                    
                  
                
              
              )
            
          
        
        ,
      
    
    
      
        (
        
          3
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          
            
              H
              ~
            
            
              (
              
                1
                ′
              
              )
            
          
          =
          
            
              (
              
                
                  
                    
                      v
                      bb
                      b
                    
                  
                  
                    
                      
                        -
                        2
                      
                      ⁢
                      
                        t
                        rb
                        b
                      
                      ⁢
                      
                        t
                        gb
                        b
                      
                    
                  
                
                
                  
                    
                      
                        -
                        2
                      
                      ⁢
                      
                        t
                        rb
                        b
                      
                      ⁢
                      
                        t
                        gb
                        b
                      
                    
                  
                  
                    
                      v
                      rg
                      b
                    
                  
                
              
              )
            
            =
            
              (
              
                
                  
                    
                      v
                      bb
                      b
                    
                  
                  
                    
                      
                        -
                        2
                      
                      ⁢
                      
                        c
                        0
                      
                      ⁢
                      
                        ε
                        2
                      
                    
                  
                
                
                  
                    
                      
                        -
                        2
                      
                      ⁢
                      
                        c
                        0
                      
                      ⁢
                      
                        ε
                        2
                      
                    
                  
                  
                    
                      v
                      rb
                      b
                    
                  
                
              
              )
            
          
        
        ,
      
    
    
      
        (
        
          3
          ⁢
          b
        
        )
      
    
  
  
    
      
        
          
            
              H
              ~
            
            
              (
              2
              )
            
          
          =
          
            
              (
              
                
                  
                    
                      v
                      bb
                      r
                    
                  
                  
                    
                      
                        -
                        2
                      
                      ⁢
                      
                        
                          (
                          
                            t
                            bb
                            r
                          
                          )
                        
                        2
                      
                    
                  
                
                
                  
                    
                      
                        -
                        2
                      
                      ⁢
                      
                        
                          (
                          
                            t
                            bb
                            r
                          
                          )
                        
                        2
                      
                    
                  
                  
                    
                      v
                      bb
                      r
                    
                  
                
              
              )
            
            =
            
              (
              
                
                  
                    
                      v
                      bb
                      r
                    
                  
                  
                    
                      
                        -
                        2
                      
                      ⁢
                      
                        ε
                        2
                      
                    
                  
                
                
                  
                    
                      
                        -
                        2
                      
                      ⁢
                      
                        ε
                        2
                      
                    
                  
                  
                    
                      v
                      bb
                      r
                    
                  
                
              
              )
            
          
        
        ,
      
    
    
      
        (
        
          3
          ⁢
          c
        
        )
      
    
  
  
    
      
        
          
            H
            ~
          
          
            (
            3
            )
          
        
        =
        
          
            (
            
              
                
                  
                    v
                    bb
                    g
                  
                
                
                  
                    
                      -
                      2
                    
                    ⁢
                    
                      
                        (
                        
                          t
                          bb
                          g
                        
                        )
                      
                      2
                    
                  
                
              
              
                
                  
                    
                      -
                      2
                    
                    ⁢
                    
                      
                        (
                        
                          t
                          bb
                          g
                        
                        )
                      
                      2
                    
                  
                
                
                  
                    v
                    bb
                    g
                  
                
              
            
            )
          
          =
          
            
              (
              
                
                  
                    
                      v
                      bb
                      g
                    
                  
                  
                    0
                  
                
                
                  
                    0
                  
                  
                    
                      v
                      bb
                      g
                    
                  
                
              
              )
            
            .
          
        
      
    
    
      
        (
        
          3
          ⁢
          d
        
        )
      
    
  




    
    
         H can now be tuned to the “small loop” value d.
    
    






  
    
      H
      ~
    
    
      (
      1
      )
    
  
  =
  
    
      
        H
        ~
      
      
        (
        
          1
          ′
        
        )
      
    
    ∝
    
      (
      
        
          
            d
          
          
            
              -
              1
            
          
        
        
          
            
              -
              1
            
          
          
            
              d
              
                -
                1
              
            
          
        
      
      )
    
  




is required because these moves change the number of small loops by one. Because a move of type 2 is an isotopy move,





  
    
      H
      ~
    
    
      (
      2
      )
    
  
  ∝
  
    (
    
      
        
          1
        
        
          
            -
            1
          
        
      
      
        
          
            -
            1
          
        
        
          1
        
      
    
    )
  




    
    
         H(3)=0 provided d>1, because it represents a “surgery” on two strands not allowed for k>1. For k=1, on the other hand,
    
    






  
    
      H
      _
    
    
      (
      3
      )
    
  
  ∝
  
    
      (
      
        
          
            1
          
          
            
              -
              1
            
          
        
        
          
            
              -
              1
            
          
          
            1
          
        
      
      )
    
    .
  




At level k=1 configurations which differ by such a surgery should have equal coefficients in any ground state vector Ψ while at levels k>1 no such relation should be imposed. Thus, for k>1, the matrix relations (3a-3d) yield equations in the model parameters:





  
    
      
        
          
            
              Types
              ⁢
              
                
              
              ⁢
              
                (
                1
                )
              
            
            ⁢
            
              
            
            &
          
          ⁢
          
            
          
          ⁢
          
            (
            
              1
              ′
            
            )
          
          ⁢
          
            :
          
          ⁢
          
            
          
          ⁢
          
            v
            gb
            b
          
        
        =
        
          
            v
            bb
            b
          
          =
          
            2
            ⁢
            d
            ⁢
            
              
            
            ⁢
            
              c
              0
            
            ⁢
            
              ε
              2
            
          
        
      
    
    
      
        (
        
          4
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          
        
        ⁢
        
          
            and
            ⁢
            
              
            
            ⁢
            
              v
              rb
              b
            
          
          =
          
            
              v
              rg
              b
            
            =
            
              2
              ⁢
              
                d
                
                  -
                  1
                
              
              ⁢
              
                c
                0
              
              ⁢
              
                ε
                2
              
            
          
        
      
    
    
      
        (
        
          4
          ⁢
          b
        
        )
      
    
  
  
    
      
        
          
            
              Types
              ⁢
              
                
              
              ⁢
              
                (
                2
                )
              
            
            ⁢
            
              
            
            &
          
          ⁢
          
            
          
          ⁢
          
            (
            3
            )
          
          ⁢
          
            :
          
          ⁢
          
            
          
          ⁢
          
            v
            bb
            r
          
        
        =
        
          
            2
            ⁢
            
              ε
              2
            
            ⁢
            
              
            
            ⁢
            and
            ⁢
            
              
            
            ⁢
            
              v
              bb
              g
            
          
          =
          0
        
      
    
    
      
        (
        
          4
          ⁢
          c
        
        )
      
    
  




    
    
        Suppose that the Hamiltonian has a bare ring exchange term, ring in Eq. (1):
    
    






  
    
      H
      ~
    
    
      (
      Ring
      )
    
  
  =
  
    (
    
      
        
          x
        
        
          
            
              -
              
                c
                3
              
            
            ⁢
            
              ε
              2
            
          
        
      
      
        
          
            
              -
              
                c
                3
              
            
            ⁢
            
              ε
              2
            
          
        
        
          y
        
      
    
    )
  




for some constants c3, x, y>0, and consider the additional equations which come from considering Ring as a fluctuation between one loop of length 8 (type 5 in FIG. 5) and four bigons. It follows from equation (2a) that





  
    
      
        
          
            
              
                H
                ~
              
              
                (
                Ring
                )
              
            
            ∝
            
              (
              
                
                  
                    
                      d
                      3
                    
                  
                  
                    
                      -
                      1
                    
                  
                
                
                  
                    
                      -
                      1
                    
                  
                  
                    
                      d
                      
                        -
                        3
                      
                    
                  
                
              
              )
            
          
          ,
          
            so
            ⁢
            
              :
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            x
            =
            
              
                d
                3
              
              ⁢
              
                c
                3
              
              ⁢
              
                ε
                2
              
            
          
          ;
          
            y
            =
            
              
                d
                
                  -
                  3
                
              
              ⁢
              
                c
                3
              
              ⁢
              
                
                  ε
                  2
                
                .
              
            
          
        
      
    
    
      
        (
        5
        )
      
    
  




are the additional equations (beyond equations (4)) to place the model at the soluble point characterized by d-isotopy. It is clear from equation (5) that the diagonal entries of

 H(Ring)x≠y—perhaps not the most natural choice (although these entries may be influenced by the local chemical environment and thus do not need to be equal). Another possibility exploits the ambiguity of whether a bigon should be considered a loop or not—which allows one to choose x=y.
This construction shows how an extended Hubbard model with the additional ring exchange term (or the equivalent quantum dimer model) can be tuned to have a d-isotopy space as its GSM. Further tuning to its special values then allows the JW projectors to reduce this manifold to the correct ground state corresponding to the unique topological phase associated with this value of d. A simple candidate for a universal quantum computer would be tuned to d=(1+sqrt 5)/2.
Although the exemplary embodiments are described with respect to a triangular lattice, it is contemplated that other lattices, such as a square or hexagonal lattice, may be used in accordance with the present invention. Moreover, the use of irregular lattices is contemplated.
Exemplary Computing Environment
FIG. 6 illustrates an example of a suitable computing system environment 100 in which the invention may be implemented. The computing system environment 100 is only one example of a suitable computing environment and is not intended to suggest any limitation as to the scope of use or functionality of the invention. Neither should the computing environment 100 be interpreted as having any dependency or requirement relating to any one or combination of components illustrated in the exemplary operating environment 100.
The invention is operational with numerous other general purpose or special purpose computing system environments or configurations. Examples of well known computing systems, environments, and/or configurations that may be suitable for use with the invention include, but are not limited to, personal computers, server computers, hand-held or laptop devices, multiprocessor systems, microprocessor-based systems, set top boxes, programmable consumer electronics, network PCs, minicomputers, mainframe computers, distributed computing environments that include any of the above systems or devices, and the like.
The invention may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. The invention may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network or other data transmission medium. In a distributed computing environment, program modules and other data may be located in both local and remote computer storage media including memory storage devices.
With reference to FIG. 6, an exemplary system for implementing the invention includes a general purpose computing device in the form of a computer 110. Components of computer 110 may include, but are not limited to, a processing unit 120, a system memory 130, and a system bus 121 that couples various system components including the system memory to the processing unit 120. The system bus 121 may be any of several types of bus structures including a memory bus or memory controller, a peripheral bus, and a local bus using any of a variety of bus architectures. By way of example, and not limitation, such architectures include Industry Standard Architecture (ISA) bus, Micro Channel Architecture (MCA) bus, Enhanced ISA (EISA) bus, Video Electronics Standards Association (VESA) local bus, and Peripheral Component Interconnect (PCI) bus (also known as Mezzanine bus).
Computer 110 typically includes a variety of computer readable media. Computer readable media can be any available media that can be accessed by computer 110 and includes both volatile and non-volatile media, removable and non-removable media. By way of example, and not limitation, computer readable media may comprise computer storage media and communication media. Computer storage media includes both volatile and non-volatile, removable and non-removable media implemented in any method or technology for storage of information such as computer readable instructions, data structures, program modules or other data. Computer storage media includes, but is not limited to, RAM, ROM, EEPROM, flash memory or other memory technology, CD-ROM, digital versatile disks (DVD) or other optical disk storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium which can be used to store the desired information and which can accessed by computer 110. Communication media typically embodies computer readable instructions, data structures, program modules or other data in a modulated data signal such as a carrier wave or other transport mechanism and includes any information delivery media. The term “modulated data signal” means a signal that has one or more of its characteristics set or changed in such a manner as to encode information in the signal. By way of example, and not limitation, communication media includes wired media such as a wired network or direct-wired connection, and wireless media such as acoustic, RF, infrared and other wireless media. Combinations of any of the above should also be included within the scope of computer readable media.
The system memory 130 includes computer storage media in the form of volatile and/or non-volatile memory such as ROM 131 and RAM 132. A basic input/output system 133 (BIOS), containing the basic routines that help to transfer information between elements within computer 110, such as during start-up, is typically stored in ROM 131. RAM 132 typically contains data and/or program modules that are immediately accessible to and/or presently being operated on by processing unit 120. By way of example, and not limitation, FIG. 6 illustrates operating system 134, application programs 135, other program modules 136, and program data 137.
The computer 110 may also include other removable/non-removable, volatile/non-volatile computer storage media. By way of example only, FIG. 6 illustrates a hard disk drive 140 that reads from or writes to non-removable, non-volatile magnetic media, a magnetic disk drive 151 that reads from or writes to a removable, non-volatile magnetic disk 152, and an optical disk drive 155 that reads from or writes to a removable, non-volatile optical disk 156, such as a CD-ROM or other optical media. Other removable/non-removable, volatile/non-volatile computer storage media that can be used in the exemplary operating environment include, but are not limited to, magnetic tape cassettes, flash memory cards, digital versatile disks, digital video tape, solid state RAM, solid state ROM, and the like. The hard disk drive 141 is typically connected to the system bus 121 through a non-removable memory interface such as interface 140, and magnetic disk drive 151 and optical disk drive 155 are typically connected to the system bus 121 by a removable memory interface, such as interface 150.
The drives and their associated computer storage media provide storage of computer readable instructions, data structures, program modules and other data for the computer 110. In FIG. 6, for example, hard disk drive 141 is illustrated as storing operating system 144, application programs 145, other program modules 146, and program data 147. Note that these components can either be the same as or different from operating system 134, application programs 135, other program modules 136, and program data 137. Operating system 144, application programs 145, other program modules 146, and program data 147 are given different numbers here to illustrate that, at a minimum, they are different copies. A user may enter commands and information into the computer 110 through input devices such as a keyboard 162 and pointing device 161, commonly referred to as a mouse, trackball or touch pad. Other input devices (not shown) may include a microphone, joystick, game pad, satellite dish, scanner, or the like. These and other input devices are often connected to the processing unit 120 through a user input interface 160 that is coupled to the system bus, but may be connected by other interface and bus structures, such as a parallel port, game port or a universal serial bus (USB). A monitor 191 or other type of display device is also connected to the system bus 121 via an interface, such as a video interface 190. In addition to the monitor, computers may also include other peripheral output devices such as speakers 197 and printer 196, which may be connected through an output peripheral interface 195.
The computer 110 may operate in a networked environment using logical connections to one or more remote computers, such as a remote computer 180. The remote computer 180 may be a personal computer, a server, a router, a network PC, a peer device or other common network node, and typically includes many or all of the elements described above relative to the computer 110, although only a memory storage device 181 has been illustrated in FIG. 6. The logical connections depicted include a LAN 171 and a WAN 173, but may also include other networks. Such networking environments are commonplace in offices, enterprise-wide computer networks, intranets and the internet.
When used in a LAN networking environment, the computer 110 is connected to the LAN 171 through a network interface or adapter 170. When used in a WAN networking environment, the computer 110 typically includes a modem 172 or other means for establishing communications over the WAN 173, such as the internet. The modem 172, which may be internal or external, may be connected to the system bus 121 via the user input interface 160, or other appropriate mechanism. In a networked environment, program modules depicted relative to the computer 110, or portions thereof, may be stored in the remote memory storage device. By way of example, and not limitation, FIG. 6 illustrates remote application programs 185 as residing on memory device 181. It will be appreciated that the network connections shown are exemplary and other means of establishing a communications link between the computers may be used.
As mentioned above, while exemplary embodiments of the present invention have been described in connection with various computing devices, the underlying concepts may be applied to any computing device or system.
The various techniques described herein may be implemented in connection with hardware or software or, where appropriate, with a combination of both. Thus, the methods and apparatus of the present invention, or certain aspects or portions thereof, may take the form of program code (i.e., instructions) embodied in tangible media, such as floppy diskettes, CD-ROMs, hard drives, or any other machine-readable storage medium, wherein, when the program code is loaded into and executed by a machine, such as a computer, the machine becomes an apparatus for practicing the invention. In the case of program code execution on programmable computers, the computing device will generally include a processor, a storage medium readable by the processor (including volatile and non-volatile memory and/or storage elements), at least one input device, and at least one output device. The program(s) can be implemented in assembly or machine language, if desired. In any case, the language may be a compiled or interpreted language, and combined with hardware implementations.
The methods and apparatus of the present invention may also be practiced via communications embodied in the form of program code that is transmitted over some transmission medium, such as over electrical wiring or cabling, through fiber optics, or via any other form of transmission, wherein, when the program code is received and loaded into and executed by a machine, such as an EPROM, a gate array, a programmable logic device (PLD), a client computer, or the like, the machine becomes an apparatus for practicing the invention. When implemented on a general-purpose processor, the program code combines with the processor to provide a unique apparatus that operates to invoke the functionality of the present invention. Additionally, any storage techniques used in connection with the present invention may invariably be a combination of hardware and software.
While the present invention has been described in connection with the preferred embodiments of the various figures, it is to be understood that other similar embodiments may be used or modifications and additions may be made to the described embodiments for performing the same functions of the present invention without deviating therefrom. Therefore, the present invention should not be limited to any single embodiment, but rather should be construed in breadth and scope in accordance with the appended claims.

## Claims

1. A method for performing quantum computation, the method comprising:
identifying a first quantum state of a lattice having a system of non-Abelian anyons disposed thereon, the lattice comprising a system of real particles, the non-Abelian anyons being excitations of least energy states of the system of real particles that form the lattice, wherein each of the non-Abelian anyons is associated with a respective initial position in the lattice;
moving the non-Abelian anyons within a 2D+1 space-time of the lattice;
identifying a second quantum state of the lattice after the non-Abelian anyons have been moved; and
determining a computational result based on the first and second quantum states of the lattice.

2. The method of claim 1, wherein the least energy state and the non-Abelian anyons are determined by a Hamiltonian operator that is defined from interactions of real particles.

3. The method of claim 2, wherein the Hamiltonian operator induces a process on the multi-loops that induces rules for creating, deforming, and annihilating loops.

4. The method of claim 1, wherein the real particles define a first dimer covering of an underlying lattice.

5. The method of claim 4, wherein the first dimer covering in combination with a second, fixed, background dimer covering define one or more multi-loops.

6. The method of claim 5, wherein the least energy state is a superposition of the multi-loops.

7. The method of claim 6, wherein the non-Abelian anyons are excitations of the superposition.

8. The method of claim 1, wherein the lattice is a Kagome lattice.

9. The method of claim 8, wherein the Kagome lattice comprises a plurality of hexagons, and wherein each hexagon contains exactly one real particle.

10. The method of claim 8, wherein the Kagome lattice arises as the edge centers of a triangular lattice.

11. The method of claim 1, wherein the lattice has a plurality of lattice sites and, in a ground state, each said lattice site is covered by one dimer, and none of said lattice sites hosts more than one dimer.

12. The method of claim 1, wherein moving the non-Abelian anyons comprises moving a first non-Abelian anyon relative to a second non-Abelian anyon.

13. The method of claim 12, wherein moving the first non-Abelian anyon relative to the second non-Abelian anyon causes a quantum braid to be formed in the space-time of the lattice.

14. The method of claim 13, wherein determining the computational result comprises determining the computational result based on the quantum braid.

15. The method of claim 13, wherein the quantum braid provides an indication as to how the first non-Abelian anyon was moved relative to the second non-Abelian anyon.

16. A method for performing quantum computation, the method comprising:
providing a lattice comprising a system of real particles, the lattice having disposed thereon a system of non-Abelian anyons, the non-Abelian anyons being excitations of least energy states of the system of real particles that form the lattice, wherein each of the non-Abelian anyons is associated with a respective initial position in the lattice;
moving at least one of a plurality of non-Abelian anyons within a 2D+1 space-time of the lattice to form a braid in the space-time of the lattice; and
determining a computational result based on the braid.

17. The method of claim 16, wherein the braid corresponds to a change of quantum state of the lattice.

18. The method of claim 16, further comprising:
moving the non-Abelian anyons within the lattice according to a rule.

19. The method of claim 16, wherein the non-Abelian anyons are determined by a Hamiltonian operator that induces a rule for cutting a first loop and reconnecting loose ends of the cut first loop to form a second loop.

20. The method of claim 16, further comprising:
moving the at least one non-Abelian anyon within the lattice via a Hamiltonian operator that defines a superposition of dimer coverings as its ground state manifold.

21. The method of claim 16, wherein moving the at least one non-Abelian anyons comprises moving the at least one non-Abelian anyons relative to another of the plurality of non-Abelian anyons.