Hahn Difference Operator Overview

• The Hahn difference operator is a family of discrete operators that replace differentiation via finite and q-lattice shifts while preserving the structure of orthogonal polynomials.
• It features a quadratic spectrum and supports Darboux transformations, exceptional and Krall deformations, which are pivotal in discrete spectral analysis.
• Its various forms, from classical finite-lattice to q-calculus and complex analytic settings, enable versatile applications in spectral theory, orthogonal systems, and value-distribution studies.

The Hahn difference operator is a context-dependent term for several discrete operators that share a common role: they replace differentiation by finite-lattice or $q$-lattice shifts, while preserving enough structure to support orthogonal polynomials, spectral theory, and Darboux-type transformations. In the finite-lattice theory of Hahn polynomials it denotes a second-order difference operator with quadratic spectrum; in $q$-calculus it denotes the affine divided-difference operator $D_{q,\omega}$ or $\mathcal D_{q,c}$, which interpolates between Jackson $q$-differences, forward differences, and the ordinary derivative. These operators underlie exceptional and Krall deformations, Heun-type extensions, Leonard-pair and meta-Hahn algebras, discrete Dirac systems, and value-distribution theory for meromorphic functions (Durán, 2021, Hıra, 2018, Wang, 23 Sep 2025, Srivastava et al., 2019).

1. Terminology and principal definitions

The same label is used in the literature for several related operators. The distinction is determined by the lattice on which the operator acts and by the analytic problem under consideration.

Setting	Operator	Defining formula
Finite uniform lattice	Classical Hahn operator	$$L^{\alpha,\beta,N}f(x)=A(x)\Delta f(x)-C(x)\nabla f(x)$$
Affine $q$-lattice	$q,\omega$-Hahn operator	$$D_{q,\omega}f(t)=\dfrac{f(qt+\omega)-f(t)}{(q-1)t+\omega}$$
Complex affine $q$-shift	Hahn operator $q$0	$q$1
Homogeneous $q$2-difference setting	$q$3	$q$4

For the classical finite-lattice operator, $q$5 and $q$6. For the $q$7 operator one sets $q$8 and defines $q$9 when $D_{q,\omega}$0 is differentiable there. For the homogeneous operator, the divided-difference kernel is

$D_{q,\omega}$1

Accordingly, the phrase “Hahn difference operator” does not designate a single universal object; it designates a family of discrete differential analogues tied to Hahn-type special functions and their generalizations (Durán, 2021, Hıra, 2018, Wang, 23 Sep 2025, Srivastava et al., 2019).

2. Classical finite-lattice Hahn operator

Let $D_{q,\omega}$2 and $D_{q,\omega}$3 be a nonnegative integer. With forward and backward shifts $D_{q,\omega}$4, $D_{q,\omega}$5, the classical Hahn difference operator is

$D_{q,\omega}$6

where

$D_{q,\omega}$7

Equivalently,

$D_{q,\omega}$8

The monic Hahn polynomials $D_{q,\omega}$9, $\mathcal D_{q,c}$0, are eigenfunctions: $\mathcal D_{q,c}$1 This quadratic spectrum is one of the defining structural features of the Hahn family. Closely related realizations appear in the algebraic literature as $\mathcal D_{q,c}$2 or $\mathcal D_{q,c}$3, written in shift form $\mathcal D_{q,c}$4, with the same quadratic eigenvalue law up to conventional sign and parameter choices (Durán, 2021, Vinet et al., 2020, Crampé et al., 19 May 2026).

In a standard hypergeometric normalization, the Hahn polynomials are orthogonal on the finite grid $\mathcal D_{q,c}$5 with Hahn weight

$\mathcal D_{q,c}$6

and satisfy

$\mathcal D_{q,c}$7

This finite-lattice self-adjoint framework is the base object from which exceptional, Krall, and Heun-type Hahn operators are constructed (Durán et al., 2014).

The same operator also has a bispectral interpretation. In the meta-Hahn and trio-Hahn formalisms, the Hahn difference operator acts in the $\mathcal D_{q,c}$8-variable while a dual operator acts in the degree index $\mathcal D_{q,c}$9, yielding the standard pair of spectral relations: a second-order difference equation in $q$0 and a three-term recurrence in $q$1 (Vinet et al., 2020, Crampé et al., 19 May 2026).

3. Exceptional Hahn operators

Exceptional Hahn polynomials are orthogonal polynomial systems that remain eigenfunctions of a second-order difference operator, but with gaps in the degree sequence. The absence of certain degrees is the main visible distinction between exceptional families and classical discrete orthogonal polynomials (Durán, 2021).

A 2021 construction introduces a finite exceptional index set $q$2, a family of real continuous parameters $q$3, and a seed polynomial

$q$4

