On the stability of topological order in tensor network states

We construct a tensor network representation of the three-dimensional (3D) toric code ground state that is stable to a generating set of uniform local tensor perturbations, including those that do not map to local operators on the physical Hilbert space. The stability is established by mapping the phase diagram of the perturbed tensor network to that of the 3D Ising gauge theory, which has a nonzero finite temperature transition. More generally, we find that the stability of a topological tensor network state is determined by the form of its virtual symmetries and the topological excitations created by virtual operators that break those symmetries. In particular, a dual representation of the 3D toric code ground state, as well as representations of the X-cube and cubic code ground states, for which pointlike excitations are created by such operators, are found to be unstable.