About the Project
18 Orthogonal PolynomialsApplications

§18.38 Mathematical Applications

Contents
  1. §18.38(i) Classical OP’s: Numerical Analysis
  2. §18.38(ii) Classical OP’s: Mathematical Developments and Applications
  3. §18.38(iii) Other OP’s

§18.38(i) Classical OP’s: Numerical Analysis

Approximation Theory

The monic Chebyshev polynomial 21nTn(x), n1, enjoys the ‘minimax’ property on the interval [1,1], that is, |21nTn(x)| has the least maximum value among all monic polynomials of degree n. In consequence, expansions of functions that are infinitely differentiable on [1,1] in series of Chebyshev polynomials usually converge extremely rapidly. For these results and applications in approximation theory see §3.11(ii) and Mason and Handscomb (2003, Chapter 3), Cheney (1982, p. 108), and Rivlin (1969, p. 31).

Differential Equations: Spectral Methods

Linear ordinary differential equations can be solved directly in series of Chebyshev polynomials (or other OP’s) by a method originated by Clenshaw (1957). This process has been generalized to spectral methods for solving partial differential equations. For further information see Mason and Handscomb (2003, Chapters 10 and 11), Gottlieb and Orszag (1977, pp. 7–19), and Guo (1998, pp. 120–151). See also the paragraph on DVRs, below.

Quadrature

Classical OP’s play a fundamental role in Gaussian quadrature. If the nodes in a quadrature formula with a positive weight function are chosen to be the zeros of the nth degree OP with the same weight function, and the interval of orthogonality is the same as the integration range, then the weights in the quadrature formula can be chosen in such a way that the formula is exact for all polynomials of degree not exceeding 2n1. See §3.5(v).

Quadrature “Extended” to Pseudo-Spectral (DVR) Representations of Operators in One and Many Dimensions

The basic ideas of Gaussian quadrature, and their extensions to non-classical weight functions, and the computation of the corresponding quadrature abscissas and weights, have led to discrete variable representations, or DVRs, of Sturm–Liouville and other differential operators. The terminology DVR arises as an otherwise continuous variable, such as the co-ordinate x, is replaced by its values at a finite set of zeros of appropriate OP’s resulting in expansions using functions localized at these points. Light and Carrington Jr. (2000) review and extend the one-dimensional analysis to solution of multi-dimensional many-particle systems, where the sparse nature of the resulting matrices is highly advantageous. Schneider et al. (2016) discuss DVR/Finite Element solutions of the time-dependent Schrödinger equation. Each of these typically require a particular non-classical weight functions and analysis of the corresponding OP’s. These are listed in §18.39(iii). These methods have become known as pseudo-spectral, and are overviewed in Cerjan (1993), and Shizgal (2015).

§18.38(ii) Classical OP’s: Mathematical Developments and Applications

Integrable Systems

The Toda equation provides an important model of a completely integrable system. It has elegant structures, including N-soliton solutions, Lax pairs, and Bäcklund transformations. While the Toda equation is an important model of nonlinear systems, the special functions of mathematical physics are usually regarded as solutions to linear equations. However, by using Hirota’s technique of bilinear formalism of soliton theory, Nakamura (1996) shows that a wide class of exact solutions of the Toda equation can be expressed in terms of various special functions, and in particular classical OP’s. For instance,

18.38.1 Vn(x)=2nHn+1(x)Hn1(x)(Hn(x))2,

with Hn(x) as in §18.3, satisfies the Toda equation

18.38.2 d2dx2lnVn(x)=Vn+1(x)+Vn1(x)2Vn(x),
n=1,2,.

Complex Function Theory

The Askey–Gasper inequality

18.38.3 m=0nPm(α,0)(x)=(α+2)nn!F23(n,n+α+2,12(α+1)α+1,12(α+3);12(1x))0,
x1, α2, n=0,1,,

