# Creation, manipulation, and detection of Abelian and non-Abelian anyons in optical lattices

## Abstract

Anyons are particle-like excitations of strongly correlated phases of matter with fractional statistics, characterized by nontrivial changes in the wavefunction, generalizing Bose and Fermi statistics, when two of them are interchanged. This can be used to perform quantum computations Kitaev (). We show how to simulate the creation and manipulation of Abelian and non-Abelian anyons in topological lattice models using trapped atoms in optical lattices. Our proposal, feasible with present technology, requires an ancilla particle which can undergo single particle gates, be moved close to each constituent of the lattice and undergo a simple quantum gate, and be detected.

###### pacs:

03.67.Lx, 03.65.Vf, 37.10.JkThe quest for physical systems where anyons WilczekOriginal () can be observed has concentrated so far in effectively 2 materials exhibiting topological order WenOrder (). Abelian anyons, whose interchange generates a nontrivial phase in the wavefunction, exist in the Fractional Quantum Hall effect. Non-Abelian anyons, whose interchange effects full unitary gates on the wavefunction, are expected at certain filling fractions MooreRead () (see recent experimental progress in Dolev ()). In spin lattice systems, anyons can appear as low-lying excitations of topologically ordered ground states (see, e.g., Kitaev (); Levin (); Fendley ()). Several implementations of lattice models with anyonic excitations have been put forward Doucot ()xue ()paredes ()duan ()micheli ()Jiang ()han (). Those involving atoms or molecules in optical lattices are especially attractive, given recent experimental progress BlochRMP (). Specifically, Kitaev’s honeycomb lattice model Kitaev:06 () can be engineered duan (); micheli (), and anyonic interferometry in its Abelian phase can be performed with cavity-mediated global string operations Jiang () or using individual addressing to braid excitations ZhangProcNatl () (but due to the perturbative nature of the effective Hamiltonian in this model, the visibility of anyonic interferometry is degraded jvidal (); Jiang ().)

Here we propose a novel scheme to create topologically ordered states, generate and braid anyons, and detect their statistics for any setup based on particles in optical lattices. We use a lattice of particles of species A to build the topological code and an ancilla of different species B that can be moved independently and brought close to any A particle to perform controlled operations on the code auxspecies (). Preparing the ancilla in superposition states, making it interact with appropriate code particles and measuring its state, the following can be achieved: creation of a topological state, or a general error-correcting code (ECC); creation, braiding, and measurement (fusion) of anyons, all operations needed to perform topological quantum computation (TQC) by braiding; and anyonic interferometry, allowing direct observation of anyonic statistics. Note that (i) by using an ancilla with different quantum states to perform the manipulation of the anyons, tasks can be carried that are not possible using classical (e.g., laser) manipulation of anyons (without single-particle addressability, all proposed methods lack the power of ours); (ii) there is no need in principle of single-particle addressability, specially to perform proof-of-principle experiments; (iii) it is based on successfully demonstrated technologies Porto:07 ()MeschedeSingleatom ()BlochMott (); (iv) it is the first realistic protocol for simulating universal TQC in an atomic, molecular, optical system system [while engineering the microscopic Hamiltonian to build topological protection may be some time off (though see Bacon ()), the method herein works independently of the existence of the background Hamiltonian.]

We consider 2 lattices loaded with atoms or molecules, e.g., Rb. The ancilla, e.g., Na, can be moved independently using a laser potential not affecting Rb atoms (see Fig. 1.) These are now routinely loaded in optical lattices, and in the Mott insulator state one can have extended regions with one particle per site JakschMott (); BlochMott (). Single particles can also be loaded in optical potentials and moved without decoherence MeschedeSingleatom (). Our scheme can be extended to layered 3 configurations Dennis (). We first consider the toric code Hamiltonian Kitaev (), with Abelian anyonic excitations only, as a toy model, but our scheme is basically model independent; later we apply it to the quantum double model Kitaev (), which has non-Abelian anyons and is universal for TQC Mochon ().)

