Topological Entanglement Entropy of Fracton Stabilizer Codes

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, but the existence of a TQFT description for these phases remains an open question. 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 topological contribution that scales linearly in subsystem size. We show via Schrieffer-Wolff transformations that the topological entanglement 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.