The cubic code model is studied in the presence of arbitrary extensive perturbations. Below a critical perturbation strength, we show that most states with finite energy are localized; the overwhelming majority of such states have energy concentrated around a finite number of defects, and remain so for a time that is near exponential in the distance between the defects. This phenomenon is due to an emergent superselection rule and does not require any disorder. Local integrals of motion for these finite energy sectors are identified as well. Our analysis extends more generally to systems with immobile topological excitations.

• Received 9 October 2015

DOI:https://doi.org/10.1103/PhysRevLett.116.027202