When the ancilla is brought close to a code atom, they experience a 2-qubit unitary between their internal levels (, , are Pauli operators.) This gate can be implemented by cold collisions JakschColdcollisions () or any other means BrennenDipole (); LukinRydberg (). Single qubit operations can be applied to the ancilla without having to address it, due to the different level structure of code atoms. Gates like can be implemented by applying appropriate gates before and after . The internal ancilla state (measurement of ) can also be detected with standard techniques without having to address it or affect the code.

We next show how to create Kitaev’s toric code Kitaev () with these tools. The code is defined as the ground level of a stabilizer Hamiltonian on a square lattice of qubits, realized as Rb atoms, at the edges of a square lattice. The Hamiltonian is the sum of mutually commuting stabilizers and , where v runs over all vertices and p over plaquettes and products involve the qubits surrounding the vertices or plaquettes. We can associate the presence or absence of particle-like excitations at the plaquettes (magnetic defects) and vertices (electric defects) with the fulfillment or not of the ground level conditions and . Plaquette and vertex excitations are thus characterised by and and appear at the ends of strings of and operators applied on a ground state. These ‘particles’ turn out to have nontrivial (anyonic) exchange statistics due to the anticommutation of the and Pauli operators, namely the wavefunction gets multiplied by when a vertex particle winds around a region containing a single plaquette excitation. Detection of this phase change is possible via interference experiments involving superpositions of states with and without anyons. Moreover, the degeneracy of the code allows to interpret it in terms of a set of logical qubits whose and operators are given in terms of chains of ’s and ’s.

We work with a rectangular surface with smooth and rough boundaries Dennis (), with appropriate 3-body vertex and plaquette operators along the boundary providing for a two-dimensional code space: One logical qubit is encoded as the eigenvalue of a chain of ’s along an edge path connecting the rough boundaries. The code space is spanned by coeigenstates of the stabilizers. To create a state in the code, we start with a well defined state ( eigenstate of each ) and measure the stabilizers sequentially, from left to right and top to bottom, using the ancilla. If all outcomes are , our goal is achieved, since VerstraetePEPS (). If is obtained, we can correct by applying to qubit at the bottom of the vertex using again the ancilla, since ; is applied to a qubit that has not been measured yet (for the last row, can be applied to the rightmost qubit.) Once we have measured all vertices and thus prepared the state, we can measure all stabilizers to detect errors and apply error-correcting ’s or ’s to the corresponding qubits by using the ancilla. We could have started at another state, measured all operators in any order, and then corrected errors in this way to prepare the desired state (this can be used to prepare the target state in models beyond Kitaev’s.) We now show how plaquette and vertex measurements, as well as ’s and ’s, are performed using the ancilla. To measure , we prepare the ancilla in state , move it to each qubit in the vertex, and apply each time. Then we apply a Hadamard gate to the ancilla and measure . If the result is , we have applied to the qubits at the vertex, thus performing the desired measurement. is measured by substituting for . To apply () to a qubit, we prepare the ancilla in state , approach it to that qubit and apply (.) Once the toric code state is prepared, operations within the code are performed by applying strings of operators, using the ancilla in state , and applying or sequentially on the desired qubits by bringing them close to the ancilla. To measure string operators, prepare the ancilla in state , follow the same sequence, and measure at the end. The toric code has two kinds of elementary excitations Kitaev (): pairs of anyons in frustrated vertices (electric defects, ) and in frustrated plaquettes (magnetic defects, ), with mutual Abelian anyonic statistics. They can be created by applying or to a given qubit, and can be moved, braided, and fused together by applying these operators along a given path using the ancilla. Superpositions of states with and without vortices, or where they are in different places (see Fig. 2), can be created, allowing the observation of fractional statistics: The simplest interference experiment is shown in Fig. 3. On how to infer anyonic statistics from interference experiments, see the Appendix.

