An Algorithm for Tambara-Yamagami Quantum Invariants of 3-Manifolds, Parameterized by the First Betti Number
Abstract
Quantum topology provides various frameworks for defining and computing invariants of manifolds inspired by quantum theory. One such framework of substantial interest in both mathematics and physics is the Turaev-Viro-Barrett-Westbury state sum construction, which uses the data of a spherical fusion category to define topological invariants of triangulated 3-manifolds via tensor network contractions. In this work we analyze the computational complexity of state sum invariants of 3-manifolds derived from Tambara-Yamagami categories. While these categories are the simplest source of state sum invariants beyond finite abelian groups (whose invariants can be computed in polynomial time) their computational complexities are yet to be fully understood. We first establish that the invariants arising from even the smallest Tambara-Yamagami categories are #P-hard to compute, so that one expects the same to be true of the whole family. Our main result is then the existence of a fixed parameter tractable algorithm to compute these 3-manifold invariants, where the parameter is the first Betti number of the 3-manifold with coefficients.
Contrary to other domains of computational topology, such as graphs on surfaces, very few hard problems in 3-manifold topology are known to admit FPT algorithms with a topological parameter. However, such algorithms are of particular interest as their complexity depends only polynomially on the combinatorial representation of the input, regardless of size or combinatorial width. Additionally, in the case of Betti numbers, the parameter itself is computable in polynomial time. Thus while one generally expects quantum invariants to be hard to compute classically, our results suggest that the hardness of computing state sum invariants from Tambara-Yamagami categories arises from classical 3-manifold topology rather than the quantum nature of the algebraic input.
Keywords and phrases:
3-manifold, quantum invariant, fixed parameter tractable algorithm, topological parameter, Gauss sums, topological quantum field theoryFunding:
Colleen Delaney: Partially supported through the Simons Collaboration on Global Categorical Symmetries (Simons Foundation Award 888988).Copyright and License:
![[Uncaptioned image]](x1.png)
2012 ACM Subject Classification:
Mathematics of computing Geometric topology ; Theory of computation Computational geometry ; Theory of computation Fixed parameter tractabilityEditors:
Oswin Aichholzer and Haitao WangSeries and Publisher:

