Para-Krawtchouk Polynomials Overview

• Para-Krawtchouk polynomials are orthogonal polynomials defined on a bilattice, characterized by the union of two interlaced arithmetic progressions.
• They satisfy explicit three-term recurrence relations and difference equations that underpin their role in inverse spectral problems and quantum spin-chain models.
• They arise naturally from a bilattice-classical classification, linking classical discrete families to modern algebraic and q-generalizations.

Para-Krawtchouk polynomials are a finite family of orthogonal polynomials defined on a bi-lattice, introduced by Vinet and Zhedanov in the analysis of quantum spin chains with perfect state transfer and later identified as a special case of a complete bilattice-classical classification of orthogonal polynomial systems. Their defining feature is that the spectral support is not a single uniform lattice but the union of two interlaced arithmetic progressions; in this sense they are a bilattice analogue of the classical Krawtchouk family, while retaining explicit recurrence relations, orthogonality, difference equations, and algebraic structure (Vinet et al., 2011, Castillo et al., 2023).

1. Bilattice setting and defining spectral data

The bilattice framework used for para-Krawtchouk polynomials is based on

$$x(s)=s+\gamma\bigl(1+(-1)^s\bigr), \qquad s\in\mathbb C,$$

with $\gamma\in\mathbb C$. This is a two-component lattice in which even and odd values of $s$ are shifted differently by the parameter $\gamma$. In the finite spin-chain realization, the spectrum is written as

$$x_s = s + \frac{\gamma}{2}(1-\delta_s), \qquad s=0,1,\ldots,N,$$

where $N$ is odd and $0<\gamma<2$. Equivalently, the support splits into the two uniform sublattices

$x_{2s}=2s,\qquad x_{2s+1}=2s+\gamma.$

When $\gamma=1$, the two sublattices merge into the uniform lattice $x_s=s$, which recovers the classical Krawtchouk situation (Castillo et al., 2023, Vinet et al., 2011).

This bilattice geometry is the structural feature that distinguishes the family from standard discrete orthogonal polynomials on a linear lattice. It also explains why para-Krawtchouk polynomials arise naturally in inverse spectral problems for Jacobi matrices: the finite spectrum is prescribed as a union of two arithmetic progressions, and the orthogonal polynomials are then determined by the associated tridiagonal matrix.

2. Recurrence relations and orthogonality

In the classification paper, the para-Krawtchouk family is presented as the monic orthogonal polynomial sequence $\gamma\in\mathbb C$0 satisfying

$\gamma\in\mathbb C$1

with $\gamma\in\mathbb C$2 and $\gamma\in\mathbb C$3. It is identified with the bilattice family $\gamma\in\mathbb C$4 through

$\gamma\in\mathbb C$5

More generally, the $\gamma\in\mathbb C$6 family has recurrence coefficients

$\gamma\in\mathbb C$7

and

$\gamma\in\mathbb C$8

under the nondegeneracy conditions

$\gamma\in\mathbb C$9

These formulas place para-Krawtchouk polynomials within a wider parametric bilattice family rather than treating them as an isolated construction (Castillo et al., 2023).

In the spin-chain formulation, the associated monic polynomials satisfy

$s$0

with $s$1, $s$2, and $s$3. Their orthogonality relation is

$s$4

with $s$5. For the finite bi-lattice spectrum, the weights are given explicitly by

$s$6

and

$s$7

where $s$8, and they satisfy the normalization

$s$9

The recurrence data are equally explicit: $\gamma$0

$\gamma$1

together with the factorization

$\gamma$2

where

$\gamma$3

These formulas are central in applications because they determine the Jacobi matrix and hence the finite orthogonality measure (Vinet et al., 2011).

3. Bilattice-classical characterization

The modern structural description of para-Krawtchouk polynomials is formulated in terms of the bilattice analogue of the classical functional equation

$\gamma$4

where $\gamma$5 is a linear functional, $\gamma$6, $\gamma$7, and

$\gamma$8

For

$\gamma$9

and

$$x_s = s + \frac{\gamma}{2}(1-\delta_s), \qquad s=0,1,\ldots,N,$$0

the functional $$x_s = s + \frac{\gamma}{2}(1-\delta_s), \qquad s=0,1,\ldots,N,$$1 is regular if and only if

$$x_s = s + \frac{\gamma}{2}(1-\delta_s), \qquad s=0,1,\ldots,N,$$2

