#### Abstract

Doubled topological phases introduced by Kitaev, Levin, and Wen supported on two-dimensional lattices are Hamiltonian versions of three-dimensional topological quantum field theories described by the Turaev-Viro state sum models. We introduce the latter with an emphasis on obtaining them from theories in the continuum. Equivalence of the previous models in the ground state is shown in case of the honeycomb lattice and the gauge group being a finite group by means of the well-known duality transformation between the group algebra and the spin network basis of lattice gauge theory. An analysis of the ribbon operators describing excitations in both types of models and the three-dimensional geometrical interpretation are given.

#### 1. Introduction

Topological quantum field theories (TQFTs) in three dimensions describe a variety of physical and toy models in many areas of modern physics. The absence of local degrees of freedom is a great simplification it often leads to complete solvability [1, 2]. Perhaps the most recent territory, where they appeared to describe real physical systems, is that of topological phases of matter, being, for example, responsible for the fractional quantum Hall effect [3]. Since the idea of fault-tolerant quantum computation appeared in the literature [4], TQFTs are also important in quantum information theory. These new applications also enhanced the mathematical research, and led to classification of the simplest models [5].

Due to their topological nature, TQFTs admit discretization yet remaining an exact description of the theory given by an action functional on a continuous manifold. One large class thereof is the so-called BF theories, whose Lagrangian density is given by the wedge product of a form and the curvature -form of a gauge field [6]. We will deal here with a special class of three-dimensional theories, which describe doubled topological phases and restrict our attention to discrete gauge groups . The context they appeared in, in recent physics literature [4, 7], is Hilbert spaces of states in two dimensions and dynamics therein, which are boundary Hilbert spaces of the relevant TQFTs. Operators acting on correspond to three-dimensional amplitudes on the thickened surface. In this paper we will explain this correspondence, which was proved for the ground-state projection recently [8], and provide the geometric interpretation of the ribbon operators, which create quasiparticle excitations from the ground state. This is a step towards extending the correspondence to identify the ribbon operators as invariants of manifolds with coloured links embedded in them in the TQFT.

The emergence of topological phases from a description of microscopic degrees of freedom is modeled by the lattice models of Kitaev [4] and Levin and Wen [7]. Since they generically have degenerate ground states and quasi-particle excitations insensitive to local disturbances, they are also investigated in the theory of quantum computation [9], their continuum limit being closely related to the spin network simulator [10, 11]. The ground states were extensively studied in the literature; their MERA (multiscale entanglement renormalization ansatz) [12, 13] and tensor network representations [14] have been constructed to study for example, their entanglement properties [15, 16]. Finding the explicit root of these structures in lattice gauge theory and TQFT can help to understand their physical properties.

Lattice gauge theories admit seemingly very different descriptions. A state can be represented by assigning elements of the gauge group to edges of the lattice. The dual description in terms of spin network states where edges are labelled by irreducible representations (irreps) of the gauge group and vertices by invariant intertwiners are also well known since the publication of [17]. To name an application, this description turned out to provide a convenient basis for most approaches to modern quantum gravity theories [18, 19]. In this paper we will show in detail how these dual descriptions give rise to Kitaev’s quantum double models in one hand and the spin net models of Levin and Wen on the other. Then the ribbon operators in both models and their identification will be discussed.

The organization of the paper is as follows. In the next section, we introduce the Turaev-Viro models via the example of BF theories. In Section 3, we briefly introduce the string net models of Levin and Wen on the honeycomb lattice in the surface and recall the proof [8] that the ground-state projection is given by the Turaev-Viro amplitude on . The boundary triangulations of are given by the dual graphs of decorated by the labels inherited from the “initial” and “final” spin nets. In Section 4 the duality between the states of the Kitaev model and the string nets will be shown by changing the basis from the group algebra to the Fourier one. By using this duality and an additional projection, we will obtain the electric constraint operators of the string net models. The matrix elements of the magnetic constraints are also recovered provided that the local rules of Levin and Wen hold. We explain that they do in all BF theories, which is a strong motivation in their favor for the case when the gauge group is finite. In Section 5, we discuss ribbon operators and give their three-dimensional geometric interpretation in terms of framed links in the Turaev-Viro picture. Finally, a summary is given with a list of questions for future research.