We now outline the preparation and manipulation of anyons in a non-Abelian setting universal for quantum computation Mochon (): the lattice quantum double model based on the group of permutations of 3 elements Kitaev () (see a brief discussion in the Appendix; the full construction is given in longversion ().) It is a lattice model generalizing the toric code, where local degrees of freedom live at the (oriented) edges of the lattice with orthonormal bases labeled by the six group elements . The Hamiltonian, also of the form , has commuting vertex and plaquette stabilizers imposing constraints on the ground states. Their violations define particle-like excitations (anyons) with topological charges (electric, magnetic and dyonic), with non-Abelian fusion and braiding rules. Creation, transport, and fusion of anyons can be achieved generalizing the controlled-NOT operations of the toric code to controlled left and right group multiplications: with , . The local degrees of freedom for the model are qudits of six dimensions, and their six basis elements can be encoded into ground electronic hyperfine states of an alkali atom with enough levels, used as code lattice A. To create, transport, and fuse pure electric and magnetic charge states (enough to simulate universal TQC Mochon ()), a 6-state ancilla species B is especially appropriate (see Fig. 4.) As in the toric code, vertex operators can be measured using ancilla-assisted operations to prepare the ground state (see the Appendix.) With one single ancilla, which need not be spatially addressed, our algorithm requires steps in an region; with an auxiliary lattice with one ancilla per face of the code lattice, assuming addressability, it can be parallelized to depth (essentially optimal Bravyi:06 (); see the Appendix.)

This scheme is independent of the method used to construct the topological state. It can be built by cooling an atomic ensemble interacting via an engineered topological Hamiltonian, providing in principle topological protection to the code except for anyonic manipulations, which should take the system to excited levels in a controlled way (as needed to perform TQC as such; on how to simulate relevant Hamiltonians, see duan (); micheli ().) But it can also be constructed by the above procedure using ancillas to impose stabilizer constraints, enough to perform proof-of-principle interference experiments, a worthy goal by itself. This also allows fault-tolerant quantum computation with general ECCs (topological codes are excellent ECCs: local operators do not mix topological sectors, string operators mixing them can be efficiently implemented.)

The arrangement of ancillas is flexible. One ancilla, individually manipulable MeschedeSingleatom (), suffices in principle. Varying degrees of parallelization are possible: the ground state can be constructed with sequential operations on one column of the sample (parallel to smooth boundaries) at a time, so as to cover the whole sample; in interferometry, commuting operations can be done simultaneously; some parallelization can be introduced with a coarser optical potential for the ancilla than for the code.

This method shares common problems of optical lattice schemes of quantum computation, in particular, spontaneous emission. Essentially, only the vertical direction of the ancilla is close to resonance; lifetimes of seconds can be reached by tuning the laser between the fine structure levels and can be enhanced by restricting manipulation of the vertical dynamics to the (short) times the ancilla is close to a code atom. The ancilla can be repumped after an operation, allowing one to repeat the tasks and detect errors. Controlled logic by cold collisions JakschColdcollisions () requires cooling the system to the physical ground state, possible for both the code (in the 1st Bloch band, the ground state of the local potential) and the ancilla. Rydberg gates BrennenDipole (); LukinRydberg () based on dipole-dipole interactions eliminate this condition. Ancilla-code interactions must break the code in a controlled way, not creating (superpositions of) stray anyons spoiling the quantum memory: Theoretical analysis Calarco () and experimental results MeschedeRegister () suggest that excellent control and small decoherence rates are achievable. Then the implementation benefits from the added protection of topological codes (or, in general, ECCs.) With a bias magnetic field, arbitrary single qudit unitaries can be realized using frequency and polarization selectivity of microwave or Raman laser pulses. Collisional gates can be realized using trap-induced shape resonances Stock:03 (), or using Raman pulses to map code and ancilla ground states ( in figure 4) to a vibrational excited state of each lattice well, evolving by a collisional phase and mapping back Strauch:07 ().

The experimental techniques required by our method are (i) independent trapping of two particle species A and B with different laser trapping potentials; (ii) diluting the population of species B so that each particle is individually addressable; (iii) bringing species A to a Mott insulator phase; (iv) initialisation of species A in a product state ; (v) single particle gates on species B; (vi) simultaneous gates on all particles of species A; (vii) independent transport of internal states of single particles of species B so as to effect cold collisions with particles of species A. Additionally, large scale simulations would require the ability to recool to vibrational states. Each one of these techniques has been demonstrated experimentally; bringing them together will pose an interesting experimental challenge.

