Braiding and fusion rules of topological excitations are indispensable topological invariants in topological quantum computation and topological orders. While excitations in two dimensions (2D) are always particlelike anyons, those in three dimensions (3D) incorporate not only particles but also loops, spatially nonlocal objects, making it novel and challenging to study topological invariants in higher dimensions. While 2D fusion rules have been well understood from bulk Chern-Simons field theory and edge conformal field theory, it is yet to be thoroughly explored for 3D fusion rules from higher-dimensional bulk topological field theory. Here, we perform a field-theoretical study on (i) how loops that carry Abelian gauge fluxes fuse and (ii) how loops are shrunk into particles in the path integral, which generates fusion rules, loop-shrinking rules, and descendent invariants, e.g., quantum dimensions. We first assign a gauge-invariant Wilson operator to each excitation and determine the number of distinct excitations through equivalence classes of Wilson operators. Then, we adiabatically shift two Wilson operators together to observe how they fuse and are split in the path integral; despite the Abelian nature of the gauge fluxes carried by loops, their fusions may be of non-Abelian nature. Meanwhile, we adiabatically deform world sheets of unknotted loops into world lines and examine the shrinking outcomes; we find that the resulting loop-shrinking rules are algebraically consistent to fusion rules. Interestingly, fusing a pair of loop and antiloop may generate multiple vacua, but fusing a pair of anyon and antianyon in 2D has one vacuum only. By establishing a field-theoretical ground for fusion and shrinking in 3D, this work leaves intriguing directions for future exploration, e.g., symmetry enrichment, quantum gates, and topological order of braided monoidal 2-category of 2-group.
Braiding and fusion rules of topological excitations are indispensable topological invariants in topological quantum computation and topological orders. While excitations in two dimensions (2D) are always particlelike anyons, those in three dimensions (3D) incorporate not only particles but also loops, spatially nonlocal objects, making it novel and challenging to study topological invariants in higher dimensions. While 2D fusion rules have been well understood from bulk Chern-Simons field theory and edge conformal field theory, it is yet to be thoroughly explored for 3D fusion rules from higher-dimensional bulk topological field theory. Here, we perform a field-theoretical study on (i) how loops that carry Abelian gauge fluxes fuse and (ii) how loops are shrunk into particles in the path integral, which generates fusion rules, loop-shrinking rules, and descendent invariants, e.g., quantum dimensions. We first assign a gauge-invariant Wilson operator to each excitation and determine the number of distinct excitations through equivalence classes of Wilson operators. Then, we adiabatically shift two Wilson operators together to observe how they fuse and are split in the path integral; despite the Abelian nature of the gauge fluxes carried by loops, their fusions may be of non-Abelian nature. Meanwhile, we adiabatically deform world sheets of unknotted loops into world lines and examine the shrinking outcomes; we find that the resulting loop-shrinking rules are algebraically consistent to fusion rules. Interestingly, fusing a pair of loop and antiloop may generate multiple vacua, but fusing a pair of anyon and antianyon in 2D has one vacuum only. By establishing a field-theoretical ground for fusion and shrinking in 3D, this work leaves intriguing directions for future exploration, e.g., symmetry enrichment, quantum gates, and topological order of braided monoidal 2-category of 2-group.
Microscopic Constructions of the BF+AAB Topological Field Theory and Borromean-Rings Braiding in (3+1) Dimensions
Topological order was first studied in fractional quantum Hall states, where anyons reflect electron fractionalization. Recent advances in continuum topological field theories have revealed braiding, fusion, shrinking, and diagrammatic consistency relations in higher dimensions [e.g., Huang, Zhang, Ye, JHEP 11 (2023) 210; JHEP 06 (2025) 238]. As the next step of research, these advances call for explicit microscopic lattice constructions that faithfully reproduce the structures and processes described by continuum field theories. Here we construct, step by step, microscopic counterparts of the (3+1)-dimensional BF field theory with an AAB twist and gauge group (Z2)3. This field theory was introduced in our earlier work on Borromean-Rings braiding [Chan, Ye, Ryu, Phys. Rev. Lett. 121, 061601 (2018)]. We start with the construction of string and thickened membrane operators that create, fuse, shrink, and braid particle and loop excitations, providing microscopic counterparts of Wilson operators in continuum field theories. We derive fusion and shrinking rules and the operators governing particle--loop and Borromean-Rings braiding. We demonstrate selective control of non-Abelian shrinking channels through the internal degrees of freedom of loop operators and establish fusion--shrinking consistency for arbitrary finite groups. An explicit correspondence between all 22 excitations of the D4 quantum double lattice model and the BF+AAB theory preserves fusion and shrinking rules and particle--loop and Borromean-Rings braiding phases. These constructions place the continuum results on a microscopic lattice footing with tensor-product local Hilbert spaces and short-range interactions.
Symmetry fractionalization on topological excitations is one of the most remarkable quantum phenomena in topological orders with symmetry, i.e., symmetry-enriched topological phases. While much progress has been theoretically and experimentally made in two dimensions (2D), the understanding on symmetry fractionalization in 3D is far from complete. A long-standing challenge is to understand symmetry fractionalization on looplike topological excitations, which are spatially extended objects. In this paper, we construct a powerful topological-field-theoretical framework approach for 3D topological orders, which leads to a systematic characterization and classification of symmetry fractionalization. For systems with Abelian gauge groups (Gg) and Abelian symmetry groups (Gs), we successfully establish equivalence classes that lead to a finite number of patterns of symmetry fractionalization, although there are notoriously infinite number of Lagrangian-descriptions of the system. Based on this, we compute topologically distinct types of fractional symmetry charges carried by particles. Then, for each type, we compute topologically distinct statistical phases of braiding processes among loop excitations and external symmetry fluxes. As a result, we are able to unambiguously list all physical observables for each pattern of symmetry fractionalization. We present detailed calculations on many concrete examples. As an example, we find that the symmetry fractionalization in an untwisted ${\mathbb{Z}}_{2}\ifmmode\times\else\texttimes\fi{}{\mathbb{Z}}_{2}$ topological order with Z2 symmetry is classified by ${({\mathbb{Z}}_{2})}^{6}\ensuremath{\bigoplus}{({\mathbb{Z}}_{2})}^{2}\ensuremath{\bigoplus}{({\mathbb{Z}}_{2})}^{2}\ensuremath{\bigoplus}{({\mathbb{Z}}_{2})}^{2}$. If the topological order is twisted, the classification reduces to (Z2)6 in which particle excitations always carry integer charge. Inspired by the field-theoretical analysis, we find that our classification of symmetry fractionalization with twist $\ensuremath{\omega}$ can be formally organized into an algebraic formalism: ${⨁}_{{\ensuremath{\nu}}_{i}}{\mathcal{H}}^{4}({G}_{g}{\ensuremath{\leftthreetimes}}_{{\ensuremath{\nu}}_{i}}{G}_{s},U(1))/{\mathrm{\ensuremath{\Gamma}}}_{\ensuremath{\omega}}({\ensuremath{\nu}}_{i}),$ where anomaly-free symmetry fractionalization on particles ${\ensuremath{\nu}}_{i}\ensuremath{\in}{\mathcal{H}}_{non}^{2}({G}_{s},{G}_{g})$. Despite the lack of detailed dependence of ${\mathrm{\ensuremath{\Gamma}}}_{\ensuremath{\omega}}({\ensuremath{\nu}}_{i})$ on $\ensuremath{\omega}$ and ${\ensuremath{\nu}}_{i}$, the present field-theoretical approach allows us to efficiently calculate and understand ${\mathrm{\ensuremath{\Gamma}}}_{\ensuremath{\omega}}({\ensuremath{\nu}}_{i})$. Several future directions are presented at the end of this paper.
Models for topological quantum computation are based on braiding and fusing anyons (quasiparticles of fractional statistics) in (2+1)-D. The anyons that can exist in a physical theory are determined by the symmetry group of the Hamiltonian. In the case that the Hamiltonian undergoes spontaneous symmetry breaking of the full symmetry group G to a finite residual gauge group H, particles are given by representations of the quantum double D(H) of the subgroup. The quasi-triangular Hopf Algebra D(H) is obtained from Drinfeld's quantum double construction applied to the algebra F(H) of functions on the finite group H. A major new contribution of this work is a program written in MAGMA to compute the particles (and their properties - including spin) that can exist in a system with an arbitrary finite residual gauge group, in addition to the braiding and fusion rules for those particles. We compute explicitly the fusion rules for two non-abelian group doubles suggested for universal quantum computation: S3 and A5, and discover some interesting results, subsystems, and symmetries in the tables. SO(3)_4 (the restriction of Chern-Simons theory SU(2)4) and its mirror image are discovered as 3-particle subsystems in the 8-particle S3 quantum double. The tables demonstrate that both S3 and A5 anyons are all Majorana, but this is not the case for all finite groups. In the appendices, the quantum doubles for the remaining nonabelian subgroups of SO(3) - S4, A4, and D4 (the second in the infinite family Dn) - are tabulated and analyzed. In addition, the probabilities of obtaining any given fusion product in quantum computation applications are determined and programmed in MAGMA. Throughout, connections to possible experiments are mentioned.
Braiding statistics of p- and q-dimensional topological excitations is conventionally defined in p+q+2 spatial dimensions. We find a novel statistical process WN(X,Y)=(Y−1X−1)N(YX)N for two order-N excitations in p+q+1 dimensions, detecting the Bockstein response A⌣β(B). This new statistics and fermionic loop statistics exhaust all loop statistics in three dimensions whose fusion rules form an Abelian group G, classified by H5(B2G,U(1)). Surprisingly, conventional three-loop braiding goes beyond this classification, so it must have non-Abelian fusion rules. We suggest viewing three-loop braiding as particle-loop braiding together with exotic fusion rules between loops and point-like defects. We also try to clarify the relationship between statistics and symmetry anomaly.
We investigate whether the braid group representation on N anyons uniquely determines the modular data (S and T matrices) of a non-Abelian topological order, and whether this framework classifies all fusion rules for Majorana zero modes (MZMs) in 2D topological superconductors. We show that the local N-anyon braid group representation constrains the modular data up to finite gauge ambiguity but does not uniquely fix the global S-matrix without specifying the underlying spin structure. For MZMs, the braid generators satisfy the Ising anyon R-matrix with the asymptotic Hecke relation σi2=1 (up to phase), which rigidly fixes the local F and R symbols of the Ising category. Consequently, the fusion rules are strictly of Ising type (σ×σ=1+ψ), precluding alternative non-Abelian fusion rules. The T-matrix is determined by topological spins, but the full S-matrix depends on the spin structure of the underlying manifold, reflecting the super-modular structure of the theory. ## 1. Introduction Topological quantum computation exploits the non-Abelian statistics of anyons—quasiparticle excitations in two-dimensional systems whose braiding operations are represented by non-commuting unitary matrices on a degenerate ground-state manifold [1]. A central question is whether the local braiding data of N anyons suffices to reconstruct the full modular data of the underlying topological order (TO), or whether additional global information is required. This question has direct physical consequences for 2D topological superconductors (TSCs) hosting Majorana zero modes (MZMs), where the anyonic excitations are vortices binding MZMs and the fermionic quasiparticle is a transparent (local) fermion [2,3]. The modular data of a TO consists of the S-matrix, which encodes the mutual braiding statistics of anyons, and the T-matrix, which encodes the topological spins. Together, these matrices form a projective representation of t
Universal topological data of topologically ordered phases can be captured by topological quantum field theory in continuous space time by taking the limit of low energies and long wavelengths. While previous continuum field-theoretical studies of topological orders in three-dimensional (3D) real space focused on either self-statistics, braiding statistics, shrinking rules, fusion rules, or quantum dimensions, it has yet to systematically put all topological data together in a unified continuum field-theoretical framework. Here, we construct the topological BF field theory with twisted terms (e.g., AAdA and AAB) as well as a K-matrix BB term to simultaneously explore all such topological data and reach anomaly-free topological orders. Following the spirit of the famous K-matrix Chern-Simons theory of two-dimensional topological orders, we present general formulas and systematically show how the K-matrix BB term confines topological excitations, and how self-statistics of particles is transmuted between bosonic one and fermionic one. To reach anomaly-free topological orders, we explore, within the present continuum field-theoretical framework, how the principle of gauge invariance fundamentally influences possible realizations of topological data. More concretely, we present the topological actions of (i) particle-loop braidings with emergent fermions, (ii) multiloop braidings with emergent fermions, and (iii) Borromean rings braidings with emergent fermions, and calculate their universal topological data. Together with previous efforts, our paper paves the way toward a more systematic and complete continuum field-theoretical analysis of exotic topological properties of 3D topological orders. Several interesting future directions are also discussed.