Quasi-Exact Solvable Systems

• Quasi-exactly solvable systems are quantum models that admit a finite invariant subspace, allowing algebraic construction of a limited set of eigenstates.
• They utilize Lie-algebraic frameworks, often based on sl(2) generators, and methods like the Bethe ansatz to derive explicit energy levels.
• These systems bridge exactly solvable and non-solvable models, impacting fields from anharmonic oscillators to PT-symmetric and non-Hermitian quantum mechanics.

Quasi-exactly solvable (QES) systems are quantum mechanical models for which a finite, nontrivial part of the spectrum—specifically, a finite number of energy levels and corresponding eigenfunctions—can be constructed algebraically in closed form, while the remainder of the spectrum is not explicitly accessible. QES systems occupy an intermediate position between exactly solvable (ES) models, where the entire spectrum is accessible, and non-solvable systems, where only numerical or perturbative techniques are available (Turbiner, 2016). QES theory provides a powerful framework to systematically uncover algebraic structures and invariant subspaces, underpinning partial solvability in both differential and difference operator settings.

1. Defining Quasi-Exact Solvability

A quantum mechanical operator (usually a Schrödinger Hamiltonian, but sometimes a discrete or matrix operator) is called quasi-exactly solvable if there exists a finite-dimensional invariant subspace of the Hilbert space, typically a space of polynomials of degree ≤ $n$, that is preserved by the operator. Explicitly, in the one-dimensional Schrödinger context,

• Exactly solvable: all eigenfunctions belong to an infinite flag of polynomial spaces $P_0 \subset P_1 \subset \cdots$, and the entire spectrum is algebraically accessible.
• QES: only a single finite-dimensional space $P_{n+1} = \mathrm{span}\{1, x, \ldots, x^n\}$ is invariant; exactly $n+1$ eigenfunctions and eigenvalues can be constructed algebraically (Turbiner, 2016, Li et al., 26 Nov 2025).
• Non-solvable: no such finite-dimensional invariant subspace exists.

The physical significance lies in explicit knowledge of a segment of the spectrum, which reveals nontrivial algebraic, analytic, and symmetry properties of the system.

2. Lie-Algebraic and Algebraic Structures

The emergence of quasi-exact solvability is deeply linked to hidden algebraic structures, most famously the embedding of the Hamiltonian (after a gauge transformation and change of variable) into the enveloping algebra of a Lie algebra—often $\mathfrak{sl}(2)$ or its deformations:

• One-dimensional cases: QES differential operators are quadratic polynomials of $\mathfrak{sl}(2)$ generators. The standard realization is

$$J^+ = x^2 \frac{d}{dx} - n x, \quad J^0 = x \frac{d}{dx} - \frac{n}{2}, \quad J^- = \frac{d}{dx}$$

ensuring that $T_2(J^+, J^0, J^-)$ preserves $P_{n+1}$ (Turbiner, 2016, Li et al., 26 Nov 2025, Agboola et al., 2011).

• Deformed algebras and higher dimensions: QES systems can arise from deformations of $\mathfrak{sl}(2)$—for example, cubic or quadratic polynomials in the commutator algebra of generalized Heun equations (Roy et al., 2013)—or from higher-rank algebras ($P_0 \subset P_1 \subset \cdots$0 for QES systems on spheres $P_0 \subset P_1 \subset \cdots$1 (Jr. et al., 2014)).
• Difference operators: Discrete QES systems utilize structures such as the Askey–Wilson algebra and involve "sinusoidal coordinates" to produce finite-dimensional polynomial invariant spaces (Sasaki, 2010).
• Non-Hermitian cases: $P_0 \subset P_1 \subset \cdots$2-QES operators, built from the Euclidean algebra in two dimensions, especially in PT-symmetric contexts, illustrate how the concept generalizes beyond $P_0 \subset P_1 \subset \cdots$3 (Fring, 2014, Fring et al., 2018).

This algebraic underpinning leads to systematic recipes for constructing QES operators and classifying their polynomial sectors.

3. Construction Methods: Polynomial Invariance and Bethe Ansatz

QES systems are identified by two main techniques: direct Lie-algebraic construction and Bethe ansatz solutions.

• Lie-algebraic construction: After an appropriate gauge transformation and change of variable, the Hamiltonian is written as a quadratic polynomial in the generators of the underlying Lie algebra; invariance of $P_0 \subset P_1 \subset \cdots$4 fixes the allowed model parameters and the dimension of the solvable subspace (Turbiner, 2016, Li et al., 26 Nov 2025).
• Bethe ansatz method: Solutions are constructed as products over roots $P_0 \subset P_1 \subset \cdots$5, $P_0 \subset P_1 \subset \cdots$6. The roots solve the algebraic Bethe ansatz equations where the residues at simple poles vanish, yielding quantization conditions for the eigenvalues and model parameters (Li et al., 2023, Agboola et al., 2011).
• Prepotential approach: An alternative, symmetry-agnostic method determines solvability in terms of two polynomials and associated Bethe ansatz equations for the roots, simultaneously constructing potential, eigenvalues, and eigenfunctions (Ho, 2023).