for each $$x_s = s + \frac{\gamma}{2}(1-\delta_s), \qquad s=0,1,\ldots,N,$$3. In this setting, regularity means that $$x_s = s + \frac{\gamma}{2}(1-\delta_s), \qquad s=0,1,\ldots,N,$$4 admits an orthogonal polynomial sequence $$x_s = s + \frac{\gamma}{2}(1-\delta_s), \qquad s=0,1,\ldots,N,$$5 such that

$$x_s = s + \frac{\gamma}{2}(1-\delta_s), \qquad s=0,1,\ldots,N,$$6

The admissibility of $$x_s = s + \frac{\gamma}{2}(1-\delta_s), \qquad s=0,1,\ldots,N,$$7 is encoded by these nonvanishing conditions (Castillo et al., 2023).

A major consequence is the bilattice Rodrigues structure. If $$x_s = s + \frac{\gamma}{2}(1-\delta_s), \qquad s=0,1,\ldots,N,$$8 is regular and satisfies the bilattice functional equation, then the monic orthogonal polynomial sequence $$x_s = s + \frac{\gamma}{2}(1-\delta_s), \qquad s=0,1,\ldots,N,$$9 associated with $N$0 satisfies

$N$1

Equivalently, for the auxiliary sequence $N$2,

$N$3

The iterative transformed functionals $N$4 remain regular for all $N$5, and the orthogonal polynomial sequence associated with $N$6 is again an orthogonal polynomial system. This places para-Krawtchouk polynomials among the families that are classical in the bilattice sense rather than merely semi-classical (Castillo et al., 2023).

A recurrent misconception is to regard para-Krawtchouk polynomials as an isolated ad hoc deformation. The bilattice functional-equation theory shows instead that they arise naturally as one admissible solution of a general classical problem on a bi-lattice.

4. Classification and relation to the Krawtchouk family

The classification theorem states that, up to an affine change of variable, the only classical orthogonal polynomial systems on the bilattice are the two families

$N$7

The standard classical discrete families on the linear lattice then appear as specializations of the bilattice theory. In particular, the paper recovers Meixner, Charlier, Krawtchouk, Hahn, and para-Krawtchouk polynomials by choosing special $N$8 and $N$9 in the bilattice functional equation. The organizing principle is:

• if $0<\gamma<2$0, one gets a Charlier-type system;
• if $0<\gamma<2$1, one gets a Meixner/Krawtchouk-type system;
• if $0<\gamma<2$2, one obtains Hahn-type and para-Krawtchouk-type systems.

The explicit identifications include

$0<\gamma<2$3

$0<\gamma<2$4

$0<\gamma<2$5

and

$0<\gamma<2$6

Within this scheme, para-Krawtchouk polynomials are not outside the classical hierarchy; they occupy a specific point in the $0<\gamma<2$7-family (Castillo et al., 2023).

The classical Krawtchouk family provides the immediate backdrop. For parameters

$0<\gamma<2$8

the Krawtchouk polynomial is defined by

$0<\gamma<2$9

where

$x_{2s}=2s,\qquad x_{2s+1}=2s+\gamma.$0

With respect to the inner product

$x_{2s}=2s,\qquad x_{2s+1}=2s+\gamma.$1

these polynomials satisfy

$x_{2s}=2s,\qquad x_{2s+1}=2s+\gamma.$2

They obey the three-term recurrence

$x_{2s}=2s,\qquad x_{2s+1}=2s+\gamma.$3

and for $x_{2s}=2s,\qquad x_{2s+1}=2s+\gamma.$4 have $x_{2s}=2s,\qquad x_{2s+1}=2s+\gamma.$5 distinct real roots in $x_{2s}=2s,\qquad x_{2s+1}=2s+\gamma.$6, with interlacing of consecutive zeros. In the binary case $x_{2s}=2s,\qquad x_{2s+1}=2s+\gamma.$7,

$x_{2s}=2s,\qquad x_{2s+1}=2s+\gamma.$8

This provides the classical finite orthogonal template from which para-Krawtchouk systems depart by replacing the single lattice with a bi-lattice (Coleman, 2011).

5. Spin-chain origin, bispectral structure, and algebraic realization

Para-Krawtchouk polynomials were introduced in the study of the XX Hamiltonian

$x_{2s}=2s,\qquad x_{2s+1}=2s+\gamma.$9

whose one-excitation restriction is a Jacobi matrix $\gamma=1$0. Perfect state transfer occurs when the spectrum satisfies the odd-spacing condition

$\gamma=1$1