also the case β=0 of (18.14.26), was used in de Branges’ proof of the long-standing Bieberbach conjecture concerning univalent functions on the unit disk in the complex plane. See de Branges (1985). For the generalized hypergeometric function F23 see (16.2.1).

Zonal Spherical Harmonics

Ultraspherical polynomials are zonal spherical harmonics. As such they have many applications. See, for example, Andrews et al. (1999, Chapter 9). See also §14.30.

Random Matrix Theory

Hermite polynomials (and their Freud-weight analogs (§18.32)) play an important role in random matrix theory. See Fyodorov (2005) and Deift (1998, Chapter 5).

Riemann–Hilbert Problems

See Deift (1998, Chapter 7) and Ismail (2009, Chapter 22).

Radon Transform

See Deans (1983, Chapters 4, 7).

§18.38(iii) Other OP’s

Group Representations

For group-theoretic interpretations of OP’s see Vilenkin and Klimyk (1991, 1992, 1993).

Coding Theory

For applications of Krawtchouk polynomials Kn(x;p,N) and q-Racah polynomials Rn(x;α,β,γ,δ|q) to coding theory see Bannai (1990, pp. 38–43), Leonard (1982), and Chihara (1987).

3j and 6j Symbols

The 3j symbol (34.2.6), with an alternative expression as a terminating F23 of unit argument, can be expressed in terms of Hahn polynomials (18.20.5) or, by (18.21.1), dual Hahn polynomials. The orthogonality relations in §34.3(iv) for the 3j symbols can be rewritten in terms of orthogonality relations for Hahn or dual Hahn polynomials as given by §§18.2(i), 18.2(iii) and Table 18.19.1 or by §18.25(iii), respectively. See Koornwinder (1981, §4) for details.

The 6j symbol (34.4.3), with an alternative expression as a terminating balanced F34 of unit argument, can be expressend in terms of Racah polynomials (18.26.3). The orthogonality relations (34.5.14) for the 6j symbols can be rewritten in terms of orthogonality relations for Racah polynomials as given by (18.25.9)–(18.25.12). See Wilson (1980, §5) for details.

Zhedanov Algebra

A symmetric Laurent polynomial is a function of the form

f(z)=c0+k=1nck(zk+zk).

Define operators K0 and K1 acting on symmetric Laurent polynomials by K0=L (L given by (18.28.6_2)) and (K1f)(z)=(z+z1)f(z). Define a further operator K2 by

18.38.4 K2=[K0,K1]q,

where the q-commutator of operators X and Y is defined by

18.38.5 [X,Y]q=q12XYq12YX.

Then

18.38.6 [K1,K2]q =BK1+C0K0+D0,
[K2,K0]q =BK0+C1K1+D1,

where

18.38.7 B =(1q1)2(e3+qe1),
C0 =(qq1)2,
C1 =q1(qq1)2e4,
D0 =q3(1q)2(1+q)(e4+qe2+q2),
D1 =q3(1q)2(1+q)(e1e4+qe3),

and e1,e2,e3,e4 are the elementary symmetric polynomials in a,b,c,d given by a+b+c+d, ab+ac++cd, abc+abd+acd+bcd, abcd, respectively. A further operator, the so-called Casimir operator

18.38.8 Q=(q12q32)K0K1K2+qK22+B(K0K1+K1K0)+qC0K02+q1C1K12+(1+q)D0K0+(1+q1)D1K1

commutes with K0,K1,K2, that is KjQ=QKj, and satisfies

18.38.9 Q=Q0,

where Q0 is a constant with explicit expression in terms of e1,e2,e3,e4 and q given in Koornwinder (2007a, (2.8)).