#### 2. Turaev-Viro Models

In three dimensions both the (the field strength) and fields of BF theory can be considered to be forms valued in the Lie algebra of the gauge group . The action can then be written as with being an invariant nondegenerate bilinear form on the Lie algebra and is a smooth, oriented, closed three-manifold. We may start from the case when is a semisimple Lie group relevant in particle physics theories and gravity, being the connection in the principal -bundle over . In three dimensions the “space-time” separated form of the Lagrangian has the structure where are spatial indices. (There is not necessarily physical time in the theory; one can do this decomposition for Euclidean signature as well.) The first term is the standard kinetic term, the second implies the (Gauss) constraint of gauge invariance, the third stands for the vanishing of the (two-dimensional) field strength ( and are Lagrange multipliers), while and stand for the covariant and the exterior derivatives, respectively.

Since locally the solution of the constraints is given by a pure gauge ( smooth function, with being the spatial hypersurface), one may discretize the theory by introducing a lattice on the spatial surface and quantize the remaining degrees of freedom: the holonomies (elements of ) describing the coordinate change between faces of the lattice. They correspond to the edges of the dual lattice, which is constructed by placing a vertex inside each face and connecting new vertices, which were put inside neighbouring faces. This dual lattice is the starting point of the models in [4], the electric constraints are the remainders of the Gauss constraint and the magnetic ones are the remainders of the flatness constraint. For a detailed exposition see, for example, [20].

The partition function of the above BF theory is formally obtained by taking the functional integral over the fields an of the phase associated to the classical action It is not so easy to give this definition a precise sense, but for the moment it is not necessary to go into further details. What does matter is that there exists a consistent way to discretize the partition function by considering an oriented triangulation of the manifold by assigning two Lie algebra elements to each edge in . The generator can be thought as the integral of along the edge , whereas as the logarithm of the group element corresponding to the holonomy around the edge . (To be more precise, one needs to introduce the dual complex by putting vertices inside every tetrahedron, connecting those vertices which were put in neighbouring tetrahedra and a vertex should be singled out on the boundary of each dual face. Then the procedure to get is the following: take the dual face corresponding to . Multiply the holonomies along the boundary edges of this dual face starting from the vertex singled out in a circular direction determined by the orientation of (say, by the right-hand rule). The logarithm of this group element is .) Then the Feynman integral in (1) can be replaced by . The integrals will yield Dirac deltas and one can now proceed with decompositions in terms of irreps of the gauge group. This way one ends up with a discrete state sum instead of the original Feynman integral, where each state is the triangulation coloured with irreps and its weight is given by the precise final form of the amplitude (examples are given below). The structure of the partition function (amplitude) for a prototypical theory, the Ponzano-Regge model [21] corresponding to , reads where is the Wigner symbol of depending on the irreps decorating the edges of the tetrahedron , is the dimension of the irreps assigned to the edge , and the sum ranges over all states, that is, all possible colourings of the edges with irreps. It turns out that this type of state sum is well-defined and independent of the chosen triangulation for a large class of models. (This is one way to define a TQFT rigorously.) For a systematic derivation of this state sum from action functionals, see [22] or [23 Section ].

