Recently, it was argued that the braiding and statistics of anyons in a two-dimensional topological phase can be extracted by studying the quantum entanglement of the degenerate ground states on the torus. This construction either required a lattice symmetry (such as $\ensuremath{\pi}/2$ rotation) or tacitly assumed that the minimum entanglement states (MESs) for two different bipartitions can be uniquely assigned quasiparticle labels. Here we describe a procedure to obtain the modular
S matrix, which encodes the braiding statistics of anyons, which does not require making any of these assumptions. Our strategy is to coherently compare MESs of three independent entanglement bipartitions of the torus, which leads to a unique modular
S. This procedure also puts strong constraints on the modular
T and
U matrices without requiring any symmetries, and in certain special cases, completely determines it. Our method applies equally to Abelian and non-Abelian topological phases.
The topological order is encoded in the pattern of long-range quantum entanglements, which cannot be measured by any local observable. Here we perform an exact diagonalization study to establish the non-Abelian topological order for topological band models through entanglement entropy measurement. We focus on the quasiparticle statistics of the non-Abelian Moore-Read and Read-Rezayi states on the lattice models with bosonic particles. We identify multiple independent minimal entangled states (MESs) in the ground state manifold on a torus. The extracted modular S matrix from MESs faithfully demonstrates the Ising anyon or Fibonacci quasiparticle statistics, including the quasiparticle quantum dimensions and the fusion rules for such systems. These findings unambiguously demonstrate the topological nature of the quantum states for these flatband models without using the knowledge of model wave functions.
Topologically ordered phases are gapped states, defined by the properties of excitations when taken around one another. Here we demonstrate a method to extract the statistics and braiding of excitations, given just the set of ground-state wave functions on a torus. This is achieved by studying the topological entanglement entropy (TEE) upon partitioning the torus into two cylinders. In this setting, general considerations dictate that the TEE generally differs from that in trivial partitions and depends on the chosen ground state. Central to our scheme is the identification of ground states with minimum entanglement entropy, which reflect the quasiparticle excitations of the topological phase. The transformation of these states allows for the determination of the modular
S and
U matrices which encode quasiparticle properties. We demonstrate our method by extracting the modular
S matrix of a chiral spin liquid phase using a Monte Carlo scheme to calculate the TEE and prove that the quasiparticles obey semionic statistics. This method offers a route to nearly complete determination of the topological order in certain cases.
Abstract An essential aspect of topological phases of matter is the existence of Wilson loop operators that keep the ground state subspace invariant. Here we present and implement an \it unbiased \rm numerical optimization scheme to systematically find the Wilson loop operators given a single ground state wave function of a gapped Hamiltonian on a disk. We then show how these Wilson loop operators can be cut and glued through further optimization to give operators that can create, move, and annihilate anyon excitations. We subsequently use these operators to determine the braiding statistics and topological twists of the anyons, yielding a way to fully extract topological order from a single wave function. We apply our method to the ground state of the perturbed toric code and doubled semion models with a magnetic field that is up to a half of the critical value. From a contemporary perspective, this can be thought of as a machine learning approach to discover emergent 1-form symmetries of a ground state wave function. From an application perspective, our approach can be relevant to find Wilson loop operators in current quantum simulators.
A topologically ordered phase on a torus possesses degenerate ground states that transform nontrivially under the modular transformations of the torus, generated by Dehn twists. Representation of modular transformations on the ground states (modular matrices) characterizes the topological order. We show that the modular matrices can be numerically measured as the non-Abelian Berry phase of adiabatic deformations of the lattice model placed on a torus. We apply this method to the example of a gauged
px+ipy superconductor, and show that the result is consistent with the topological quantum field theory descriptions.
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