A method based on the availability of an ancilla species thus allows one to perform all TQC tasks on a given topological state, independently of the way in which this state is constructed, in optical lattices; e.g., universal TQC based on braiding can be performed on top of the ground state of the model. This scheme is likely to prove the most practical and general way to perform TQC in optical lattices. Large computations will face a steep scaling (a problem not exclusive of topological settings), but observing interference phenomena and applying gates by anyon braiding is feasible with today’s technology. This method can also be used advantageously in general ECCs.

Acknowledgements: Financial support from the EU (SCALA) and DFG (MAP and NIM Excellence Clusters) is gratefully acknowledged. M. A. thanks S. Iblisdir, D. Pérez-García and J. Pachos for discussions on the quantum double models.

## Appendix A Appendix: Supplementary material

In this Appendix we describe in some detail how to simulate non-Abelian anyonic models. Our construction is based on a generalization of the toric code which is spin lattice realization of the quantum double of a discrete group . To realize this we work on a connected two complex which is a cellulation of a two dimensional surface described by vertex, edge, and face sets. Particles with levels (qudits) are placed on the edges and physical states live in a Hilbert space where . As in the toric code, particles on edges that meet at a vertex all interact via a vertex operator . Similarly, all particles on edges that are on the boundary of a face interact via . We pick an orientation for each edge with denoting an edge with arrow pointing from vertex to . The choice of edge orientations is not important as long as a consistent convention is used. The face orientations are all chosen consistent with the surface orientation. Whereas in the toric code model, the ground states are invariant under local gauge transformations generated by , here we demand gauge invariance under transformations:

(1) |

where , the dimensional unitary group, are the permutation representations of the left and right action of multiplication by the group element on the system particle located at edge . For the particle states we make the identification , where by convention , with the identity element. The action of left and right group multiplication on the basis states is then , . A suitable spin lattice model was introduced by Kitaev Kitaev ():

(2) |

where

(3) |

The operator is a projector onto gauge symmetrized states. Similarly, the operator , is a projector onto states with flux at face . Here the sum is taken over all products of group elements acting on a counterclockwise cycle of edges on the boundary of such that the product is . The function according to whether the orientation of the edge is the same as (opposite to) the face orientation. Assigning the origin from which the product is taken is important as failure to do so is equivalent to only specifying the conjugacy class or magnetic charge. Since the trivial flux is its own conjugacy class we write it , which is the projection we seek for ground states of .

By construction . The ground states of are then manifestly gauge invariant states defined as coeigenstates of . Excited states are described by violations of the local constraints and are particle-like corresponding to irreducible representations (irreps) where labels the magnetic charge denoting a conjugacy class of , and labels the electric charge which denotes a unitary irrep of the normalizer of an element in the conjugacy class . Notice that for the group , and , and is precisely the toric code Hamiltonian.

Henceforth, we fix the symmetric group of three objects and let be an square lattice with boundary. We place qudits on edges labeling each by it’s location and vertex and face ancillas labeled by their locations and . For a two complex with boundary the ground state of is unique Bullock:07 () and can be constructed by measuring the vertex projectors in an analogous way to the toric code. We begin with all system particles in state which satisfies the zero flux condition . All vertex and face ancillae are prepared in the state where . We then apply a sequence of operations columnwise beginning on the left column of vertices and ending on the right. For column we apply the controlled operation between each vertex ancilla and its edge neighbors. These operations can be done in parallel since . Next we measure each vertex ancilla in the basis and given the outcome apply the single qudit correction gate , where , on the right edge. Finally when we reach the rightmost column of vertices we apply a sequence of operations from bottom to top involving applying followed by measurement of the ancilla and the correction gate on the top edge. When we reach vertex we are done since the operator is not independent but can be written as a product of the others meaning that already.