The Ponzano-Regge partition function is formally independent of , but is divergent. However, the Turaev-Viro (TV) model [24], a regularized version thereof, has a well defined partition function, given by where the underlying algebraic structure is the quasitriangular Hopf algebra with fixed, denotes the irreps of (with nonzero trace), is the so-called quantum dimension of , the constant is defined by , the quantity in the brackets is the quantum symbol, and is the number of vertices of the triangulation. One finds the precise definitions of all the quantities along with the algebraic properties assuring consistency and triangulation independence [21] of the amplitudes in [24]. We will briefly mention the origin of the latter property in the next section (Note, that in the case of the Turaev-Viro model based on these properties hold, as in the case of the Ponzano-Regge model, for arbitrary , but it is only for being a root of unity, when the partition function is finite.). The final fact for this introductory section is about the form of the amplitude for manifolds with nonempty boundary. The associated boundary triangulation, whose edges are decorated by labels , derives from a given triangulation in the 3D bulk and is kept fixed. The amplitude reads where is the number of internal vertices, the index ranges over internal edges, over boundary edges, over boundary vertices (each boundary vertex is the endpoint of an internal edge; is its colour), and over all tetrahedra, and the summation is done for internal edge labels only, while those on the boundary are kept fixed.

Note that there is a quasi-triangular Hopf algebra associated to finite groups as well, the so-called Drinfeld (quantum) double [25]. There, the dimensions as well as the symbols can be obtained from the representation theory of the group .

#### 3. String Nets

Levin and Wen [7] started off from the algebraic structure underlying the above models (consistent set of symbols and quantum dimensions), which serves as the algebraic data in defining TQFT’s. Taking these data for granted, they constructed a two-dimensional lattice model, which we will now introduce briefly. Consider a surface with a fixed oriented honeycomb lattice embedded in it. The Hilbert space is spanned by all possible decorations of the edges with labels ; we will refer to them as irreps (of or the finite group ) as we will not need to treat the most general TQFT’s. The Hamiltonian is a sum of two families of mutually commuting constraint operators: where the first sum is over all vertices and the second is over all plaquettes of the lattice. for being the irreps decorating the edges adjacent to the vertex (the numbers are referred to as fusion coefficients between the irreps: ). The magnetic constraints are written as a sum over irreps and the action of the individual terms is while its action on the rest of the state supported on the honeycomb lattice is trivial. The numbers are the symbols, part of the algebraic data of a TQFT; denotes the irreps dual to . Changing the orientation of an edge is equivalent to changing its label to its dual . Levin and Wen use a different normalisation that of , : Before proceeding, let us write down an important algebraic property of the symbols: This identity is called the Biedernharn-Elliot identity or pentagon equation, which holds in every TQFT. In the concrete examples mentioned above, they can be proved by the definition of the symbols as connecting the two different fusion channels of recoupling irreps (graphically encoded by in Section 4.2) (by means of using two different ways of coupling five irreps.)

##### 3.1. Reconstructing 3D Geometry

In our work [8] we recovered a three-dimensional Turaev-Viro invariant [26, 27] from the algebra of Levin and Wen. We associated geometric tetrahedra to the algebraic symbols, where the edges are decorated with irreps from the symbols. In that a convention needs to be adopted; for example, the upper row should correspond to a (triangular) face of the tetrahedron and labels in the same column should correspond to opposite edges. In the examples we are looking at there is always a normalization of the ’s such that they have the same symmetry as the tetrahedron. Orientation of edges can also be taken care of in a consistent manner; we will however omit them for most of what follows. Now we can translate the Biedernharn-Elliot identity to geometry. Then Figure 1 arises: where the two configurations (three tetrahedra joined at the edge and two tetrahedra glued along the triangle ) correspond to the left-hand and right-hand sides of (9), respectively. This is a cornerstone in proving triangulation independence of the amplitudes (3) and (4) and shows that whenever tetrahedra are glued labels corresponding to internal edges have to be summed over.

**(a)**

**(b)**

