The Elliptic Gaudin Model: a Numerical Study

The elliptic Gaudin model describes completely anisotropic spin systems with long range interactions. The model was proven to be quantum integrable by Gaudin and latter the exact solution was found by means of the algebraic Bethe ansatz. In spite of the appealing properties of the model, it has not yet been applied to any physical problem. We here generalize the exact solution to systems with arbitrary spins, and study numerically the behavior of the Bethe roots for a system with three different spins. Then, we propose an integrable anisotropic central spin model that we study numerically for very large systems.