and the chain is mirror-symmetric: $\gamma=1$2 For the para-Krawtchouk family, the recurrence coefficients supply explicit couplings and fields, and the perfect-transfer amplitude is

$\gamma=1$3

Under an affine rescaling

$\gamma=1$4

the recurrence coefficients transform as

$\gamma=1$5

and the perfect-state-transfer normalization imposes

$\gamma=1$6

where $\gamma=1$7 are positive coprime integers and $\gamma=1$8 is odd (Vinet et al., 2011).

The family also satisfies a second-order difference equation

$\gamma=1$9

with

$x_s=s$0

Because the shifts $x_s=s$1 preserve each sublattice, this is interpreted as a fourth-order relation on the underlying bi-lattice. When $x_s=s$2, one recovers the classical symmetric Krawtchouk case,

$x_s=s$3

and the difference equation becomes the square of the Krawtchouk eigenvalue equation. The explicit hypergeometric realization is obtained through complementary Bannai–Ito polynomials, and the characteristic polynomial factorizes as

$x_s=s$4

so that its roots are exactly the bi-lattice points (Vinet et al., 2011).

The algebraic structure behind the family is the quadratic Hahn algebra. With operators

$x_s=s$5

$x_s=s$6

and

$x_s=s$7

the commutation relations are

$x_s=s$8

with

$x_s=s$9

$\gamma\in\mathbb C$00

The Casimir operator acts as a scalar, and the eigenbasis of $\gamma\in\mathbb C$01 is precisely the para-Krawtchouk basis: $\gamma\in\mathbb C$02 This gives a direct algebraic explanation of the recurrence structure (Vinet et al., 2011).

6. $\gamma\in\mathbb C$03-generalizations and later algebraic developments

A $\gamma\in\mathbb C$04-generalization of the para-Krawtchouk family is obtained from a singular truncation of the Askey–Wilson polynomials. The truncation condition is

$\gamma\in\mathbb C$05

After regularization and passage to a finite family $\gamma\in\mathbb C$06 of $\gamma\in\mathbb C$07-para-Racah polynomials, a scaling limit produces the $\gamma\in\mathbb C$08-para-Krawtchouk polynomials. Writing

$\gamma\in\mathbb C$09

and rescaling by

$\gamma\in\mathbb C$10

the support becomes the exponential bi-lattice

$\gamma\in\mathbb C$11

The resulting family satisfies

$\gamma\in\mathbb C$12

and is orthogonal on that exponential bi-lattice. In this construction, the family is bispectral, but the eigenvalues on the $\gamma\in\mathbb C$13-difference side are doubly degenerate after finite truncation, so it is not classical in the strict Leonard-pair sense (Lemay et al., 2017).

A further development gives an algebraic interpretation through S-Heun operators on linear and $\gamma\in\mathbb C$14-linear grids. In the discrete linear case, truncation of the Continuous Hahn family occurs when

$\gamma\in\mathbb C$15

and yields the para-Krawtchouk polynomials. In the $\gamma\in\mathbb C$16-linear case, truncation of Big $\gamma\in\mathbb C$17-Jacobi occurs when

$\gamma\in\mathbb C$18

and yields the $\gamma\in\mathbb C$19-para-Krawtchouk family. The central claim is that these polynomial systems are finite-dimensional representation bases of Sklyanin-like algebras: $\gamma\in\mathbb C$20 in the discrete linear case and $\gamma\in\mathbb C$21 in the $\gamma\in\mathbb C$22-linear case. This resolves what the paper describes as a missing algebraic interpretation of para-Krawtchouk and $\gamma\in\mathbb C$23-para-Krawtchouk polynomials. The same framework also explains the subtle dimension count: the degree-$\gamma\in\mathbb C$24 polynomial annihilated by the raising operator is the characteristic polynomial of the upper block of the truncated Jacobi matrix, vanishes on the orthogonality grid, and therefore represents a null vector; the true representation dimension is $\gamma\in\mathbb C$25 (Bergeron et al., 2020).

Taken together, these developments place para-Krawtchouk polynomials at the intersection of finite orthogonal polynomial theory, inverse spectral problems for Jacobi matrices, exact spin-chain design, and nontrivial algebraic structures on discrete and $\gamma\in\mathbb C$26-discrete grids. The classification on bi-lattices shows that they belong to the classical orthogonal-polynomial landscape in a precise bilattice sense, while the spin-chain and Sklyanin-like perspectives explain why they continue to recur in finite-dimensional spectral models.