Now, one can try to find the geometric counterpart of the operator (7). Constructing the dual (triangle) graph of the honeycomb lattice such that a dual edge inherits the label of the original edge it corresponds to (recall that edges of the original and the dual graph are in 1-1 correspondence in two dimensions), we proved the equality [8] The left-hand side. means the matrix element of the operator between two spin nets, that is, the honeycomb lattice decorated by labels and another copy of decorated by . The right-hand side. is the Turaev-Viro amplitude of the three-dimensional manifold with fixed triangulations on the two boundaries given by the dual graph with the labels inherited from . In Figure 2(a) below, the dashed lines show a part of . Let us concentrate on the middle vertex in Figure 2 corresponding to the dual triangle . There is an symbol corresponding to that vertex from all three operators of the three plaquettes sharing that vertex. To each we associate a tetrahedron and they induce the internal triangulation of depending on the order how the operators are multiplied one after the other. These different orders of multiplication correspond to different decompositions of the prism (built from translating the triangle in to the corresponding one in ) into three tetrahedra. The fact that they commute is nicely reflected by the independence of the TV amplitude on the internal triangulation (Note that for proving this equality it is necessary that the coefficients of be given by .) Figure 2(b) shows a part of the amplitude corresponding to . The labels with one bar are to be summed over when the full amplitude is written in accordance with the fact that the underlying edges are the internal edges of the triangulation of . Summation of the labels comes from the sum in .

**(a)**

**(b)**

Note that if we naturally define the symbols to be zero whenever a triple corresponding to a dual triangle of a tetrahedron has . This way, what we found is that the TV amplitude gives the ground-state projection. For the precise matching of all the weights to those of and the consistency of the full amplitude, see [8, the last Section].

#### 4. Duality and the Quantum Double Lattice Models

Now, we will restrict our attention to the case when the TQFT is given by the structure of the double of a finite group and show the equivalence between the lattice models of Kitaev [4] based on that structure and the corresponding string net models. Our method relies on the duality in the underlying lattice gauge theory [17]. Essentially the same idea was employed also in the very recent paper in [28].

The Hilbert space of the Kitaev model (and that of a lattice gauge theory) is spanned by the group algebra basis , supported on an oriented lattice with edges, vertices, and faces. The scalar product is given by The Hamiltonian consists of two families of sums of constraint operators, which are projections and mutually commuting. We will follow the strategy of imposing the electric constraints first and find the corresponding operators in the string net model of Levin and Wen. Then we will study the action of the magnetic constraints in the range of the set of electric constraints, and determine their matrix elements in the dual basis, recovering the magnetic operators in the string net model this way.

The basic idea is the well-known expansion of any function with G being any compact Lie group in terms of irreps ( are the representation matrices and are coefficients). The statement is known as the Peter-Weyl theorem. We now define a new basis by means of the scalar product The () denote the target (source) index of the oriented edge and they range over the dimensions of the irreps . We will need a linear combination of this basis defined in the following way. Consider all elements with fixed irreps . For every vertex of take a three-index tensor , where the indices range over the dimension of the irreps associated to the three edges (the honeycomb lattice is trivalent) incident to . Then contract all indices with the corresponding ones in . A simple example corresponding to the theta-graph is given in Figure 3. For these states (Figure 3)associated to the graph we use the notation

##### 4.1. The Electric Constraints

Let us recall the electric constraints of the Kitaev model. They are written in terms of the following operators:The local gauge transformation acting at vertex reads (the product is over edges incident to the vertex ) and the electric constraint is the projection defined as the average of the latter over the group Note that the range of the set of electric constraints are gauge invariant states, that is, they are invariant under with arbitrary tuple , as shown by the following calculation: Hence, at each vertex, the projection implements gauge invariance. Let us see how this is done in the general set of states also called spin networks. The action of a gauge transformation on a spin network can be determined by rewriting the scalar product as Let us write down this action explicitly for a vertex whose incident edges are oriented outwards and use a simpler notation for a generic state supported on with being the labels of the edges incident to . Let us also use a similar abbreviation for as the remaining parts are not important for the case at hand. Since with denoting the opposite orientation for the edge , we can write In the above the dots stand for the remaining part of the spin network, which is not affected. In the fourth equality the group homomorphism property of the matrices () wass used. The last equality is the definition of as the quantity in the big parenthesis.

