  The exact solution for the two particle spectrum of the Heisenberg-Peierls
one dimensional spin chain is given by working in the fermionic representation.
The resulting equations for the eigenvalues are, in some sense, similar to
those of the Richardson's solution of the BCS model and must be solved
numerically.