Our algorithm has a computational depth of as measured by the number of parallel elementary two qudit gates. We find contructions of the operators or in controlled phase gates and roughly twice that number of local pairwise basis state couplings and suspect this count is optimal. One might wonder if a faster ground state preparation procedure is possible. The answer is no if the initial state is uncorrelated and the available set of operations is quasi-local. The reason is that the final state has global correlations that are created by quasi-local operations. In our algorithm these operations are measurements but they could also be adiabatic turn on of the summands of . The time scale to perform the quasi-local operations (here the measurement time of ) establishes a light cone for the flow of correlations. In Bravyi:06 () it was shown by an application of the Lieb-Robinson bound that the minimal time to prepare a topologically ordered state beginning in a completely uncorrelated state is of the order of the length of the correlations. Since the correlation length scales as the linear dimension of our lattice, our algorithm is essentially optimal.

The group has three conjugacy classes: the identity , transpositions and cyclic permutations , and three irreps: the two one dimensional irreps , the signature representation

and the two-dimensional irrep

For there are irreps: the vacuum state , pure magnetic charges , pure electric charges , and dyonic combinations . A complete derivation of the fusion rules and braid relations for this model is given in Propitius:95 (). We focus here pure electric or pure magnetic charge states, which in fact are sufficient for universal topological quantum computation Mochon ().

Pure electric charges are labeled by basis states where is the projection operator onto a subspace belonging to the unirep . Since we are considering pure electric charges, the unireps of the normalizer of are equivalent to the unireps of itself. In the context of the spin lattice model this is interpreted as electric charge created by applying the projection operator onto the system with the group action being local gauge transformations . Magnetic fluxes are labeled by basis states: and are understood as the result of a projection .

Excitations created in the bulk always come in particle anti-particle pairs. A generic state of a magnetic charge pair is

with . The unique magnetic charge pair state invariant under conjugation by fluxes is the state with , denoted . Note the vacuum state with neighboring magnetic charge pairs can be written

(3) |

We can prepare this by beginning with the ancilla in the state where , with and we have labeled the group elements in ), and applying the two qudit unitary

(4) |

Next we measure the face ancilla in the basis . For the outcome the target magnetic charge state is created. Otherwise for outcome , we need a correction step. To do this prepare the ancilla in the state and apply the controlled operation where

(5) |

which maps the ancilla to state when the flux at evaluated with base point is . Such a controlled operation can be decomposed into elementary two qudit controlled rotation operators with each edge surrounding as a control and the ancilla as the target.

In the toric code and in fact for discrete gauge theories for all finite Abelian groups Bullock:07 (), excitations can be propagated by applying one local operator which simultaneously annihilates a charge at one location and create one at a neighboring location. This is not true for non-Abelian theories as doing so violates a face or vertex constraint. To propagate magnetic charges from one face to an adjacent face essentially involves coherently mapping the value of flux at to the face ancilla using and applying a controlled operation on the shared edge of the faces . After this controlled operation the face ancilla is disentangled from the system by mapping the flux at face to ancilla and performing a controlled operation between ancillae and finally reversing the mapping on . In this protocol we are careful to demand only single qudit and nearest neighbor two qudit interactions. This entire process respects superpositions over flux states and can therefore be used to propagate magnetic charges around the lattice longversion (). Fusion of a magnetic charge pair can be measured by using controlled operations to bring the constituent charges in conjugacy class adjacent to one another at faces and applying the operator followed by measurement of the ancilla in the basis . The probability to obtain the outcome equals the probability for the pair to fuse into the vacuum.

A generic state of an electric charge pair at vertices will be represented as an matrix

with . There is a unique electric charge state which is invariant under conjugation by fluxes: . The state can be prepared as follows. First prepare the vertex ancilla in state . and apply the conditional unitary defined by

(6) |

Depending on the representation, , apply a single qubit operation on the ancilla :

(the latter is not unitary but can be constructed using adaptive measurements as detailed in longversion ()) and finally apply to disentangle the ancilla from the system. One can prepare more distant vacuum charge pairs, e.g. using the same procedure but after each application of , swapping the vertex qudit with the next relevant vertex qudit on the path (via an intermediary swap with the in between edge qudit) then applying the operator at the end, and finally inverting the steps targeting the ancilla. Fusion of electric charge pair is measured by beginning at vertex and applying a sequence of controlled operations acting on a path of vertices from to . At vertex , the operator is applied and then the set of controlled operations along the vertex path is inverted. Finally, the first ancilla is measured in any unitary extension of the basis . The probability to measure the outcome state is the probability for fusion in the vacuum or into the signed irrep.