The abstract associative algebra with generators K0,K1,K2 and relations (18.38.4), (18.38.6) and with the constants B,C0,C1,D0,D1 in (18.38.6) not yet specified, is called the Zhedanov algebra or Askey–Wilson algebra AW(3). If we consider this abstract algebra with additional relation (18.38.9) and with dependence on a,b,c,d according to (18.38.7) then it is isomorphic with the algebra generated by K0=L given by (18.28.6_2), (K1f)(z)=(z+z1)f(z) and K2 given by (18.38.4), and K0,K1,K2 act on the linear span of the Askey–Wilson polynomials (18.28.1). See Zhedanov (1991), Granovskiĭ et al. (1992, §3), Koornwinder (2007a, §2) and Terwilliger (2011). Similar algebras can be associated with all families of OP’s in the q-Askey scheme and the Askey scheme.

Dunkl Type Operators and Nonsymmetric Orthogonal Polynomials

The Dunkl operator, introduced by Dunkl (1989), is an operator associated with reflection groups or root systems which has terms involving first order partial derivatives and reflection terms. Analogues of the original Dunkl operator (the rational case) were introduced by Heckman and Cherednik for the trigonometric case, and by Cherednik for the q-case. Algebraic structures were built of which special representations involve Dunkl type operators. In the q-case this algebraic structure is called the double affine Hecke algebra (DAHA), introduced by Cherednik. Eigenvalue equations involving Dunkl type operators have as eigenfunctions nonsymmetric analogues of multivariable special functions associated with root systems. This gives also new structures and results in the one-variable case, but the obtained nonsymmetric special functions can now usually be written as a linear combination of two known special functions.

In the one-variable case the Dunkl operator eigenvalue equation

18.38.10 f(x)+(α+12)f(x)f(x)x=iλf(x)

has a solution

18.38.11 f(x)=Γ(α+1)(2λx)α(Jα(λx)+iJα+1(λx)),

where the Bessel function Jν(z) is defined in (10.2.2).

For the one-variable q-case see Noumi and Stokman (2004), Koornwinder (2007a, §§3,4), Koornwinder and Bouzeffour (2011, §§4,5) and Terwilliger (2013). The Dunkl type operator is a q-difference-reflection operator acting on Laurent polynomials and its eigenfunctions, the nonsymmetric Askey–Wilson polynomials, are linear combinations of the symmetric Laurent polynomial Rn(z;a,b,c,d|q) and the ‘anti-symmetric’ Laurent polynomial z1(1az)(1bz)Rn1(z;qa,qb,c,d|q), where Rn(z) is given in (18.28.1_5). See Koornwinder (2007a, (3.13), (4.9), (4.10)) for explicit formulas.

Dunkl type operators and nonsymmetric polynomials have been associated with various other families in the Askey scheme and q-Askey scheme, in particular with Wilson polynomials, see Groenevelt (2007), and with Jacobi polynomials, see Koornwinder and Bouzeffour (2011, §7).

Supersymmetric Quantum Mechanics (SUSY)

The solved Schrödinger equations of §18.39(i) involve shape invariant potentials, and thus are in the family of supersymmetric or SUSY potentials. SUSY leads to algebraic simplifications in generating excited states, and partner potentials with closely related energy spectra, from knowledge of a single ground state wave function. These generalize the ladder operators, as reviewed and extended by Infeld and Hull (1951), and also called creation and annilhilation operators. Overviews appear in Dutt et al. (1988), Cooper et al. (1995), and Quesne (2011).

Exceptional OP’s

Exceptional OP’s (EOP’s) are those which are ‘missing’ a finite number of lower order polynomials, but yet form complete sets with respect to suitable measures. Detailed examples appear in the Section on Miscellaneous OP’s §§18.36(v) and 18.36(vi). A review is Milson (2017). Hermite EOP’s appear in solutions of a rationally modified Schrödinger equation in §18.39. Much of the exploration of the EOP’s is based on the operator algebra as developed in SUSY, above.

EOP’s, Painlevé Transcendents, and Quantum Mechanics

EOP’s are the subject of recent work on rational solutions to the fourth Painlevé equation, see Clarkson (2003a) and Marquette and Quesne (2016),where use of Hermite EOP’s makes a connection to quantum mechanics.

Non-Classical Weight Functions

See §§18.32, 18.39(iii), and 32.15, for recent developments and applications.