The transformation rule of means that . The Clebsch-Gordan series shows that (after a suitable unitary transformation) is block diagonal with and each block transforms according to an irreps of . The case for all corresponds to the trivial representation, which appears in the block decomposition if and only if thus but . We see now that gauge invariance at vertex is achieved by acting with the projection that projects into the invariant subspace of the decomposition. This should correspond to the projection of Levin and Wen. Assuming that for every triple of irreps ensures that there is one unique gauge invariant tensor (for coupling three irreps), an intertwiner, so that the invariant subspace of is spanned by and is understood to be at every vertex contracting all indices of (13). For these states associated to the honeycomb graph and only for these, we will use the notation .

A shorter way to arrive at invariant spin network states is to consider a generic gauge invariant state supported on in the group algebra basis. These are the so-called cylindrical functions with the invariance property for every where () denotes the target (source) vertex of the edge . It can be shown that the spin network states constitute an orthonormal basis in this Hilbert space [29].

##### 4.2. The Magnetic Constraints

To recall the construction of the magnetic operators of the Kitaev model, we define auxiliary operators associated to pairs where is a face (plaquette) of and is an edge on the boundary of : where we have () in the argument of the Dirac delta, when is to the right (left) of the edge oriented forward. The magnetic constraint is the special case of the operator To adapt to the string net model we took to be the honeycomb lattice. After straightforward calculation one finds the action of to be given by whenever all edges bounding the hexagonal face consecutively point to the counterclockwise direction. Should a boundary edge point to the opposite direction, needs to be replaced by in the above expression. To proceed we write down the Plancherel decomposition of the Dirac delta function, which reads Each term in the above sum is the scalar product of a spin network based on a two-valent graph, the hexagon with . This is a spin network of only bivalent vertices (the only intertwiner for bivalent vertex is the trivial Dirac delta connecting identical representations) “evaluated” on the same group elements that appear in the bounding plaquette of the spin network . So we may write the action of , where the state is a generalized spin network with double lines inside the plaquette . The notion, used also in [7], nonetheless, still requires proper definition. Were the use of the local rules of [7] allowed, we could just refer to the calculation given by formula (C1) in that article, which gives the expansion of in terms of bona fide spin networks . We could then just take it as the definition and we would be done. However, to argue in favour of these local rules in the Kitaev model, we need to get back to the theory in the continuum. It has been mentioned that in the case when the group is a Lie group the electric constraints are the lattice versions of the Gauss constraint that imposes local gauge invariance. This was explicitly justified in the previous section. Turning to the flatness constraint, any flat connection has trivial holonomy along a closed curve that is contractible (otherwise we could contract the curve to the point, whose curvature would be proportional to the generator of the holonomy). The converse is also true, to every decoration of with group elements satisfying the constraint (21) for all plaquettes; there exist smooth flat connection(s) in the manifold is embedded into. Suppose that we have constructed one for the embedding surface of the honeycomb lattice. Then a spin network state with any graph makes sense and it is invariant of the homotopy class of the graph . This justifies (23). The connection is flat, so the holonomy along a contractible curve is , Tr, which gives (24) There is no nontrivial intertwiner between two different irreps, whereas the left-hand side. of (25) is a composition of invariant maps with included, so that rule also holds. Finally, we can smoothly contract the edge with label in (26) without changing the value of (27), we have where the middle equality is a property of intertwiners and the rightmost formula coincides with the inner part of the rhs. of (26), when its middle edge with label is contracted. Note that we have omitted also the representation matrices for the irreps as they are not affected by the above, as well as the other parts of the spin networks.

