Symmetry-enriched topological order and quasi-fractonic behavior in $\mathbb{Z}_N$ stabilizer codes

Abstract: We study a broad class of qudit stabilizer codes, termed $\mathbb{Z}_N$ bivariate-bicycle (BB) codes, arising either as two-dimensional realizations of modulated gauge theories or as $\mathbb{Z}_N$ generalizations of binary BB codes. Our central finding, derived from the polynomial representation, is that the essential topological properties of these $\mathbb{Z}_N$ codes can be determined by the properties of their $\mathbb{Z}_p$ counterparts, where $p$ are the prime factors of $N$, even when $N$ contains prime powers ($N = \prod_i p_i^{k_i}$). This result yields a significant simplification by leveraging the well-studied framework of codes with prime qudit dimensions. In particular, this insight directly enables the generalization of the algebraic-geometric methods (e.g., the Bernstein-Khovanskii-Kushnirenko theorem) to determine anyon fusion rules in the general qudit situation. Moreover, we analyze the model's symmetry-enriched topological order (SET) to reveal a quasi-fractonic behavior, resolving the anyon mobility puzzle in this class of models. We also present a computational algebraic method using Gr\"{o}bner bases over the ring of integers to efficiently calculate the topological order and its SET properties.