Entropy-Like Properties and ${\mathfrak{L}}_q$-Norms of Hypergeometric Orthogonal Polynomials: Degree Asymptotics

Abstract

In this work, the spread of hypergeometric orthogonal polynomials (HOPs) along their orthogonality interval is examined by means of the main entropy-like measures of their associated Rakhmanov’s probability density—so, far beyond the standard deviation and its generalizations, the ordinary moments. The Fisher information, the Rényi and Shannon entropies, and their corresponding spreading lengths are analytically expressed in terms of the degree and the parameter(s) of the orthogonality weight function. These entropic quantities are closely related to the gradient functional (Fisher) and the ${\mathfrak{L}}_q$-norms (Rényi, Shannon) of the polynomials. In addition, the degree asymptotics for these entropy-like functionals of the three canonical families of HPOs (i.e., Hermite, Laguerre, and Jacobi polynomials) are given and briefly discussed. Finally, a number of open related issues are identified whose solutions are both physico-mathematically and computationally relevant.

1. Introduction

Hypergeometric orthogonal polynomials (HOPs), also called Shohat–Favard polynomials and classical orthogonal polynomials, play a key role in the development of the theory of special functions, and they are instrumental in numerous scientific problems, ranging from approximation theory to quantum theory and mathematical physics [1,2,3,4,5,6,7,8]. Here, we examine and review the knowledge of the spreading measures of the HOPs in a real continuous variable, $\left\{p_n\left(x\right)\right\}$, orthogonal with respect to the weight function $h\left(x\right)$ on the interval $\Lambda \subseteq \mathbb{R}$. These quantities measure the different spreading facets of their associated probability density as follows:

$${\rho }_n\left(x\right)={\widehat{p}}_n^2\left(x\right)h\left(x\right)=\frac{1}{{\kappa }_n}{p}_n^2\left(x\right)h\left(x\right), \mathrm{with} {\kappa }_n={\int }_{\Lambda }{\left|p_n\left(x\right)\right|}^{2}h\left(x\right)dx,$$

(1)

where the symbol ${\widehat{p}}_n\left(x\right)=p_n\left(x\right)/{\kappa }_n^{\frac{1}{2}}$ denotes the orthonormal polynomial. This quantity is called Rakhmanov’s density [9] because he first discovered it to govern the asymptotical ($n\to +\infty$) behavior of the ratio of two polynomials with consecutive orders. Later on, it was realized that this normalized-to-unity probability density characterizes the Born’s probability density of the bound stationary states of a great deal of quantum systems in one and many dimensions (harmonic oscillator, hydrogenic atoms, etc.) [3,10,11]. This is because the physical solutions (wave functions) of the Schrödinger equation for a broad family of quantum potentials, such as the spherically symmetric ones [3,12,13], are controlled by the three canonical HOPs families of Hermite $H_n\left(x\right)$, Laguerre $L_n^{\left(\alpha \right)}\left(x\right)$ and Jacobi $P_n^{(\alpha ,\beta )}\left(x\right)$ types, respectively. So, the spreading measures of the HOPs quantify numerous physical and chemical properties of quantum systems, such as the physical entropy-like measures, which are the basic variables of the classical and quantum information theories [14,15,16]. These spreading measures include standard deviation (and its generalizations, the ordinary moments) and the entropy-like measures (Fisher information [17,18], Shannon entropy [14,19] and Rényi entropies [20,21]) and ${\mathfrak{L}}_q$-norms of the density ${\rho }_n\left(x\right)$. In order for all these quantifiers to be mutually compared on the same footing, we have examined the corresponding information-theoretic lengths of Fisher, Shannon and Rényi types [22,23,24,25]. These quantities, together with the standard deviation, are direct spreading measures [26] of Rakhmanov’s density because they have some common properties: namely, the same units as the variable x, invariance under translations and reflections, linear scaling with x and vanishing in the limit as ${\rho }_n\left(x\right)$ approaches a delta function. Note that, contrary to the standard deviation, which measures the concentration of the probability density around a particular point (the centroid) of the orthogonality interval, the information-theoretic lengths do not refer to any specific point of that interval.

The Fisher information $F\left(\rho \right)$ is a functional of the derivative of the density $\rho \left(x\right)$; so, it is sensitive to the fluctuations of the density. The Fisher information controls the localization of the density around its nodes, appropriately grasping the oscillatory nature of the density; this confers it a relevant role in the characterization of a great diversity of scientific phenomena [18]. It is a local spreading measure of the density, which quantifies the pointwise concentration of $\rho$. Moreover, the higher this quantity, the more localized the density, and the higher the accuracy in predicting the localization of the particle.

The Rényi entropies $R_q\left[\rho \right],0<q<\infty ,q\ne 1,$ conform a family of spreading measures of the density $\rho \left(x\right)$, depending on a real parameter q. They are q-power functionals of the density, closely related to the ${\mathfrak{L}}_q$-norms, so they have a global character, contrary to the previous Fisher information. According to the values of the parameter q, these quantities supply complementary ways to quantify the different facets of the extent/shape of $\rho \left(x\right)$ all over the orthogonality interval, including the Onicescu functional [27] or disequilibrium $\mathcal{D}\left[\rho \right]=exp(-R_2\left[\rho \right])$ and the Shannon entropy $S\left[\rho \right]={lim}_{q\to 1}{R}_q\left[\rho \right]$. These quantities satisfy a large number of interesting physico-mathematical properties [14,28,29,30,31]. In particular, the Shannon entropy is the only one that satisfies all the hypotheses of Shannon’s theorem of information theory [14,19] as well as some other important criteria [32]; for instance, Shannon entropy becomes thermodynamical entropy in the case of a thermal ensemble. It is worth noting that $S\left[\rho \right]$ can have any values in $\left[-\infty ,\infty \right]$, contrary to differential Shannon entropy, $-{\sum }_i{p}_{i}ln{p}_i$, of a probability on a discrete sample space, which is always positive. Moreover, any sharp peaks in $\rho \left(x\right)$ tends to make $S\left[\rho \right]$ negative, whereas a slowly decaying density’s tail provokes positive values for $S\left[\rho \right]$; hence, the Shannon entropy $S\left[\rho \right]$ estimates the total extent of the density $\rho$.

In this work, we update the analytical determination of the spreading measures of the HOPs with emphasis on the entropy-like quantities and ${\mathfrak{L}}_q$-norms because of their relevance in the information theory of special functions and quantum systems, and to facilitate its numerical and symbolic computation. We should keep in mind that the naive numerical evaluation of these quantities using quadratures is often not convenient due to the increasing number of integrable singularities when the polynomial degree n is increasing, which spoils any attempt to achieve reasonable accuracy, even for rather small n. In addition, the asymptotics ($n\to \infty$) of these polynomial functionals has also been considered, although briefly discussed. The degree asympotics of HOPs and its generalizations was initiated in the middle of the nineties by A.I. Aptekarev, J.S. Dehesa, W. van Assche and their collaborators [33,34,35] and reviewed with some physical and mathematical applications in 2001 [11] and 2010 [36], respectively. The q-asymptotics and the weight-function’s parameter asymptotics for unweighted and weighted ${\mathfrak{L}}_q$-norms are not considered here. The degree asymptotics was recently used to evaluate physical Rényi and Shannon entropies for the highly-excited (Rydberg) states of quantum systems of harmonic (oscillator-like) [37,38] and Coulombic (hydrogenic-like) types [39,40,41] as well as for Rydberg atoms, which are used as building elements of logical gates in quantum computation.

Here, we do not consider the entropy-like measures of polynomials with varying weights (i.e., polynomials whose weight-function’s parameter does depend on the polynomial degree), which are also of great mathematical and physical interest [42,43], nor do we consider another kind of entropy, the discrete Shannon entropy of HOPs [44] introduced in 2009, which was explicitly calculated for Chebyshev polynomials and whose degree asymptotics for Jacobi polynomials was determined in 2015 [45] (see also [46] for the discrete Tsallis and Rényi entropies).

This paper has the following structure. In Section 2, we give the spreading measures for an arbitrary probability density $\rho \left(x\right)$ of dispersion (standard deviation and ordinary moments) and information-theoretic (Fisher, Shannon and Rényi) types, which are considered here. In Section 3, the moments around the origin and the standard deviation for the three canonical families of real HOPs (namely, Hermite, Laguerre and Jacobi) are explicitly calculated. In Section 4 and Section 5, we determine the Fisher and Rényi (and the associated weighted ${\mathfrak{L}}_q$-norms) spreading measures of the HOPs for all degree n, respectively. In Section 6, Section 7 and Section 8, we calculate the degree asymptotics of the weighted and generalized weighted ${\mathfrak{L}}_q$-norms for Hermite, Laguerre and general (i.e., Szegö type) orthonormal polynomials and Jacobi polynomials, respectively. In Section 9 and Section 10 we show the Shannon spreading measures of the HOPs and their degree asymptotics via the appropriate (unweighted, weighted) ${\mathfrak{L}}_q$-norms, respectively. Finally, some concluding remarks are pointed out and a number of open related issues are identified.

2. Spreading Measures of a Probability Density

In this section, we describe the spreading measures of dispersion (ordinary moments, standard deviation) and information-theoretic (Fisher information, Shannon entropy, Rényi entropies) types for a random variable X characterized by the continuous probability density $\rho \left(x\right)$, $x\in \Lambda \subseteq \mathbb{R}$.

To estimate the spread of X over the interval $\Lambda$, we often take a familiar dispersion measure [47] related to the moments around a particular point of $\Lambda$; usually, one chooses the origin (ordinary moments or moments around the origin $\langle x^k\rangle$) or the centroid $\langle x\rangle$ (central moments or moments around the centroid $\langle {(x-\langle x\rangle )}^k\rangle$), or we use some information-theoretic quantities connected with the frequency or entropic moments $W_q\left[\rho \right]=\langle {\left[\rho \left(x\right)\right]}^{q-1}\rangle ={\parallel {\rho }_n\parallel }_q^q$ [47,48,49], which do not depend on any particular point of the interval and are closely related to the ${\mathfrak{L}}_q$-norms of the density. Each set of moments determines the probability density $\rho \left(x\right)$ under certain conditions as the associated ordinary [50] and entropic moment problems state. The expectation value of $f\left(x\right)$ with respect to the (normalized-to-unity) density $\rho \left(x\right)$ is given by $\langle f\left(x\right)\rangle ={\int }_{\Lambda }f\left(x\right)\rho \left(x\right)dx$. The most familiar dispersion measure is the statistical root-mean-square or standard deviation $\Delta x$, which is the square root of the variance:

$$V\left[\rho \right]={\left(\Delta x\right)}^2=\langle x^2\rangle -{\langle x\rangle }^2,$$

The entropic moments quantify the extent/shape to which the probability is, in fact, distributed. They are, at times, much better probability quantifiers than the ordinary moments [48]; moreover, they are fairly efficient in the range where the ordinary moments are fairly inefficient [51]. Two relevant spreading measures related to the entropic moments are the entropy-like measures of Rényi and Shannon types. The Rényi entropies $R_q\left[\rho \right]$ of $\rho \left(x\right)$ are defined as follows [20,21]:

$$R_q\left[\rho \right]=\frac{1}{1-q}ln{W}_q\left[\rho \right]=\frac{1}{1-q}ln{\int }_{\Lambda }{\left[\rho \left(x\right)\right]}^{q}dx,q>0,q\ne 1,$$

(2)

The Shannon entropy [14,19] is given by the limiting value $q\to 1$, taking into account the normalization condition $W_1\left[\rho \right]=1$ as follows:

$$S\left[\rho \right]=\underset{q\to 1}{lim}{R}_q\left[\rho \right]=-{\int }_{\Lambda }\rho \left(x\right)ln\rho \left(x\right)dx.$$

(3)

The Rényi parameter q allows to enhance or diminish the contribution of the integrand over different regions to the whole integral in (2). Higher values of q make the function ${\left[\rho \left(x\right)\right]}^q$ to concentrate around the local maxima of the distribution, while the lower values have the effect of smoothing that function over its whole domain. It is in this sense that the parameter q provides a powerful tool in order to obtain information on the configuration shape of the probability density by means of the Rényi entropies. Moreover, let us mention here the monotonicity relations given by the following:

$$R_p\left[\rho \right]\ge R_q\left[\rho \right],\mathrm{if}p\le q;\mathrm{and}\frac{p-1}{p}{R}_p\left[\rho \right]\ge \frac{q-1}{q}{R}_q\left[\rho \right],\mathrm{if}p\ge q>1$$

(4)

which, in particular, allow one to lowerbound all the Rényi entropies by means of the second-order entropy as $R_q\left[\rho \right]\ge \frac{1}{2}{R}_2\left[\rho \right],\mathrm{for}q>0$.

A third remarkable, qualitatively different, entropy-like spreading measure of $\rho \left(x\right)$ is the Fisher information [17,18] which is defined as follows:

$$F\left[\rho \right]:={\int }_{\Lambda }\frac{{\left[{\rho }^{\prime }\left(x\right)\right]}^2}{\rho \left(x\right)}dx.$$

(5)

Opposite to the previous dispersion and entropy-like entropies, which have a global character because they are power and logarithmic functionals of the density, the Fisher information has a locality property because it is a functional of the derivative of $\rho \left(x\right)$. Moreover, these three information-theoretic spreading measures (Shannon, Rényi, Fisher) do not depend on any particular point of the interval $\Lambda$, contrary to the standard deviation. However, they have different units from the standard deviation (i.e., the units of X) so that they can not be mutually compared. To overcome this difficulty, the following information-theoretic lengths were introduced [26]:

$${\mathcal{L}}_q^R\left[\rho \right]=e^{R_q\left[\rho \right]}={\left(W_q\left[\rho \right]\right)}^{\frac{1}{1-q}}={\left\{{\int }_{\Lambda }{\left[\rho \left(x\right)\right]}^{q}dx\right\}}^{\frac{1}{1-q}},\mathrm{R}\acute{\mathrm{e}}\mathrm{nyi}\mathrm{length},$$

(6)

