  We consider two-dimensional lattice models that support Ising anyonic
excitations and are coupled to a thermal bath. We propose a phenomenological
model for the resulting short-time dynamics that includes pair-creation,
hopping, braiding, and fusion of anyons. By explicitly constructing topological
quantum error-correcting codes for this class of system, we use our
thermalization model to estimate the lifetime of the quantum information stored
in the encoded spaces. To decode and correct errors in these codes, we adapt
several existing topological decoders to the non-Abelian setting. We perform
large-scale numerical simulations of these two-dimensional Ising anyon systems
and find that the thresholds of these models range between 13% to 25%. To our
knowledge, these are the first numerical threshold estimates for quantum codes
without explicit additive structure.
