Computation of entanglement entropy in inhomogeneous free fermions chains by algebraic Bethe ansatz

The computation of the entanglement entropy for inhomogeneous free fermions chains based on q-Racah polynomials is considered. The eigenvalues of the truncated correlation matrix are obtained from the diagonalization of the associated Heun operator via the algebraic Bethe ansatz. In the special case of chains based on dual q-Hahn polynomials, the eigenvectors and eigenvalues are expressed in terms of symmetric polynomials evaluated on the Bethe roots.