Phase diagram of a symmetry-enriched toric code

We analyze the effect of perturbations on a toric code-like Hamiltonian enriched with global symmetries. The ground states of this Hamiltonian are superpositions of loops decorated with cluster states, resulting in anyonic excitations that fractionalize a global on-site $\mathbb{Z}_2\times\mathbb{Z}_2$, inversion and time reversal symmetry. We use a modified version of the string order parameters proposed in [New Journal of Physics 21, 113016 (2019)], adapted to work beyond the framework of tensor networks, to capture the symmetry fractionalization pattern of the anyons. We show that this order parameter accurately captures the symmetry fractionalization in a direct Hamiltonian interpolation between the decorated and trivial toric codes, and in the presence of symmetry-preserving local magnetic fields and Ising interactions. We also observe how the condensation of an anyon that fractionalizes the symmetry forces the symmetry to spontaneously break at the phase transition. This notably changes the phase diagram of this Hamiltonian under parallel magnetic fields compared to the regular toric code since the condensation of charge and flux anyons drives the system to two distinct trivial phases, a symmetry broken phase and symmetric phase. Moreover, we provide a proof of spontaneous symmetry breaking from anyon condensation in the framework of projected entangled pair states.