of degree $q$5. The exceptional Hahn operator is then obtained by gauging the classical operator with $q$6: $q$7 where

$q$8

$q$9

$$L^{\alpha,\beta,N}f(x)=A(x)\Delta f(x)-C(x)\nabla f(x)$$0

with $$L^{\alpha,\beta,N}f(x)=A(x)\Delta f(x)-C(x)\nabla f(x)$$1 and $$L^{\alpha,\beta,N}f(x)=A(x)\Delta f(x)-C(x)\nabla f(x)$$2. This is a Darboux or gauge transformation of the classical Hahn operator, since the coefficients differ by rational factors $$L^{\alpha,\beta,N}f(x)=A(x)\Delta f(x)-C(x)\nabla f(x)$$3 (Durán, 2021).

The spectrum retains the same quadratic eigenvalues,

$$L^{\alpha,\beta,N}f(x)=A(x)\Delta f(x)-C(x)\nabla f(x)$$4

but only for the admissible index set

$$L^{\alpha,\beta,N}f(x)=A(x)\Delta f(x)-C(x)\nabla f(x)$$5

Hence a finite set of degrees is missing. Under an admissibility condition on $$L^{\alpha,\beta,N}f(x)=A(x)\Delta f(x)-C(x)\nabla f(x)$$6, the operator $$L^{\alpha,\beta,N}f(x)=A(x)\Delta f(x)-C(x)\nabla f(x)$$7 is symmetric for the discrete inner product

$$L^{\alpha,\beta,N}f(x)=A(x)\Delta f(x)-C(x)\nabla f(x)$$8

with positive weight $$L^{\alpha,\beta,N}f(x)=A(x)\Delta f(x)-C(x)\nabla f(x)$$9, and the exceptional polynomials $q$0, $q$1, are orthogonal and form a basis of $q$2 (Durán, 2021).

An earlier construction uses Casorati determinants indexed by a pair $q$3. The allowed degrees are

$q$4

and the exceptional polynomials $q$5 satisfy a second-order difference equation

$q$6

with rational coefficients expressed through Casorati determinants $q$7 and $q$8. Under Hahn-admissibility they are orthogonal and complete for a positive discrete measure. In the limit $q$9, the Casorati construction passes to a Wronskian construction for exceptional Jacobi polynomials (Durán, 2015).

4. Higher-order, Heun, and algebraic extensions

Beyond second-order operators, the Hahn setting supports higher-order difference operators whose eigenfunctions remain orthogonal. In the $q,\omega$0-operator method, one starts with a polynomial sequence $q,\omega$1, an algebra $q,\omega$2 of finite-order difference operators, and an operator $q,\omega$3 with $q,\omega$4. One then constructs

$q,\omega$5

and arranges the coefficients so that the new family is again an eigenbasis: $q,\omega$6 For Hahn polynomials, this yields Krall-Hahn families orthogonal with respect to Christoffel-transformed Hahn weights and satisfying higher-order difference equations. The operator has the form

$q,\omega$7

and, under nondegeneracy of the relevant Casorati determinants, the eigenvalues $q,\omega$8 become explicit polynomials of degree $q,\omega$9 in $$D_{q,\omega}f(t)=\dfrac{f(qt+\omega)-f(t)}{(q-1)t+\omega}$$0 (Durán et al., 2014).

A standard Krall-Hahn weight is

$$D_{q,\omega}f(t)=\dfrac{f(qt+\omega)-f(t)}{(q-1)t+\omega}$$1

For this weight there exists a family $$D_{q,\omega}f(t)=\dfrac{f(qt+\omega)-f(t)}{(q-1)t+\omega}$$2 of orthogonal polynomials and a single difference operator of order $$D_{q,\omega}f(t)=\dfrac{f(qt+\omega)-f(t)}{(q-1)t+\omega}$$3 such that

$$D_{q,\omega}f(t)=\dfrac{f(qt+\omega)-f(t)}{(q-1)t+\omega}$$4

In the symmetric case

$$D_{q,\omega}f(t)=\dfrac{f(qt+\omega)-f(t)}{(q-1)t+\omega}$$5

the operator order drops to

$$D_{q,\omega}f(t)=\dfrac{f(qt+\omega)-f(t)}{(q-1)t+\omega}$$6

which is smaller than the generic order (Durán et al., 2018).

A different higher-order deformation is the Heun-Hahn operator on the uniform grid $$D_{q,\omega}f(t)=\dfrac{f(qt+\omega)-f(t)}{(q-1)t+\omega}$$7. It is the most general second-order difference operator that maps polynomials of degree $$D_{q,\omega}f(t)=\dfrac{f(qt+\omega)-f(t)}{(q-1)t+\omega}$$8 to polynomials of degree $$D_{q,\omega}f(t)=\dfrac{f(qt+\omega)-f(t)}{(q-1)t+\omega}$$9. It is tridiagonal both in the Pochhammer basis and in the Hahn-polynomial basis, and it is bilinear in the generators $q$0 and $q$1 of the Hahn algebra: $q$2 Adjoining $q$3 to the Hahn algebra yields an extended Heun-Racah algebra with cubic commutation relations (Vinet et al., 2018).

