E2-quasi-exact solvability for non-Hermitian models

We propose the notion of $E_{2}$-quasi-exact solvability and apply this idea to find explicit solutions to the eigenvalue problem for a non-Hermitian Hamiltonian system depending on two parameters. The model considered reduces to the complex Mathieu Hamiltonian in a double scaling limit, which enables us to compute the exceptional points in the energy spectrum of the latter as a limiting process of the zeros for some algebraic equations. The coefficient functions in the quasi-exact eigenfunctions are univariate polynomials in the energy obeying a three-term recurrence relation. The latter property guarantees the existence of a linear functional such that the polynomials become orthogonal. The polynomials are shown to factorize for all levels above the quantization condition leading to vanishing norms rendering them to be weakly orthogonal. In two concrete examples we compute the explicit expressions for the Stieltjes measure.