An integrable model with parafermion zero energy modes

Parafermion zero energy modes are a vital element of fault-tolerant topological quantum computation. Although it is believed that such modes form on the border between topological and normal phases, this has been demonstrated only for Z$_2$ (Majorana) and Z$_3$ parafermions. I consider an integrable model where such demonstration is possible for Z$_N$ parafermions with any $N$. The procedure is easely generalizable for more complicated symmetry groups.