[Paper Review] Dual -1 Hahn polynomials: "classical" polynomials beyond the Leonard duality

Satoshi Tsujimoto, Luc Vinet+1 more|arXiv (Cornell University)|Jul 31, 2011

Molecular spectroscopy and chirality3 citations

TL;DR

This paper introduces dual -1 Hahn polynomials as a q→-1 limit of dual q-Hahn polynomials, establishing them as orthogonal polynomials on a finite discrete grid that do not satisfy Leonard duality. Unlike classical families, they satisfy a 4th-order difference eigenvalue equation, leading to a generalized Leonard pair where one matrix is 5-diagonal in the eigenbasis of the other, revealing a new class of 'classical' orthogonal polynomials beyond the standard duality framework.

ABSTRACT

We introduce the -1 dual Hahn polynomials through an appropriate $q o -1$ limit of the dual q-Hahn polynomials. These polynomials are orthogonal on a finite set of discrete points on the real axis, but in contrast to the classical orthogonal polynomials of the Askey scheme, the -1 dual Hahn polynomials do not exhibit the Leonard duality property. Instead, these polynomials satisfy a 4-th order difference eigenvalue equation and thus possess a bispectrality property. The corresponding generalized Leonard pair consists of two matrices $A,B$ each of size $N+1 imes N+1$. In the eigenbasis where the matrix $A$ is diagonal, the matrix $B$ is 3-diagonal; but in the eigenbasis where the matrix $B$ is diagonal, the matrix $A$ is 5-diagonal.

Motivation & Objective

• To define and characterize a new family of orthogonal polynomials arising from the q→-1 limit of dual q-Hahn polynomials.
• To investigate whether these polynomials retain the Leonard duality property seen in classical orthogonal polynomials like Bannai-Ito polynomials.
• To derive a difference eigenvalue equation for the new polynomials and analyze its structure.
• To establish the bispectrality of the dual -1 Hahn polynomials through a 5-term recurrence and a 4th-order difference operator.

Proposed method

• The dual -1 Hahn polynomials are constructed as a limiting case of dual q-Hahn polynomials in the q→-1 regime.
• The polynomials are shown to be orthogonal on a finite discrete set of points on the real line, with explicit weight functions derived.
• A second-order Dunkl-type shift operator H is derived that acts as a 4th-order difference operator on the polynomial eigenfunctions.
• The operator H is expressed in terms of shift operators T^±j and the reflection operator R, preserving polynomial degree.
• The eigenvalue equation HR_n^{(-1)}(x) = 2n R_n^{(-1)}(x) is established, with H explicitly given in terms of rational functions of x.
• The structure of the generalized Leonard pair is analyzed: in the A-eigenbasis, B is 3-diagonal; in the B-eigenbasis, A is 5-diagonal.

Experimental results

Research questions

• RQ1 Do the dual -1 Hahn polynomials inherit the Leonard duality property from the dual q-Hahn polynomials in the q→-1 limit?
• RQ2 What is the structure of the difference eigenvalue equation satisfied by the dual -1 Hahn polynomials?
• RQ3 How does the matrix representation of the generalized Leonard pair differ from classical cases in terms of diagonal structure?
• RQ4 Can the dual -1 Hahn polynomials be described via a second-order Dunkl-type shift operator?
• RQ5 What is the explicit form of the 5-term recurrence or difference equation governing the dual -1 Hahn polynomials?

Key findings

• The dual -1 Hahn polynomials are orthogonal on a finite discrete set of N+1 points on the real line, with explicit weight functions derived from the q→-1 limit.
• These polynomials do not satisfy the Leonard duality property, as their dual eigenvalue equation is of 4th order rather than 2nd order.
• They satisfy a 4th-order difference eigenvalue equation of the form HR_n^{(-1)}(x) = 2n R_n^{(-1)}(x), where H is a second-order Dunkl shift operator.
• The generalized Leonard pair consists of two (N+1)×(N+1) matrices A and B: A is 3-diagonal in its eigenbasis, while B is 5-diagonal in A's eigenbasis.
• The operator H is explicitly constructed as H = E1(x)T^4 + E2(x)T^{-4} + G1(x)T^2R + G2(x)T^{-2}R - (E1+E2+G1+G2)I, with rational coefficients depending on parameters α, β, and N.
• For even and odd N, the coefficients E1, E2, G1, G2 are given in closed-form rational functions of x, α, β, and N, confirming the operator's polynomial-preserving property.