Orthogonal Polynomials with Singularly Perturbed Freud Weights

In this paper, we are concerned with polynomials that are orthogonal with respect to the singularly perturbed Freud weight functions. By using Chen and Ismail’s ladder operator approach, we derive the difference equations and differential-difference equations satisfied by the recurrence coefficients. We also obtain the differential-difference equations and the second-order differential equations for the orthogonal polynomials, with the coefficients all expressed in terms of the recurrence coefficients.


Introduction
Orthogonal polynomials are of great importance in Random Matrix Theory (RMT), integrable systems, numerical analysis, representation theory, etc. It is well known that classical orthogonal polynomials (Hermite, Laguerre and Jacobi) are orthogonal with respect to a weight function w(x) that satisfies the Pearson equation where σ(x) is a polynomial of degree ≤ 2 and τ(x) is a polynomial of degree 1. Semiclassical orthogonal polynomials have a weight w(x) that satisfies the Pearson Equation (1), where σ(x) and τ(x) are polynomials with deg σ(x) > 2 or deg τ(x) = 1 (see, e.g., ([1], Section 1.1.1)). A motivation of this paper is the fact that the recurrence coefficients of semi-classical orthogonal polynomials are usually related to the solutions of the Painlevé equations. For example, Chen and Its [2] proved that the recurrence coefficients of orthogonal polynomials with a singularly perturbed Laguerre weight are expressed in terms of a particular Painlevé III equation. Filipuk, Van Assche and Zhang [3] showed that the recurrence coefficients of a class of semi-classical Laguerre polynomials are related to the Painlevé IV equation. Basor, Chen and Ehrhardt [4] established the relation between the recurrence coefficients of time-dependent Jacobi polynomials and the Painlevé V equation. See [5][6][7][8][9][10][11] and also the recent monograph of Van Assche [1] for more information.
A Freud weight is a weight function of the form ( [12], Section 18.32) where Q(x) is real, even, nonnegative and continuously differentiable. Of special interest are the cases Q(x) = x 2m , m = 1, 2, 3, . . .. In a seminal paper [9], Magnus studied the relations between the Painlevé equations and many semi-classical orthogonal polynomials, in which there are two examples of the one-parameter Freud weights: with t ∈ R a parameter, and Magnus obtained a series of differential and difference equations satisfied by the recurrence coefficients of the corresponding orthogonal polynomials. See also [3,7,13,14].
In this paper, we consider the singularly perturbed Freud weight with t ≥ 0, m = 1, 2, 3, . . . . When t > 0, the factor e − t x 2 induces an infinitely strong zero at the origin. This is a semi-classical weight since it satisfies the Pearson Equation (1) with Orthogonal polynomials with singularly perturbed Gaussian, Laguerre and Jacobi weights have been studied in [2,5,15]. We mention that the weights with an essential singularity at the origin, such as (2), play an important role in many mathematical and physical problems, such as the study of statistics for zeros of the Riemann zeta function [16], the calculation of finite temperature expectation values in integrable quantum field theory [17], the study of the Wigner time-delay distribution [18][19][20], etc.
Note that we do not display the t-dependence of many quantities for simplicity. By the definitions (9) and (10) and with the aid of the three-term recurrence relation (5), it can be shown that A n (x) and B n (x) satisfy the compatibility conditions Multiplying by A n (x) on both sides of (12) and using (11) to eliminate xA n (x), we have the following identity after taking a telescopic sum: Finally, eliminating P n−1 (x) from (7) and (8), we obtain the second-order differential equation satisfied by the orthogonal polynomials: where use has been made of (13) to simplify the coefficient of P n (x). From the point of view of RMT, the weight (2) can be used to define a singularly perturbed Freud unitary ensemble with probability distribution where Z n is the partition function and can be expressed as a Hankel determinant generated by the weight (2) by using Andréief or Heine's identity. See [26][27][28] for more information about this topic. In this respect, there are many problems to be considered, including the large n asymptoics of the partition function and the gap probabilities of the ensemble. We will leave these problems to a future investigation. We mention that Claeys, Krasovsky and Minakov [29] recently studied the asymptotics of the partition function and gap probabilities of a certain Freud random matrix ensemble. In addition, construction of the quadrature formulas related to the weight (2) as in [30] may be an interesting problem.
The main purpose of this paper is to derive the differential and difference equations for the recurrence coefficient β n and also the orthogonal polynomials with respect to the singularly perturbed Freud weight (2) by using Chen and Ismail's method [24,25].