$${\mathcal{L}}_1^S\left[\rho \right]=\underset{q\to 1}{lim}{\mathcal{L}}_q^R\left[\rho \right]=exp\left(S\left[\rho \right]\right)=exp\left\{-{\int }_{\Lambda }\rho \left(x\right)ln\rho \left(x\right)dx\right\},\mathrm{Shannon}\mathrm{length},$$

(7)

$$\delta x=\frac{1}{\sqrt{F\left[\rho \right]}}={\left\{{\int }_{\Lambda }\frac{{\left[{\rho }^{\prime }\left(x\right)\right]}^2}{\rho \left(x\right)}dx\right\}}^{-\frac{1}{2}},\mathrm{Fisher}\mathrm{length}.$$

(8)

We remark that the quantities ($V\left[\rho \right]$, $R_q\left[\rho \right]$, $S\left[\rho \right]$, $F\left[\rho \right]$), and its related spreading lengths ($\Delta x$, ${\mathcal{L}}_q^R\left[\rho \right]$, ${\mathcal{L}}_1^S\left[\rho \right]$, $\delta x$), are complementary since each of them grasps a single different facet of the probability density $\rho \left(x\right)$. So, the variance measures the concentration of the density around the centroid, while the Rényi and Shannon entropies are measures of the configurational aspects of the extent/shape to which the density is, in fact, concentrated, and the Fisher information is a quantifier of the oscillatory character of the density since it estimates the pointwise concentration of the probability over its support $\Lambda$. All spreading lengths share the following properties [26,52,53]: same units as the random variable, translation, reflection invariance and linear scaling under adequate boundary conditions, vanishing when the density tends to a delta function. Moreover, they fulfill an uncertainty property [24,54,55,56] and the Cramér-Rao [52] and Shannon [19] inequalities given by the following:

$$\delta x\le \Delta x,\mathrm{and}{\mathcal{L}}_1^S\left[\rho \right]\le {\left(2\pi e\right)}^{\frac{1}{2}}\Delta x,$$

(9)

respectively.

In the next sections, we determine the previous spreading measures for the Rakhmanov density (1) of the real hypergeometric polynomials $\left\{p_n\left(x\right)\right\}$, orthogonal with respect to the weight function $h\left(x\right)$ on the interval $\Lambda$ so that the following holds:

$${\int }_{\Lambda }{p}_n\left(x\right)p_m\left(x\right)h\left(x\right)dx={\kappa }_n{\delta }_{n,m},\mathrm{deg}{p}_n=n$$

(10)

where the weight function has the following expressions:

$$\begin{array}{cc}\hfill h_H\left(x\right)& =e^{-x^2},x\in (-\infty ,+\infty )\hfill \\ \hfill h_{\alpha }^L\left(x\right)& =x^{\alpha }{e}^{-x},(\alpha >-1),x\in [0,+\infty )\hfill \\ \hfill h_{\alpha ,\beta }^J\left(x\right)& ={(1-x)}^{\alpha }{(1+x)}^{\beta },(\alpha ,\beta >-1),x\in [-1,+1]\hfill \end{array}$$

(11)

for the HOPs families of Hermite $H_n\left(x\right)$, Laguerre $L_n^{\left(\alpha \right)}\left(x\right)$ and Jacobi $P_n^{(\alpha ,\beta )}\left(x\right)$ types, respectively. The corresponding normalization constants are as follows:

$${\kappa }_n^H\equiv d_n^2=\sqrt{\pi }n!2^n,$$

(12)

$${\kappa }_n^L\equiv d_{n,\alpha }^2=\Gamma (n+\alpha +1)/n!,$$

(13)

$${\kappa }_n^J\equiv d_{n,\alpha ,\beta }^2=\frac{2^{\alpha +\beta +1}\Gamma (\alpha +n+1)\Gamma (\beta +n+1)}{n!(\alpha +\beta +2n+1)\Gamma (\alpha +\beta +n+1)},$$

(14)

respectively. Note that ${\kappa }_n=1$ for the orthonormal polynomials ${\widehat{p}}_n\left(x\right)$ of Hermite ${\widehat{H}}_n\left(x\right)$, Laguerre ${\widehat{L}}_n^{\left(\alpha \right)}\left(x\right)$ and Jacobi ${\widehat{P}}_n^{(\alpha ,\beta )}\left(x\right)$ types.

3. Ordinary Moments and Standard Deviation of HOPs

In this section, we show the ordinary moments and the standard deviation for the three canonical families of the hypergeometric orthonormal polynomials defined by the expressions (10) and (11). They are given by the corresponding quantities of the associated Rakhmanov density given by (1):

$${\langle x^k\rangle }_n={\int }_a^b{x}^k{\rho }_n\left(x\right)dx=\frac{1}{{\kappa }_n}{\int }_a^b{x}^k{p}_n^2\left(x\right)h\left(x\right)dxk=0,1,2,...$$

(15)

for the ordinary moments, and

$${(\Delta x)}_n={\left({\langle x^2\rangle }_n-{\langle x\rangle }_n^2\right)}^{\frac{1}{2}},$$

(16)

for the standard deviation. These quantities can be obtained from the different characterizations [3,4,6] of the HOPs, such as the explicit expression [22,23,24], the orthogonality relation (10), the three-term recurrence relation [57] and the second-order differential equation [3,58]. Briefly, for Hermite polynomials these quantities have the following values:

