Entanglement Hamiltonian and orthogonal polynomials

We study the entanglement Hamiltonian for free-fermion chains with a particular form of inhomogeneity. The hopping amplitudes and chemical potentials are chosen such that the single-particle eigenstates are related to discrete orthogonal polynomials of the Askey scheme. Due to the bispectral properties of these functions, one can construct an operator which commutes exactly with the entanglement Hamiltonian and corresponds to a linear or parabolic deformation of the physical one. We show that this deformation is interpreted as a local inverse temperature and can be obtained in the continuum limit via methods of conformal field theory. Using this prediction, the properly rescaled eigenvalues of the commuting operator are found to provide a very good approximation of the entanglement spectrum and entropy.