The bulk-boundary correspondence is a hallmark feature of topological phases of matter. Nonetheless, our understanding of the correspondence remains incomplete for phases with intrinsic topological order, and is nearly entirely lacking for more exotic phases, such as fractons. Intriguingly, for the former, recent work suggests that bulk topological order manifests in a non-local structure in the boundary Hilbert space; however, a concrete understanding of how and where this perspective applies remains limited. Here, we provide an explicit and general framework for understanding the bulk-boundary correspondence in Pauli topological stabilizer codes. We show -- for any boundary termination of any two-dimensional topological stabilizer code -- that the boundary Hilbert space cannot be realized via local degrees of freedom, in a manner precisely determined by the anyon data of the bulk topological order. We provide a simple method to compute this "obstruction" using a well-known mapping to polynomials over finite fields. Leveraging this mapping, we generalize our framework to fracton models in three-dimensions, including both the X-Cube model and Haah's code. An important consequence of our results is that the boundaries of topological phases can exhibit emergent symmetries that are impossible to otherwise achieve without an unrealistic degree of fine tuning. For instance, we show how linear and fractal subsystem symmetries naturally arise at the boundaries of fracton phases.
The bulk-boundary correspondence is a hallmark feature of topological phases of matter. Nonetheless, our understanding of the correspondence remains incomplete for phases with intrinsic topological order, and is nearly entirely lacking for more exotic phases, such as fractons. Intriguingly, for the former, recent work suggests that bulk topological order manifests in a non-local structure in the boundary Hilbert space; however, a concrete understanding of how and where this perspective applies remains limited. Here, we provide an explicit and general framework for understanding the bulk-boundary correspondence in Pauli topological stabilizer codes. We show -- for any boundary termination of any two-dimensional topological stabilizer code -- that the boundary Hilbert space cannot be realized via local degrees of freedom, in a manner precisely determined by the anyon data of the bulk topological order. We provide a simple method to compute this "obstruction" using a well-known mapping to polynomials over finite fields. Leveraging this mapping, we generalize our framework to fracton models in three-dimensions, including both the X-Cube model and Haah's code. An important consequence of our results is that the boundaries of topological phases can exhibit emergent symmetries that are impossible to otherwise achieve without an unrealistic degree of fine tuning. For instance, we show how linear and fractal subsystem symmetries naturally arise at the boundaries of fracton phases.
Topological holography (TH), or SymTFT, realizes symmetries and dualities of a quantum system as boundary data of a topological bulk in one higher dimension. We formulate fracton topological holography (FTH), extending this mechanism from liquid topological orders to fracton stabilizer codes. The construction is organized as a general four-stage framework: prepare the bulk model and compute its excitations, determine boundary data and admissible gapped top boundaries, identify the low-energy preserving operator algebra together with its symmetry, relation, and twist data, and then switch among top boundaries to compare the induced boundary descriptions. As a type-I example, we develop FTH for the X-cube model with smooth and rough top boundaries; for a minimal effective Hamiltonian, both yield transverse-field plaquette Ising models, with exchanged subsystem symmetry and twist data, and the boundary switch is implemented by a linear-depth local unitary sequential quantum circuit (SQC). As a type-II example, we formulate FTH for Haah's cubic code in the Laurent-polynomial stabilizer formalism and analyze the natural (Z) and (X) top boundaries, which induce two two-dimensional qubit systems related locally by exchanging generalized plaquette Ising and transverse-field terms and nonlocally by a symmetry--relation duality. These results show that FTH is a genuine extension of TH to both type-I and type-II fracton orders. FTH therefore provides a concrete framework for organizing and understanding duality, with the prospect of offering a systematic route to new dualities.
Entanglement entropy provides a powerful characterization of two-dimensional gapped topological phases of quantum matter, intimately tied to their description by topological quantum field theories (TQFTs). Fracton topological orders are three-dimensional gapped topologically ordered states of matter that lack a TQFT description. We show that three-dimensional fracton phases are nevertheless characterized, at least partially, by universal structure in the entanglement entropy of their ground-state wave functions. We explicitly compute the entanglement entropy for two archetypal fracton models, the ``X-cube model'' and ``Haah's code,'' and demonstrate the existence of a nonlocal contribution that scales linearly in subsystem size. We show via Schrieffer-Wolff transformations that this piece of the entanglement entropy of fracton models is robust against arbitrary local perturbations of the Hamiltonian. Finally, we argue that these results may be extended to characterize localization-protected fracton topological order in excited states of disordered fracton models.
Quantum low-density parity-check (qLDPC) codes promise low-overhead fault-tolerant quantum memory, but many leading constructions are naturally defined on periodic lattices, whereas scalable quantum hardware favors planar architectures with open boundaries. In this work, we develop a general operator-algebraic framework and computational algorithms for constructing all gapped boundaries and defects of two-dimensional topological Pauli stabilizer codes. Starting from a truncated bulk code, our method identifies boundary gauge operators, boundary anyons, and their string operators, and then constructs microscopic boundary Hamiltonians through anyon condensation and topological-order completion, thereby eliminating local logical operators near the boundary. We establish a one-to-one bulk-boundary correspondence between bulk anyons and excitations of the anomalous boundary subsystem code, enabling fusion rules, topological spins, and braiding statistics to be extracted directly from the boundary operator algebra. For translation-invariant codes, Laurent-polynomial representations reduce these tasks to finite matrix computations based on Hermite and Smith normal forms, enabling automated constructions for Z d qudits of both prime and composite dimension. We demonstrate the framework on toric, color, double-semion, six-semion, and anomalous three-fermion codes, and apply it to bivariate bicycle codes, revealing large logical dimensions and unusually long translation periods for anyon strings. By systematically converting bulk topological data into explicit boundary and defect Hamiltonians, our framework can serve as a foundation for constructing low-overhead planar qLDPC codes with open boundaries and for engineering fault-tolerant logical operations.
We present a general approach to the bulk-boundary correspondence of noninvertible topological phases, including both topological and fracton orders. This is achieved by a novel bulk construction protocol where solvable (d+1) ( d + 1 ) -dimensional bulk models with noninvertible topology are constructed from the so-called generalized Ising (GI) models in d d dimensions. The GI models can then terminate on the boundaries of the bulk models. The construction generates abundant examples, including not only prototype ones such as \mathbb{Z}_2 ℤ 2 toric code models in any dimensions no less than two, and the X-cube fracton model, but also more diverse ones such as the \mathbb{Z}_2× \mathbb{Z}_2 ℤ 2 × ℤ 2 topological order, the 4d \mathbb{Z}_2 ℤ 2 topological order with pure-loop excitations, etc. The boundary of the solvable model is potentially anomalous and corresponds to precisely only sectors of the GI model that host certain total symmetry charges and/or satisfy certain boundary conditions. We derive a concrete condition for such bulk-boundary correspondence. The condition is violated only when the bulk model is either trivial or fracton ordered. A generalized notion of Kramers-Wannier duality plays an important role in the construction. Also, utilizing the duality, we find an example where a single anomalous theory can be realized on the boundaries of two distinct bulk fracton models, a phenomenon not expected in the case of topological orders. More generally, topological orders may also be generated starting with lattice models beyond the GI models, such as those with symmetry protected topological orders, through a variant bulk construction, which we provide in an appendix.
ZN generalizations of three-dimensional stabilizer codes
In this work, we generalize several three-dimensional Z2 stabilizer models--including the X-cube model, the three-dimensional toric code, and Haah's code--to their ZN counterparts. Under periodic boundary conditions, we analyze their ground state degeneracies and topological excitations, and uncover behaviors that strongly depend on system size. For the X-cube model, we identify excitations with mobility restricted under local operations but relaxed under nonlocal ones derived from global topology. These excitations, previously confined to open boundaries in the Z2 model, now appear even under periodic boundaries. In the toric code, we observe nontrivial braiding between string and point excitations despite the absence of ground state degeneracy, indicating long-range entanglement independent of topological degeneracy. Again, this effect extends from open to periodic boundaries in the generalized models. For Haah's code, we find new excitations--fracton tripoles and monopoles--that remain globally constrained, along with a relaxation of immobility giving rise to lineons and planons. These results reveal new forms of topological order and suggest a broader framework for understanding fracton phases beyond the conventional Z2 setting.
Stabilizer codes provide exact lattice realizations of bosonic topological orders. In contrast, systematic stabilizer descriptions of intrinsically fermionic topological phases remain much less developed. In this work, we introduce Majorana-Pauli stabilizer codes, a class of exactly solvable fermionic lattice models whose stabilizers are built from both generalized Pauli operators and Majorana operators. As a main example, we construct an exactly solvable stabilizer realization of the fermionic toric code: an intrinsically fermionic Z2 topological order in (2+1) dimensions, using Z8 Pauli operators coupled to Majorana modes. Within this stabilizer framework, the anyons, string operators, fusion rules, and braiding statistics all follow naturally from the stabilizer algebra. More broadly, we show that the fermionic toric code belongs to a duality web generated by anyon condensation and by gauging or ungauging bosonic or fermionic symmetries. This web connects bosonic topological orders, symmetry-enriched topological (SET) phases, and both bosonic and fermionic symmetry-protected topological (SPT) phases, all within a common stabilizer description. We extend this construction to all Abelian fermionic topological orders with gapped boundaries and to fermionic SPT phases whose fully gauged theories are Abelian. Going beyond Majorana operators, we introduce fermionic versions of the clock and shift operators and use them to construct an exact bosonization map for ZDF symmetries for D even. Using this, we realize a stabilizer model for a nontrivial Z8F fermionic SPT phase with no free-fermion analog. Altogether, these results extend the stabilizer-code paradigm to a broad class of intrinsically fermionic phases bridging fermionic quantum many-body physics to quantum error correction.