Non-Abelian fusion rules from an Abelian system

We demonstrate the emergence of non-Abelian fusion rules for excitations of a two dimensional lattice model built out of Abelian degrees of freedom. It can be considered as an extension of the usual toric code model on a two dimensional lattice augmented with matter fields. It consists of the usual C(Zp)C(Zp) gauge degrees of freedom living on the links together with matter degrees of freedom living on the vertices. The matter part is described by a nn dimensional vector space which we call HnHn. The ZpZp gauge particles act on the vertex particles and thus HnHn can be thought of as a C(Zp)C(Zp) module. An exactly solvable model is built with operators acting in this Hilbert space. The vertex excitations for this model are studied and shown to obey non-Abelian fusion rules. We will show this for specific values of nn and pp, though we believe this feature holds for all n>pn>p. We will see that non-Abelian anyons of the quantum double of C(S3)C(S3) are obtained as part of the vertex excitations of the model with n=6n=6 and p=3p=3. Ising anyons are obtained in the model with n=4n=4 and p=2p=2. The n=3n=3 and p=2p=2 case is also worked out as this is the simplest model exhibiting non-Abelian fusion rules. Another common feature shared by these models is that the ground states have a higher symmetry than ZpZp. This makes them possible candidates for realizing quantum computation.