We construct a previously unknown E 2 -quasi-exactly solvable non-Hermitian model whose eigenfunctions involve weakly orthogonal polynomials obeying three-term recurrence relations that factorize beyond the quantization level. The model becomes Hermitian when one of its two parameters is fixed to a specific value. We analyze the double scaling limit of this model leading to the complex Mathieu equation. The norms, Stieltjes measures and moment functionals are evaluated for some concrete values of one of the two parameters.In [1] we introduced E 2 -quasi-exactly solvable models in analogy to the notion of sl 2 (C)-quasi-exactly solvability originally proposed by Turbiner [2,3]. The different setting is motivated mathematically by the fact that solutions for E 2 -quasi-exactly solvable models do not belong to the general class of hypergeometric functions which emerge as solutions from an sl 2 (C)-setting. The physical motivation results from the current interest in extending the study of solvable models [4,5,6,7, 1] to non-Hermitian quantum mechanical systems [8,9,10,11]. The E 2 -quasi-exactly solvable models are especially interesting in optical settings [12,13,14,15,16,17,18,19,20] where the fact is exploited that the Helmholtz equation results as a reduction from the Schrödinger equation. Solvable models are rare exceptions in the study of quantum mechanical systems and the model presented here should be added to that list.The starting point for the construction of E 2 -quasi-exactly systems consists of expressing the Hamiltonian operator H of the model in terms of the E 2 -basis operators u, v and J obeying the commutation relationsinstead of the standard sl 2 (C)-generators. We now use the particular realization [21]and demand a specific anti-linear symmetry [22], PT 3 : J → J, u → v, v → u, i → −i, as defined in [18], which for (2) becomes PT 3 : θ → π/2 − θ, i → −i. The operators in (2)