We now have all the steps to create and braid and fuse anyons. One can verify that all the usual rules for exchanging and braiding anyons are satisfied longversion (). In particular it is easy to see that braiding charges around each other acts trivially since the electric charge pair creation operators are diagonal in the logical basis. Furthermore, the gauge transformations can be viewed as creating a magnetic flux pair, braiding them around the vertex and annihilating. To verify non trivial braiding, one can look for imperfect fusion of electric or magnetic charge pairs into the vacuum. Let’s describe a simplified interference experiment that could be performed in principle on a single face of the lattice. We begin by preparing the state of adjacent electric charge pairs and the vertex ancilla in where for some . Next we apply the controlled operation which creates the state followed by measurement of the ancilla in the basis with outcome . The outcome distribution satisfies

which is precisely the fusion probability for .

It is important to clearly define what we mean by inferring anyonic statistics from interference measurements. For any physical theory, the phase accumulated when braiding one particle around another has contributions from the statistics as well as possible dynamical and geometric (Berry’s phase) contributions Levin:03 (). Considering the interferometry experiment depicted in Fig. 3, including the other phases in the process gives the transformations . If, e.g., there were a background Hamiltonian present, then where is the time to perform the braiding. The statistical phase can be extracted by comparing this measured phase to a proper reference phase. To do so we can perform a similar set of operations but starting in the state such that the created electric defects neighboring A0 are not enclosed by the braiding path of the magnetic defects. Since the braiding paths and particle numbers in each component are the same in both experiments, the difference in the measured phases accumulated on A2 is .

Our scheme is flexible in how the ancillary particles are introduced and one could work entirely within a single two dimensional lattice with in place face and vertex ancilla. For example, the retroreflected square lattice of Porto:07 () provides a bipartite square lattice with addressable systems and ancillary particles. By dynamically changing the trapping fields strengths and phases, the face and vertex ancilla can be moved along any cardinal direction to interact pairwise with the system particles on the edges. However, to build the ground states and manipulate anyons, individual addressability of the ancillary particles using, e.g., the techniques in Cho:07 (); Gorshkov:07 () seems necessary.

### References