$${\langle x^k\rangle }_n=\left\{\begin{array}{cc}\frac{k!}{2^k\Gamma \left(\frac{k}{2}+1\right)}{}_2{F}_1\left(\left.\begin{array}{c}-n,-\frac{k}{2}\\ 1\end{array}\right|2\right),\hfill & \mathrm{even}k\hfill \\ 0,\hfill & \mathrm{odd}k\hfill \end{array}\right.,$$

(17)

where ${}_2{F}_1(a,b;c;2)$ denotes the Gauss hypergeometric function ${}_2{F}_1(a,b;c;x)$ [6] evaluated at $x=2$ for the ordinary moments, and

$${(\Delta x)}_n=\sqrt{n+\frac{1}{2}},$$

(18)

for the standard deviation of the orthonormal Hermite polynomials. For Laguerre polynomials the following values

$${\langle x^k\rangle }_{n,\alpha }=\frac{n!\Gamma (k+\alpha +1)}{\Gamma (n+\alpha +1)}\sum _{r=0}^n{\left(\genfrac{}{}{0pt}{}{k}{n-r}\right)}^2\left(\genfrac{}{}{0pt}{}{k+\alpha +r}{r}\right),$$

(19)

$${(\Delta x)}_{n,\alpha }=\sqrt{2n^2+2(\alpha +1)n+\alpha +1},$$

(20)

were found for the ordinary moments and for the standard deviation of the orthonormal Laguerre polynomials, respectively. For Jacobi polynomials, one has the following expressions:

$$\begin{array}{cc}\hfill {\langle x^k\rangle }_{n,\alpha ,\beta }& =\sum _{i=n-k}^n\left(\sum _{t=i}^n{(-1)}^{t-i}\left(\genfrac{}{}{0pt}{}{n}{t}\right)\frac{\Gamma (\alpha +n+1)\Gamma (\alpha +\beta +n+t+1)}{n!2^t\Gamma (\alpha +\beta +n+1)\Gamma (\alpha +t+1)}\left(\genfrac{}{}{0pt}{}{t}{i}\right)\right)\hfill \\ \hfill & \times \frac{{(-1)}^{k+i-n}{2}^n(k+i)!\Gamma (n+\alpha +\beta +1)}{(k+i-n)!\Gamma (2n+\alpha +\beta +1)}\hfill \\ \hfill & \times {}_2{F}_1\left(\left.\begin{array}{c}n-k-i,n+\beta +1\\ 2n+\alpha +\beta +2\end{array}\right|2\right).\hfill \end{array}$$

(21)

for the ordinary moments of the orthonormal Jacobi polynomials, and

$$\begin{array}{c}{\left(\Delta x\right)}_{n,\alpha ,\beta }=\left[\frac{4(n+1)(n+\alpha +1)(n+\beta +1)(n+\alpha +\beta +1)}{(2n+\alpha +\beta +1){(2n+\alpha +\beta +2)}^2(2n+\alpha +\beta +3)}\right.\hfill \\ \hfill +{\left.\frac{4n(n+\alpha )(n+\beta )(n+\alpha +\beta )}{(2n+\alpha +\beta -1){(2n+\alpha +\beta )}^2(2n+\alpha +\beta +1)}\right]}^{1/2},\end{array}$$

(22)

for the standard deviation. Note, in particular, that the asymptotical behavior of the standard deviation of the HOPs goes as $\sqrt{n}$ (Hermite), $\sqrt{2}n$ (Laguerre) and $1/\sqrt{2}$ (Jacobi) when the weight’s parameters are fixed. Finally, it is worth mentioning here that other kinds of moments of orthogonal polynomials were recently considered, such as the exponential and logarithmic moments as well as the generalized Krein-like moments of HOPs [59].

4. Fisher’s Spreading Length of HOPs

In this section, we give the values of the Fisher spreading length (8) for the three canonical families of the HOPs defined by the expressions (10) and (11). They are given by the corresponding quantities of the associated Rakhmanov probability density ${\rho }_n\left(x\right)$ given by (1):

$${\left(\delta x\right)}_n\equiv \frac{1}{\sqrt{F\left[{\rho }_n\right]}}\equiv {〈{\left[\frac{d}{dx}ln{\rho }_n\left(x\right)\right]}^2〉}^{-\frac{1}{2}}={\left\{{\int }_{\Delta }dx\frac{{\left[{\rho }_n^{\prime }\left(x\right)\right]}^2}{{\rho }_n\left(x\right)}\right\}}^{-\frac{1}{2}},$$

(23)

where $F\left[{\rho }_n\right]\equiv F\left[p_n\right]$ denotes the Fisher information of the polynomial $p_n\left(x\right)$. This quantity was derived for the first time from the second-order differential equation of HOPS [60] (see also [25,61]), obtaining the following values:

$${\left(\delta x\right)}_n=\frac{1}{\sqrt{4n+2}},$$

(24)

for Hermite polynomials $H_n\left(x\right)$,

$$\begin{array}{c}\hfill {\left(\delta x\right)}_{n,\alpha }=\left\{\begin{array}{cc}\frac{1}{\sqrt{4n+1}};& \alpha =0,\\ \sqrt{\frac{{\alpha }^2-1}{(2n+1)\alpha +1}};& \alpha >1,\\ 0;& \alpha \in (-1,+1],\alpha \ne 0.\end{array}\right.\end{array}$$

(25)

for Laguerre polynomials $L_n^{\left(\alpha \right)}\left(x\right)$, $\alpha >-1$, and

$${\left(\delta x\right)}_{n,\alpha ,\beta }=\frac{1}{\sqrt{F\left[{\rho }_{n,\alpha ,\beta }\right]}},$$

(26)

where the Fisher information $F\left[{\rho }_{n,\alpha ,\beta }\right]$ is given by the following:

$$\begin{array}{c}\hfill F\left[{\rho }_{n,\alpha ,\beta }\right]=\left\{\begin{array}{cc}2n(n+1)(2n+1),\hfill & \alpha ,\beta =0,\hfill \\ \frac{2n+\beta +1}{4}\left[\frac{n^2}{\beta +1}+n+(4n+1)(n+\beta +1)+\frac{{(n+1)}^2}{\beta -1}\right],\hfill & \alpha =0,\beta >1,\hfill \\ \frac{2n+\alpha +\beta +1}{4(n+\alpha +\beta -1)}\left[n(n+\alpha +\beta -1)\left(\frac{n+\alpha }{\beta +1}+2+\frac{n+\beta }{\alpha +1}\right)\right.\hfill & \\ +\left.(n+1)(n+\alpha +\beta )\left(\frac{n+\alpha }{\beta -1}+2+\frac{n+\beta }{\alpha -1}\right)\right],\hfill & \alpha ,\beta >1,\hfill \\ \infty ,\hfill & \mathrm{otherwise},\hfill \end{array}\right.\end{array}$$

for Jacobi polynomials $P_n^{(\alpha ,\beta )}\left(x\right)$, with $\alpha ,\beta >-1$. We can easily note from these values and the ones for the standard deviation of the previous section that the Cramér-Rao inequality (9) is duly fulfilled for the three families of HOPs. Moreover, we observe that the asymptotical behavior of the Fisher spreading length $\delta \left(x\right)$ for the HOPs is $n^{-1/2}$ (Hermite), $n^{-1/2}$(Laguerre) and $n^{-3/2}$ (Jacobi) when the weight’s parameters are fixed.

5. Rényi’s Spreading Lengths and Weighted ${\mathcal{L}}_{\mathbf{q}}$-Norms of HOPs

In this section, we give the values of the Rényi spreading length (6) for the three canonical families of HOPs $\left\{p_n\left(x\right)\right\}$, defined by the expressions (10) and (11). These quantities, denoted as ${\mathcal{L}}_q^R\left[p_n\right]$ for convenience, are given by the corresponding quantities of the associated Rakhmanov density ${\rho }_n\left(x\right)$ given by (1):

$${\mathcal{L}}_q^R\left[{\widehat{p}}_n\right]=e^{R_q\left[{\widehat{p}}_n\right]}={\left(W_q\left[{\widehat{p}}_n\right]\right)}^{\frac{1}{1-q}}={\left\{{\int }_{\Lambda }{\left[{\rho }_n\left(x\right)\right]}^{q}dx\right\}}^{\frac{1}{1-q}}={\left\{{\int }_{\Lambda }{\left[{\widehat{p}}_n^2\left(x\right)h\left(x\right)\right]}^{q}dx\right\}}^{\frac{1}{1-q}},q>0,q\ne 1,$$

(27)

where the symbols $R_q\left[{\widehat{p}}_n\right]$ and $W_q\left[{\widehat{p}}_n\right]$ for the Rényi entropy and the weighted ${\mathcal{L}}_q$-norm of the orthonormal polynomials ${\widehat{p}}_n\left(x\right)$ denote the qth-order Rényi entropy $R_q\left[{\rho }_n\right]$ and the entropic moment $W_q\left[{\rho }_n\right]$ of ${\rho }_n\left(x\right)$, respectively, so that the following holds:

$$R_q\left[{\widehat{p}}_n\right]=\frac{1}{1-q}ln{\int }_{\Lambda }{\left[{\widehat{p}}_n^2\left(x\right)h\left(x\right)\right]}^{q}dx,q>0,q\ne 1,$$

(28)

and

$$W_q\left[{\widehat{p}}_n\right]={\int }_{\Lambda }{\left[{\widehat{p}}_n^2\left(x\right)h\left(x\right)\right]}^{q}dx,q>0,q\ne 1,$$

(29)

Then, we have in particular that the weighted ${\mathcal{L}}_q$-norm and the Rényi spreading length of the orthogonal (HOPs) and orthonormal polynomials are mutually related as follows:

$$W_q\left[p_n\right]={\int }_{\Lambda }{\left[p_n^2\left(x\right)h\left(x\right)\right]}^{q}dx={\kappa }_n^q{W}_q\left[{\widehat{p}}_n\right]$$

(30)

and

$${\mathcal{L}}_q^R\left[p_n\right]={\left(W_q\left[p_n\right]\right)}^{\frac{1}{1-q}}={\kappa }_n^{\frac{q}{1-q}}{\left(W_q\left[{\widehat{p}}_n\right]\right)}^{\frac{1}{1-q}}.$$

(31)

The analytical determination of these norms has been a long standing problem in the theory of special functions and extremal polynomials itself since the times of Bernstein and Steklov (see [62,63,64]) and in the theory of trigonometric series [65]. More recently, these quantities were obtained by using either (a) the series expansion of the powers of the HOPs $p_n\left(x\right)$ by means of the combinatorial Bell polynomials [66], or (b) the Srivastava–Niukkanen’s linearization method [67,68] of the positive integer powers ${\left[p_n\left(x\right)\right]}^{2q}$ of HOPs by means of the multivariate Lauricella function $F_A^{\left(2q\right)}\left(\frac{1}{q},\dots ,\frac{1}{q}\right)$ [69] (Hermite and Laguerre cases) and in terms of the Srivastava–Daoust generalized function $F_{1:1;1\dots ;1}^{1:2;\dots ;2}(1,\dots ,1)$ (Jacobi case), as Sánchez-Moreno et al. showed in [22,70], respectively. The combinatorial approach is based on the expansion p-th power of the arbitrary polynomial as follows:

$$y_n\left(x\right)=\sum _{t=0}^n{c}_t{x}^t$$

(32)

given by the following [22]:

$${\left[y_n\left(x\right)\right]}^p=\sum _{t=0}^{np}\frac{p!}{(t+p)!}{B}_{t+p,p}(c_0,2!c_1,\dots ,(t+1)!c_t)x^t,$$

(33)

with $c_i=0$ for $i>n$; the remaining expansion coefficients are given by Equations (44), (49) and (56) in the Hermite, Laguerre and Jacobi cases, respectively. The B-symbols denote the renowned multivariate Bell polynomials of combinatorics [66] which are given by the following:

$$B_{m,l}(c_1,c_2,\dots ,c_{m-l+1})=\sum _{\widehat{\pi }(m,l)}\frac{m!}{j_1!j_2!\cdots j_{m-l+1}!}{\left(\frac{c_1}{1!}\right)}^{j_1}{\left(\frac{c_2}{2!}\right)}^{j_2}\cdots {\left(\frac{c_{m-l+1}}{(m-l+1)!}\right)}^{j_{m-l+1}}$$

(34)

where the sum runs over all partitions $\widehat{\pi }(m,l)$ such that the following holds:

$$j_1+j_2+\cdots +j_{m-l+1}=l,\mathrm{and}{j}_1+2j_2+\cdots +(m-l+1)j_{m-l+1}=m.$$

(35)

On the other hand, the algebraic approach [70] uses the Srivastava–Niukanen linearization formulas of the HOPS. In particular, for the Laguerre polynomials we have the following [70]:

$$x^{\mu }{\left[L_n^{\left(\alpha \right)}\left(tx\right)\right]}^r=\sum _{i=0}^{\infty }{c}_i(\mu ,r,t,n,\alpha ,\gamma )L_i^{\left(\gamma \right)}\left(x\right)$$

(36)

with the following linearization coefficients:

$$c_i(\mu ,r,t,n,\alpha ,\gamma )={(\gamma +1)}_{\mu }{\left(\genfrac{}{}{0pt}{}{n+\alpha }{n}\right)}^r{F}_A^{(r+1)}\left(\begin{array}{c}\gamma +\mu +1;\stackrel{r}{\overbrace{-n,\dots ,-n}},-i\\ \underset{r}{\underbrace{\alpha +1,\dots ,\alpha +1}},\gamma +1\end{array};\stackrel{r}{\overbrace{t,\dots ,t}},1\right),$$

(37)

where the Pochhammer symbol ${\left(z\right)}_a=\frac{\Gamma (z+a)}{\Gamma \left(z\right)}$ [6] and $F_A^{(r+1)}(x_1,\dots ,x_r)$ denote the Lauricella function of type A of $r+1$ variables and $2r+3$ parameters defined as follows: [69]

$$F_A^{\left(s\right)}\left(\begin{array}{c}a;b_1,\dots ,b_s\\ c_1,\dots ,c_s\end{array};x_1,\dots ,x_s\right)=\sum _{j_1,\dots ,j_s=0}^{\infty }\frac{{\left(a\right)}_{j_1+\cdots +j_s}{\left(b_1\right)}_{j_1}\cdots {\left(b_s\right)}_{j_s}}{{\left(c_1\right)}_{j_1}\cdots {\left(c_s\right)}_{j_s}}\frac{x_1^{j_1}\cdots x_s^{j_s}}{j_1!\cdots j_s!}.$$

(38)

Similarly, we have the following linearization formula for the Jacobi polynomials [70]:

$${\left(P_n^{(\alpha ,\beta )}\left(x\right)\right)}^r=\sum _{i=0}^{\infty }{\tilde{c}}_i(r,n,\alpha ,\beta ,\gamma ,\delta )P_i^{(\gamma ,\delta )}\left(x\right),$$

(39)

where the coefficients are given by the following:

$$\begin{array}{c}{\tilde{c}}_i(r,n,\alpha ,\beta ,\gamma ,\delta )={\left(\genfrac{}{}{0pt}{}{n+\alpha }{n}\right)}^r\frac{\gamma +\delta +2i+1}{\gamma +\delta +i+1}\sum _{j_1,\dots ,j_r=0}^n\sum _{j_{r+1}=0}^i\frac{{(\gamma +1)}_{j_1+\cdots +j_r+j_{r+1}}}{{(\gamma +\delta +i+2)}_{j_1+\cdots +j_r}}\hfill \\ \hfill \times \frac{{(-n)}_{j_1}{(\alpha +\beta +n+1)}_{j_1}\cdots {(-n)}_{j_r}{(\alpha +\beta +n+1)}_{j_r}{(-i)}_{j_{r+1}}}{{(\alpha +1)}_{j_1}\cdots {(\alpha +1)}_{j_r}{(\gamma +1)}_{j_{r+1}}{j}_1!\cdots j_r!j_{r+1}!},\end{array}$$

(40)

To obtain the linearization formula for the Hermite polynomials within the framework of this algebraic approach, we use Equation (36) together with the known relation of the Hermite and Laguerre polynomials [6]:

$$H_{2n}\left(x\right)={(-1)}^n{2}^{2n}n!L_n\left(x^2\right);H_{2n+1}\left(x\right)={(-1)}^n{2}^{2n+1}n!x{L}_n\left(x^2\right).$$

(41)

5.1. Hermite Polynomials

From (27) and the explicit expressions (32) of the Hermite polynomials $H_n\left(x\right)$, the combinatorial approach allows one to obtain the following values of the Rényi spreading lengths ${\mathcal{L}}_q^R\left[H_n\right]$ for $2q\in \mathbb{N},q>2$:

$${\mathcal{L}}_q^R\left[H_n\right]={\left\{W_q\left[H_n\right]\right\}}^{\frac{1}{1-q}},$$

(42)

with the weighted norm $W_q\left[H_n\right]$ of the Hermite polynomials given by the following:

$$W_q\left[H_n\right]={\int }_{-\infty }^{\infty }{\left[H_n^2\left(x\right)h_H\left(x\right)\right]}^{q}dx=\sum _{j=0}^{nq}\frac{\Gamma \left(j+\frac{1}{2}\right)}{q^{j+\frac{1}{2}}}\frac{\left(2q\right)!}{(2j+2q)!}{B}_{2j+2q,2q}(c_0^{\left(n\right)},2!c_1^{\left(n\right)},\dots ,(2j+1)!c_{2j}^{\left(n\right)}),$$

(43)

with the following values:

$$c_t^{\left(n\right)}=\frac{{(-1)}^{\frac{3n-t}{2}}n!}{{\left(2^{n}n!\sqrt{\pi }\right)}^{\frac{1}{2}}}\frac{2^t}{\left(\frac{n-t}{2}\right)!t!}\frac{{(-1)}^t+{(-1)}^n}{2}.$$

(44)

for the expansion coefficients $c_t$ of the orthonormal Hermite polynomials (see [3,6,22]). Then, the expressions (43) and (42) together with (34) and (44) provide an algorithmic procedure to find the qth-weighted norm and the qth-Rényi spreading length of the Hermite polynomial, respectively, in terms of its degree n and the parameter q. See Section $3.2$ of [22] for further details.

Alternatively, from Equations (36) and (41) and the orthogonality relation (10) of the Hermite polynomials, the algebraic approach [70] allows one to find the following values:

$$\begin{array}{c}\hfill W_q\left[H_n\right]=q^{-\nu q-\frac{1}{2}}{\left(\frac{\Gamma \left(\frac{n+\nu +1}{2}\right)}{\Gamma \left(\frac{n-\nu }{2}+1\right)}\right)}^q{\left(\Gamma \left(\nu +\frac{1}{2}\right)\right)}^{-2q}\Gamma \left(q\nu +\frac{1}{2}\right)F_A^{\left(2q\right)}\left(\begin{array}{c}q\nu +\frac{1}{2};\stackrel{2q}{\overbrace{\frac{\nu -n}{2},\dots ,\frac{\nu -n}{2}}}\\ \underset{2q}{\underbrace{\nu +\frac{1}{2},\dots ,\nu +\frac{1}{2}}}\end{array};\stackrel{2q}{\overbrace{\frac{1}{q},\dots ,\frac{1}{q}}}\right)\end{array}$$

(45)

for the weighted norm of Hermite polynomials. Note that this quantity is completely determined by n and q. Let us recall here that $F_A^{(r+1)}(x_1,\dots ,x_r)$ denotes the Lauricella function of type A of $r+1$ variables and $2r+3$ parameters defined as follows [69]:

$$F_A^{\left(s\right)}\left(\begin{array}{c}a;b_1,\dots ,b_s\\ c_1,\dots ,c_s\end{array};x_1,\dots ,x_s\right)=\sum _{j_1,\dots ,j_s=0}^{\infty }\frac{{\left(a\right)}_{j_1+\cdots +j_s}{\left(b_1\right)}_{j_1}\cdots {\left(b_s\right)}_{j_s}}{{\left(c_1\right)}_{j_1}\cdots {\left(c_s\right)}_{j_s}}\frac{x_1^{j_1}\cdots x_s^{j_s}}{j_1!\cdots j_s!}.$$

(46)

which depends on $N={\sum }_{i=0}^r(p_i+q_i)$ parameters.

Then, from expressions Equations (42), (43) and (45), one has the Rényi spreading length of the Hermite polynomials in the two combinatorial and algebraic approaches. In particular, for $q=2$, we obtain the following values for the Onicescu information-theoretic length ${\mathcal{L}}_2^R\left[H_n\right]={\mathcal{L}}_2^R\left[{\rho }_n\right]$

$${\mathcal{L}}_2^R\left[H_0\right]=\sqrt{2\pi },{\mathcal{L}}_2^R\left[H_1\right]=\frac{4}{3}\sqrt{2\pi },{\mathcal{L}}_2^R\left[H_2\right]=\frac{64}{41}\sqrt{2\pi }$$

of the first few Hermite polynomials with degrees $n=0,1,2$. See [22,70,71] for further details and application to the determination of the total Rényi entropies of the one- and multidimensional harmonic systems in position and momentum spaces.

5.2. Laguerre Polynomials

From (27) and the explicit expressions (32) of the orthonormal Laguerre polynomials ${\widehat{L}}_n^{\left(\alpha \right)}\left(x\right)$, the combinatorial approach allows one to obtain the following values of the Rényi spreading lengths ${\mathcal{L}}_q^R[{\widehat{L}}_n^{\left(\alpha \right)}]={\mathcal{L}}_q^R\left[{\rho }_{n,\alpha }\right]$ for $2q\in \mathbb{N},q>2$ and being ${\rho }_{n,\alpha }={\widehat{L}}_n^{\left(\alpha \right)}{\left(x\right)}^2{h}_{\alpha }^L\left(x\right)$:

$${\mathcal{L}}_q^R[{\widehat{L}}_n^{\left(\alpha \right)}]={\left\{W_q[{\widehat{L}}_n^{\left(\alpha \right)}]\right\}}^{\frac{1}{1-q}},$$

(47)

with the weighted norm $W_q\left[L_n^{\left(\alpha \right)}\right]$ of the Laguerre polynomials given by the following:

$$W_q\left[{\widehat{L}}_n^{\left(\alpha \right)}\right]={\int }_0^{\infty }{\left[{\widehat{L}}_n^{\left(\alpha \right)}{\left(x\right)}^2{h}_{\alpha }^L\left(x\right)\right]}^{q}dx=\sum _{k=0}^{2nq}\frac{\Gamma (\alpha q+k+1)}{q^{\alpha q+k+1}}\frac{\left(2q\right)!}{(k+2q)!}{B}_{k+2q,2q}\left(c_0^{(n,\alpha )},2!c_1^{(n,\alpha )},...,(k+1)!c_k^{(n,\alpha )}\right),$$

(48)

where B-symbol describes the multivariate Bell polynomials to be calculated at the following values:

$$c_t^{(n,\alpha )}=\sqrt{\frac{\Gamma (n+\alpha +1)}{n!}}\frac{{(-1)}^t}{\Gamma (\alpha +t+1)}\left(\genfrac{}{}{0pt}{}{n}{t}\right).$$

(49)

for the expansion coefficients of the orthonormal Laguerre polynomials (see [6,24]). Then, the expressions (48) and (47) together with (34) and (49) provide an algorithmic procedure to find the qth-weighted norm and the qth-Rényi spreading length of the Laguerre polynomial, respectively, in terms of its degree n and the parameters $\alpha$ and q. In particular, for $q=2$, we obtain the following values for the Onicescu information-theoretic length ${\mathcal{L}}_2^R\left[{\widehat{L}}_n^{\left(\alpha \right)}\right]$:

$${\mathcal{L}}_2^R\left[{\widehat{L}}_0^{\left(\alpha \right)}\right]={\left(\frac{2^{2\alpha +1}{\left(\Gamma (\alpha +1)\right)}^2}{\Gamma (2\alpha +1)}\right)}^{\frac{1}{2}},{\mathcal{L}}_2^R\left[{\widehat{L}}_1^{\left(\alpha \right)}\right]={\left(\frac{2^{2\alpha +3}{\left(\Gamma (\alpha +2)\right)}^2}{(1+\alpha )(2+3\alpha )\Gamma (2\alpha +1)}\right)}^{\frac{1}{2}}$$

(50)

of the Laguerre polynomials $L_n^{\left(\alpha \right)}\left(x\right)$ with $n=0,1$.

Alternatively, from the expression (36) and the orthogonality relation (10) of the Laguerre polynomials, the algebraic approach [70] allows one to find the following values:

$$\begin{array}{c}\hfill W_q^R\left[L_n^{\left(\alpha \right)}\right]=\frac{\Gamma (\alpha q+1)}{q^{\alpha q+1}}\frac{{\left(\Gamma (n+\alpha +1)\right)}^q}{{(n!)}^q{\left(\Gamma (\alpha +1)\right)}^{2q}}{F}_A^{\left(2q\right)}\left(\begin{array}{c}\alpha q+1;\stackrel{2q}{\overbrace{-n,\dots ,-n}}\\ \underset{2q}{\underbrace{\alpha +1,\dots ,\alpha +1}}\end{array};\stackrel{2q}{\overbrace{\frac{1}{q},\dots ,\frac{1}{q}}}\right),\end{array}$$

(51)

for the weighted norm of the Laguerre polynomials $L_n^{\left(\alpha \right)}\left(x\right)$ for all integer orders (i.e., when $q\in \mathbb{N}$). Note that these quantities only depend on the order q, the polynomial degree n and the polynomial parameter $\alpha$, where $F_A^{\left(2q\right)}(x_1,\dots ,x_r)$ denotes the Lauricella function of type A given by Equation (46).

Then, from expressions Equations (47) and (51) one has the Rényi spreading length of the Laguerre polynomials in the two combinatorial and algebraic approaches. In particular, we have the following values:

$${\mathcal{L}}_q^R\left[{\widehat{L}}_0^{\left(\alpha \right)}\right]={\left[\frac{1}{\Gamma {(\alpha +1)}^q}\frac{\Gamma (\alpha q+1)}{q^{\alpha q+1}}\right]}^{\frac{1}{1-q}},$$

(52)

which recalls the Mehta–Selberg integral, and

$${\mathcal{L}}_q^R\left[{\widehat{L}}_1^{\left(\alpha \right)}\right]={\left[\frac{\Gamma (\alpha q+1){(1+\alpha )}^{2q}}{{(\Gamma (\alpha +2))}^q{q}^{\alpha q+1}}{}_2{F}_0\left(\begin{array}{c}-2q,\alpha q+1\\ -\end{array};\frac{1}{q(\alpha +1)}\right)\right]}^{\frac{1}{1-q}}.$$

(53)

for the Laguerre polynomials $L_n^{\left(\alpha \right)}\left(x\right)$ with $n=0,1$. Note also that Equations (52) and (53) with $q=2$ boil down to the previous values (50), respectively, of the Onicescu lengths. See [70,72] for further details and application to the determination of the radial Rényi entropies of the multidimensional hydrogenic systems in position and momentum spaces.

5.3. Jacobi Polynomials

From (27) and the explicit expressions (32) of the Jacobi polynomials $P_n^{(\alpha ,\beta )}\left(x\right)$, the combinatorial approach allows one to obtain the following values of the Rényi spreading lengths ${\mathcal{L}}_q^R[P_n^{(\alpha ,\beta )}]={\mathcal{L}}_q^R\left[{\rho }_{n,\alpha ,\beta )}\right]$ for $2q\in \mathbb{N},q>2$:

$${\mathcal{L}}_q^R[P_n^{(\alpha ,\beta )}]={\left\{W_q\left[P_n^{(\alpha ,\beta )}\right]\right\}}^{\frac{1}{1-q}},$$

(54)

with the weighted norm $W_q[P_n^{(\alpha ,\beta )}]$ of the Jacobi polynomials given by the following:

$$W_q[P_n^{(\alpha ,\beta )}]={\int }_{-1}^{+1}{\left[P_n^{(\alpha ,\beta )}{\left(x\right)}^2{h}_{\alpha ,\beta }^J\left(x\right)\right]}^{q}dx=\sum _{k=0}^{2nq}\frac{\left(2q\right)!}{(k+2q)!}{B}_{k+2q,2q}\left(c_0^{(n,\alpha ,\beta )},2!c_1^{(n,\alpha ,\beta )},\dots ,(k+1)!c_k^{(n,\alpha ,\beta )}\right)\mathcal{I}(k,q,\alpha ,\beta ),$$

(55)

with the parametric function

$$\mathcal{I}(k,q,\alpha ,\beta )=\frac{{(-1)}^k{2}^{1+\alpha q+\beta q}\Gamma (\alpha q+1)\Gamma (\beta q+1)}{\Gamma (\alpha q+\beta q+2)}{}_2{F}_1\left(\begin{array}{c}-k,1+\beta q\hfill \\ 2+(\alpha +\beta )q\hfill \end{array};2\right),$$

with the Jacobi expansion coefficients (see [3,6,23])

$$\begin{array}{c}\hfill c_t^{(n,\alpha ,\beta )}=\sqrt{\frac{\Gamma (\alpha +n+1)(2n+\alpha +\beta +1)}{n!2(\alpha +\beta +1)\Gamma (\alpha +\beta +n+1)\Gamma (n+\beta +1)}}\sum _{i=t}^n{(-1)}^{i-t}\left(\genfrac{}{}{0pt}{}{n}{i}\right)\left(\genfrac{}{}{0pt}{}{i}{t}\right)\frac{\Gamma (\alpha +\beta +n+i+1)}{2^i\Gamma (\alpha +i+1)}\end{array}$$

(56)

Then, the expressions (55) and (54) together with (34) and (56) provide an algorithmic procedure to find the qth-weighted norm and the qth-Rényi spreading length of the Jacobi polynomial, respectively, in terms of its degree n and the parameters $(\alpha ,\beta )$ and q. See Section 4 of [23] for further details.

Alternatively, from the expression (39) and the orthogonality relation (10) of the Jacobi polynomials, the algebraic approach [70] allows one to find the following values:

$$\begin{array}{c}{W}_q\left[P_n^{(\alpha ,\beta )}\right]=\left[\frac{2^{\alpha q+\beta q+1}\Gamma (\alpha q+1)\Gamma (\beta q+1)}{\Gamma (\alpha q+\beta q+2)}{\left(\frac{n!(\alpha +\beta +2n+1)\Gamma (\alpha +\beta +n+1)}{2^{\alpha +\beta +1}\Gamma (\alpha +n+1)\Gamma (\beta +n+1)}\right)}^q\right.\hfill \\ \hfill \left.\times {\left(\genfrac{}{}{0pt}{}{n+\alpha }{n}\right)}^{2q}{F}_{1:1;\dots ;1}^{1:2;\dots ;2}\left(\begin{array}{c}\alpha q+1:\stackrel{2q}{\overbrace{-n,\alpha +\beta +n+1;\dots ;-n,\alpha +\beta +n+1}}\\ \alpha q+\beta q+2:\underset{2q}{\underbrace{\alpha +1;\dots ;\alpha +1}}\end{array};\stackrel{2q}{\overbrace{1,\dots ,1}}\right)\right]\end{array}$$

(57)

for the weighted norms of the Jacobi polynomials $P_n^{(\alpha ,\beta )}\left(x\right)$ with order $q=2,3,4,\dots$. Note that these quantities depend only on the order q, the polynomial degree n and parameters $\alpha$ and $\beta$. Let us recall here that the symbol $F_{1:1;\dots ;1}^{1:2;\dots ;2}(x_1,\dots ,x_r)$ denotes the r-variate Srivastava–Daoust function defined as follows [67,70]:

$$\begin{array}{ccc}\hfill F_{1:1;\dots ;1}^{1:2;\dots ;2}\left(\begin{array}{cc}{a}_0^{\left(1\right)}:a_1^{\left(1\right)},a_1^{\left(2\right)};\dots ;a_r^{\left(1\right)},a_r^{\left(2\right)}& \\ & ;x_1,\dots ,x_r\\ b_0^{\left(1\right)}:b_1^{\left(1\right)};\dots ;b_r^{\left(1\right)}& \end{array}\right)& =& \hfill \\ \hfill =\sum _{j_1,\dots ,j_r=0}^{\infty }\frac{{\left(a_0^{\left(1\right)}\right)}_{j_1+\dots +j_r}}{{\left(b_0^{\left(1\right)}\right)}_{j_1+\dots +j_r}}\frac{{\left(a_1^{\left(1\right)}\right)}_{j_1}{\left(a_1^{\left(2\right)}\right)}_{j_1}\cdots {\left(a_r^{\left(1\right)}\right)}_{j_r}{\left(a_r^{\left(2\right)}\right)}_{j_r}}{{\left(b_1^{\left(1\right)}\right)}_{j_1}{\left(b_r^{\left(1\right)}\right)}_{j_r}}\frac{x_1^{j_1}{x}_2^{j_2}\cdots x_r^{j_r}}{j_1!j_2!\cdots j_r!},& \end{array}$$

(58)

which depends on $N={\sum }_{i=0}^r(p_i+q_i)$ parameters.

Then, from the expressions of Equations (54), (55) and (57), one has the Rényi spreading length of the Laguerre polynomials in the two combinatorial and algebraic approaches [23]. See [67,70,72] for further details and application to the determination of the angular Rényi entropies of the multidimensional hydrogenic systems in position and momentum spaces.

Finally, since the expressions obtained for the weighted norms $W\left[p_n\right]$ and, consequently, for the Rényi entropies $R\left[p_n\right]$ and the Rényi spreading lengths ${\mathcal{L}}_1^R\left[p_n\right]$ of the HOPs are not easily manipulated analytically, it would be interesting at least to have the asymptotical values ($n\to \infty$) and to obtain simple, compact and accurate upper bounds to these quantities in terms of the degree and the parameters of their weight function. The latter problem has not yet been determined, but the former one related to the asymptotics ($n\to \infty$) is extensively shown in the next three sections for the Hermite, Laguerre and Jacobi polynomials and its generalizations.

6. Degree Asymptotics for the Weighted ${\mathfrak{L}}_{\mathbf{q}}$-Norms of Hermite Polynomials

The aim of this section is the strong asymptotic ($n\to \infty$) determination of the entropic moments or weighted ${\mathfrak{L}}_q$-norms of Hermite polynomials, i.e., the following:

$$W_q\left[H_n\right]={\int }_{-\infty }^{+\infty }{\left[H_n^2\left(x\right)e^{-x^2}\right]}^{q}dx;q>0.$$

(59)

This problem was initiated by Aptekarev et al. [34,36] who solved it for $q\in \left[0,\frac{4}{3}\right]$ by use of an extension of the Plancherel and Rotach asymptotics of the orthonormal Hermite polynomials [1]. They obtained [34] the following expression:

$${\int }_{-\infty }^{+\infty }{\left[{\widehat{H}}_{n-1}^2\left(x\right)e^{-x^2}\right]}^{q}dx=c_q{\left(2n\right)}^{\frac{1-q}{2}}\left(1+o\left(1\right)\right),forq\le \frac{4}{3},$$

(60)

for the main term of the weighted norm of the orthonormal polynomials ${\tilde{H}}_{n-1}\left(x\right)$, with the following constant:

$${\displaystyle c_q={\left({\displaystyle \frac{2}{\pi }}\right)}^q\frac{\Gamma (q+\frac{1}{2})}{\Gamma (q+1)}\frac{\Gamma (1-\frac{q}{2})}{\Gamma (\frac{3}{2}-\frac{q}{2})}.}$$

(61)

Now, taking into account the normalization constant (12) and the Stirling formula for the gamma function [6], this expression gives rise to the following asymptotics for the weighted norm of the orthogonal Hermite polynomial $H_{n-1}\left(x\right)$:

$$W_q\left[H_{n-1}\right]=c_q{\kappa }_{n-1}^q{\left(2n\right)}^{\frac{1-q}{2}}\left(1+o\left(1\right)\right)=c_q{\pi }^q{\left(2n\right)}^{q(n-1)+1/2}{e}^{-qn}\left(1+o\left(1\right)\right).$$

(62)

To improve this expression and extends its validity for all values of q, we have to go beyond the Plancherel–Rotach asymptotics since it cannot give the asymptotics of the Hermite polynomials on the whole real line. This can be achieved either by means of the powerful Riemann–Hilbert method of Deift et al. [73] (see also [74]) or by the Tulyakov approach [75], whose starting point is the weight of orthogonality and the recurrence relation, which characterizes these orthogonal polynomials, respectively. The application of the latter approach to Hermite polynomials has allowed us to find [76] the following asymptotical result:

Theorem 1.

Let $H_n\left(x\right)$ be the Hermite polynomials with the standard orthogonality (10)–(12). Then, the entropic moments or weighted ${\mathfrak{L}}_q$-norms $W_q\left[{\rho }_{n-1}\right]$, given by Equation (59), have the following asymptotic ($n\to \infty$) behavior:

$$W_q\left[H_{n-1}\right]=\left\{\begin{array}{cc}{c}_q{\pi }^q{\left(2n\right)}^{q(n-1)+1/2}{e}^{-qn}\left(1+o\left(1\right)\right),\hfill & q<2,\hfill \\ 2{\left(2n\right)}^{2n-\frac{3}{2}}{e}^{-2n}\left(ln\left(n\right)+O\left(1\right)\right),\hfill & q=2,\hfill \\ 2C_q{2}^{-q}{\left(2n\right)}^{q\left(n-\frac{2}{3}\right)-\frac{1}{6}}{e}^{-qn}\left(1+o\left(1\right)\right),\hfill & q>2.\hfill \end{array}\right.$$

(63)

where the constant $c_q$ is defined in (61) and the constant $C_q$ is equal to the following:

$$C_q={\int }_{-\infty }^{+\infty }{\left[\frac{2\pi }{\sqrt[3]{2}}{\mathrm{Ai}}^2\left(-\frac{z\sqrt[3]{2}}{2}\right)\right]}^{q}dz.$$

with$z:=\frac{2n}{x^{\frac{2}{3}}}-x^{\frac{4}{3}},$and the Airy function$\mathrm{Ai}$(see [6]).

We note that the first asymptotic formula in the right hand side of (63) coincides with (62), but now it holds true in the maximal range of p (when $p=2$, then $c_p=\infty$); let us also highlight that the main term of the asymptotics is growing. Moreover, the smaller terms contain a constant which depends on p, and when $p\to 0$, this constant tends to infinity; however, our formula is correct for any small fixed $p>0$. We also note that the leading term of all three formulae in the right hand side of (63) match each other when $p\to 2$.

Let us highlight the great simplicity of the main term of asymptotics ($n\to \infty$) of the weighted norm of the Hermite polynomials, $W_q\left[H_n\right]$, given by expressions (63). Indeed, this result is impressively simple, given the complexity of the expressions (43) and (45) obtained in Section 5 for the corresponding quantities of a Hermite polynomial with fixed degree n, which require the evaluation of either the multivariate Bell polynomials (34) at the Hermite expansion coefficients (44) or the multivariate Lauricella function of type A of the type (46) at $1/q$, respectively.

This theorem facilitates the evaluation of the Rényi entropies of the highly-excited (i.e., Rydberg) states of the one-dimensional oscillator-like systems [76], basically because the radial wave functions of these quantum systems are controlled by Hermite’s polynomials. Briefly, the various physical quantities described by the Rényi entropies with different orders q have a clearly increasing dependence on the principal quantum number n characterizing the Rydberg state of the system.

7. Degree Asymptotics for the Generalized Weighted ${\mathfrak{L}}_{\mathbf{q}}$-Norms of Laguerre Polynomials

In this section, we study and solve the asymptotics ($n\to \infty$) of the generalized weighted ${\mathfrak{L}}_q$-norms of the Laguerre polynomials defined as follows:

$$N_n^{\left(q\right)}(\alpha ,\beta )=\underset{0}{\overset{\infty }{\int }}{\left({\left[{\widehat{L}}_n^{\left(\alpha \right)}\left(x\right)\right]}^2{h}_{\alpha }^L\left(x\right)\right)}^q{x}^{\beta }dx,q>0,\alpha >-1.$$

(64)

with $\beta$ so that $\beta +p\alpha >-1$ to get convergence, and the weight function $h_{\alpha }^L\left(x\right)=x^{\alpha }{e}^{-x}$. Note that for $\beta =0$ these quantities denote the weighted ${\mathfrak{L}}_q$-norms (29) of the orthonormal Laguerre polynomials, $W_q\left[{\widehat{L}}_n^{\left(\alpha \right)}\right]$. The orthonormal polynomials ${\left[{\widehat{L}}_n^{\left(\alpha \right)}\left(x\right)\right]}^2={\left[L_n^{\left(\alpha \right)}\left(x\right)\right]}^2/{\kappa }_n^L$, with the norm ${\kappa }_n^L={\parallel L_n^{\left(\alpha \right)}\parallel }^2=d_{n,\alpha }^2$ given by (13) and the orthogonal polynomials:

$$L_n^{\left(\alpha \right)}\left(x\right)=\sum _{\nu =0}^n\left(\genfrac{}{}{0pt}{}{n+\alpha }{n-\nu }\right)\frac{{(-x)}^{\nu }}{\nu !}.$$

(65)

This problem, which essentially depends on the values of the parameters $\alpha ,\beta$ and q, was tackled by Aptekarev et al. [77] in 2016. They realized that the dominant contribution to the integral $N_n^{\left(q\right)}(\alpha ,\beta )$ is given by different regions of integration in (64) depending on the parameters’ values, which control the behavior of the polynomials along its orthogonality interval. Accordingly, we have to use different asymptotical representations for the Laguerre polynomials in the various integration’s regions of the orthogonality interval $(0,\infty )$; these representations appear to follow different scales.

First, at the extreme left of the orthogonality interval (i.e., in the neighborhood of zero) we have the Bessel regime because the Laguerre polynomials can be asymptotically represented by means of Bessel functions. Then, to the right, in the bulk region of zeros’ location, the oscillatory behavior of the polynomials is modeled asymptotically by means of the trigonometric functions. At the neighborhood of the extreme right zeros, asymptotics of the polynomials is given by Airy functions: it is the Airy regime. Finally, in the neighborhood of ∞ the polynomials have growing asymptotics. Furthermore, in some contiguous integration’s regions, the asymptotical behaviors appear to match each other. Namely, asymptotics of the Bessel functions for large arguments match the trigonometric function, and asymptotics of the Airy functions do the same.

There are in total five asymptotical regimes that can give (depending on $\alpha$ and q) the dominant contribution in the asymptotics of $N_n^{\left(q\right)}(\alpha ,\beta )$. Three of them show a growing dependence on n, according to a power law with an exponent, which depends on $\alpha$ and q; we refer to them as Bessel, Airy and cosine (or oscillatory) regimes. In addition, there are two more asymptotical regimes associated to the transition regions: cosine-Bessel and cosine-Airy. If these regimes dominate in integral (64), then the asymptotics of $N_n^{\left(q\right)}(\alpha ,\beta )$ has a factor $lnn$ besides the power law in n. Summarizing, a detailed and highbrow analysis led Aptekarev et al. [77] to the following result.

Theorem 2.

The asymptotics ($n\to \infty$) of the generalized weighted ${\mathfrak{L}}_q$-norms $N_n^{\left(q\right)}(\alpha ,\beta )$ of the Laguerre polynomials defined by (64), is given by the following expressions:

A. For $\beta >0$,

$$N_n^{\left(q\right)}(\alpha ,\beta )=\left\{\begin{array}{cc}C(\beta ,q){\left(2n\right)}^{1+\beta -q}(1+o\left(1\right)),\hfill & q\in (0,2)\hfill \\ \\ {\displaystyle \frac{lnn+O\left(1\right)}{{\pi }^2{\left(4n\right)}^{1-\beta }},}\hfill & q=2\hfill \\ \\ {\displaystyle \frac{C_A\left(q\right)}{{\pi }^q}{\left(4n\right)}^{(\frac{1-2q}{3}+\beta )}(1+o\left(1\right)),}\hfill & q\in (2,2+3\beta )\hfill \\ \\ \left({\displaystyle \frac{C_A\left(q\right)}{{\pi }^q{4}^{\beta +1}}+C_B(\alpha ,\beta ,q)}\right)n^{-\beta -1},\hfill & q=2+3\beta \hfill \\ \\ C_B(\alpha ,\beta ,q)n^{-\beta -1},\hfill & q\in (2+3\beta ,\infty )\hfill \end{array}\right.,$$

(66)

B. For $\beta =0$,

$$N_n^{\left(q\right)}(\alpha ,0)\equiv W_q\left[{\widehat{L}}_n^{\left(\alpha \right)}\right]=\left\{\begin{array}{cc}C(0,q){\left(2n\right)}^{(1-q)}(1+o\left(1\right)),\hfill & q\in (0,2)\hfill \\ \\ {\displaystyle \frac{lnn+O\left(1\right)}{{\pi }^{2}n},}\hfill & q=2\hfill \\ \\ {\displaystyle \frac{C_B(\alpha ,0,q)}{n}(1+o\left(1\right)),}\hfill & q>2\hfill \end{array}\right.,$$

(67)

C. For $\beta <0$ (and $p>0$),

$$N_n^{\left(q\right)}(\alpha ,\beta )=\left\{\begin{array}{cc}C(\beta ,q){\left(2n\right)}^{1+\beta -q}(1+o\left(1\right)),\hfill & q\in (0,2+2\beta )\hfill \\ \\ {\displaystyle \frac{2\Gamma (q+1/2)(lnn+O(1\left)\right)}{{\pi }^{q+1/2}\Gamma (q+1){\left(4n\right)}^{1+\beta }},}\hfill & q=2+2\beta \hfill \\ \\ C_B(\alpha ,\beta ,q)n^{-(1+\beta )}(1+o\left(1\right)),\hfill & q>2+2\beta \hfill \end{array}\right.,$$

(68)

respectively, where the characteristics constants $C,C_B,C_A$ are given by the following:

$$C_B(\alpha ,\beta ,q):=2\underset{0}{\overset{\infty }{\int }}{t}^{2\beta +1}{\left|J_{\alpha }\left(2t\right)\right|}^{2q}dt.$$

(69)

for the Bessel regime,

$$C_A\left(q\right):={\int }_{-\infty }^{+\infty }{\left[\frac{2\pi }{\sqrt[3]{2}}{\mathrm{Ai}}^2\left(-\frac{t\sqrt[3]{2}}{2}\right)\right]}^{q}dt.$$

(70)

for the Airy regime, and

$${\displaystyle C(\beta ,q):=\frac{2^{\beta +1}}{{\pi }^{q+1/2}}\frac{\Gamma (\beta +1-q/2)\Gamma (1-q/2)\Gamma (q+1/2)}{\Gamma (\beta +2-q)\Gamma (1+q)}.}$$

(71)

for the cosine regime, respectively. The symbols $J_{\alpha }\left(z\right)$ and $\mathrm{Ai}(-z)$ denote the known Bessel and Airy functions [6,77], respectively.

Note that in all cases, the great simplicity of the main term of asymptotics ($n\to \infty$) for the the generalized weighted ${\mathfrak{L}}_q$-norms $N_n^{\left(q\right)}(\alpha ,\beta )$. This is more impressive for the weighted ${\mathfrak{L}}_q$-norms of the Laguerre polynomials, $N_n^{\left(q\right)}(\alpha ,0)=W_q\left[{\widehat{L}}_n^{\left(\alpha \right)}\right]$, given by expressions (67). Indeed, this result is amazingly simple, given the complexity of the expressions (48) and (51) obtained in Section 5 for the corresponding quantities of a Laguerre polynomial with arbitrary, fixed degree n, which require the evaluation of either the multivariate Bell polynomials (34) at the Laguerre expansion coefficients (49), or the multivariate Lauricella function of type A of the type (46) at $1/q$, respectively.

These simple, compact asymptotical results were used to calculate the Rényi and Shannon entropies of the highly-excited (Rydberg) states of the three and multidimensional harmonic [77] and hydrogenic [40] systems. This is basically because the radial wave functions of these quantum systems are controlled by Laguerre’s polynomials.

8. Degree Asymptotics for Generalized Weighted ${\mathfrak{L}}_{\mathbf{q}}$-Norms of OPs with General and Jacobi Weights

In this section, we investigate and solve the asymptotics of the generalized weighted ${\mathfrak{L}}_q$-norms $N_n^{\left(q\right)}(h,F)$ of the polynomials ${\widehat{p}}_n\left(x\right)$ orthonormal with respect to a general weight function (defined almost everywhere on $[-1,1]$) given by the following:

$$N_n^{\left(q\right)}(h,F):=\underset{-1}{\overset{1}{\int }}{\left[{\widehat{p}}_n^2\left(x\right)h\left(x\right)\right]}^{q}F\left(x\right)dx$$

(72)

with

$$F\left(x\right):={(1-x)}^{\alpha }{(1+x)}^{\beta }f\left(x\right),\alpha \ge \beta >-1,f\left(x\right)\in {\mathfrak{L}}_1(-1,1),$$

(73)

Then, we use this general asymptotics result to obtain the generalized weighted ${\mathfrak{L}}_q$-norms of the Jacobi polynomials $P_n^{(\alpha ,\beta )}\left(x\right)$ and, as instances, the (simple) weighted ${\mathfrak{L}}_q$-norms denoted by $W_q\left[P_n^{(\alpha ,\beta )}\right]$.

Let us begin by considering the sequence ${\left\{{\widehat{p}}_n\left(x\right)\right\}}_{n=0}^{\infty }$ of the orthonormal polynomials as follows:

$$\underset{-1}{\overset{1}{\int }}{\widehat{p}}_n\left(x\right){\widehat{p}}_m\left(x\right)h\left(x\right)dx={\delta }_{n,m},\mathrm{deg}{p}_n=n$$

(74)

with respect to a general weight function

$$h\left(x\right)>0,\mathrm{a}.\mathrm{e}.x\in [-1,1].$$

(75)

To obtain the dominant term in the asymptotics of the norms $N_n^{\left(q\right)}(h,F)$, we first realize that the weight function $h\left(x\right)$ in (74) satisfies the Szegö condition as follows:

$$\underset{-1}{\overset{1}{\int }}\frac{h\left(x\right)}{\sqrt{1-x^2}}dx>-\infty ,$$

(76)

(which is more restrictive than (75)), then the orthonormal polynomials are known to have [1] the following asymptotical behavior on $[-1,1]$:

$$\underset{-1}{\overset{1}{\int }}{\left|{\widehat{p}}_n\left(x\right)-{\overline{p}}_n\left(x\right)\right|}^{2}h\left(x\right)dx\to 0,\mathrm{when}n\to \infty ,$$

(77)

where the polynomials

$${\overline{p}}_n\left(x\right):=\sqrt{\frac{2}{h_0}}cos\left(narccosx+\overline{\gamma }\left(x\right)\right)$$

(78)

form a set of polynomials orthonormal with respect to the weight $h_0\left(x\right):=\pi \sqrt{1-x^2}h\left(x\right)$, and where $\overline{\gamma }\left(x\right)$ is the Hilbert transform of the function $ln{h}_0\left(x\right)$. Second, we take into account the following Lemma [34].

Lemma 1.

Let g be a continuous, periodic function on $\mathbb{R}$, f be an integrable and $\tilde{\gamma }$ be a measurable function on $(0,\pi )$:

$$g\in C\left(\mathbb{R}\right),g(\theta +\pi )=g\left(\theta \right),f\in L^1(0,\pi ),\tilde{\gamma }<\mathrm{const}a.e.on(0,\pi ).$$

Then, the following holds:

$$\underset{0}{\overset{\pi }{\int }}g\left(n\theta +\tilde{\gamma }\left(\theta \right)\right)f\left(\theta \right)d\theta \to \frac{1}{\pi }\underset{0}{\overset{\pi }{\int }}g\left(\theta \right)d\theta \underset{0}{\overset{\pi }{\int }}f\left(\theta \right)d\theta ,n\to \infty .$$

(79)

Now, with these two previous results, Aptekarev et al. [78] proved in 2021 the following powerful result.

Theorem 3.

Let $\left\{{\widehat{p}}_n\left(x\right)\right\}$ be the polynomials (74) orthonormal on $[-1,1]$ with respect to the weight function $h\left(x\right)$ for which (a) has singularities at the end points $\pm 1$ of the order ${(1-x)}^{{\alpha }_h}$ and ${(1+x)}^{{\beta }_h}$, respectively, with ${\alpha }_h\ge {\beta }_h>-1$, and (b) satisfies the Szegö condition (76). Let also the following hold:

$$F\left(x\right):={(1-x)}^{\alpha }{(1+x)}^{\beta }f\left(x\right),\alpha \ge \beta >-1,f\left(x\right)\in {\mathfrak{L}}_1(-1,1),$$

(80)

where $f\left(x\right)$ may have a possible growth at the end points $\pm 1$ but not faster than logarithmic. Then, we have the following asymptotical behavior:

$$N_n^{\left(q\right)}(h,F):=\underset{-1}{\overset{1}{\int }}{\left[{\widehat{p}}_n^2\left(x\right)h\left(x\right)\right]}^{q}F\left(x\right)dx\underset{n\to \infty }{⟶}\frac{2^q\overline{\Gamma }\left(q\right)}{{\pi }^{q+1}}\underset{-1}{\overset{1}{\int }}{(1-x^2)}^{-\frac{q}{2}}F\left(x\right)dx,$$

(81)

for

$$-\frac{\beta +1}{{\beta }_h}<q<\beta +1$$

(82)

and where

$$\overline{\Gamma }\left(q\right):=\frac{\Gamma (q+1/2)\Gamma (1/2)}{\Gamma (q+1)}.$$

(83)

Moreover, when $h\left(x\right)$ and $F\left(x\right)$ are the Jacobi weights, we can find, using the methodology of [34], the following asymptotical result [78] in a straightforward manner.

Theorem 4.

Let $\left\{{\widehat{p}}_n\left(x\right)\right\}$ be a system of Jacobi polynomials, orthonormal with respect to an arbitrary Jacobi weight function $h\left(x\right):={(1-x)}^{{\alpha }_h}{(1+x)}^{{\beta }_h},{\alpha }_h\ge {\beta }_h>-1$, and $F\left(x\right)$ be the Jacobi weight:

$$F\left(x\right):={(1-x)}^{\alpha }{(1+x)}^{\beta },\alpha \ge \beta >-1.$$

Then, the generalized weighted norm of the Jacobi polynomials

$$N_n^{\left(q\right)}(h,F):=\underset{-1}{\overset{1}{\int }}{\left[{\widehat{p}}_n^2\left(x\right)h\left(x\right)\right]}^{q}F\left(x\right)dx,$$

has the following limiting ($n\to \infty$) behavior.

A. For $-\frac{\beta +1}{{\beta }_h}<q<2(\beta +1)$

$$\begin{array}{cc}\hfill N_n^{\left(q\right)}(h,F)& ⟶\frac{2^q}{{\pi }^{q+1}}\overline{\Gamma }\left(q\right)\underset{-1}{\overset{1}{\int }}{(1-x^2)}^{-\frac{q}{2}}F\left(x\right)dx=\hfill \\ \hfill & =\frac{2^{\alpha +\beta +1}}{{\pi }^{q+1}}\frac{\Gamma (q+\frac{1}{2})\Gamma \left(\frac{1}{2}\right)}{\Gamma (q+1)}\frac{\Gamma (\alpha -\frac{q}{2}+1)\Gamma (\beta -\frac{q}{2}+1)}{\Gamma (\alpha +\beta -q+2)};\hfill \end{array}$$

(84)

B. For $q=2(\beta +1)$

$$N_n^{\left(q\right)}(h,F)\sim lnn,i.e.\exists c_0:N_n^{\left(q\right)}/lnn\to c_0;$$

(85)

C. For $q>2(\beta +1)$

$$N_n^{\left(q\right)}(h,F)\asymp n^{q-2(\beta +1)},i.e.\exists c_1,c_2:c_1\le N_n^{\left(q\right)}/n^{q-2(\beta +1)}\le c_2.$$

(86)

Note the diversity of the asymptotics of $N_n^{\left(q\right)}$ when $n\to \infty$ for various q. For q relatively small to the parameter $\beta$, the integrals remain bounded for big n. When q approaches its critical value $2(\beta +1)$, the right hand side of the bound (84) explodes since the argument of the multiplier $\Gamma (\beta -\frac{q}{2}+1)$ is close to the pole of gamma function. Then, for the critical $q=2(\beta +1)$, the bound (84) does not work any more and according to (65), a new asymptotical regime takes place where $N_n^{\left(q\right)}$ grows as $lnn$. Finally, when $q>2(\beta +1)$, the norms shows a n-degree exponential scale governed by the positive difference between q and its critical value.

It is remarkable the great simplicity of the main term of asymptotics ($n\to \infty$) of the weighted ${\mathfrak{L}}_q$-norms of the Jacobi polynomials, $W_q\left[{\widehat{P}}_n^{(\alpha ,\beta )}\right]$, given as particular instances of the expressions (84)–(86). Indeed, this result is impressively simple, given the complexity of the expressions (55) and (57) obtained in Section 5 for the corresponding quantities of a Jacobi polynomial with fixed degree n, which require the evaluation of either the multivariate Bell polynomials (34) at the Jacobi expansion coefficients (56), or the multivariate Srivastava–Daoust of the type (58) at unity, respectively.

These mathematical results were recently used [78] to evaluate in a compact and elegant way the Rényi and Shannon entropies of the Rydberg hydrogenic three and multidimensional states. This was possible because the radial and angular wave functions of these quantum states are controlled by the Gegenbauer polynomials, which are instances of the great family of Jacobi polynomials.

9. Shannon’s Spreading Length and Logarithmic Potential of HOPs

In this section, we examine the Shannon spreading length ${\mathcal{L}}_1^S\left[p_n\right]={\mathcal{L}}_1^S\left[{\rho }_n\right]$ of the three canonical families of the HOPs $\left\{p_n\left(x\right)\right\}$ for all n by means of the logarithmic-potential-based approach, which is the only existing one with general validity, although it requires the evaluation of the logarithmic potentials of the polynomials at their zeros on the orthogonality support. The latter, however, is a non-trivial task which makes necessary, almost mandatory, the consideration and analysis of the asymptotics ($n\to \infty$) for the Shannon entropy of the HOPs; this is shown in the next section. Likewise, it is helpful to obtain analytical upper bounds to ${\mathcal{L}}_1^S\left[p_n\right]$ as simply and accurately as possible; some results, which are briefly pointed out at the end of this section, were found for Hermite, Laguerre and Jacobi polynomials [22,23,24] in terms of the ordinary moments of their associated Rakhmanov density (1).

According to Equations (7), (10) and (11), this quantity is given by the Shannon spreading length of the associated Rakhmanov probability density ${\rho }_n\left(x\right)$ given by (1), that is, the following:

$${\mathcal{L}}_1^S\left[{\rho }_n\right]=\underset{q\to 1}{lim}{\mathcal{L}}_q^R\left[{\rho }_n\right]=exp\left(S\left[{\rho }_n\right]\right)=exp\left\{-{\int }_{\Delta }{\rho }_n\left(x\right)ln{\rho }_n\left(x\right)dx\right\},$$

(87)

where $S\left[{\rho }_n\right]$ denotes the Shannon-like integral functional $S\left[p_n\right]$ of the polynomials given by the following:

$$S\left[p_n\right]=S\left[{\rho }_n\right]=-{\int }_{\Lambda }\frac{1}{{\kappa }_n}{p}_n^2\left(x\right)h\left(x\right)ln\left[\frac{1}{{\kappa }_n}{p}_n^2\left(x\right)h\left(x\right)\right]dx=ln{\kappa }_n+\frac{1}{{\kappa }_n}\left(E\left[p_n\right]+I\left[p_n\right]\right),$$

(88)

with the polynomial functionals

$$I\left[p_n\right]:=-{\int }_{\Lambda }{p}_n^2\left(x\right)h\left(x\right)lnh\left(x\right)dx$$

(89)

and

$$E\left[p_n\right]:=-{\int }_{\Lambda }{p}_n^2\left(x\right)h\left(x\right)ln{p}_n^2\left(x\right)dx.$$

(90)

The functional $I\left[p_n\right]$ was determined [79] for all HOPs in terms of the coefficients of its second-order differential equation. However, the explicit computation of the functional $E\left[p_n\right]$ for all n in terms of the degree and the parameters of the weight function $h\left(x\right)$ is a formidable task, not yet solved for the HOPs despite the efforts of many authors, except (a) for a few specific families of the Jacobi polynomials, such as the Chebyshev polynomials of the first and second type, and (b) in some asymptotical cases, such as when ($n\to \infty$) or when the parameters of the weight $h\left(x\right)$ go toward ∞. This functional $E\left[p_n\right]$ is called Shannon entropy of the HOPs$p_n\left(x\right)$. Note that, to fix the notation, the Shannon-like integral functional $S\left[{\widehat{p}}_n\right]$ of the orthonormal polynomial ${\widehat{p}}_n$ is the following:

$$S\left[{\widehat{p}}_n\right]=-{\int }_{\Lambda }{\widehat{p}}_n^2\left(x\right)h\left(x\right)ln\left[{\widehat{p}}_n^2\left(x\right)h\left(x\right)\right]dx=E\left[{\widehat{p}}_n\right]+I\left[{\widehat{p}}_n\right],$$

(91)

and the Shannon entropy of ${\widehat{p}}_n\left(x\right)$ is the following:

$$E\left[{\widehat{p}}_n\right]:=-{\int }_{\Lambda }{\widehat{p}}_n^2\left(x\right)h\left(x\right)ln{\widehat{p}}_n^2\left(x\right)dx,$$

(92)

Then, for illustration, the Shannon spreading length for the orthonormal Hermite polynomials ${\widehat{H}}_n\left(x\right)$ is given by the following:

$${\mathcal{L}}_1^S\left[{\widehat{H}}_n\right]=exp\left\{S\left[{\widehat{H}}_n\right]\right\},$$

(93)

where the Shannon-like integral functional of ${\widehat{H}}_n\left(x\right)$ is the following:

$$S\left[{\widehat{H}}_n\right]=-{\int }_{\Lambda }{\widehat{H}}_n^2\left(x\right)h_H\left(x\right)ln\left[{\widehat{H}}_n^2\left(x\right)h_H\left(x\right)\right]dx=E\left[{\widehat{H}}_n\right]+n+\frac{1}{2},$$

(94)

where $E\left[{\widehat{H}}_n\right]$ denotes the Shannon entropy of ${\widehat{H}}_n\left(x\right)$ defined as follows:

$$E\left[{\widehat{H}}_n\right]:=-{\int }_{\Lambda }{\widehat{H}}_n^2\left(x\right)h_H\left(x\right)ln{\widehat{H}}_n^2\left(x\right)dx,$$

(95)

whose explicit value for all n has not yet been found despite serious attempts [10,34,80]. Only the asymptotic case ($n\to \infty$) has been determined in a simple compact way as discussed in the next section.

Now, the Shannon spreading length of the orthonormal Laguerre polynomials is given by the following:

$${\mathcal{L}}_1^S\left[{\widehat{L}}_n^{\left(\alpha \right)}\right]=exp\left\{S\left[{\widehat{L}}_n^{\left(\alpha \right)}\right]\right\},$$

(96)

where the Shanon-like integral functional of the orthonormal Laguerre polynomials ${\widehat{L}}_n^{(\alpha }\left(x\right)$ is the following:

$$\begin{array}{ccc}\hfill S\left[{\widehat{L}}_n^{\left(\alpha \right)}\right]& =& -{\int }_{-1}^{+1}{\left[{\widehat{L}}_n^{\left(\alpha \right)}\left(x\right)\right]}^2{h}_{\alpha }^L\left(x\right)ln\left\{{\left[{\widehat{L}}_n^{\left(\alpha \right)}\left(x\right)\right]}^2{h}_{\alpha }^L\left(x\right)\right\}dx=E\left[{\widehat{L}}_n^{\left(\alpha \right)}\right]+I\left[{\widehat{L}}_n^{\left(\alpha \right)}\right],\hfill \end{array}$$

with the integral functional as follows [24,79]:

$$\begin{array}{ccc}\hfill I\left[{\widehat{L}}_n^{\left(\alpha \right)}\right]& =& -{\int }_0^{\infty }{\left[{\widehat{L}}_n^{\left(\alpha \right)}\left(x\right)\right]}^2{h}_{\alpha }^L\left(x\right)ln{h}_{\alpha }^L\left(x\right)=2n+\alpha +1-\alpha \psi (\alpha +n+1)\hfill \end{array}$$

(97)

and the Shannon entropy of ${\widehat{L}}_n^{(\alpha }\left(x\right)$ is given by the following:

$$\begin{array}{ccc}\hfill E\left[{\widehat{L}}_n^{\left(\alpha \right)}\right]& =& -{\int }_0^{\infty }{\left[{\widehat{L}}_n^{\left(\alpha \right)}\left(x\right)\right]}^2{h}_{\alpha }^L\left(x\right)ln{\left[{\widehat{L}}_n^{\left(\alpha \right)}\left(x\right)\right]}^{2}dx,\hfill \end{array}$$

(98)

whose explicit value for all n has not yet been found in spite of different attempts [10,34,80]; only the asymptotic case ($n\to \infty$) has been determined in a simple compact way as is discussed in the next section.

Finally, the Shannon spreading length of the orthonormal Jacobi polynomials ${\widehat{P}}_n^{(\alpha ,\beta )}\left(x\right)$, which is given by the following:

$${\mathcal{L}}_1^S\left[{\widehat{P}}_n^{(\alpha ,\beta )}\right]=exp\left\{S\left[{\widehat{P}}_n^{(\alpha ,\beta )}\right]\right\},$$

(99)

where the Shannon-like integral functional of ${\widehat{P}}_n^{(\alpha ,\beta )}\left(x\right)$ is the following:

$$\begin{array}{ccc}\hfill S\left[{\widehat{P}}_n^{(\alpha ,\beta )}\right]& =& -{\int }_{-1}^{+1}{\left[{\widehat{P}}_n^{(\alpha ,\beta )}\left(x\right)\right]}^2{h}_{\alpha ,\beta }^J\left(x\right)ln\left\{{\left[{\widehat{P}}_n^{(\alpha ,\beta )}\left(x\right)\right]}^2{h}_{\alpha ,\beta }^J\left(x\right)\right\}dx\hfill \\ & =& E\left[{\widehat{P}}_n^{(\alpha ,\beta )}\right]+I\left[{\widehat{P}}_n^{(\alpha ,\beta )}\right],\hfill \end{array}$$

with the integral functional as follows [57,79]:

$$\begin{array}{ccc}\hfill I\left[{\widehat{P}}_n^{(\alpha ,\beta )}\right]& =& -{\int }_{-1}^{+1}{\left[{\widehat{P}}_n^{(\alpha ,\beta )}\left(x\right)\right]}^2{h}_{\alpha ,\beta }^J\left(x\right)ln{h}_{\alpha ,\beta }^J\left(x\right)\hfill \\ \hfill & =& -\alpha \psi (n+\alpha +1)-\beta \psi (n+\beta +1)+(\alpha +\beta )\hfill \end{array}$$

(100)

$$\begin{array}{ccc}& \times & \left[-ln2+\frac{1}{2n+\alpha +\beta +1}+2\psi (2n+\alpha +\beta +1)-\psi (n+\alpha +\beta +1)\right],\hfill \end{array}$$

(101)

and the Shannon entropy of ${\widehat{P}}_n^{(\alpha ,\beta )}\left(x\right)$ is given by the following:

$$\begin{array}{ccc}\hfill E\left[{\widehat{P}}_n^{(\alpha ,\beta )}\right]& =& -{\int }_{-1}^{+1}{\left[{\widehat{P}}_n^{(\alpha ,\beta )}\left(x\right)\right]}^2{h}_{\alpha ,\beta }^J\left(x\right)ln{\left[{\widehat{P}}_n^{(\alpha ,\beta )}\left(x\right)\right]}^{2}dx\hfill \end{array}$$

(102)

The explicit determination of this entropic functional for the Jacobi polynomials was not known up until now, except in the asymptotic case (whose value is given in the next section) and when the parameters $(\alpha ,\beta )$ are equal to $(-1/2,-1/2)$ and $(1/2,1/2)$, which correspond to the families of Chebyshev polynomials of the first and second kind, denoted by $T_n\left(x\right)$ and $U_n\left(x\right)$, respectively. The following values were found [10,80]:

$$\begin{array}{cc}\hfill E\left[T_n\right]& =-\frac{1}{\pi }{\int }_{-1}^{+1}{T}_n^2\left(x\right)ln{T}_n^2\left(x\right)\frac{dx}{\sqrt{1-x^2}}=ln2-1;\mathrm{for}n\ge 1\hfill \end{array}$$

for the Shannon entropy of the orthogonal Chebyshev polynomials of the first kind, and the values [80]

$$\begin{array}{cc}\hfill E\left[U_n\right]& =-\frac{2}{\pi }{\int }_{-1}^{+1}{U}_n^2\left(x\right)ln{U}_n^2\left(x\right)\sqrt{1-x^2}dx=\frac{1}{n+1}-1;\mathrm{for}n\ge 0\hfill \end{array}$$

for the Shannon entropy of the orthogonal Chebyshev polynomials of the second kind. Moreover, when the parameters $(\alpha ,\beta )$ are equal to $(\lambda -1/2,\lambda -1/2)$ the corresponding Jacobi polynomials conform to the family of Gegenbauer polynomials $C_n^{\left(\lambda \right)}\left(x\right)$. Then, it is known the value for the Shannon entropy of the modified Gegenbauer polynomials $G_n^{(\lambda =2)}\left(x\right)$, which are the following polynomials:

$$\begin{array}{cc}\hfill G_k^{\left(\lambda \right)}\left(x\right)& ={\left(\frac{k!(k+\lambda )\Gamma \left(2\lambda \right)}{\lambda \Gamma (k+2\lambda )}\right)}^{\frac{1}{2}}{C}_k^{\left(\lambda \right)}\left(x\right)={\gamma }_k^{\lambda }{x}^k+\mathrm{lower}\mathrm{degree}\mathrm{terms}\hfill \end{array}$$

orthogonal with respect to the positive unit weight on $[-1,+1]$,

$$h_{\lambda }\left(x\right)=\frac{\Gamma (\lambda +1)}{\sqrt{\pi }\Gamma (\lambda +1/2)}{(1-x^2)}^{\lambda -\frac{1}{2}}.$$

Buyarov et al. [81] (see also [11,82]) found the following:

$$\begin{array}{c}\hfill E\left[G_k^{\left(\lambda \right)}\right]=-{\int }_{-1}^{+1}{\left[G_k^{\left(\lambda \right)}\left(x\right)\right]}^{2}ln{\left[G_k^{\left(\lambda \right)}\left(x\right)\right]}^2{h}_{\lambda }\left(x\right)dx\\ \hfill =-ln\left(\frac{3(n+1)}{n+3}\right)-\frac{n^3-5n^2-29n-27}{(n+1)(n+2)(n+3)}-\frac{1}{n+2}{\left(\frac{n+3}{n+1}\right)}^{n+2},\mathrm{for}\lambda =2.\end{array}$$

This expression was extended for Gegenbauer polynomials with integer order $\lambda$ (see [81,83] and the references therein).

In the general case, the Shannon entropy of the orthonormal polynomials ${\widehat{p}}_n\left(x\right)$ on the real line can be expressed in terms of the logarithmic potential of its Rakhmanov density ${\rho }_n\left(x\right)={\widehat{p}}_n^2\left(x\right)h\left(x\right)$ and, equivalently, by means of the mutual logarithmic energy between the normalized zero counting distribution ${\lambda }_n\left(x\right)=\frac{1}{n}{\displaystyle \sum _{j=1}^n}\delta \left(x-x_{j,n}\right)$ and the Rakhmanov density ${\rho }_n\left(x\right)$. Indeed, let $x_{j,n}(j=1,...,n)$ be the zeros of ${\widehat{p}}_n\left(x\right)$ and ${\gamma }_n$ the leading coefficient, so that the following holds:

$${\widehat{p}}_n\left(x\right))={\gamma }_n\prod _{j=1}^n\left(x-x_{j,n}\right),{\gamma }_n>0.$$

Then, it can be proved that the Shannon entropy of the orthonormal polynomials ${\widehat{p}}_n\left(x\right)$ can be expressed [80] as follows:

$$E\left[{\widehat{p}}_n\right]:=-2ln{\gamma }_n+2n\int \mathcal{V}(x;{\lambda }_n){\widehat{p}}_n^2\left(x\right)h\left(x\right)dx,$$

(103)

where

$$\mathcal{V}(z;\mu )=\int ln\frac{1}{|z-x|}d\mu \left(x\right)$$

is known as the logarithmic potential of the probability measure $\mu$ [84]. Moreover, taking into account the mutual logarithmic energy of the two measures $\mu$ and $\nu$, given [84] by

$$I[\nu ,\mu ]=\int \mathcal{V}(z;\nu )d\mu (x)=-\int \int ln|z-t|d\nu (t)d\mu (t),$$

one has that the Shannon entropy of the orthonormal polynomials can be written [80] as follows:

$$E\left[{\widehat{p}}_n\right]=-2ln{\gamma }_n+2nI\left[{\lambda }_n,{\rho }_n\right]$$

When the support interval of $\mu \left(x\right)$ is compact, the potential $\mathcal{V}(z;\mu )$ oscillates around its Robin (or extremal) constant. The zeros $x_{j,n}$ turn out to be [80] points of local minima for the potential $\mathcal{V}(x;{\rho }_n)$; this implies that to calculate $E\left[{\widehat{p}}_n\right]$, we have to sum up the values of the logarithmic potential $\mathcal{V}(x;{\rho }_n)$ at its local minima. Then, by applying this expression to the Jacobi polynomials and taking into account Fubini’s theorem, one finally has the Shannon entropy of the orthonormal Jacobi polynomials ${\widehat{P}}_n^{(\alpha ,\beta )}\left(x\right)$:

$$E\left[{\widehat{P}}_n^{(\alpha ,\beta )}\right]:=-2ln{\gamma }_n-2\sum _{j=1}^n\mathcal{V}(x_{j,n};{\rho }_{n,\alpha ,\beta }),$$

(104)

in terms of the following logarithmic potential evaluated at the zeros of the polynomial:

$$\mathcal{V}(z;{\rho }_{n,\alpha ,\beta })=\int {\widehat{P}}_n^{(\alpha ,\beta )}{\left(x\right)}^2{h}_{\alpha ,\beta }^J\left(x\right)ln\frac{1}{|z-x|}dx$$

Similarly, one can obtain the Shannon entropy of the remaining HOPs in terms of the values of their logarithmic potentials evaluated at their zeros, such as the Hermite, Laguerre and Gegenbauer polynomials, as explained in detail in [11,36,80,85]. Additionally, the Shannon entropy of the HOPs was calculated [58] in terms of the eigenvalues of their associated Jacobi matrix and the ordinary moments (15) of their corresponding Rakhmanov density (1). In practice, however, we have a nontrivial task in both cases, especially when the degree of the polynomials is not sufficiently low. Then, it appears that the determination of the asymptotical ($n\to \infty$) values of the Shannon entropy of the HOPs is necessary, which is shown in the next section.

Since the expressions obtained for the Shannon entropy $E\left[p_n\right]$ and, consequently, for the Shannon spreading lengths ${\mathcal{L}}_1^S\left[p_n\right]$ of the HOPs are not easily manipulated analytically, it is also interesting to obtain simple, compact and accurate upper bounds to these quantities in terms of the degree and the parameters of their weight function. Some upper bounds were found for Hermite [22], Laguerre [24] and Jacobi [23] polynomials by means of the ordinary moments of their associated Rakhmanov densities (1). In the Hermite and Laguerre cases, we used an optimized information-theoretical technique to obtain the upper bounds based on the non-negativity of the Kullback–Leibler functional of the Rakhmanov density of the polynomial and a probability density of exponentially decreasing type of the form $exp\left(-x^k)\right)$. It is found that the following holds:

$${\mathcal{L}}_1^S\left[{\widehat{H}}_n\right]\le \frac{2{\left(ek\right)}^{\frac{1}{k}}}{k}\Gamma \left(\frac{1}{k}\right){\left({\langle x^k\rangle }_n\right)}^{1/k};k=2,4,\dots$$

and

$${\mathcal{L}}_1^S\left[{\widehat{L}}_n^{\left(\alpha \right)}\right]\le \frac{{\left(ke\right)}^{\frac{1}{k}}}{k}\Gamma \left(\frac{1}{k}\right){\left({\langle x^k\rangle }_{n,\alpha }\right)}^{1/k};k>0,$$

for the Shannon spreading length of Hermite and Laguerre polynomials, respectively, where the ordinary moments of such polynomials are given by (17) and (19). In the Jacobi case, beyond the general upper bound ${\mathcal{L}}_1^S\left[p_n\right]\le 2$ valid for any probability density with the support interval $[-1,+1]$, some variational bounds were found by means of the following ordinary and logarithmic moments: $\langle x\rangle$, $\langle x^2\rangle$, $\langle ln{x}^2\rangle$, $\langle ln(1\pm x)\rangle$, $\langle ln(1-x^2)\rangle$. In particular, it is known that

$${\mathcal{L}}_1^S\left[{\widehat{P}}_n^{(\alpha ,\beta )}\right]\le \frac{2sinh\left({\lambda }_2\right)}{{\lambda }_2}exp(-1-{\lambda }_{2}coth\left({\lambda }_2\right)),$$

(105)

where ${\lambda }_2$ is given in terms of ${\langle x\rangle }_{n,\alpha ,\beta }$ by means of the following implicit equation:

$$\frac{1}{{\lambda }_2}-coth\left({\lambda }_2\right)={\langle x\rangle }_{n,\alpha ,\beta }.$$

(106)

where, according to (21), the first order moment ${\langle x\rangle }_{n,\alpha ,\beta }$ of the Rakhmanov density ${\rho }_{n,\alpha ,\beta }\left(x\right)$ of the Jacobi polynomials is given explicitly by the following:

$${\langle x\rangle }_{n,\alpha ,\beta }=\frac{{\beta }^2-{\alpha }^2}{(2n+\alpha +\beta )(2n+\alpha +\beta +2)}.$$

(107)

These variational bounds to ${\mathcal{L}}_1^S\left[{\widehat{P}}_n^{(\alpha ,\beta )}\right]$ have been analytically and numerically discussed in detail [23].

Finally, let us also remark that Beckermann et al. [86] found some asymptotic upper bounds to the Shannon entropy $E\left[{\widehat{p}}_n\right]$ for polynomials ${\widehat{p}}_n\left(x\right)$ orthonormal with respect to a weight function $h\left(x\right)$ on the support $[-1,+1]$, which belongs to the Szegö class, although the determination of the entropy itself remains an open problem. See [43,87,88,89] for recent extensions of these results.

10. Shannon Entropy $\mathit{E}\left[{\mathit{p}}_{\mathit{n}}\right]$ and ${\mathfrak{L}}_{\mathbf{q}}$-Norms ${\mathcal{N}}_{\mathbf{q}}\left[{\mathit{p}}_{\mathit{n}}\right]$ of HOPs Degree Asymptotics

In this section, we tackle and solve the asymptotics ($n\to \infty$) for the Shannon entropy of the three canonical families of the HOPs and the Freud polynomials (which is a generalization of the Hermite polynomials). This study was initiated in 1994 [33,34,35] and reviewed in 2001 [11] and 2010 [36].

The following approach was used. Taking into account the limiting relation (3), we can convert the calculation of this Shannon entropy of the HOPs to the evaluation of the (unweighted) ${\mathfrak{L}}_q$-norms ${\mathcal{N}}_q\left[{\widehat{p}}_n\right]$ since the following holds:

$$\begin{array}{cc}\hfill E\left[{\widehat{p}}_n\right]& =-\underset{q\to 1}{lim}\frac{1}{q-1}ln\int {\left|{\widehat{p}}_n\left(x\right)\right|}^{2q}h\left(x\right)dx=-\underset{q\to 1}{lim}\frac{\partial }{\partial q}{\mathcal{N}}_{2q}\left[{\widehat{p}}_n\right],\hfill \end{array}$$

(108)

where

$${\mathcal{N}}_q\left[{\widehat{p}}_n\right]:=\int {\left|{\widehat{p}}_n\left(x\right)\right|}^{q}h\left(x\right)dx$$

(109)

denotes the (unweighted) ${\mathfrak{L}}_q\left(h\right)$ norms of ${\widehat{p}}_n\left(x\right)$. Then, the asymptotical value for the Shannon entropy $E\left[{\widehat{p}}_n\right]$ is as follows [34,36]:

$$\begin{array}{cc}\hfill E_{\infty }& =\underset{n\to \infty }{lim}E\left[{\widehat{p}}_n\right]=-\underset{n\to \infty }{lim}{\left.\frac{\partial {\mathcal{N}}_{2q}\left[{\widehat{p}}_n\right]}{\partial q}\right|}_{q=1}={\left.-{\mathcal{N}}_{2q}^{\prime }\right|}_{q=1}\hfill \end{array}$$

(110)

This approach based on the (unweighted) ${\mathfrak{L}}_q$-norms can be used for the HOPs with both bounded and unbounded supports. However, for an absolutely continuous weight function $h\left(x\right)$ on the real line, it is more convenient to calculate the Shannon entropy $E\left[{\widehat{p}}_n\right]$ from the weighted ${\mathfrak{L}}_q$-norms $W_q\left[{\widehat{p}}_n\right]$ given by (30) so that, according to (91) and (89), one has the following [34]:

$$E\left[{\widehat{p}}_n\right]=S\left[{\widehat{p}}_n\right]-I\left[{\widehat{p}}_n\right]=-{\left.\frac{\partial W_q\left[{\widehat{p}}_n\right]}{\partial q}\right|}_{q=1}-I\left[{\widehat{p}}_n\right].$$

(111)

Thus, the asymptotical value $E_{\infty }$ for the Shannon entropy is basically controlled by the asymptotics of $W_q\left[{\widehat{p}}_n\right]$. Let us now apply this second approach to the Freud polynomials $F_n^{\left(\kappa \right)}\left(x\right)$, which are orthogonal with respect to the weight function $h_{\kappa }\left(x\right)=e^{-{\left|x\right|}^{\kappa }}$, $\kappa >1$, $\infty <x<+\infty$. In this case we have that the Shannon entropy is the following:

$$E\left[F_n^{\left(\kappa \right)}\right]=-{\left.\frac{\partial W_q\left[F_n^{\left(\kappa \right)}\right]}{\partial q}\right|}_{q=1}-I\left[F_n^{\left(\kappa \right)}\right]=-{\left.\frac{\partial W_q\left[F_n^{\left(\kappa \right)}\right]}{\partial q}\right|}_{q=1}-\frac{2n+1}{\kappa }$$

(112)

Now, the asymptotical behavior ($n\to \infty$) of the weighted ${\mathfrak{L}}_q$-norms $W_q\left[F_n\right]$ is found to be the following:

$$W_q\left[{\widehat{F}}_n^{\left(\kappa \right)}\right]={\left(\frac{2}{\pi }\right)}^q\frac{\Gamma (q+1/2)}{\Gamma (q+1)}\frac{\Gamma (1-q/2)}{\Gamma (3/2-q/2)}{\alpha }_n^{1-q}(1+o\left(1\right))$$

(113)

with

$${\alpha }_n={\left(\frac{2n+1}{2\beta }\right)}^{1/\kappa }\mathrm{and}\beta =\frac{\Gamma (\kappa /2+1/2)}{\Gamma (1/2)\Gamma (\kappa /2)}.$$

(114)

Then, from the last three expressions, we have that the ($n\to \infty$)-asymptotics of the Shannon entropy of the Freud polynomials [35] is given by the following:

$$\begin{array}{cc}\hfill E\left[{\widehat{F}}_n^{\left(\kappa \right)}\right]=& =-\frac{2n+1}{\kappa }+\frac{1}{\kappa }ln\left(2n\right)-\frac{1}{\kappa }ln\left(\frac{\sqrt{\pi }\Gamma (\kappa /2)}{2\Gamma (\kappa /2+1/2)}\right)+ln\pi -1+o\left(1\right)\hfill \end{array}$$

Consequently, for $\kappa =2$ we have that the asymptotics ($n\to \infty$) for the Shannon entropy of the Hermite polynomials $H_n\left(x\right)$ is given by the following:

$$E[{\widehat{H}}_n]=-n+ln\sqrt{2n}-\frac{3}{2}+ln\pi +o\left(1\right)$$

(115)

The same technique can be used to find the asymptotic behaviour ($n\to \infty$) for the Shannon entropy of the orthonormal Laguerre polynomials ${\widehat{L}}_n^{\left(\alpha \right)}\left(x\right)$, obtaining [90] for fixed real $\alpha >-1$ such that the following holds:

$$E\left[{\widehat{L}}_n^{\left(\alpha \right)}\right]=-2n+(\alpha +1)lnn-\alpha -2+ln\pi +o\left(1\right)$$

(116)

Similar operations with the orthonormal Jacobi polynomials ${\widehat{P}}_n^{(\alpha ,\beta )}\left(x\right)$ by using the (unweighted) ${\mathfrak{L}}_q$-norms ${\mathcal{N}}_q\left[{\widehat{p}}_n\right]$ given by (109) allowed us, according to (108), to find the following expression:

$$E\left[{\widehat{P}}_n^{(\alpha ,\beta )}\right]=-2n+(\alpha +1)lnn-\alpha -2+ln\pi +o\left(1\right)$$

(117)

for the asymptotic behavior of the Shannon entropy of the orthonormal Jacobi polynomials. A particularly relevant family of the Jacobi polynomials are the Gegenbauer polynomials $C_n^{\left(\lambda \right)}\left(x\right)$, which are orthogonal with respect to the weight function $h_{\lambda }\left(x\right)={(1-x^2)}^{\lambda -\frac{1}{2}}$ on the interval $[-1,+1]$; then, the orthonormal Gegenbauer polynomial ${\widehat{C}}_n^{\left(\lambda \right)}\left(x\right)$ are as follows:

$${\tilde{C}}_n^{\left(\lambda \right)}\left(x\right)=\frac{C_n^{\left(\lambda \right)}}{{\kappa }_n};\mathrm{with}{\kappa }_n^2=\frac{2^{1-2\lambda }\pi \Gamma (n+2\lambda )}{{\left[\Gamma \left(\lambda \right)\right]}^2(n+\lambda )n!},$$

(118)

so that, according to (109), the Shannon entropy of ${\widehat{C}}_n^{\left(\lambda \right)}\left(x\right)$ has the following asymptotic expression [34]:

$$\begin{array}{cc}\hfill E\left({\widehat{C}}_m^{\alpha }\right)& \equiv -{\int }_{-1}^{+1}{h}_{\alpha }\left(x\right){\left[{\widehat{C}}_m^{\alpha }\left(x\right)\right]}^{2}ln{\left[{\widehat{C}}_m^{\alpha }\left(x\right)\right]}^{2}dx=ln\pi +(1-2\alpha )ln2-1+o\left(1\right),\hfill \end{array}$$

(119)

for fixed $\lambda$ and large degree n. This asymptotic expression was improved by obtaining further terms as shown in [81,83]. Further details about the proof of these results and their extensions to more general orthonormal polynomials belonging to the Bernstein and Szegö classes can be seen in [34,36,86,87,88,89].

Moreover, from the previous expressions (115)–(117) for the asymptotics of the Shannon entropy $E\left[{\widehat{p}}_n\right]$ of HOPS and the exact expressions (94), (97) and (100) for the integrals $I\left[{\widehat{p}}_n\right]$ given by (89), we have the following asymptotics ($n\to \infty$) for the Shannon-like integral functionals $S\left[{\widehat{p}}_n\right]$ of the HOPS [34]:

$$S[{\widehat{H}}_n]=ln\sqrt{2n}+ln\pi -1+o\left(1\right),S[{\widehat{L}}_n^{\left(\alpha \right)}]=(\alpha +1)lnn-\alpha \psi (\alpha +n+1)-1+ln\left(2\pi \right)+o\left(1\right),$$

and

$$S[{\widehat{P}}_n^{(\alpha ,\beta )}]=ln\pi -1+o\left(1\right)$$

for the Hermite, Laguerre and Jacobi polynomials, respectively. To obtain the last expression, we have also taken into account that ${lim}_{n\to \infty }I\left[{\tilde{P}}_n^{(\alpha ,\beta )}\right]=(\alpha +\beta )ln2$. Then, the Shannon spreading length ${\mathcal{L}}_1^S\left[{\widehat{p}}_n\right]={\mathcal{L}}_1^S\left[{\rho }_n\right]$ of the three canonical families of the HOPs $\left\{{\widehat{p}}_n\left(x\right)\right\}$ have the following asymptotical values:

$${\mathcal{L}}_1^S[{\widehat{H}}_n]\simeq \frac{\pi }{e}\sqrt{2n},{\mathcal{L}}_1^S[{\widehat{L}}_n^{\left(\alpha \right)}]\simeq \frac{\pi }{e}2n,\mathrm{and}{\mathcal{L}}_1^S[{\widehat{P}}_n^{(\alpha ,\beta )}]\simeq \frac{\pi }{e}$$

(120)

for Hermite, Laguerre and Jacobi polynomials, respectively. By comparing these expressions with the asymptotical values of the standard deviation $\Delta x$ given by Equations (18), (20) and (22) for Hermite, Laguerre and Jacobi polynomials, respectively, one is led to the following linear relation:

$${\mathcal{L}}_1^S\left[{\widehat{p}}_n\right]\simeq \frac{\pi \sqrt{2}}{e}\Delta x\left[{\widehat{p}}_n\right]\simeq 1.6389\Delta x\left[{\widehat{p}}_n\right],n\to \infty$$

(121)

which holds for the three canonical families of HOPs [91] (see also [22,23,24]); note that it fulfills the general Cramér-Rao relation (9). In fact, the validity of this relation for HOPs extends to the whole class of Bernstein-Szegö polynomials, although it fails for arbitrary orthogonal polynomials. For example, it is violated for Freud polynomials for which the Shannon length and the standard deviation have a quadratic relation [91].

11. Conclusions

Hypergeometric orthogonal polynomials (HOPs) are the simplest non-elementary functions in the theory of special functions of applied mathematics and mathematical physics. The three canonical families of HOPs (Hermite, Laguerre and Jacobi) and their generalizations are not only interesting per se, but also because of their physico-mathematical applications in science and technology. In particular, they have been shown to control the physical solutions (wave functions) of the Schrödinger equation of numerous classical and quantum systems, which characterize, for example, the bound stationary states of systems subject to spherically-symmetric quantum-mechanical potentials (harmonic systems of oscillator type, Coulomb systems of hydrogenic type, etc.).

Keeping this in mind, we have examined and reviewed in this work the present knowledge of the spreading of the HOPs along the orthogonality support of their weight function by means of the entropy-like measures of Fisher, Shannon and Rényi types (and their associated ${\mathfrak{L}}_q$-norms)—far beyond the ordinary moments and the standard deviation of their associated Rakhmanov probability density. We have shown the analytical expressions that allow one to compute such spreading measures of the HOPs in a numerical and symbolic manner by means of the degree and the weight function’s parameters of the polynomials. In addition, we have shown the asymptotic behavior ($n\to \infty$) for the ${\mathfrak{L}}_q$-norms of the HOPs $\left\{p_n\left(x\right)\right\}$ and some generalizations, which provides simple and compact expressions for the first dominant terms of these quantities in terms of q and the weight parameters. The latter controls the physical entropies of the highly excited states of atomic systems and the corresponding Rydberg logical gates of quantum computation.

As always, some related open issues deserve to be mentioned, including the following: the calculation of the upper bounds to the Rényi spreading lengths ${\mathcal{L}}_q^R\left[p_n\right]$ for all q and n, extending the corresponding quantities obtained for the Shannon spreading length as mentioned above; the algebraic determination of the Shannon entropies for the HOPs at all n by means of a technique that does not require the evaluation of the logarithmic potential evaluated at the zeros of the polynomials, perhaps at the cost of evaluating some generalized hypergeometric functions at unity (see e.g., [92]); and the asymptotics of the (unweighted) ${\mathfrak{L}}_q$-norms of HOPs when $n\to \infty$ and when the parameters of the weight function become large or very large. In addition, the asymptotics of the Shannon entropy of orthogonal polynomials in the whole Szegö class remains open, although some recent relevant efforts relative to asymptotic upper bounds and universality limits were carried out (see [43,86,89]).