The finite-lattice Hahn operator also sits naturally inside algebraic bispectral frameworks. In the meta-Hahn algebra $q$4, one embeds the Hahn algebra by

$q$5

and $q$6 is represented by the Hahn difference operator up to a scalar shift. The pair $q$7 forms a Leonard pair: one operator is diagonal in the basis in which the other is tridiagonal, and vice versa (Vinet et al., 2020). In the trio-Hahn formulation, the generators $q$8, $q$9, and $q$00 satisfy

$q$01

and the overlaps of the associated eigenbases recover Hahn polynomials and biorthogonal rational functions together with their bispectrality and orthogonality (Crampé et al., 19 May 2026).

5. The $q$02-Hahn operator and its calculus

For $q$03, $q$04, and $q$05, the $q$06-Hahn operator is

$q$07

If $q$08 with $q$09 fixed, then $q$10 tends to the forward finite difference $q$11. If $q$12, then

$q$13

the Jackson $q$14-difference operator. As $q$15 and $q$16, the operator tends to the ordinary derivative. In this precise sense, $q$17 unifies the forward difference, Jackson $q$18-difference, and differentiation (Hıra, 2018, Filipuk et al., 2020).

The operator satisfies linearity and a shifted product rule,

$q$19

It also supports a Jackson-Nörlund integral and a fundamental theorem of $q$20-calculus: $q$21 These identities make $q$22 a viable derivative substitute in spectral and variational constructions (Hıra, 2018).

One application is the $q$23-Dirac system

$q$24

with boundary conditions at $q$25 and $q$26, where $q$27. In the corresponding Hilbert space $q$28, existence and uniqueness hold for the initial-value problem, the $q$29-Wronskian is constant on solutions, and all eigenvalues are real and simple. A characteristic entire function $q$30 determines the spectrum, and in model cases one obtains explicit eigenvalue asymptotics through the zeros of the Hahn-cosine and Hahn-sine functions (Hıra, 2018).

A second application appears in discrete Sobolev orthogonality. For

$q$31

the associated monic Sobolev orthogonal polynomials $q$32 admit ladder operators

$q$33

satisfying lowering and raising relations. Combining them yields a second-order equation

$q$34

This framework includes the cases $q$35, $q$36, and $q$37 as special or limiting cases (Filipuk et al., 2020).

6. Complex-analytic and homogeneous $q$38-difference extensions

For meromorphic function theory in the complex plane, the Hahn operator is

$q$39

Its limits are

$q$40

For nonconstant meromorphic functions of zero order and $q$41, a Hahn-version logarithmic derivative lemma holds: $q$42 outside an exceptional set of logarithmic density zero. The same setting supports a Hahn second fundamental theorem, modified counting functions $q$43, a deficiency relation, a Hahn-Picard theorem stating that at most two values can be omitted in the Hahn-Picard sense, and a five-value theorem asserting uniqueness for two zero-order meromorphic functions that Hahn-share five distinct values (Wang, 23 Sep 2025).

The operator also governs growth questions for complex difference equations. For entire solutions of

$q$44

if

$q$45

then every nontrivial entire solution satisfies

$q$46

For the Hahn-Fermat equation

$q$47

there are no nonconstant meromorphic solutions of finite order (Wang, 23 Sep 2025).

A separate homogeneous $q$48-difference development defines

$q$49

At $q$50, it reduces to the earlier one-parameter homogeneous $q$51-difference operator, and at $q$52 it reproduces the classical Hahn difference operator in homogeneous form. The generalized Hahn polynomials are given by the Rodrigues-type formula

$q$53

with $q$54 the Cauchy polynomials. This operator yields ordinary generating functions, extended generating functions, Rogers-type bilinear identities, and Mehler’s formula for generalized Hahn polynomials (Srivastava et al., 2019). A trivariate extension,

$q$55

is governed by a common homogeneous $q$56-difference equation in $q$57, from which extended generating and Rogers-type formulas follow (Arjika et al., 2021).

Across these finite-lattice, affine $q$58-lattice, complex-analytic, and homogeneous $q$59-difference settings, the Hahn difference operator functions less as a single formula than as a structural archetype: a shift-based operator that preserves polynomial or meromorphic structure, carries a tractable spectral theory, and admits systematic deformation by Darboux, Christoffel-Geronimus, and algebraic methods (Durán, 2021, Wang, 23 Sep 2025).