- A. Yu. Kitaev, Annals Phys. 303, 2 (2003), arXiv:quant-ph/ 9707021v1.
- F. Wilczek, Phys. Rev. Lett. 48, 1144 (1982); Phys. Rev. Lett. 49, 957 (1982).
- X.-G. Wen, Adv. Phys. 44, 405 (1995), arXiv:cond-mat/ 9506066v2.
- G. Moore and N. Read, Nucl. Phys. B360, 362 (1991).
- M. Dolev et al., Nature 452, 829 (2008), arXiv:0802.0930v1 [cond-mat.mes-hall].
- M. A. Levin and X.-G. Wen, Phys. Rev. B71, 045110 (2005), arXiv:cond-mat/0404617v2 [cond-mat.str-el].
- P. Fendley and E. Fradkin, Phys. Rev. B72, 024412 (2005), arXiv:cond-mat/0502071v2 [cond-mat.str-el].
- B. Douçot, L. B. Ioffe, and J. Vidal, Phys. Rev. B69, 214501 (2004), arXiv:quant-ph/0302104v2.
- Zh.-Y. Xue, Z. D. Wang, and Sh.-L. Zhu, Phys. Rev. A77, 024301 (2008), arXiv:0710.2407v4 [quant-ph].
- B. Paredes and I. Bloch, Phys. Rev. A77, 023603 (2008), arXiv:0711.3796v1 [cond-mat.other].
- L. M. Duan, E. Demler, and M. D. Lukin, Phys. Rev. Lett. 91 090402 (2003), arXiv:cond-mat/0210564v2.
- A. Micheli, G. K. Brennen, and P. Zoller, Nat. Phys. 2, 341 (2006), arXiv:quant-ph/0512222v2.
- L. Jiang et al., Nature Phys. 4, 482 (2008), arXiv: 0711.1365v1 [quant-ph].
- Y. J. Han, R. Raussendorf, and L.-M. Duan, Phys. Rev. Lett. 98 (2007) 150404, arXiv:quant-ph/0702031v1.
- I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 (2008), arXiv:0704.3011v2 [cond-mat.other].
- A. Yu. Kitaev, Annals Phys. 321, 2 (2006), arXiv:cond-mat/ 0506438v3 [cond-mat.mes-hall].
- Ch. Zhang et al., Proc. Natl. Acad. Sci. USA 104 18415 (2007), arXiv:quant-ph/0609101v4.
- S. Dusuel, K. P. Schmidt, and J. Vidal, Phys. Rev. Lett. 100, 177204 (2008), arXiv:0802.0379v1 [cond-mat.other].
- Using auxiliary particles as “moving heads” was already proposed in J. I. Cirac and P. Zoller, Nature 404, 579 (2000); K. G. H. Vollbrecht, E. Solano, and J. I. Cirac, Phys. Rev. Lett. 93, 220502 (2004), arXiv:quant-ph/0405014v2; T. Calarco et al., Phys. Rev. A70, 012306 (2004), arXiv:quant-ph/0403197v1.
- M. Anderlini et al., Nature 448, 452 (2007), arXiv: 0708.2073v1 [quant-ph].
- S. Kuhr et al., Phys. Rev. Lett. 91, 213002 (2003), arXiv: quant-ph/0304081v1.
- M. Greiner et al., Nature 415, 39 (2002).
- D. Bacon, Phys. Rev. A78, 042324 (2008), arXiv: 0806.2160v2 [quant-ph].
- D. Jaksch et al., Phys. Rev. Lett. 81, 3108 (1998), arXiv:cond-mat/9805329v3.
- E. Dennis et al., J. Math. Phys. 43, 4452 (2002), arXiv:quant-ph/0110143v1.
- C. Mochon, Phys. Rev. A69, 032306 (2004), arXiv:quant-ph/0306063v2.
- D. Jaksch et al., Phys. Rev. Lett. 82, 1975 (1999), arXiv:quant-ph/9810087v1.
- G. K. Brennen et al., Phys. Rev. Lett. 82, 1060 (1999), arXiv:quant-ph/9806021v4.
- M. D. Lukin et al., Phys. Rev. Lett. 87, 037901 (2001), arXiv:quant-ph/0011028v1.
- F. Verstraete et al., Phys. Rev. Lett. 96, 220601 (2006), arXiv:quant-ph/0601075v2.
- G. K. Brennen, M. Aguado, and J. I. Cirac, arXiv: 0901.1345v1 [quant-ph].
- S. Bravyi, M. B. Hastings, and F. Verstraete, Phys. Rev. Lett. 97, 050401 (2006), arXiv:quant-ph/0603121v1.
- E. Charron et al., Phys. Rev. A74, 012308 (2006), arXiv: quant-ph/0603138v1.
- D. Schrader et al., Phys. Rev. Lett. 93, 150501 (2004), arXiv:quant-ph/0409037v1.
- R. Stock, E. L. Bolda, and I. H. Deutsch, Phys. Rev. Lett. 91, 183201 (2003), arXiv:quant-ph/0304093v2.
- F. W. Strauch et al., Phys. Rev. A77, 050304(R) (2008), arXiv:0712.1856v1 [quant-ph].
- S. S. Bullock and G. K. Brennen, J. Phys. A: Math. Theor. 40, 3481 (2007), arXiv:quant-ph/0609070v2.
- M. de Wild Propitius, Topological interactions in broken gauge theories, Ph. D. thesis, University of Amsterdam, 1995, arXiv:hep-th/9511195v1.
- M. Levin and X.-G. Wen, Phys. Rev. B67, 245316 (2003), arXiv:cond-mat/0302460v2 [cond-mat.str-el].
- J. Cho, Phys. Rev. Lett. 99, 020502 (2007), arXiv:quant-ph/ 0703185v2.
- A. V. Gorshkov et al., Phys. Rev. Lett. 100, 093005 (2008), arXiv:0706.3879v2 [quant-ph].