The m = 1 Case
In this section, we consider the simplest case (m = 1) and the weight function now reads It is also called the singularly perturbed Gaussian weight and has been studied by Min, Lyu and Chen [15] (see also [31]). One of the main results in [15] is that the authors establish the relation between the recurrence coefficient and the Painlevé III equation (see Theorem 3 below). The following results in Lemma 1 and Lemma 2 are obtained in [15].

Lemma 1. For this problem, we have
where R n (t) and r n (t) are the auxiliary quantities defined by Substituting (15) and (16) into (11) and (13), we have the following. Lemma 2. The auxiliary quantities R n (t), r n (t) and the recurrence coefficient β n satisfy the following relations: Theorem 1. The recurrence coefficient β n satisfies the nonlinear second-order difference equation Proof. From (18) and (17), we have Substituting (21) and (22) into (19), we establish the theorem.
Theorem 2. The recurrence coefficient β n satisfies the differential-difference equation Proof. Taking a derivative with respect to t in the orthogonality condition Using (6), we have 2tβ n (t) = β n (R n−1 (t) − R n (t)).
Substituting (22) into the above, we obtain the desired result.
Theorem 3. The auxiliary quantity R n (t), related to the recurrence coefficient β n by satisfies a particular Painlevé III equation ( [32], (2)): Proof. This is obtained by eliminating r n (t) from the coupled Riccati equations satisfied by R n (t) and r n (t). See [15] for details.

The m = 2 Case
In this section, we consider the m = 2 case and the weight function is Lemma 3. We have where R n (t) and r n (t) are the auxiliary quantities given by Proof. Substituting (25) into the definition of A n (x) in (9) and using the parity of the integrand, we obtain By the three-term recurrence relation (5), we have y 2 P 2 n (y) = (P n+1 (y) + β n P n−1 (y)) 2 = P 2 n+1 (y) + β 2 n P 2 n−1 (y) + 2β n P n+1 (y)P n−1 (y).
where use has been made of (6). Inserting (29) into (28) gives (26). Similarly, substituting (25) into the definition of B n (x) in (10), we have It is easy to see that Then, we arrive at (27). The proof is complete.
Theorem 6. The recurrence coefficient β n satisfies the differential-difference equation

Proof. Taking a derivative with respect to t in the equality
and then 2tβ n (t) = β n (R n−1 (t) − R n (t)).
Substituting (38) into the above gives the desired result.
In the m = 2 case, we find that the expressions of A n (x) and B n (x) include the terms of the recurrence coefficient β n , which is different from the m = 1 case. This leads to the result that β n and r n (t) do not have a simple relation as in the m = 1 case by using the compatibility conditions. Thus, we cannot derive the coupled Riccati equations for the auxiliary quantities R n (t) and r n (t), and the second-order differential equations for R n (t) and r n (t). Therefore, the relation between our problem and the Painlevé equations is not clear in the m = 2 case. A similar phenomenon arises in the next m = 3 case; see also [33].

The m = 3 Case
In this section, we consider the weight function Lemma 4. We have where R * n (t), R n (t), r * n (t) and r n (t) are the auxiliary quantities given by y 3 P n (y)P n−1 (y)w(y)dy.
Proof. Taking a derivative with respect to t in the equality and then 2tβ n (t) = β n (R n−1 (t) − R n (t)).
Substituting (60) into the above gives the desired result.

Conclusions
In this paper, we study the singularly perturbed Freud weights (2), the associated orthogonal polynomials and the recurrence coefficients in the m = 1, 2, 3 cases. We derive the differential and difference equations for the recurrence coefficients and the orthogonal polynomials in the three cases. It can be seen that for increasing m, the order of the difference equations increases by 2 in each case. We can use our method to consider the higher-order singularly perturbed Freud weights. However, one can imagine that the results for the differential and difference equations will be more and more complicated. Finally, we find that the recurrence coefficient is related to the Painlevé III equation when m = 1, but it is not clear whether there is any connection to the Painlevé equations when m ≥ 2.