An archetypal example is the sextic anharmonic oscillator, with quasi-exactly solvable energies and wavefunctions constructed via both approaches and analyzed via truncated series and Bender–Dunne polynomials (Dorey et al., 2012, Handy et al., 2014).

4. Invariant Subspaces and Polynomial Sectors

A central feature of QES systems is the existence of a single finite-dimensional invariant subspace under the action of the transformed Hamiltonian:

• For QES models, only the first $P_0 \subset P_1 \subset \cdots$7 eigenvalues and eigenfunctions can be found algebraically. The truncation in the polynomial ansatz or moment equation (as a "kink" in the recurrence order) reflects this non-extendability to higher degrees (Handy et al., 2014, Dorey et al., 2012).
• In difference operator settings (discrete QES), the invariant subspace is spanned by $P_0 \subset P_1 \subset \cdots$8 for a suitable sinusoidal coordinate $P_0 \subset P_1 \subset \cdots$9, and invariance is achieved only for degrees $P_{n+1} = \mathrm{span}\{1, x, \ldots, x^n\}$0 (Sasaki, 2010).

Beyond the algebraic sector, the system admits no further closed-form solutions; higher eigenvalues require numerical or perturbative computation.

5. Extensions: Deformations, Supersymmetry, Superintegrability

• Algebraic deformations: QES systems are not restricted to conventional $P_{n+1} = \mathrm{span}\{1, x, \ldots, x^n\}$1 embeddings. Heun-type equations, for instance, yield cubic deformations, and certain QES models cannot be represented with standard $P_{n+1} = \mathrm{span}\{1, x, \ldots, x^n\}$2 generators (Roy et al., 2013, Leon et al., 2014).
• Supersymmetric QES models: Supersymmetric quantum mechanics techniques, such as the Krein–Adler construction, extend QES theory to the systematic state addition/deletion beyond shape-invariant potentials. This allows for the creation of finite-dimensional chains of QES partner Hamiltonians, often entirely classified via their Bethe ansatz roots (Li et al., 26 Nov 2025).
• Multi-dimensional QES systems: QES theory extends to higher dimensions (e.g., on $P_{n+1} = \mathrm{span}\{1, x, \ldots, x^n\}$3 or the plane), with invariance under $P_{n+1} = \mathrm{span}\{1, x, \ldots, x^n\}$4 or other higher-rank algebras, allowing for superintegrable or non-maximally superintegrable structures (Jr. et al., 2014, Tremblay et al., 2009).

6. Physical Realizations and Applications

• Potentials: QES structure is exhibited in families such as the sextic and symmetrized quartic oscillators (Znojil, 2016), singular and non-polynomial oscillators, isotonic and soft-core Coulomb potentials (Agboola et al., 2011), and magnetic two-body problems (Kreshchuk, 2015).
• Non-Hermitian and PT-symmetric models: QES concepts extend to non-Hermitian, PT-symmetric models with logical extensions of the Lie-algebraic machinery and occurrence of truncated orthogonal polynomial recurrences (Fring, 2014, Fring et al., 2018).
• Higher-order ODEs: QES methods generalize to third- and higher-order ordinary differential equations under specific algebraic constraints detectable via projective triviality criteria (Dorey et al., 2012).
• Time-dependence: Explicitly time-dependent QES Hamiltonians can be constructed using dynamical invariants and appropriate metrics, extending solvable subspaces and spectral results to non-stationary scenarios (Fring et al., 2018).

7. Moment Methods, Orthogonal Polynomials, and Further Developments

• Moment equation and OPPQ: The order-changing "kink" in finite-difference moment recursions precisely signals the existence and boundaries of the QES subspace. Orthogonal Polynomial Projection Quantization (OPPQ) utilizes this for robust and unified computation of both QES and non-QES spectra (Handy et al., 2014).
• Canonical and difference-operator correspondence: There exists an isospectral mapping between continuous QES operators and their difference- or $P_{n+1} = \mathrm{span}\{1, x, \ldots, x^n\}$5-difference analogs, leveraging Fock-space structures and quantum canonical transformations (Turbiner, 2016).
• Stäckel transform and coupling constant metamorphosis: These tools systematically relate QES systems with distinct parameterizations, enabling the algebraization and solution of otherwise inaccessible models (notably Hooke's atoms in magnetic fields, cosmological models) via transformations that interchange energy and coupling constants (Li et al., 19 Feb 2025).

---

For comprehensive technical expositions, see foundational works and systematic classifications in (Li et al., 26 Nov 2025, Turbiner, 2016, Li et al., 2023, Jr. et al., 2014, Sasaki, 2010), and references therein. Current frontiers include the extension to higher ranks, multi-variable systems, rational extensions, non-Hermitian and time-dependent models, and deeper analysis of root distributions and spectral transitions across parameter space.