1 Introduction
Quantum physics – especially quantum field theory – is a profuse source of invariants in low-dimensional topology, and these invariants often provide information about manifolds beyond the reach of classical topology. For example, on one hand, the Alexander polynomial of a knot is determined “classically” (from the first homology of the knot’s universal abelian cover) and infinitely many different knots have trivial Alexander polynomial; on the other hand, the Jones polynomial (the prototypical example of a quantum invariant) has never been shown to be trivial on any non-trivial knot. However, the native paradigm for (approximately) computing quantum invariants is manifestly quantum computational, so that in general one should expect them to be hard to compute on a classical computer [19]. Naturally then we seek to understand the full landscape of the computational complexities of quantum invariants. From the point of view of computational topology we are particularly interested in a middle ground where there may be interesting invariants that are not too hard to compute, and to develop algorithms to compute them as efficiently as possible.
In this work, we investigate the complexity landscape of a specific class of quantum invariants of 3-manifolds, namely, the Turaev-Viro-Barrett-Westbury (TVBW) state sum construction [30, 2]. Unlike other quantum invariants (e.g. Donaldson invariants), the TVBW construction is a priori algorithmically computable, combining a finite amount of algebraic data in the form of a (skeletalized) spherical fusion category with a triangulation of a closed, oriented 3-manifold to generate a complex algebraic number that is an orientation-preserving homeomorphism invariant of . The invariant is defined as the contraction of a certain tensor network built by decorating with the data of , and can be formally described via a combinatorial state sum formula. We refer the reader to the full version of this article [8] for a brief review of the TVBW construction (in the case of multiplicity-free, pseudo-unitary spherical categories) and [32] for more background.
As we review at the end of this introduction, it is desirable to understand the computational complexity of evaluating the TVBW invariant on a given triangulated 3-manifold
More precisely, we would like to understand how the computation of these invariants depends on the choice of spherical fusion category , since there are countably infinitely many equivalence classes of such categories, and their classification is at least as hard as the classification of finite groups. As usual, we aim to understand the complexity of computing as a function of the size of a triangulation used to encode .
We approach this complexity-theoretic classification problem from the perspective that fusion categories are “quantum” generalizations of (the representation theory of) finite groups. While this perspective has not been made explicitly overt in the computational topology literature, it is common in quantum algebra, where the extent to which a fusion category deviates from a finite group can be made precise in various ways. For the sake of space, we will not review these notions here, but will note that anyway one cuts it, the simplest spherical fusion categories are of the type – where is an abelian group – consisting of -graded finite-dimensional vector spaces, with a tensor product that comes from the group operation in . In fact, for any finite group , the category of -graded vector spaces forms a spherical fusion category, but if is non-abelian, we consider this category to be more “complicated” than . One reason for this is that is always polynomial-time computable, whereas is often -hard for non-abelian groups [21].
Arguably, after , the next simplest family of fusion categories consists of the Tambara-Yamagami categories . Such a category is determined by a finite abelian group , together with a small amount of additional data: a bicharacter (i.e., a non-degenerate symmetric bilinear pairing) , and a choice of square root [28]. (See [8] for a complete description of these categories using this data.) And yet, in the full version of this article [8], we show that already in the simplest non-trivial case with , the invariants are -hard to compute. This is more-or-less a restatement of known results about the Turaev-Viro invariant [14, 23], since . However, as we shall explain, our approach offers a conceptual clarification by placing within the broader taxonomy of spherical fusion categories. In particular, we conjecture that the TVBW state sum invariant of every non-trivial Tambara-Yamagami category is -hard.
Our main result is a parameterized algorithm to compute the TVBW state sum invariants associated to Tambara-Yamagami categories, where our parameter is the first Betti number .
Main Result (Informal; see Theorem 4).
Fix an integer and a Tambara-Yamagami category . Then there exists a polynomial time algorithm to evaluate for triangulated 3-manifolds with .
While the hardness of general Tambara-Yamagami invariants is still conjectural, we interpret our main result as showing that any such hardness is explainable by a simple classical topological fact: there exist triangulations of 3-manifolds with simplices where can be as large as .
Our algorithm succeeds by converting the natively quantum topological computational problem of computing a TVBW state sum invariant of a triangulated 3-manifold – which a priori requires an expensive tensor network contraction – into a sequence of other much more efficient classical computational problems in algebraic number theory and combinatorial-algebraic topology. Prior to our work, to the best of our knowledge, the only published algorithm in 3-manifold topology directly parameterized by an efficiently computable topological quantity (whose value is independent of the presentation of the input) is the algorithm of Maria and Spreer to compute the Turaev-Viro quantum invariant of 3-manifolds [30] at a 4th root of unity in operations [23].
The technical underpinnings of our algorithm are explained in Section 3; the algorithm itself is precisely described and analyzed in Section 4. At the heart of our algorithm is a (non-obvious!) observation: can be identified with a (normalized) sum of many Gauss sums of -valued quadratic forms on certain finite abelian groups derived from and the given triangulation of .
The second and third equalities above follow for a rather general reason: possesses the structure of a -grading and so by general facts related to homotopy quantum field theory [31], we are able to reduce to a sum of terms in a manner similar in spirit to [23]. Notably, however, the concept of graded fusion categories is absent from [23], and this structure enables us to analyze the state sum invariant of directly from any simplicial triangulation, without having to pass to a one-vertex triangulation via the crushing algorithm of Jaco and Rubinstein [12] (as in [23]).
The bulk of Section 3 is then concerned with explaining and justifying , and our methods here diverge substantially from [23]. The efficient algorithmic evaluation of quadratic Gauss sums appears to be well-known to experts in algebraic number theory, although we were unable to find specific references for this fact, the closest apparently being [3]; see Section 2.1 and [8] for further discussion.
The general nature of our main result allows us to establish new conceptual insight into the complexity-theoretic landscape of TVBW invariants, a portion of which we summarize in Figure 1. It is known that TVBW 3-manifold invariants are often #P-hard to compute exactly, and sometimes hard even just to approximate [13, 14, 18, 19, 20, 21]; moreover, in many of these cases, it is known that there can be no tractable algorithm in . Thus the existence of a parametrized algorithm in a topological and efficiently computable parameter for the Tambara-Yamagami invariants is an interesting counterpoint to their expected hardness. This suggests that the hardness of computing the invariants is due to some intrinsic hardness in 3-manifold topology, rather than the Tambara-Yamagami category alone, whose construction is only a small modification of .
spherical category | complexity of | FPT in ? | Reference |
---|---|---|---|
for abelian | in | Yes | well-known, e.g. [20] |
-hard | Yes | [14, 23], long version [8] | |
unknown in general | Yes | this work | |
for non-abelian solvable | unknown | unknown | |
for non-abelian simple | -complete | No (unless | [20] |
-mod ( a non-lattice root of 1) | -hard | No (unless ) | [13, 19] |
After scrutinizing Figure 1, it seems sensible to expect that there should be a dichotomy theorem for TVBW invariants. A naïve conjecture would be that is in if and only if , and otherwise is -hard. However, this cannot be correct, as there are “twisted” versions of -graded vector spaces , , for which is also known to be polynomial-time computable. More subtly, if is any 2-step nilpotent group, then should also be in .
A systematic understanding of the computational complexity of state sum invariants is significant because they fit into the larger framework of a fully-extended 3D TQFT. Indeed, this connection endows the invariants with important applications to fault-tolerant quantum computation and the characterization of topologically ordered phases of matter, e.g. via topological quantum computing [16, 17, 26, 34], as well as to the developing paradigm of “non-invertible” or “categorical” symmetries and dualities present in quantum field theories and lattice models [1, 6, 9, 10, 29]. Of course, there is plenty of motivation from low-dimensional topology as well, where these invariants can be used in practice to distinguish non-homeomorphic 3-manifolds. TVBW invariants are heuristically strong, with a typical category capable of distinguishing “most” pairs of non-homeomorphic manifolds.
2 Background
In this section we introduce the TVBW invariants of 3-manifolds derived from Tambara-Yamagami categories, as well as some of the basic components necessary for our algorithm.
2.1 Elements of algorithmic number theory
Finite abelian groups.
Let be a finite abelian group with elements, denoted additively. The fundamental theorem of finitely generated abelian groups shows that admits a unique primary decomposition of the form:
for finitely many primes , and positive integers . For a finite abelian group , we denote by the number of non-trivial summands in its primary decomposition. We use to denote the subgroup elements whose order is a power of . Hence,
The following is a standard algorithmic result for finite abelian groups, formulated according to our needs in this article; see for example [25, Chap. 1.11].
Theorem 1.
Let be a finite abelian group presented via a primary decomposition with generators for each of its prime summands. If is a subgroup of specified via a set of generators (each of which is encoded by a list of integers such that ), then there exists a polynomial time algorithm using operations in that finds a presentation of the primary decomposition of – that is, a new minimal set of generators of denoted (each of which is, as before, represented by a list of integers such that ) together with relations of the form for primes and positive integers .
Pairings and forms on finite abelian groups.
We now quickly review the basics of quadratic forms on finite abelian groups. We refer the reader to [27] for more background.
Let be a finite abelian group, denoted additively. A symmetric bilinear pairing into the (not necessarily finite) abelian group is a map satisfying
A pairing is non-degenerate if for every , , there exists some such that is not trivial in . When the symmetric bilinear pairing takes value in , where is the additive group of rational numbers modulo the integers, we call the pair a discriminant pairing. For any real number , denote by the quantity . Let be the multiplicative group of complex numbers on the unit circle. There is a natural group isomorphism . When a symmetric bilinear pairing is non-degenerate and takes value in , we call it a bicharacter. Note that is a non-degenerate discriminant pairing if and only if is a bicharacter. If is a discriminant pairing and , then the Gram matrix of is the matrix .
For a map , define by for any . The map is a (homogeneous) quadratic form if, for any and , and the map is bilinear. For a quadratic form , we call the bilinear pairing associated to . A quadratic form is non-degenerate if the bilinear pairing is non-degenerate. If takes value in , the pair is called a pre-metric group, and if moreover is non-degenerate, is a metric group.
Gauss sums on pre-metric groups.
If , let . To every pre-metric group we associate a complex number called the Gauss sum:
Our next theorem is one of the central tools in our algorithm for computing the TVBW invariants of Tambara-Yamagami categories.
Theorem 2.
Let be a finite abelian group presented via a primary decomposition with summands. If oracle access to a quadratic form is provided, then the Gauss sum can be evaluated in polynomial time using operations, oracle accesses, and memory words.
Theorem 2 is seemingly folklore in the algebraic number theory literature, but we were unable to find any references proving it. We understand it as a consequence of a scattered collection of works on pre-metric groups and algorithms. Wall [33] abstractly classifies metric groups by showing how to decompose them into direct sums of elementary pieces; Basak and Johnson [3] later provide an algorithmic formulation of the decomposition theorem and explain that the Gauss sum of a metric group can be computed readily from the decomposition. We have however not found a complexity analysis of this construction anywhere in the literature. The closest we are aware of is [5, Thm. 1.1], which shows that Gauss sums of (not-necessarily-homogeneous) quadratic forms on cyclic groups are computable in polynomial time; however, even if itself is cyclic, the Gauss sums we will need to compute in our algorithm for are on non-cyclic pre-metric groups of the form . In order to fill this gap, we prove Theorem 2 in the full version of the paper [8] by supplementing the algorithm of [3] with a detailed complexity analysis.
2.2 Combinatorial and algebraic topology
(Co)homology.
Throughout this article we employ the simplicial cohomology of simplicial complexes with coefficients. We assume the reader is familiar with the basic definitions (see [11] for an introduction to the topic). The following well-known fact can be understood as an application of Theorem 1.
Theorem 3.
Let be a simplicial complex with simplices. For any dimension , the following can be computed in polynomial time using operations: a basis of ; a basis of ; a set of elements in that represent a basis of . In particular, the dimension of each of these spaces can be computed in polynomial time.
Triangulated 3-manifolds.
The most common data type in computational 3-manifold topology is that of a generalized triangulation; e.g., see [24]. However, strictly speaking, the TVBW invariants are only defined for simplicial (or PL) triangulations. It is possible to extend the definition to “polytope decompositions” (including generalized triangulations) [15], but we will not need to work at this level of generality.
Thus, throughout this work, we make the following conventions: a 3-manifold is always closed and oriented, and presented by a simplicial triangulation, denoted , together with an orientation of . The orientability of is reflected in the fact that the simplicial homology group is isomorphic to (this can be checked quickly from ), and then an orientation of is encoded by labeling each tetrahedron in with a sign or in a way that gives rise to a -valued cellular 0-cocycle on the cellulation dual to .
Later, it will also be convenient to decorate each edge of an oriented triangulation with an arbitrarily chosen orientation. We stress that these edge orientations are arbitrary, and have nothing to do with the orientation of the manifold’s tetrahedra. Going forward, by a small abuse of language, if we say triangulation, then we will usually mean an oriented triangulation equipped with such an edge orientation. We denote the set of vertices of by , the set of edges by , the set of triangular faces by , and the set of tetrahedra by .
2.3 Tambara-Yamagami quantum invariants
The quantum invariants we consider in this article are -valued topological invariants of closed 3-manifolds derived from algebraic data (a spherical fusion category) attached to a combinatorial presentation of the 3-manifold (a simplicial triangulation), together with evaluation rules described in the works of Turaev and Viro [30] and Barrett and Westbury [2]. We now give an extremely brief introduction to Tambara-Yamagami quantum invariants, which are the TVBW invariants produced by a particular infinite family of spherical fusion categories called Tambara-Yamagami categories [28]. For the sake of space, here we provide only the details necessary to follow the proofs of our main results. We refer the reader to [8] for the full details of the combinatorial description of spherical fusion categories in general and for Tambara-Yamagami categories in particular, as well as the general construction of TVBW invariants (for unitary and multiplicity-free categories).
Let be a finite abelian group with elements, a bicharacter on , and a choice of square root . The Tambara-Yamagami category is an algebraic object which gives rise to the following data:
-
a set of colors (where is some symbol not in ), together with a duality operation ∗ defined by and .
-
scalars and associated to every color in called quantum dimensions, as well as the (global) quantum dimension scalar . These scalars are defined as follows:
-
a set of admissible triples of colors
-
and a system of -valued tetrahedron weights defined below.
A coloring of a triangulation is an assignment of a color to each edge of . Given an edge on the boundary of a face , the orientation of determines a direction along which to traverse the three boundary edges of ; call these edges and , according to the order in which we first see them while doing this traversal. Now associate a triple of colors , where means we apply the duality operation ∗ to the label if the direction of traversal around disagrees with the orientation of . The coloring is called admissible if every such triple associated to an edge-face pair is admissible. The set of admissible colorings is denoted by .
We now describe the tetrahedral weights involved in the definition of the Tambara-Yamagami invariants. Let be a coloring of and let denote the label of an edge . For each tetrahedron with edges and orientation , we assign a tetrahedral weight . If is inadmissible, then . Otherwise, is admissible and so, as described in [8], takes on one of four possible types; taking into account we then assign the tetrahedral weights as follows:
where are arbitrary elements in .
Finally, the Tambara-Yamagami quantum invariant of a 3-manifold with triangulation is defined via the following “state sum formula:”
We note that for a fixed coloring and tetrahedron , testing its admissibility on and determining its type can be carried out in constant time. Thus, there is a naïve algorithm to compute directly from this definition with running time exponential in the size of . In fact, this is true of for any spherical fusion category . We now turn to showing that we can do much better when if we bound the Betti number .
3 The role of the Betti number in Tambara-Yamagami invariants
We are interested in the following algorithmic problem:
-invariant computation |
---|
Input: a triangulation of an oriented closed 3-manifold |
Output: |
Since we are not working uniformly in the data , there is no harm in assuming we have full “explicit” access to it. Our main result is
Theorem 4.
For a fixed Tambara-Yamagami category , there is a deterministic algorithm that, given as input a triangulation of an orientable closed 3-manifold with tetrahedra, computes the quantum invariant in operations and using memory words.
In this section we expound the ideas introduced in Section 1 that lead to this algorithm; its proof is deferred until Section 4. Hereafter and are fixed.
3.1 Admissible colorings induce simplicial 1-cocycles
We begin by noting that each admissible coloring determines a -valued simplicial 1-cocycle. More precisely, consider the function
(1) |
that assigns to an admissible coloring of a -valued simplicial 1-cochain of the triangulation (that is, a function ).
Lemma 5.
For any admissible coloring of a triangulation , is a 1-cocycle.
In the case of , this lemma was proved in [4, 22] and used in the algorithm of [23]. It can be proved easily in the general case by examining the definition of admissible triples in – the calculation boils down to the simple observation that in any admissible coloring, the number of edges around any face of colored by is either 0 or 2. More conceptually, this lemma follows from the fact that is a spherical -graded fusion category, and a similar result holds for any spherical -graded fusion category, cf. [31].
3.2 The partial state sum at a 1-cocycle
Fix a 1-cocycle and consider the partial state sum at , defined by only summing over the admissible labelings associated to :
By definition, all admissible colorings of have the same set of edges colored by the object . We partition , , and , such that for every coloring in :
-
1.
where is the set of edges colored by an element of , and is the set of edges colored by ,
-
2.
where is the set of faces whose edges are colored exclusively by elements of , and is the set of faces with exactly two edges colored by ,
-
3.
, where is the set of -empty faces, is the set of -triangles, and are respectively the set of positively and negatively oriented flat -quads, and and are respectively the set of positively and negatively oriented crooked -quads (see end of Section 2.3 for conventions).
Replacing by the appropriate values and factorizing, we now rewrite the partial state sum at the 1-cocycle .
For an admissible coloring and an -quad tetrahedron in (either flat or crooked), define to be the value where are the colors in assigned to the two opposite edges not colored by in the (flat or crooked) -quad tetrahedron . One can then check:
(2) |
where .
In Section 3.4 we prove:
Theorem 6.
Fix a category . Let be an arbitrary ordering of the edges of . For any collections of pairs of indices , , with potentially repeated pairs, such that , there is an algorithm to evaluate the sum:
(3) |
in operations, using memory words.
Using the theorem, we arrive at a key step towards our main result.
Corollary 7.
Fix a category . Given as input a 3-manifold triangulation with tetrahedra, and a 1-cocycle in , the sum (2) can be evaluated in operations, using memory words.
Proof.
By definition of the projection (Equation (1)), the 1-cocycle gives a description of the edges colored by the object in all admissible colorings of . Additionally, the multiplicative constant in front of the sum (Equation (2)) depends exclusively on the set of edges colored by . In consequence, the multiplicative constant can be evaluated in linear time by checking the set of -colored edges of each tetrahedron. Using the notations of Theorem 6, we consider the pairs (resp. ) of indices of opposite edges not colored by in tetrahedra of type (resp. ), in order to evaluate the sum in cubic time.
3.3 Efficiently characterizing colorings over a 1-cocycle
For a given 1-cocycle , we now characterize the set . Recall that is the set of edges not colored by and is the set of faces with no edge colored by .
Lemma 8.
If the finite abelian group is given by its primary decomposition with summands, and the triangulation has tetrahedra, then the set is in bijection with the elements of a finite abelian group with generators.
Proof.
Fix an arbitrary orientation of the edges of the triangulation. The admissibility constraints for colorings in a Tambara-Yamagami category imply that the set is in bijection with the set of assignments satisfying the following linear equations: for any -empty face with edges ,
where depending whether the orientations of edges agree or not with an arbitrarily chosen orientation of the triangular face. Consider the group homomorphism where the -many values are the relations above. The set is in bijection with the kernel of this application. Because is a finite abelian group with summands in its primary decomposition, with , then is a finite abelian group with generators.
Additionally, a minimal family of generators , for can be computed in polynomial time.
Lemma 9.
For a fixed abelian group , with the notations above, a minimal family of generators for can be computed in operations, using memory words.
Proof.
The proof is standard linear algebra over using Smith normal form, see the proof of [8, Lemma 3.0] for details.
In conclusion, any generator of the kernel can be interpreted as an assignment of colors to the empty edges of the triangulation . For a fixed 1-cocycle of the triangulation, any admissible coloring in can be expressed as a sum (with integer coefficients) of elements of the kernel (where the sum of colors is edge-wise) following the rule of addition in .
3.4 Partial state sums are (normalized) Gauss sums
The goal of this subsection is to prove Theorem 6. Consider a finite abelian group with generators for the summands of its primary decomposition, its product with generators , and a subgroup with .
For a group element , we express , for integer coefficients , in the set of generators for . We also write for the projection of into the -th component of the product .
Let be arbitrary sets of pairs of indices , with , and possible multiplicities for the pair , and for the pair . Consider the product:
(4) |
and define the complex valued sum:
(5) |
Lemma 10.
The sum in Equation (5) is a quadratic Gauss sum over , i.e., there exists a quadratic form s.t.,
Proof.
Define the pairing satisfying, for any , the equality . The map is a bilinear pairing. This follows from the definition of a bicharacter. For any ,
and in consequence . By symmetry of , bilinearity follows. Define by:
such that . The map is a quadratic form. Indeed, first notice that, for any integer , follows from the fact that . Second, the form is bilinear, as it can be expressed as a sum of evaluations of the bilinear pairing on for some .
Proof of Theorem 6.
Using the notation of Theorem 6, we identify the following sum on the triangulation:
(6) |
with the quantity in Equation (5), where the group is the subgroup of the product , whose presentation we identify in polynomial time via Lemma 9. The pairs of indices and represent indices of edges in Equation (6), and particular copies of in the product in Equation (4), with a natural correspondence induced by the fact that copies of in are in bijection with edges not colored by (once the 1-cocycle is fixed). The order of the finite group can be evaluated by finding the primary decomposition of , which can be done in polynomial time (Theorem 1). Let be the quadratic form defined in the proof of Theorem 10. We can compute the Gram matrix of a minimal family of generators of (computed in Lemma 9) for the discriminant pairing by evaluating , and compute the Gauss sum in polynomial time (by virtue of Theorem 2), with the expected complexity.
3.5 Cohomologous 1-cocycles have equal partial state sums
The tools introduced so far in this section would only be enough to construct a parameterized algorithm in terms of the number of 1-cocycles . The final ingredient necessary to further reduce the parameter to follows from a general statement about TVBW invariants of spherical -graded fusion categories (for any finite group ) first described by Turaev and Virelizier [31, Thm. 7.1], specialized to our situation: is a spherical -graded fusion category.
Theorem 11 (special case of Thm. 7.1 in [31]).
If and are cohomologous elements of the group of cocycles , then .
4 Fixed parameter tractable algorithm in the first Betti number
Fix the data of . Given as input an arbitrary triangulation of a closed oriented 3-manifold , our algorithm to compute the TVBW invariant is as follows:
-
1.
Compute a set of 1-cocycles such that forms a basis of the 1-cohomology vector space .
Compute a basis for ; the cardinality of the set of 1-coboundaries is equal to . [Theorem 3]
-
2.
For each of the distinct 1-cohomology classes of , construct a representative by enumerating all -combinations of elements in , i.e., all , for .
Now, for each such , apply Corollary 7 to compute the partial state sum , with the following steps:
-
(a)
Compute a minimal family of generators for , [Lemma 9]
-
(b)
Compute the primary decomposition of with generators for the summands, [Theorem 1]
-
(c)
Compute the Gram matrix of for the pre-metric form , by evaluating on all pairs . [defined in Lemma 10]
-
(d)
Evaluate the Gauss sum of , [Theorem 2]
and multiply by . [proof of Theorem 6]
-
(a)
-
3.
Sum over all :
Our main Theorem 4 follows readily from the following:
Theorem 12.
Fix a Tambara-Yamagami category . Given a triangulation of a closed 3-manifold with tetrahedra, the above algorithm computes in operations and using memory words.
Proof.
Computing a basis for , as well as a set of generators whose cohomology classes form a basis of , in step (1) is a standard procedure in computer algebra (Theorem 3), done by normalizing a matrix with coefficients. This can be done in operations in with Gaussian eliminations, where an arithmetic operation in has constant complexity.
Given the 1-cocycles , we can enumerate every combination in (amortized) operations per . Indeed, enumerating all is akin to incrementing a -bits binary counter from to , which induces an (amortized) change of bits per increment [7, Chapter 16]. The 1-cocycles are represented by formal sums of edges, and computing the sum of two 1-cocycles takes operations.
For each as above, we compute the partial state sum using Corollary 7 in polynomial time. Specifically, step (2a) is solved by computing a Smith normal form, together with transformation matrices, of a particular -integer matrix, which has complexity (Lemma 9) with (modular) Gaussian eliminations, once is fixed.
We can read-off the generators of from the transformation matrices, and compute the generators of step (2b) by computing the primary decomposition of the subgroup (Theorem 1). Computing the Gram matrix in step (2c) is done by evaluating the bilinear pairing , defined in Lemma 10, on all pairs . The quadratic form is defined as a sum of terms of the form , and the overall Gram matrix can be computed in operations. Recall that we have an oracle access to the function (see beginning of Section 3). By Theorem 2, the Gauss sum of can be evaluated from the Gram matrix of in operations in , where the arithmetic complexity of operations in depend solely on the maximal order of an element of (independently of the size of the Gram matrix). See the full version [8] for details.
Finally, the invariant computed in step (3) is the sum of complex numbers represented with a constant number of memory words each for a fixed group (with ). In consequence, the overall complexity of the algorithm, knowing and considering constant for a fixed group , is . The memory complexity is bounded by the size of the matrices, which is , since we have a fixed (coefficients are represented with a constant number of memory words, and is a constant).
References
- [1] David Aasen, Roger S. K. Mong, and Paul Fendley. Topological defects on the lattice: I. The Ising model. J. Phys. A, 49(35):354001, 46, 2016. doi:10.1088/1751-8113/49/35/354001.
- [2] John W. Barrett and Bruce W. Westbury. Invariants of piecewise-linear -manifolds. Trans. Amer. Math. Soc., 348(10):3997–4022, 1996. doi:10.1090/S0002-9947-96-01660-1.
- [3] Tathagata Basak and Ryan Johnson. Indicators of Tambara-Yamagami categories and Gauss sums. Algebra & Number theory, 9(8):1793–1823, 2015.
- [4] Benjamin A. Burton, Clément Maria, and Jonathan Spreer. Algorithms and complexity for Turaev-Viro invariants. J. Appl. Comput. Topol., 2(1-2):33–53, 2018. doi:10.1007/S41468-018-0016-2.
- [5] Jin-Yi Cai, Xi Chen, Richard Lipton, and Pinyan Lu. On tractable exponential sums. In Frontiers in Algorithmics: 4th International Workshop, FAW 2010, Wuhan, China, August 11-13, 2010. Proceedings 4, pages 148–159. Springer, 2010. doi:10.1007/978-3-642-14553-7_16.
- [6] Yichul Choi, Ho Tat Lam, and Shu-Heng Shao. Noninvertible Global Symmetries in the Standard Model. Phys. Rev. Lett., 129:161601, October 2022. doi:10.1103/PhysRevLett.129.161601.
- [7] Thomas H. Cormen, Clifford Stein, Ronald L. Rivest, and Charles E. Leiserson. Introduction to Algorithms. McGraw-Hill Higher Education, 4th edition, 2001.
- [8] Colleen Delaney, Clément Maria, and Eric Samperton. An algorithm for Tambara-Yamagami quantum invariants of 3-manifolds, parameterized by the first Betti number. CoRR, abs/2311.08514, 2023. doi:10.48550/arXiv.2311.08514.
- [9] Daniel S. Freed, Gregory W. Moore, and Constantin Teleman. Topological symmetry in quantum field theory, 2023. arXiv:2209.07471.
- [10] Daniel S. Freed and Constantin Teleman. Topological dualities in the Ising model. Geom. Topol., 26(5):1907–1984, 2022. doi:10.2140/gt.2022.26.1907.
- [11] Allan Hatcher. Algebraic Topology. Cambridge University Press, 1 edition, 2002.
- [12] William Jaco and J. Hyam Rubinstein. 0-efficient triangulations of 3-manifolds. J. Differential Geom., 65(1):61–168, 2003.
- [13] F. Jaeger, D. L. Vertigan, and D. J. A. Welsh. On the computational complexity of the Jones and Tutte polynomials. Math. Proc. Cambridge Philos. Soc., 108(1):35–53, 1990. doi:10.1017/S0305004100068936.
- [14] Robion Kirby and Paul Melvin. Local surgery formulas for quantum invariants and the Arf invariant. Geometry & Topology Monographs, 7:213–233, 2004.
- [15] Alexander Kirillov Jr and Benjamin Balsam. Turaev-Viro invariants as an extended TQFT. arXiv preprint, 2010. arXiv:1004.1533.
- [16] A. Yu. Kitaev. Fault-tolerant quantum computation by anyons. Ann. Physics, 303(1):2–30, 2003. doi:10.1016/S0003-4916(02)00018-0.
- [17] Robert Koenig, Greg Kuperberg, and Ben W. Reichardt. Quantum computation with Turaev-Viro codes. Ann. Physics, 325(12):2707–2749, 2010. doi:10.1016/j.aop.2010.08.001.
- [18] Hari Krovi and Alexander Russell. Quantum Fourier transforms and the complexity of link invariants for quantum doubles of finite groups. Comm. Math. Phys., 334(2):743–777, 2015. doi:10.1007/s00220-014-2285-5.
- [19] Greg Kuperberg. How hard is it to approximate the Jones polynomial? Theory Comput., 11:183–219, 2015. doi:10.4086/toc.2015.v011a006.
- [20] Greg Kuperberg and Eric Samperton. Computational complexity and 3-manifolds and zombies. Geom. Topol., 22(6):3623–3670, 2018. doi:10.2140/gt.2018.22.3623.
- [21] Greg Kuperberg and Eric Samperton. Coloring invariants of knots and links are often intractable. Algebr. Geom. Topol., 21(3):1479–1510, 2021. doi:10.2140/agt.2021.21.1479.
- [22] Clément Maria and Jonathan Spreer. Admissible colourings of 3-manifold triangulations for Turaev-Viro type invariants. In Piotr Sankowski and Christos D. Zaroliagis, editors, 24th Annual European Symposium on Algorithms, ESA 2016, August 22-24, 2016, Aarhus, Denmark, volume 57 of LIPIcs, pages 64:1–64:16. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2016. doi:10.4230/LIPICS.ESA.2016.64.
- [23] Clément Maria and Jonathan Spreer. A polynomial-time algorithm to compute Turaev-Viro invariants TV4,q of 3-manifolds with bounded first Betti number. Found. Comput. Math., 20(5):1013–1034, 2020. doi:10.1007/s10208-019-09438-8.
- [24] Sergei Matveev. Algorithmic Topology and Classification of 3-Manifolds. Springer, 9 edition, 2003.
- [25] James R. Munkres. Elements of algebraic topology. CRC Press, 1 edition, 1984.
- [26] Eric C. Rowell and Zhenghan Wang. Mathematics of topological quantum computing. Bull. Amer. Math. Soc. (N.S.), 55(2):183–238, 2018. doi:10.1090/bull/1605.
- [27] Winfried Scharlau. Quadratic and Hermitian Forms. Grundlehren der mathematischen Wissenschaften. Springer, 1985.
- [28] Daisuke Tambara and Shigeru Yamagami. Tensor Categories with Fusion Rules of Self-Duality for Finite Abelian Groups. Journal of Algebra, 209:692–707, 1998.
- [29] Ryan Thorngren and Yifan Wang. Fusion category symmetry I: anomaly in-flow and gapped phases. arXiv preprint, 2019. arXiv:1912.02817.
- [30] V. G. Turaev and O. Ya. Viro. State sum invariants of -manifolds and quantum -symbols. Topology, 31(4):865–902, 1992. doi:10.1016/0040-9383(92)90015-A.
- [31] Vladimir Turaev and Alexis Virelizier. On 3-dimensional homotopy quantum field theory, i. International Journal of Mathematics, 23(09):1250094, 2012.
- [32] Vladimir G. Turaev. Quantum invariants of knots and 3-manifolds, volume 18 of De Gruyter Studies in Mathematics. De Gruyter, Berlin, third edition, 2016. doi:10.1515/9783110435221.
- [33] Charles T C Wall. Quadratic forms on finite groups and related topics. Topology, 2:281–298, 1963.
- [34] Zhenghan Wang. Topological quantum computation, volume 112 of CBMS Regional Conference Series in Mathematics. Conference Board of the Mathematical Sciences, Washington, DC; by the American Mathematical Society, Providence, RI, 2010. doi:10.1090/cbms/112.