Heun operator of Lie type and the modified algebraic Bethe ansatz

The generic Heun operator of Lie type is identified as a certain $BC$-Gaudin magnet Hamiltonian in a magnetic field. By using the modified algebraic Bethe ansatz introduced to diagonalize such Gaudin models, we obtain the spectrum of the generic Heun operator of Lie type in terms of the Bethe roots of inhomogeneous Bethe equations. We show also that these Bethe roots are intimately associated to the roots of polynomial solutions of the differential Heun equation. We illustrate the use of this approach in two contexts: the representation theory of $O(3)$ and the computation of the entanglement entropy for free Fermions on the Krawtchouk chain.