  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.