Let us summarize what we have achieved. If we have a Lie group and impose gauge invariance on the honeycomb lattice , the matrix elements of the magnetic operators in the Kitaev model in the spin network basis can be done in two steps. First, one constructs a smooth connection in the manifold in which is embedded. Then one uses the local rules for transforming the spin network (27) as in [7, formula (C1)]. During this process, the group elements also change as we deform the edges, whose holonomies are these group elements, but in the end, we can deform all edges to their original location. This way we find a linear combination of the spin network states corresponding to the magnetic constraint given by the expression (7).

Nevertheless, for finite groups the local rules are, even if well motivated, postulates. The magnetic operator has been derived in a more direct way by introducing some auxiliary degrees of freedom in the very recent paper in [28].

#### 5. Ribbon Operators

In Section 3.1 we have been studying the ground state, the constraints that it stabilizes and the projection from the Hilbert space into the ground state as a three-dimensional TV amplitude. One of the main physical interests, however, is the string-like excitations, the ribbon operators, which correspond to quasiparticles. We are going to sketch the corresponding preliminary results to illustrate that the logic which worked for the ground-state projection, provides us with the three-dimensional interpretation of these quantities as well.

A general ribbon in the spin net model is a string running along a certain path in the honeycomb lattice. The corresponding operator has the following structure: where runs through the vertices of the string, and is the label of the third edge adjacent to the th vertex, which is not part of the string. The label is the “type” of the string. The index structure of the symbols is given by where is the th vertex of the string. The matrices in a string operator have the index structure and generically are matrices (so they have two more indices, which are suppressed above). We would like to proceed as in Section 3 and find a TV amplitude that a string operator describes. Before asking what the matrices correspond to, let us see what geometry we find by passing the description to the dual graph and gluing a tetrahedron whenever there is an symbol.

In Figure 4(a) we have depicted a part of a string, indicating the dual graph along. In the following we will mean this line when referring to the string, and we will mean the collection of dual triangles (shown by green dashed lines in Figure 4(a)) when referring to the ribbon. In Figure 4(b) we took the ribbon and drew a tetrahedron over each triangle it consists of, as dictated by (29). The decoration of the edges follow the index structure of the operator. The edges with the same label belong to the same edge of the spin net, so they are to be glued. This results in the Figure 4(c).

**(a)**

**(b)**

**(c)**

We may interpret the above in the following way. There is the path in the dual graph , which is a continuous line of dual edges of the ribbon that correspond to edges of , which connect vertices with different turning directions of the string (Figures 4(a) and 4(c)). There is an analogous path in . The gluing dictated by the algebraic structure of the operator is such that the line in winds around the one in exactly once during each segment of the line. For each such segment an matrix is present in the form of the operator and the notation ) refers to the initial and the final string nets.

The observables in the TV model are typically ribbon graphs, fat graphs or links embedded in a manifold, over the labels of which, there is no summation in the amplitude [26, 27]. They are invariant under isotopy transformations. This property is ensured by the precise form of the braiding matrices, which then satisfy the Yang-Baxter equation. The latter equation seems to be related to equation of [7] in the spin nets. However, the precise relation and the identification of the braiding matrix in the TV models with an expression of as, for example, the work in [30] suggests should be found for a complete equivalence.

##### 5.1. Kitaev’s Ribbons

The ribbon operators are present in the Kitaev model as well [4].

