Asymptotics for Recurrence Coefficients of X1-Jacobi Polynomials and Christoffel Function
\'A. P. Horv\'ath

TL;DR
This paper derives asymptotic behaviors of recurrence coefficients for X1-Jacobi polynomials and explores their impact on the Christoffel function, addressing challenges posed by their non-standard recurrence relations.
Contribution
It provides new asymptotic formulas for recurrence coefficients of X1-Jacobi polynomials and analyzes their influence on Christoffel functions, overcoming the limitations of classical methods.
Findings
Asymptotic formulas for recurrence coefficients derived
Limit behavior of Christoffel function analyzed
Handling of five-term recurrence relations in exceptional polynomials
Abstract
Computing asymptotics of the recurrence coefficients of X1-Jacobi polynomials we investigate the limit of Christoffel function. We also study the relation between the normalized counting measure based on the zeros of the modified average characteristic polynomial and the Christoffel function in limit. The proofs of corresponding theorems with respect to ordinary orthogonal polynomials are based on the three-term recurrence relation. The main point is that exceptional orthogonal polynomials possess at least five-term formulae and so the Christoffel-Darboux formula also fails. It seems that these difficulties can be handled in combinatorial way.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsRandom Matrices and Applications · Quantum Mechanics and Non-Hermitian Physics · Quantum chaos and dynamical systems
Asymptotics for Recurrence Coefficients of X1-Jacobi Polynomials and Christoffel Function
Á. P. Horváth
Abstract.
Computing asymptotics of the recurrence coefficients of X1-Jacobi polynomials we investigate the limit of Christoffel function. We also study the relation between the normalized counting measure based on the zeros of the modified average characteristic polynomial and the Christoffel function in limit. The proofs of corresponding theorems with respect to ordinary orthogonal polynomials are based on the three-term recurrence relation. The main point is that exceptional orthogonal polynomials possess at least five-term formulae and so the Christoffel-Darboux formula also fails. It seems that these difficulties can be handled in combinatorial way.
Key words and phrases:
exceptional orthogonal polynomials, recurrence relation, Chistoffel function
2010 Mathematics Subject Classification:
33C47,33C45
Supported by the NKFIH-OTKA Grant K128922.
1. Introduction
Given a (non-trivial) measure supported on a set of the complex plane, say with finite moments, by Gram-Schmidt orthogonalization the sequence of orthogonal polynomials can be generated, . The kernel of the projection , where is the subspace of polynomials of dimension is
[TABLE]
There is an extensive literature of the following two questions investigated in parallel, the asymptotic distribution of zeros of the orthogonal polynomials, , where is the normalized counting measure based at the zeros, and the limit of the weighted reciprocal of the Christoffel function that is the measure
[TABLE]
These questions are in close contact with each other and are related to the recurrence relation of orthogonal polynomials.
For instance on the real line if the recurrence coefficients and have limits and respectively, the polynomials have asymptotic zero distribution with density on the interval , see [17] and the references therein.
On the other hand for a supported on the unite circle or on the interval Máté, Nevai and Totik proved that provided the so-called Szegő condition, tends to the equilibrium measures, to the arc- or to the arcsin measure respectively, cf. [21]. The asymptotics of is investigated in rather general circumstances, we mention here a result due to Totik on arcs and curves, see [34]. Finally it can be shown that the zeros of orthogonal polynomials are the eigenvalues of the operator , where is the projection mentioned above and , which ensures that in some sense, see [32]. This observation ensures a direct relation between the two measures derived by ordinary orthogonal polynomials and this technique leads to determinantal point processes.
The aim of the paper is to make small contribution to examination of the problem sketched above in case of exceptional orthogonal polynomials. We focus onto the Christoffel function part.
The notion of exceptional orthogonal polynomials is motivated by problems in quantum mechanics. It turned up in the last twenty years; one of the earliest papers is [7], so the topic is fairly new. Owing to the widespread investigation of exceptional systems its literature became rather rich. We mention here just some examples: [5], [30], [8], etc. and see also the references therein.
The study of Christoffel functions of ordinary orthogonal systems depend on the three-term recurrence formula. Exceptional orthogonal polynomials fulfil -term recurrence relations, where is at least 2. The other important difference is while in ordinary case can be expressed by a (three-term) linear combination of orthogonal polynomials, in exceptional case instead of there is a polynomial of degree at least two. This is why most of the methods developed to ordinary orthogonal polynomials are not applicable to the exceptional systems. This investigation based on combinatorial methods inspired by [11] and [12]. Below we deal with Jacobi polynomials in general (generated by one Darboux transform from Jacobi polynomials) and specially Jacobi polynomials of codimension 1.
The paper is organized as it follows. In the next section we define exceptional orthogonal polynomials and list up some properties which are important for further investigation. In Section 3 we deal with recurrence relation and give the asymptotics of the recurrence coefficients. In the fourth section we study the asymptotic behavior of the weighted reciprocal of the Christoffel function for codimension 1 exceptional Jacobi case, and in the last section we introduce the modified average characteristic polynomial and study the connection between its zero distribution and .
2. Exceptional orthogonal polynomials
In this section we summarize some properties of exceptional orthogonal polynomials which are important hereinafter. Finally we deal with -Jacobi polynomials, but our investigation is based on the general characterization of exceptional systems rather than the known examples. The definition is as it follows, cf. [5, Definition 7.4]
Exceptional polynomials means a co-finite real-valued polynomial sequence , where from the sequence of degrees finite many are missing provided
(1) The polynomials are eigenfunctions of a differential operator of second order with rational coefficients.
(2) There is an interval and a positive weight function on with finite moments and at the endpoints of , where is the coefficient of the second derivative in the differential operator.
(3) The vector space spanned by the elements of the sequence is dense in the weighted space .
Second property ensures orthogonality of the exceptional polynomials on with respect to . Exceptional orthogonal polynomial systems (XPS) are characterized by their construction cf. [5, Theorem 1.2]: all XPS can be obtained by applying finite sequence of Darboux transformations to a classical orthogonal polynomial system.
Below we need the construction thus we summarize the -step Darboux transformation case in short.
Classical orthogonal polynomials are eigenfunctions of the second order linear differential operator
[TABLE]
with eigenvalues . can be decomposed as
[TABLE]
where , are rational functions and
[TABLE]
Then the exceptional polynomials are the eigenfunctions of that is the partner operator of which is
[TABLE]
where
[TABLE]
and fulfils the Ricatti equation
[TABLE]
cf. [5, Propositions 3.5 and 3.6].
[TABLE]
so exceptional polynomials can be obtained from the classical ones by application of (finite many) appropriate first order differential operator(s) to the classical polynomials. This observation motivates the notation below
[TABLE]
(and recursively in case of Darboux transforms.) The degree of is usually greater than , so subsequently is an index not a degree any more; it shows that is generated from . Since this notation does not follows where the gaps are in the sequence of degrees, it proved to be useful in handling recurrence formulae, cf. [15].
is an orthogonal system on with respect to the weight , where is one of the classical weights. Since at the endpoints of (if there is any) may have a zero, can be zero as well here, but it has no zeros inside , otherwise would not have finite moments. Canceling the zeros at the endpoints we introduce the following notation.
Notation.
Let , such that and let .
Remark.
(R1) As -s are polynomials for all , applying operator to it can be seen that must be a polynomial and to shows that itself is a polynomial too. Thus we can write as , where is also polynomial.
(R2) According to (4) satisfies the following differential equation:
[TABLE]
[TABLE]
Since all the other terms are polynomials is also polynomial for all , thus the coefficient of is polynomial. Furthermore denoting by
[TABLE]
(8) ensures that is divisible by for all .
(R3) Observe that if then .
3. Recurrence relations
In this section after introducing recurrence relation with respect to exceptional orthogonal polynomials, we investigate the asymptotics of the recurrence coefficients.
Ordinary orthogonal polynomial systems contain polynomials of all degrees which fulfil the three-term recurrence relation
[TABLE]
Recurrence formulae can be derived for exceptional orthogonal polynomials as well, but these are -term formulae, where . As (10) can be rearranged in two ways, for exceptional polynomials there are two different kinds of recurrence relations: with constant- and with variable dependent coefficients. Recalling the Darboux decomposition of the differential operator which generates exceptional polynomials, formulae with dependent coefficients can be derived by operators (cf. (7)) (see e.g. [31], [6], [15], [24], [25], [26]) and with constant coefficient ones by operators (see e.g. [25], [26], [9], [27], [3], [22], [2]). The length of the formula, in the variable dependent coefficient case depends on the number of Darboux transforms, that is
[TABLE]
In contrast in the constant coefficient case the length depends on the codimension independently of how many Darboux transforms resulted the codimension of the space of exceptional polynomials. In variable dependent coefficient case it was more convenient to work with monic polynomials, but in constant coefficient case we use the orthonormal form of the polynomials.
Subsequently we deal with exceptional polynomials generated by one-step Darboux transform, so for sake of simplicity further properties are written in this case. We denote by the orthonormal classical () or exceptional () polynomials. Below we summarize the results of Odake, G omez-Ullate, Kasman, Kuijlaars, and Milson (see [25], [26], [9]). Because we need some steps of the proof too, below the sketch of the proofs are also given in brief and for sake of simplicity in one-step Darboux transform case.
Let us introduce the linear space
[TABLE]
and the stabilizer ring
[TABLE]
where denotes the set of polynomials.
Theorem A. With the notation above Let . If is divisible by , then .
Proof.
By (R2) and by comparison of the codimensions of the two sets in question, it can be seen that a polynomial belongs to if and only if is divisible by . According to (R3) . Thus using (R2) if is divisible by then is also divisible by and so it belongs to . By orthogonality properties it can be expand as it is given below
[TABLE]
Corollary B. If
[TABLE]
then
[TABLE]
Proof.
[TABLE]
[TABLE]
That is
[TABLE]
[TABLE]
which ensures that for all .
3.1. Asymptotics of recurrence coefficients
The purpose of this section is determination of the limits of recurrence coefficients i.e. , cf. (13). To this end we return to the expression in Corollary B.
Subsequently we use the simplest element of :
[TABLE]
where . Applying (R3), (2) and (1)
[TABLE]
[TABLE]
Let
[TABLE]
In Theorem A the assumption refers to thus we can take with an arbitrary constant term. Here the constant is chosen to be zero.
Classical orthonormal polynomials satisfy the following relation (cf. e.g. [33]).
[TABLE]
where , and are constants and is the coefficient of the second derivative in the differential equation. Thus
[TABLE]
[TABLE]
[TABLE]
That is knowing the limits of the recurrence coefficients of classical orthonormal polynomials similar limits can be derived to the corresponding exceptional ones.
Wide class of ordinary orthogonal polynomials supported on have recurrence coefficients with limit ; , see [23, Theorem 4.5.7]. In classical cases, the situation is easier, after a normalization if it is necessary the limit of the recurrence coefficients; ; can be got by direct computations. That is (10) can be written as
[TABLE]
where , tend to zero. Taking the product of and , where the error term can be estimated by , say, where is a constant depending on . The symmetry of (10) ensures the same symmetry of the iterated recurrence relation. It has the following form:
[TABLE]
where that is it tends to zero when tends to infinity. Substituting it into (16) we have
[TABLE]
[TABLE]
[TABLE]
[TABLE]
In view of (19) and Corollary B the coefficients in the expansion of are
[TABLE]
[TABLE]
and if
[TABLE]
[TABLE]
[TABLE]
[TABLE]
To finish this general part we compute .
Lemma 1**.**
In (18) the iterated coefficients in limit are
[TABLE]
**Remark. **Certainly similar results can be found in the literature, see e.g. [1] with certain Freud weights. We introduce here the method we use mainly in the next section to get the result in the required form.
Proof.
Recalling that is the main part of the coefficient of in the expansion of we can compute it on an oriented weighted graph , where the vertices and the edges
[TABLE]
The weight on an edge is given by
[TABLE]
We start from zero level, say and we have to reach level (or ) in steps. As we can step forward (on the same level) with weight , up and down with weight , on the -long path we can take steps up or down which contains steps upwards and downwards, and the remainder ones forwards. (For we get the same.) That is
[TABLE]
which is just the statement.
In the rest of this section we compute the asymptotic recurrence coefficients .
3.2. X1 Jacobi polynomials
Let the Jacobi polynomial defined by Rodrigues’ formula:
[TABLE]
where .
[TABLE]
the orthonormal Jacobi polynomials which fulfil the following differential equation (cf. [33, (4.2.1)])
[TABLE]
With
[TABLE]
so taking into consideration that , etc. (for see [33, (4.3.4)]), according to [33, (4.5.5)]
[TABLE]
, etc. - for sake of simplicity we omit and and means that is between two positive constants. Again after normalization in view of [33, (4.5.1)]
[TABLE]
Applying Lemma 1 with the values of (26) we immediately have
Lemma 2**.**
[TABLE]
Now we are in position to compute the limit coefficients. Since does not depend on , by (23) . Taking into consideration this and substituting (26), (25) and Lemma 2 into (20) and (21) we arrive to
Proposition 1**.**
If is even or odd respectively
[TABLE]
4. Jacobi polynomials of codimension 1
For ordinary orthogonal polynomials it is proved in rather general circumstances that tends to the equilibrium measure of the (essential) support of in some sense, where is the Radon-Nikodym derivative of and the Christoffel-Darboux kernel, see e.g. [32] and the references therein. All he proofs of these types of theorems use the three-term recurrence relation of orthogonal polynomials. In this section we show in the simplest, one codimensional case that is with five-term recurrence formula, that the statement above still fulfils.
With the notation above let us choose
[TABLE]
Here we choose the constant term of to [math] again. To define monic polynomials are also appropriate. Since in the known examples is not monic (cf. e.g. [10]) we take it in the form given above and it does not cause any further difficulty. Let
[TABLE]
where are the exceptional Jacobi polynomials orthonormal with respect to on . We introduce the following notation:
[TABLE]
[TABLE]
The equilibrium measure of , is defined by
[TABLE]
Let us recall the Wallis’ formula:
[TABLE]
Theorem 1**.**
If is defined as in (28) by -Jacobi polynomials of codimension 1, on
[TABLE]
in weak-star sense.
To prove Theorem 1 we need the following Lemma.
Lemma 3**.**
For all
[TABLE]
Proof.
First we compute the effect of to , . Binomial theorem and Wallis’ formula give
[TABLE]
[TABLE]
On the other hand, according to Theorem A for all (fixed) there exist coefficients such that
[TABLE]
thus
[TABLE]
where the limit exists with regard to Proposition 1. We can proceed as previously, that is recalling the oriented weighted graph with , we define edges as
[TABLE]
and weights on the edges as
[TABLE]
where, according to Proposition 1
[TABLE]
Thus can be expressed as the collection of all possible weighted -long paths from to ,
[TABLE]
That is in each step we can move forwards, denoted by (with weight ); one level up or down denoted by (with weight ) and two levels up or down denoted by (with weight ). Since finally we have to arrive to the starting level it is clear that we need even number of . Let the number of is . Then we can choose pieces of (that is two levels big steps upwards) and . The rest of the steps are zeros (i.e. forwards). That is
[TABLE]
[TABLE]
Indeed we can choose the -s in ways, and to neutralize the pieces of we need and we have further , that is the number of the choices of -s is . Then we choose from the remainder elements the elements with absolute value , and from these the pieces of . If , we can do the same with -s, etc. and it gives a factor in these cases. Now taking into consideration the weights we arrive to the formula above.
Notice that (32) can be rearranged as
[TABLE]
where
[TABLE]
In order to prove that , it is enough to show that for all ,
[TABLE]
First we observe that if and is arbitrary, then
[TABLE]
[TABLE]
If is even and , then on both of the sides of (33) stands .
To show (33) we use induction from to . Notice that replacing instead of the right-hand side of (33) doubles up. We show the same for the left-hand side of (33).
We can consider as all the weighted paths from zero level to zero level in steps with fixed -s and with the following new weights: have weight and , have weight . In we have one more step and two more -s. Thus from a paths of we can get a path of if we replace a by two -s (or a by two -s), or a by and and there are no other ways. Observe that all these steps double up the weighted sum, because if we replace a by -s we replace by and when we replace a , we can change the places of and , so we get two different paths. Proceeding this way and then omitting the paths we have got more than one times, we get all the paths of . Conversely, taking two paths from such that one can be derived from the other by a change of a pair of and replacing them by a , or two -s by a , etc., omitting the repetitions again we get all the paths of , thus and together with (34) it gives the statement.
Proof.
(of Theorem 1) By arithmetic mean law Lemma 3 ensures for all
[TABLE]
The linear space
[TABLE]
is also an algebra with unit element. Recalling that has finite moments, and is in the denominator of , the zero of is out of the interval of orthogonality, that is has constant sign here and so is monotone on which ensures that separates the points of . So by Stone-Weierstrass theorem is dense in . As for all and are probability measures, the theorem follows from density.
5. ”Average characteristic polynomial”
In this section we look at (see (28)) from another perspective. It shows up in connection with determinantal point processes. Without going into this theory deep, we sketch how it relates to average characteristic polynomial, .
In ordinary case this correspondence ensures results about limit zero distribution of orthogonal polynomials. Although the known examples show that normalized counting measures based on the zeros of exceptional polynomials tend to the equilibrium measure too (see e.g. [10], [13], [18], [14], [19]), the situation in general case is a little bit more complicated.
To get relationship between and similar to the ordinary case, we have to modify the average characteristic polynomial - we denote the modified version with again.
Let be an ordinary or exceptional orthonormal system of polynomials on a real interval with respect to an appropriate measure or weight function (which is the Radon-Nikodym derivative of ), that is orthonormal in and . Let
[TABLE]
and as above
[TABLE]
We define a polynomial and such that
[TABLE]
with some fixed depending only on . We also define two mappings acting on
[TABLE]
Note that in ordinary case , as in [11] and [12].
are random variables, the joint probability distribution on is
[TABLE]
where is a normalization factor. The expectation refers to .
Let be real valued, bounded and measurable function with bounded support contained by and assume that
[TABLE]
Then
[TABLE]
[TABLE]
see e.g. [16, Proposition 2.2].
The empirical distribution
[TABLE]
Specially if (38) becomes
[TABLE]
We also need and the modified average characteristic polynomial
[TABLE]
Denote by the zeros of and define
[TABLE]
After introducing notation above we are in position to state the main result of this section.
Theorem 2**.**
If there is a such that for all , then for all
[TABLE]
Let , , , that is . Proposition 1 and Theorem 2 immediately ensures the next corollary.
**Corollary. ** In the Jacobi case for all
[TABLE]
Lemma 4**.**
[TABLE]
[TABLE]
Proof.
To apply (38) first we have to check the conditions assumed above in our case. We take , which has the required properties. Orthogonality together with Dyson’s theorem (see [20, Theorem 5.14]) ensures the convergence in (37).
[TABLE]
[TABLE]
[TABLE]
[TABLE]
Considering the Fredholm’s expansion of integral operator acting on with kernel we have
[TABLE]
cf. [4], from which follows the next equality
[TABLE]
That is the roots of are the eigenvalues of , so .
Let
[TABLE]
The proof of the next lemma is independent of the form of . The weight of a path, .
Lemma 5**.**
[11, Lemma 2.5, Lemma 2.6]** For all
[TABLE]
[TABLE]
To make the discussion self-contained we insert the proof of Theorem 2 in brief, cf. [12].
Proof.
(of Theorem) In view of Lemma 4 and 5 we estimate
[TABLE]
According to (35), if then . By the assumption on the weight of an -long path is at most and the vertices of admissible paths are contained by the set which has at most elements. That is
[TABLE]
**Remark. **By Theorem 1 and Theorem 2 with an appropriately defined in Jacobi case of codimension 1 the following relation can be derived
[TABLE]
Indeed, Theorem 2 ensures that there is a compact interval which contains eigenvalues for all . Define such that . Thus the main part of is supported on . Applying the binomial theorem to both of the sides of Lemma 3 we get the same result with as well. So the density argument mentioned above leads to the weak star convergence of the restriction of .
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] S. S. Bonan, D. S. Clark, Estimates of the Hermite and the Freud polynomials, J. of Approximation Theory , 63 (1990), 210-224.
- 2[2] P. A. Clarkson, D. Gómez-Ullate, Y. Grandati, R. Milson, Rational solutions of higher order Painlevé systems I., (2018) ar Xiv:1811.09274.
- 3[3] A. J. Durán, Higher order recurrence relation for exceptional Charlier, Meixner, Hermite and Laguerre orthogonal polynomials, Integral Transforms and Special Functions , 26 (2015), 357 376.
- 4[4] I. Gohberg, S. Goldberg and N. Krupnik, Traces and Determinants of Linear Operators Birkh user, Basel, 2000.
- 5[5] M. Á. García-Ferrero, D. Gómez-Ullate, R. Milson, A Bochner type classification theorem for exceptional orthogonal polynomials, J. Math. Analysis and Appl. , 472 (1), (2019) 584-626.
- 6[6] D. G omez-Ullate, Y. Grandati and R. Milson, Rational extensions of the quantum Harmonic oscillator and exceptional Hermite polynomials, J. Phys. A , 47 (2014), 015203.
- 7[7] D. Gómez-Ullate, N. Kamran, R. Milson, An extended class of orthogonal polynomials defined by a Sturm-Liouville problem, J. Math. Analysis and Appl. , 359 (2009), 352-367.
- 8[8] D. Gómez-Ullate, N. Kamran, R. Milson, On Orthogonal Polynomials Spanning a Non-standard Flag, Algebraic Aspects of Darboux Transformations, Quantum Integrable Systems and Supersymmetric Quantum Mechanics, Contemp. Math. , 563 (2012), 51-72.