A prototypical example shown in the Figure 5 is given by a strip between a path along the edges of the original lattice (thick lines) and a neighboring path in the dual lattice (dashed lines). It can be composed of elementary operators associated to triangles, which connect sites, that is, pairs of a plaquette and a vertex on its boundary. In Figure 5 sites are indicated by green dotted lines. One elementary building block is a triangle, which is composed of an edge and a dual vertex (which corresponds to a plaquette ). The other elementary building block is also a triangle composed of a dual edge (which corresponds to the original edge ) and a vertex . The associated operators depend on elements of the double , which can be represented by pairs of group elements . The two types of elementary ribbon operators read Recall that the lattice is assumed to be oriented, so these formulae make sense. The composition of these elementary operators into a long ribbon is done by the comultiplication, which is given by It is desirable to express these ribbons in the spin network basis to recover their corresponding matrix elements in the spin net model. However, there are several obstacles, which should be overcome to accomplish this task. In [7], there are additional local rules to reduce a generalized spin net containing ribbons, to the basis ; see the beginning of Section 4. In order for this to work, one should find a generalized spin network representation of the above operators. Another difficulty comes about when the dual string crosses the original one. In this case, the elementary triangles overlap and the corresponding comultiplication operations do not commute. One needs to find a consistent rule to define their comultiplication in a non-ambiguous way. Note that the simplest ribbon operator in the spin net model, which is the one that winds around one hexagon, is easily found to correspond to . We find the following equality: The procedure to get it is doing the comultiplication for the six elementary operators, all corresponding to the second type in to arrive at . Here the are the group elements corresponding to the edges in the group algebra basis . Then one draws a generalized spin network representation corresponding to the Plancherel decomposition of the Dirac delta as shown in Figure 6 (similarly to those for the magnetic constraints) and resolves it to the spin network basis by using the local moves. It is, however, not straightforward to generalize it (Furthermore, the argument given in Section 4.2 in favour of the local rules is also lost, since the underlying connection here is not flat.)

#### 6. Summary and Outlook

In this paper we have been studying the lattice models of Levin, Wen, and Kitaev from two perspectives. On one hand we identified the ground states and the constraint operators of these models in case the underlying lattice is the honeycomb and the gauge group is a finite group. This has been achieved by changing the basis from that of the group algebra, that is, when edges are decorated by group elements, to the Fourier basis. This basis is spanned by the matrix elements of the irreps. A special linear combination by means of invariant intertwiners at the vertices has been shown to provide the range of all electric constraints and the projection at individual vertices has been identified with the projection to the invariant subspace. Then, the magnetic operators in the group algebra basis have been shown to correspond to those in the spin net model once the local rules postulated in the latter are satisfied. We gave an argument in favour of them from lattice gauge theory with continuous gauge group.

A second focus of the paper was on mapping the spin net to the Turaev-Viro state sum. We have used the idea of building up simplicial manifolds by tetrahedra with edges decorated with irreps corresponding to symbols in the algebraic expressions of operators in the spin net model. This provided the three-dimensional geometric interpretation for the ribbon operator. Having a precise TV amplitude identified with the ribbon operator in the spin net needs further investigation.

One would also like to match these ribbon operators also in the model of Kitaev and the spin net of Levin and Wen. However, finding generalized spin network representations of the previous so that one could reduce them to the spin network basis is not straightforward.

In a series of papers [31–33], families of -deformed “spin network automata” were implemented for processing efficiently classes of computationally-hard problems in geometric topology in particular, approximate calculations of topological invariants of links (collections of knots) and of closed -manifolds. A prominent role was played there by “universal” unitary braiding operators associated with suitable representations of the braid group in the tensor algebra of . Traces of matrices of these representations provide polynomial invariants of -colored links (actually framed links), while weighted sums of the latter give topological invariants of -manifolds presented as complements of framed knots in the -sphere. These invariants are in turn recognized as partition functions and vacuum expectation values of physical observables (Wilson loop operators) in -dimensional Chern-Simons-Witten (CSW) Topological Quantum Field Theory [1]. As is well known (see, e.g., [6], the review in [34], and the original references therein), any 3D TQFT of BF type can be presented as a “double” CSW model, on one hand; and the square modulus of the Witten invariant for a closed oriented 3-manifold equals the TV invariant for the same manifold, on the other.

The remarks above make it manifest that the efficient (approximate) quantum algorithms proposed in [31–33] could be extended in a quite straightforward way to the string-net ground states and ribbon-like excitations framed in the “naturally discretized” double CSW environment given by the TV approach, as we have done in the present paper. Work is in progress in this direction.

#### Acknowledgments

Z. Kádár would like to thank Dirk Schlingemann and Zoltán Zimborás for helpful discussions.