Observe that a naive evaluation of these functionals by means of quadratures encounters the difficulty of the zeros of $p_n$  , that belong to the interval of orthogonality. Some of the contributions in this sense are [5, 7, 8, 9, 22] . A second problem is the study of the asymptotic behavior of $\left\{E_n\right\}$  , $\left\{F_n\right\}$  , and $\left\{G_n\right\}$  when $n\to \infty$  , which has a special interest in the analysis of the highly-excited (Rydberg) states of numerous quantum-mechanical systems [22] . In this sense there have been important contributions in the last few years [2, 6, 7, 10, 11, 14, 20] . In a recent paper [3] the authors have studied the asymptotic behavior of these functionals under the assumption that the weight of orthogonality satisfies the Szego condition,

$$\begin{array}{c}{\int }_{-1}^1\frac{log\left(w\left(x\right)\right)}{\sqrt{1-x^2}}dx>-\infty .\end{array}$$

(4)

Under an additional assumption on the growth of the polynomials on the interval of orthogonality they proved that both $E_n$  and $F_n$  (and in consequence, also $G_n$  ) converge, and

$$\begin{array}{c}{lim}_n{F}_n\left(w\right)=log\left(\pi \right)-1\end{array}$$

(5)

(notice that $F_n$  is taken here with a slightly different normalization than in [3] ). The authors of [3] conjectured that the limit in  5 is valid for a larger class of weights; from their work it follows also that if $w>0$  on $(-1,1)$  does not satisfy  4 , then $E_n\left(w\right)$  and $G_n\left(w\right)$  diverge to $-\infty$  . The Pollaczek polynomials constitute the first and the best known example of a family of orthogonal polynomials on $[-1,1]$  with respect to a weight not satisfying the Szego condition  4 . In this paper we deal with the symmetric Pollaczek polynomials, $p_n^{\lambda }(x;a)$  , that depend on two real parameters, $\lambda >0$  , $a\ge 0$  , and that may be defined by the recurrence relation

$$\begin{array}{c}x{p}_n^{\lambda }(x;a)=a_{n+1}{p}_{n+1}^{\lambda }(x;a)+a_n{p}_{n-1}^{\lambda }(x;a),p_{-1}^{\lambda }(x;a)=0,p_0^{\lambda }(x;a)=1,\end{array}$$

(6)

with the coefficients

$$\begin{array}{c}{a}_n=\frac{1}{2}\sqrt{\frac{n(n+2\lambda -1)}{(n+\lambda +a)(n+\lambda +a-1)}}.\end{array}$$

(7)

It is known (see [21,Appendix] ) that these polynomials are orthonormal on $[-1,1]$  with respect to the unitary weight function

$$\begin{array}{c}{w}_{\lambda }(x;a)=\frac{2^{2\lambda }(\lambda +a)}{2\pi \Gamma \left(2\lambda \right)}(1-x^2{)}^{\lambda -1/2}{e}^{(2arccosx-\pi )\frac{ax}{\sqrt{1-x^2}}}{\left|\Gamma \left(\lambda +i\frac{ax}{\sqrt{1-x^2}}\right)\right|}^2,\end{array}$$

(8)

where $\Gamma \left(x\right)$  denotes the gamma function. From  8 it is clear that Pollaczek polynomials $p_n^{\lambda }(x;0)$  (that is, for $a=0$  ) reduce to orthonormal Gegenbauer polynomials with parameter $\lambda$  . In the sequel, whenever it cannot lead us into confusion, we omit the explicit reference to the parameters $\lambda$  and $a$  from the notation of the polynomials.

Our main goal is to study the asymptotic behavior of the sequences $E_n\left(w\right)$  , $F_n\left(w\right)$  and $G_n\left(w\right)$  as $n\to \infty$  , when $w$  is the symmetric Pollaczek weight. We can summarize our main results saying that limit  5 is proved to be valid also for $w=w_{\lambda }(\cdot ;a)$  , with the restriction $\lambda \ge 1$  (so, this fact extends beyond the Szego class, as conjectured), and we find the main part of the asymptotic expansion of $E_n\left(w\right)$  and $G_n\left(w\right)$  , up to the $o\left(1\right)$  terms.

Namely, we establish the following: for $F_n$  we prove that indeed, also for the symmetric Pollaczek weight, limit  5 is still valid:

Theorem 1.1For the symmetric Pollaczek weight $w=w_{\lambda }(\cdot ;a)$  , with $a\ge 0$  and $\lambda \ge 1$  , $$F_n\left(w\right)=log\left(\pi \right)-1+o\left(1\right),n\to \infty .$$ 

Remark 1.2The restriction $\lambda \ge 1$  comes from the method of proof; we believe that Theorem  1.1 is valid for the whole range of the parameter $\lambda$  , that is, for $\lambda >0$  .

For the divergent sequence $\left\{G_n\right\}$  we find the first two terms of its asymptotic expansion:

Theorem 1.3For the symmetric Pollaczek weight $w=w_{\lambda }(\cdot ;a)$  , with $a\ge 0$  and $\lambda >0$  ,

$$\begin{array}{cc}{G}_n\left(w\right)& =-2alog\left(n\right)+2a+log\left(\frac{\Gamma (\lambda +a)\Gamma (\lambda +a+1)}{\pi \Gamma \left(2\lambda \right)}\right)+o\left(1\right),n\to \infty .\end{array}$$

As a straightforward corollary we obtain

Corollary 1.4For the symmetric Pollaczek weight $w=w_{\lambda }(\cdot ;a)$  , with $a\ge 0$  and $\lambda \ge 1$  ,

$$\begin{array}{c}{E}_n\left(w\right)=-2alog\left(n\right)+\tau (\lambda ;a)+o\left(1\right),n\to \infty ,\end{array}$$

(9)

where

$$\begin{array}{c}\tau (\lambda ;a)\stackrel{\text{def}}{=}2a-1+log\left(\frac{\Gamma (\lambda +a)\Gamma (\lambda +a+1)}{\Gamma \left(2\lambda \right)}\right).\end{array}$$

(10)

Remark 1.5The value $\tau (\lambda ;0)=-1+log\left(\Gamma \left(\lambda \right)\Gamma (\lambda +1)/\Gamma \left(2\lambda \right)\right)$  matches ${lim}_n{E}_n$  for orthonormal Gegenbauer polynomials, found in [2] .

For illustration, we have computed the entropy $E_n\left(w_{\lambda }\right(\cdot ;a\left)\right)$  for $n=1,2,...,500$  , $\lambda =5,15$  , and $a=5,10,15$  (Fig.  1 ), using the numerical algorithm from [8] , which admits as the only input data the expression of the recurrence coefficients $a_n$  in  7 . For comparison, in Fig.  2 we plot the difference $E_n\left(w_{\lambda }\right(\cdot ;a\left)\right)-{\tilde{E}}_n\left(w_{\lambda }\right(\cdot ;a\left)\right)$  , where $${\tilde{E}}_n\left(w_{\lambda }\right(\cdot ;a\left)\right)\stackrel{\text{def}}{=}-2alog\left(n\right)+\tau (\lambda ;a).$$ 

Figure 1 : Entropy $E_n\left(w_{\lambda }\right(\cdot ;a\left)\right)$  for $n=1,2,...,500$  , with $\lambda =5$  (upper) and $\lambda =15$  (lower). We use the values $a=5$  (‘*’), $a=10$  (‘x’) and $a=15$  (‘o’).

Figure 2 : $E_n\left(w_{\lambda }\right(\cdot ;a\left)\right)-{\tilde{E}}_n\left(w_{\lambda }\right(\cdot ;a\left)\right)$  for $n=1,2,...,500$  , with $\lambda =5$  (upper) and $\lambda =15$  (lower). We use the values $a=5$  (‘*’), $a=10$  (‘x’) and $a=15$  (‘o’).

In [14] , a rather general result about the leading term of the asymptotics of $E_n$  has been established, that we state here for the symmetric case:

for even weights functions $w$  on $[-1,1]$  which belong to the class $\mathcal{ℱ}(C^2+)$  introduced in [13] (and whose definition we recall in Section  2 ), if an additional assumption on the behavior of $w$  at $\pm 1$  (see Eq. (1) in [14] ) is satisfied, then

$$\begin{array}{c}{E}_n\left(w\right)=-\frac{2}{\pi }{\int }_{-{\alpha }_n}^{{\alpha }_n}\frac{Q\left(x\right)}{\sqrt{{\alpha }_n^2-x^2}}dx\left(1+o\left(1\right)\right),n\to \infty ,\end{array}$$

(11)

where

$$\begin{array}{c}Q\left(x\right)\stackrel{\text{def}}{=}-\frac{1}{2}log\frac{w\left(x\right)}{w\left(0\right)}\end{array}$$

(12)

is the “external field”, associated with the weight $w$  , and ${\alpha }_n$  is the Mahskar-Rakhmanov-Saff number (or MRS number), defined as the unique solution of the integral equation

$$\begin{array}{c}\frac{1}{\pi }{\int }_{-{\alpha }_n}^{{\alpha }_n}\frac{x{Q}^{\prime }\left(x\right)}{\sqrt{{\alpha }_n^2-x^2}}dx=n\end{array}$$

(13)

(see e.g. [19] for details). Unfortunately, the assumptions from [14] on the behavior of $w$  at $\pm 1$  are not fulfilled by the symmetric Pollaczek weights.

Nevertheless, the result of Corollary  1.4 above shows that the assertion in [14] is still valid:

Corollary 1.6For the symmetric Pollaczek weight $w=w_{\lambda }(\cdot ;a)$  , with $a\ge 0$  and $\lambda \ge 1$  , formula  11 holds.

Remark 1.7We should observe however that the result in [14] is sharp: there the $o\left(1\right)$  term of  11 has a power decay, fact which is not true for the Pollaczek weight, where the decay is logarithmic.

Finally, asymptotic formula  9 has a curious interpretation in terms of the behavior of the mutual energy of two relevant sequences of probability measures on $[-1,1]$  associated with $p_n$  ’s:

$$\begin{array}{c}{\rho }_n=\frac{1}{n}{\sum }_{j=1}^n{\delta }_{{\zeta }_j^{\left(n\right)}}\text{and}d{\nu }_n\left(x\right)=p_n^2\left(x\right)w\left(x\right)dx,\end{array}$$

(14)

where $-1<{\zeta }_1^{\left(n\right)}<\cdot \cdot \cdot <{\zeta }_n^{\left(n\right)}<1$  are the zeros of the polynomial $p_n$  . Both measures are standard objects of study in the analytic theory of orthogonal polynomials. For instance, the normalized zero counting measure ${\rho }_n$  is closely connected with the $n$  -th root asymptotics of $p_n$  , while ${\nu }_n$  is associated with the behavior of the ratio $p_{n+1}/p_n$  as $n\to \infty$  (see [16, 18] ).

If $\rho$  and $\nu$  are Borel (generally speaking, real signed) measures on $\mathbb{C}$  , we denote by

$$\begin{array}{c}{V}^{\rho }\left(z\right)\stackrel{\text{def}}{=}-\int log|z-t|d\rho \left(t\right)\end{array}$$

(15)

the logarithmic potential of $\rho$  , and by $$I[\nu ,\rho ]\stackrel{\text{def}}{=}\int V^{\nu }\left(z\right)d\rho \left(z\right)=-\int \int log|z-t|d\nu \left(t\right)d\rho \left(z\right),$$  the mutual energy of $\nu$  and $\rho$  . The latter is connected with the entropy  1 by the formula

$$\begin{array}{c}{E}_n=-2log{\gamma }_n+2{\sum }_{j=1}^n{V}^{{\nu }_n}\left({\zeta }_j^{\left(n\right)}\right)=-2log{\gamma }_n+2nI[{\rho }_n,{\nu }_n],\end{array}$$

(16)

where

$$\begin{array}{c}{\gamma }_n=2^n{\left(\frac{(\lambda +a+1{)}_n(\lambda +a{)}_n}{n!(2\lambda {)}_n}\right)}^{1/2}>0\end{array}$$

(17)

is the leading coefficient of $p_n$  , and $(z{)}_n=\Gamma (z+n)/\Gamma (z)$  denotes as usual the Pochhammer’s symbol. It is well known that as long as the orthogonality weight $w>0$  a.e. on $[-1,1]$  , both ${\rho }_n$  and ${\nu }_n$  tend (as $n\to \infty$  ) in the weak-* sense to the Chebyshev (equilibrium) distribution of the interval, which implies that ${lim}_{n\to \infty }I[{\rho }_n,{\nu }_n]=log\left(2\right)$  . In [3] a rather surprising “universal” behavior of the next term of the asymptotic expansion of $I[{\rho }_n,{\nu }_n]$  was observed. Namely, if the orthogonality weight satisfies the Szego condition  4 with an additional assumption on the growth of the sequence of $\left\{p_n\right\}$  on $[-1,1]$  , then $I[{\rho }_n,{\nu }_n]=log\left(2\right)-1/\left(2n\right)+o(1/n)$  , $n\to \infty$  . However, this result is no longer valid for the Pollaczek polynomials, as it follows from formula  9 :

Corollary 1.8For the symmetric Pollaczek weight $w=w_{\lambda }(\cdot ;a)$  , with $a\ge 0$  and $\lambda \ge 1$  , $$I[{\rho }_n,{\nu }_n]=log\left(2\right)-\frac{1-2a}{2n}+o\left(\frac{1}{n}\right),n\to \infty .$$ 

Observe that the second term of asymptotics is still independent of the main parameter $\lambda$  , and matches the result in [3] for $a=0$  .

The structure of this paper is as follows. In Section  2 we gather some technical facts about the weight function $w_{\lambda }(\cdot ;a)$  ; in particular, we show that this weight does not satisfy  4 , and for $\lambda \ge 1$  , it belongs to the class $\mathcal{ℱ}(C^2+)$  . Section  3 contains some results about the equilibrium measure of total mass $n$  in the external field $Q$  ; it is needed for the proof of Theorem  1.1 and Corollary  1.4 (Section  4 ). We defer the proof of the asymptotics of the sequence $\left\{G_n\right\}$  (Theorem  1.3 ) to Section  5 . Finally, corollaries  1.6 and  1.8 are established in Section  6 .

2 The weight function

As a first step in our analysis we study the behavior of the symmetric Pollaczek weight function $w_{\lambda }(\cdot ;a)$  defined in  8 , which we denote simply by $w$  whenever it cannot lead us into confusion. Using the notation $$t\stackrel{\text{def}}{=}\frac{ax}{\sqrt{1-x^2}},x\in \Delta \stackrel{\text{def}}{=}[-1,1],$$  we rewrite its definition as

$$\begin{array}{c}w\left(x\right)=\frac{2^{2\lambda }(\lambda +a)}{2\pi \Gamma \left(2\lambda \right)}(1-x^2{)}^{\lambda -1/2}{e}^{(2arccosx-\pi )t}{\left|\Gamma (\lambda +it)\right|}^2,x\in \Delta .\end{array}$$

(18)

This is an even function on $\Delta$  , strictly positive in $(-1,1)$  , but vanishing at the end points. The fast (exponential) decay at $\pm 1$  is precisely the reason why $w$  does not satisfy  4 . Indeed, using the asymptotic formula (6.1.40) from [1] (see also [15,§2.11] ), we can easily obtain that $$2log\left|\Gamma \left(\lambda +ti\right)\right|=-\pi t+2(\lambda -1/2)logt+2log\sqrt{2\pi }+\frac{\lambda (\lambda -1)(2\lambda -1)}{6t^2}+\mathcal{O}\left(t^{-4}\right),$$  when $x\to 1^{-}$  (which denotes in what follows the one-sided limit from the left). Thus, ( 18 ) yields that

$$\begin{array}{cc}logw\left(x\right)& =log\left(\frac{2^{2\lambda }(\lambda +a)}{2\pi \Gamma \left(2\lambda \right)}\right)+(\lambda -1/2)log(1-x^2)+(2arccosx-\pi )t\end{array}$$

$$\begin{array}{cc}& +log{\left|\Gamma (\lambda +it)\right|}^2\end{array}$$

$$\begin{array}{cc}& =-2\pi t+(\lambda -1/2)log\left(x^2\right)+2tarccosx+log\left(\frac{2^{2\lambda }(\lambda +a)a^{2\lambda -1}}{2\pi \Gamma \left(2\lambda \right)}\right)\end{array}$$

$$\begin{array}{cc}& +2log\sqrt{2\pi }+\frac{\lambda (\lambda -1)(2\lambda -1)}{6t^2}+\mathcal{O}\left(t^{-4}\right),x\to 1^{-}.\end{array}$$

(19)

Since

$$\begin{array}{c}2tarccosx=2a-\frac{2a}{3}(1-x^2)+\mathcal{O}(1-x^2{)}^2,x\to 1^{-},\end{array}$$

(20)

we obtain that $$logw\left(x\right)=-2\pi t+\mathcal{O}\left(1\right),x\to 1^{-}.$$  In consequence,

$$\begin{array}{c}w\left(x\right)=exp(-2\pi |t|+\mathcal{O}(1\left)\right)=exp\left(\frac{-2\pi a\left|x\right|}{\sqrt{1-x^2}}+\mathcal{O}\left(1\right)\right),\left|x\right|\to 1^{-},\end{array}$$

(21)

showing that for this weight the integral in  4 is divergent. The previous analysis motivates the introduction of functions $w_0$  and $s$  on $\Delta$  , such that

$$\begin{array}{c}{w}_0\left(x\right)\stackrel{\text{def}}{=}{e}^{-2\pi \left|t\right|}=exp\left(-\frac{2\pi a\left|x\right|}{\sqrt{1-x^2}}\right),\text{and}w\left(x\right)=w_0\left(x\right)e^{s\left(x\right)}.\end{array}$$

(22)

Lemma 2.1Function $s\in C^{\infty }(-1,1)$  is even and continuous in $[-1,1]$  .

Proof. ] We need to check only the existence of finite limits of $s$  at $\pm 1$  (the rest is trivial). From ( 19 ) it is clear that

$$\begin{array}{cc}s\left(x\right)=log\frac{w\left(x\right)}{w_0\left(x\right)}=& log\left(\frac{2^{2\lambda }(\lambda +a)}{2\pi \Gamma \left(2\lambda \right)}{a}^{2\lambda -1}\right)+(\lambda -1/2)log\left(x^2\right)\end{array}$$

$$\begin{array}{cc}& +2tarccosx+\frac{\lambda (\lambda -1)(2\lambda -1)}{6t^2}+\mathcal{O}\left(t^{-4}\right),\end{array}$$

and using ( 20 ) we get

$$\begin{array}{cc}s\left(x\right)=log\left(\frac{2^{2\lambda }(\lambda +a)}{\Gamma \left(2\lambda \right)}{a}^{2\lambda -1}{e}^{2a}\right)+\mathcal{O}(x-1),x\to 1^{-},& \end{array}$$

which concludes the proof.

Proposition 2.2For $\lambda \ge 1$  , the weight $w=w_{\lambda }(\cdot ;a)$  belongs to the class $\mathcal{ℱ}(C^2+)$  .

Recall that $w\in \mathcal{ℱ}(C^2+)$  (see [13] ) if the corresponding external field $Q$  introduced in  12 is positive and verifies the following conditions:

• a) $Q^{\prime }$  is continuous in $\Delta$  .
• b) ${Q^{\prime }}^{\prime }$  exists and is positive in $\Delta \backslash \left\{0\right\}$  .
• c) ${lim}_{\left|x\right|\to 1^{-}}Q\left(x\right)=\infty .$ 
• d) Function $$T\left(x\right)\stackrel{\text{def}}{=}\frac{x{Q}^{\prime }\left(x\right)}{Q\left(x\right)},$$  is cuasi-increasing in $(0,1)$  cuasi-decreasing in $(-1,0)$  , and

  $$\begin{array}{c}T\left(x\right)\ge \Lambda >1,x\in (-1,1).\end{array}$$

  (23)

• e) There exists $C_1>0$  such that $$\frac{{Q^{\prime }}^{\prime }\left(x\right)Q\left(x\right)}{\left(Q^{\prime }\right(x){)}^2}\le C_1,x\in (-1,1).$$ 
• f ) There exist a compact subinterval $J$  , contained in $(-1,1)$  , and $C_2>0$  such that $$\frac{{Q^{\prime }}^{\prime }\left(x\right)Q\left(x\right)}{\left(Q^{\prime }\right(x){)}^2}\ge C_2,$$  for all $x\in (-1,1)\backslash J$  , except a subset with zero measure.

A function $f:I\to [0,+\infty )$  is cuasi-increasing if $\forall x<y\in I,\exists C>0$  , such that $f\left(x\right)<Cf\left(y\right).$  We prove Proposition  2.2 in several steps.

Lemma 2.3If a function $f\in C[0,+\infty )$  is positive and decreasing, then for $t>0$  ,

$$\begin{array}{c}I\left(t\right)\stackrel{\text{def}}{=}{\int }_0^{+\infty }f\left(u\right)sin\left(ut\right)du\in (0,+\infty ].\end{array}$$

(24)

Proof. ] With the change of variable $v=ut$  we can write $$I\left(t\right)=\frac{1}{t}{\int }_0^{\infty }f(v/t)sin\left(v\right)dv.$$  We denote $g\left(v\right)=f(v/t)$  ; then

$$\begin{array}{cc}tI\left(t\right)& ={\sum }_{k=0}^{\infty }{\int }_{k\pi }^{(k+1)\pi }g\left(v\right)sinvdv={\sum }_{k=0}^{\infty }{\int }_0^{\pi }g(k\pi +v)sin(k\pi +v)dv\end{array}$$

$$\begin{array}{cc}& ={\sum }_{k=0}^{\infty }{\int }_0^{\pi }g(2k\pi +v)sin(2k\pi +v)dv\end{array}$$

$$\begin{array}{cc}& +{\sum }_{k=0}^{\infty }{\int }_0^{\pi }g\left(\right(2k+1)\pi +v)sin\left(\right(2k+1)\pi +v)dv\end{array}$$

$$\begin{array}{cc}& ={\sum }_{k=0}^{\infty }{\int }_0^{\pi }g(2k\pi +v)sin\left(v\right)dv-{\sum }_{k=0}^{\infty }{\int }_0^{\pi }g\left(\right(2k+1)\pi +v)sin\left(v\right)dv\end{array}$$

$$\begin{array}{cc}& ={\sum }_{k=0}^{\infty }{\int }_0^{\pi }\left(g(2k\pi +v)-g\left(\right(2k+1)\pi +v)\right)sin\left(v\right)dv.\end{array}$$

Since each integral in the series is strictly positive, it proves ( 24 ).

Now we gather some properties of the digamma and trigamma functions in the following technical lemma:

Lemma 2.4For the digamma function $\psi \left(x\right)={\Gamma }^{\prime }\left(x\right)/\Gamma \left(x\right)$  the following statements hold: for $\lambda \ge 1$  ,

• i) $Re{\psi }^{\prime }(\lambda +it)$  is a strictly positive even function of $t\in \mathbb{R}$  , and

  $$\begin{array}{c}{lim}_{t\to \pm \infty }tRe{\psi }^{\prime }(\lambda +it)=0.\end{array}$$

  (25)

• ii) $Im\psi (\lambda +it)$  is an odd function of $t\in \mathbb{R}$  , strictly positive in $(0,+\infty )$  , and

  $$\begin{array}{c}Im\psi (\lambda \pm it)=\pm \frac{\pi }{2}+\mathcal{O}\left(t^{-1}\right),t\to +\infty .\end{array}$$

  (26)

Remark 2.5Limits  25 – 26 are valid in fact for $\lambda >0$  .

Proof. ] The symmetry of both the real and the imaginary parts of $\psi (\lambda +it)$  follows from the well known property $$\psi \left(\overline{z}\right)=\overline{\psi \left(z\right)};$$  so, we restrict our attention to $t>0$  .

For i) we consider the integral representation of the trigamma function,

$$\begin{array}{c}{\psi }^{\prime }\left(z\right)={\int }_0^{+\infty }{e}^{-uz}\frac{u}{1-e^{-u}}du,\end{array}$$

(27)

(see e.g. [12,formula3.41.371.6] ), from where $$Re{\psi }^{\prime }(\lambda +it)={\int }_0^{\infty }{e}^{-u\lambda }\frac{u}{1-e^{-u}}cos\left(ut\right)du.$$  Integrating by parts it can be reduced to $$Re{\psi }^{\prime }(\lambda +it)=\frac{1}{t}{\int }_0^{\infty }{e}^{-\lambda u}\frac{u{e}^{-u}-(1-\lambda u)(1-e^{-u})}{(1-e^{-u}{)}^2}sin\left(ut\right)du.$$  It is easy to check that for $\lambda \ge 1$  , function $$f\left(u\right)=e^{-\lambda u}\frac{u{e}^{-u}-(1-\lambda u)(1-e^{-u})}{(1-e^{-u}{)}^2},$$  is positive and decreasing on $(0,+\infty )$  . Hence, the first part of i) follows from Lemma  2.3 . On the other hand, by [1,formula6.4.12] , $${\psi }^{\prime }\left(z\right)=\frac{1}{z}+\mathcal{O}\left(\frac{1}{z^2}\right),z\to \infty ,|argz|<\pi ,$$  so that $$(\lambda +it){\psi }^{\prime }(\lambda +it)=1+\mathcal{O}\left(\frac{1}{t}\right),t\to \pm \infty .$$  In particular,

$$\begin{array}{cc}0& ={lim}_{t\to \pm \infty }Im\left((\lambda +it){\psi }^{\prime }(\lambda +it)\right)=\lambda {lim}_{t\to \pm \infty }Im\left({\psi }^{\prime }(\lambda +it)\right)+{lim}_{t\to \pm \infty }tRe{\psi }^{\prime }(\lambda +it)\end{array}$$

$$\begin{array}{cc}& ={lim}_{t\to \pm \infty }tRe{\psi }^{\prime }(\lambda +it),\end{array}$$

which proves  25 . For ii) we may use the series expansion [1,formula6.3.16] , $$\psi (1+z)=-\gamma +{\sum }_{k=1}^{\infty }\frac{z}{k(k+z)},-z/\in \mathbb{N},$$  according to which $$Im\psi (1+z)=-{\sum }_{k=1}^{\infty }Im\left(\frac{1}{k+z}\right)={\sum }_{k=1}^{\infty }\frac{Im\left(z\right)}{|k+z{|}^2}>0\text{if}Im\left(z\right)>0.$$  Finally, by the asymptotic formula [1,formula6.3.18] , $$\psi \left(z\right)=logz+\mathcal{O}\left(\frac{1}{z}\right),z\to \infty ,|argz|<\pi ,$$  from which  26 is immediate.

Now we are ready to analyze whether the weight belongs to $\mathcal{ℱ}(C^2+)$  . Conditions a) and b) are a straightforward consequence of the following statement:

Lemma 2.6The even function $Q\in C^{\infty }(-1,1)$  satisfies $\frac{d^{k}Q\left(x\right)}{d{x}^k}>0\text{for}x\in (0,1)\text{and}k=0,1,2\text{.}$ 

Proof. ] By definition, $Q\left(0\right)=0$  , and by symmetry, $Q^{\prime }\left(0\right)=0$  , so $$Q^{\prime \prime }\left(x\right)>0\Rightarrow Q^{\prime }\left(x\right)>0\Rightarrow Q\left(x\right)>0,x\in (0,1).$$  But for $x\in (0,1)$  we have $t>0$  , and

$$\begin{array}{c}\begin{array}{cccc}{Q}^{\prime \prime }\left(x\right)& =\frac{(\lambda -1/2)(1+x^2)}{(1-x^2{)}^2}+\frac{a{x}^2+2a}{(1-x^2{)}^2}+\frac{3ax}{(1-x^2{)}^{5/2}}\frac{1}{2}(\pi -2arccosx)& & +\frac{3ax}{(1-x^2{)}^{5/2}}Im\psi (\lambda +it)+\frac{a^2}{(1-x^2{)}^3}Re{\psi }^{\prime }(\lambda +it)>0,\end{array}\end{array}$$

(28)

where we have used Lemma  2.4 .

Since by ( 21 ),

$$\begin{array}{c}Q\left(x\right)=-\frac{1}{2}logw\left(x\right)+\frac{1}{2}logw\left(0\right)=\frac{1}{2}\frac{2\pi a\left|x\right|}{\sqrt{1-x^2}}+\mathcal{O}\left(1\right),\left|x\right|\to 1^{-},\end{array}$$

(29)

condition c) also trivially holds.

We turn now to the even function $$T\left(x\right)=\frac{x{Q}^{\prime }\left(x\right)}{Q\left(x\right)};$$  let us show that it is cuasi-increasing in $(0,1)$  . Since $T$  is continuous and positive on the bounded interval $(0,1)$  , it is sufficient to show that it does not blow up at the left end point, nor it vanishes at the right one. Recall that $$Q\left(0\right)=Q^{\prime }\left(0\right)=0,{Q^{\prime }}^{\prime }\left(0\right)=(\lambda -1/2)+2a+a^2{\psi }^{\prime }\left(\lambda \right)>0;$$  in particular,

$$\begin{array}{c}Q\left(x\right)=\frac{(\lambda -1/2)+2a+a^2{\psi }^{\prime }\left(\lambda \right)}{2}{x}^2+\mathcal{O}\left(x^3\right),x\to 0.\end{array}$$

(30)

Hence,

$$\begin{array}{c}{lim}_{x\to 0}T\left(x\right)={lim}_{x\to 0}\frac{x{Q}^{\prime }\left(x\right)}{Q\left(x\right)}={lim}_{x\to 0}\frac{Q^{\prime }\left(x\right)/x}{Q\left(x\right)/x^2}=\frac{{Q^{\prime }}^{\prime }\left(0\right)}{\frac{1}{2}{Q^{\prime }}^{\prime }\left(0\right)}=2>0.\end{array}$$

(31)

On the other hand,

$$\begin{array}{c}{Q}^{\prime }\left(x\right)=\frac{(\lambda -1/2)x}{1-x^2}+\frac{ax}{1-x^2}+\frac{1}{2}\frac{(\pi -2arccosx)a}{(1-x^2{)}^{3/2}}+\frac{a}{(1-x^2{)}^{3/2}}Im\psi (\lambda +it),\end{array}$$

(32)

and by Lemma  2.4 ,

$$\begin{array}{c}{lim}_{x\to 1^{-}}(1-x^2{)}^{3/2}{Q}^{\prime }(x)=a\pi .\end{array}$$

(33)

Together with ( 29 ) it shows that

$$\begin{array}{c}{lim}_{x\to 1^{-}}T\left(x\right)={lim}_{x\to 1^{-}}\frac{x{Q}^{\prime }\left(x\right)}{Q\left(x\right)}={lim}_{x\to 1^{-}}\frac{1}{1-x^2}\frac{\pi a}{\pi a}=+\infty .\end{array}$$

(34)

In conclusion, $T$  is cuasi-increasing in $(0,1)$  .

On the other hand, if $\zeta \in (0,1)$  is a local minimum of $T$  , then $$\frac{dlog\left(T\right(x)}{dx}{|}_{x=\zeta }=0\Rightarrow \frac{1}{\zeta }+\frac{Q^{\prime \prime }}{Q^{\prime }}\left(\zeta \right)-\frac{Q^{\prime }}{Q}\left(\zeta \right)=0,$$  or equivalently, $$T\left(\zeta \right)=1+\frac{\zeta Q^{\prime \prime }\left(\zeta \right)}{Q^{\prime }\left(\zeta \right)}>1.$$  Taking into account also the behavior at $x=0$  and $x=1$  (see  31 and  34 ), we obtain  23 .

Finally, let us check conditions e) and f ). Denote $$H\left(x\right)\stackrel{\text{def}}{=}\frac{Q\left(x\right)Q^{\prime \prime }\left(x\right)}{\left(Q^{\prime }\right(x){)}^2}.$$  This is an even, continuous and positive function on $(-1,1)$  , with $H\left(0\right)=1/2$  , where we have used  30 . On the other hand, $${lim}_{x\to 1^{-}}(1-x^2{)}^{1/2}Q(x)={lim}_{x\to 1^{-}}(1-x^2{)}^{3/2}{Q}^{\prime }\left(x\right)=a\pi ,$$  where we have used ( 29 ) and ( 33 ). Also from  28 it follows that $${lim}_{x\to 1^{-}}(1-x^2{)}^{5/2}{Q}^{\prime \prime }(x)=\frac{3}{2}a\pi +\frac{3}{2}a\pi +a{lim}_{x\to 1^{-}}tRe{\psi }^{\prime }(\lambda +it),$$  and by  25 , ${lim}_{x\to 1^{-}}(1-x^2{)}^{5/2}{Q^{\prime }}^{\prime }(x)=3a\pi$  . Gathering these identities we obtain that $${lim}_{x\to 1^{-}}H\left(x\right)={lim}_{x\to 1^{-}}\frac{(1-x^2{)}^{1/2}Q(x\left)\right(1-x^2{)}^{5/2}{Q}^{\prime \prime }\left(x\right)}{(1-x^2{)}^3(Q^{\prime }\left(x\right){)}^2}=3.$$  Hence, $H$  can be extended as a strictly positive and continuous (and thus, uniformly continuous) function on $[-1,1]$  ; from this fact conditions e) and f ) follow automatically. This concludes the proof of Proposition  2.2 .

3 Equilibrium measure

The equilibrium measure ${\mu }_n$  on $[-1,1]$  of total weight $n$  in the external field $Q$  plays a prominent role in the asymptotics of the orthogonal polynomials, and we gather in this section some of its properties needed further. By Lemma  2.6 , function $Q\in C^{\infty }(-1,1)$  is strictly convex, and $Q(-1^{+})=Q\left(1^{-}\right)=+\infty$  . In consequence, (see e.g. [19] ), ${\mu }_n$  is absolutely continuous and supported on the interval $[-{\alpha }_n,{\alpha }_n]$  , where ${\alpha }_n$  is the MRS number, defined by  13 . If we denote by ${\sigma }_n\left(x\right)$  the density ( ${\mu }_n^{\prime }$  ) of ${\mu }_n$  , then $${\int }_{-{\alpha }_n}^{{\alpha }_n}{\sigma }_n\left(x\right)dx=n,$$  and the characterizing property of the equilibrium is $$V^{{\mu }_n}\left(x\right)+Q\left(x\right)\{\begin{array}{cc}=b_n(=const),& x\in [-{\alpha }_n,{\alpha }_n],\\ >b_n,& {\alpha }_n<\left|x\right|\le 1,\end{array}$$  where $V^{{\mu }_n}$  is the logarithmic potential of ${\mu }_n$  (cf.  15 ). We analyze first the asymptotic behavior of $\left\{{\alpha }_n\right\}$  , needed in the proof of Corollary  1.6 . It is known that ${\alpha }_n\to 1^{-}$  , and even more, that $1-{\alpha }_n=\mathcal{O}(1/n)$  (see [13,§1.6] ), but we are looking for a more precise information. The following technical lemma is useful for the estimation of the behavior of the integral in  13 :

Lemma 3.1Let $f(u,x)$  be defined for $u,x\in [0,1]$  with $f$  and $\partial f/\partial x$  continuous in $[0,1{]}^2$  , $f(1,1)\ne 0$  ; then when $u\to 1^{-}$  ,

$$\begin{array}{cc}{\int }_0^1\frac{f(u,x)}{\sqrt{1-u^2{x}^2}}\frac{dx}{\sqrt{1-x^2}}& =-\frac{f(1,1)}{2}log(1-u)\left(1+o\left(1\right)\right),\end{array}$$

(35)

$$\begin{array}{cc}{\int }_0^1\frac{f(u,x)}{(1-u^2{x}^2)}\frac{dx}{\sqrt{1-x^2}}& =\frac{\pi f(1,1)}{2\sqrt{2}}\frac{1}{\sqrt{1-u}}\left(1+o\left(1\right)\right),\end{array}$$

(36)

$$\begin{array}{cc}{\int }_0^1\frac{f(u,x)}{(1-u^2{x}^2{)}^{3/2}}\frac{dx}{\sqrt{1-x^2}}& =\frac{f(1,1)}{2}\frac{1}{1-u}\left(1+o\left(1\right)\right).\end{array}$$

(37)

Formula  35 appears in [12,formula8.113.3] for $f(u,x)=1$  ; the proof is standard, and we omit it here for the sake of brevity.

Now we can obtain the first two terms of the asymptotics of ${\alpha }_n$  :

Proposition 3.2The MRS numbers ${\alpha }_n$  satisfy $${\alpha }_n=1-\frac{a}{n}+o\left(\frac{1}{n}\right),n\to \infty .$$ 

Proof. ] Formula  13 defining the MRS numbers may be rewritten as

$$\begin{array}{c}\frac{2}{\pi }{\int }_0^1\frac{{\alpha }_{n}x{Q}^{\prime }\left({\alpha }_{n}x\right)}{\sqrt{1-x^2}}dx=n,\end{array}$$

(38)

which motivates the study of the asymptotics (as $u\to 1^{-}$  ) of the integral

$$\begin{array}{cc}\frac{2}{\pi }& {\int }_0^1\frac{ux{Q}^{\prime }\left(ux\right)}{\sqrt{1-x^2}}dx\end{array}$$

$$\begin{array}{cc}& =\frac{1}{\pi }{\int }_0^1\frac{ux\left((\lambda -1/2)2ux+2aux\right)}{(1-u^2{x}^2)\sqrt{1-x^2}}dx+\frac{1}{\pi }{\int }_0^1\frac{ux(-2)aarccos\left(ux\right)}{(1-u^2{x}^2{)}^{3/2}\sqrt{1-x^2}}dx\end{array}$$

$$\begin{array}{cc}& +\frac{1}{\pi }{\int }_0^1\frac{ux\left(\pi a+2aIm\psi \left(\lambda +i\frac{aux}{\sqrt{1-u^2{x}^2}}\right)\right)}{(1-u^2{x}^2{)}^{3/2}\sqrt{1-x^2}}dx=I_1\left(u\right)+I_2\left(u\right)+I_3\left(u\right).\end{array}$$

Using Lemma  3.1 we have that for $u\to 1^{-}$  ,

$$\begin{array}{cc}{I}_1\left(u\right)& =\frac{1}{\pi }\left(\right(2\lambda -1)+2a)\frac{\pi }{2\sqrt{2}}\frac{1}{\sqrt{1-u}}\left(1+o\left(1\right)\right)=o\left(\frac{1}{1-u}\right),\end{array}$$

$$\begin{array}{cc}{I}_3\left(u\right)& =\frac{1}{\pi }(\pi a+2a\pi /2)\frac{1}{2}\frac{1}{1-u}\left(1+o\left(1\right)\right)=\frac{a}{1-u}\left(1+o\left(1\right)\right).\end{array}$$

Moreover, since $$f(u,x)=\frac{arccos\left(ux\right)}{\sqrt{1-u^2{x}^2}},$$  satisfies the conditions of Lemma  3.1 , we have that $$I_2\left(u\right)=-\frac{1}{\pi }{\int }_0^1\frac{ux2af(u,x)}{(1-u^2{x}^2)\sqrt{1-x^2}}=o\left(\frac{1}{1-u}\right),u\to 1^{-}.$$  Summarizing, $$\frac{2}{\pi }{\int }_0^1\frac{ux{Q}^{\prime }\left(ux\right)}{\sqrt{1-x^2}}dx=\frac{a}{1-u}\left(1+o\left(1\right)\right),$$  and equation  38 for the MRS numbers can be rewritten as $$\frac{a}{1-{\alpha }_n}\left(1+o\left(1\right)\right)=n,$$  which proves the statement.

We turn now to the analysis of the density ${\sigma }_n$  . It is convenient to introduce the normalized translation of ${\sigma }_n$  to $\Delta$  ,

$$\begin{array}{c}{\sigma }_n^{*}\left(u\right)=\frac{{\alpha }_n}{n}{\sigma }_n\left({\alpha }_{n}u\right),u\in \Delta ,\end{array}$$

(39)

as well as the cumulative distribution

$$\begin{array}{c}{\Phi }_n\left(\theta \right)\stackrel{\text{def}}{=}\pi {\int }_{cos\theta }^1{\sigma }_n^{*}\left(t\right)dt,\end{array}$$

(40)

which is obviously a smooth and strictly increasing function on $[0,\pi ]$  ; moreover, ${\Phi }_n:[0,\pi ]\to [0,\pi ]$  is a bijection, and the inverse function ${\Phi }_n^{[-1]}$  exists (observe that our definition differs in normalization from that used in [13, 14] ). We summarize some properties of ${\Phi }_n$  in the following lemma.

Lemma 3.3For ${\Phi }_n$  defined in  40 ,

• i) ${\Phi }_n^{\prime }\left(\theta \right)\to 1$  pointwise in $(0,\pi )$  .
• ii) $${\int }_0^{\pi }\left|\frac{1}{{\Phi }_n^{\prime }\left({\Phi }_n^{[-1]}\right(\eta \left)\right)}-1\right|d\eta =o\left(1\right),n\to \infty .$$ 
• iii) ${\Phi }_n\left(\theta \right)\to \theta$  and ${\Phi }_n^{[-1]}\left(\theta \right)\to \theta$  as $n\to \infty$  uniformly in $[0,\pi ]$  .

Proof. ] Let $0<\varepsilon <1/2$  , and $x=cos\theta \in (-1+2\varepsilon ,1-2\varepsilon )$  . By Lemma 6.5 of [13] , $${\Phi }_n^{\prime }\left(\theta \right)-1=\pi {\sigma }_n^{*}(cos\theta )sin\theta -1=\frac{{\alpha }_n}{\pi n}PV{\int }_{-1}^1\frac{Q^{\prime }\left({\alpha }_{n}u\right)\sqrt{1-u^2}}{u-x}du,$$  where $PV$  means the principal value of integral. Hence,

$$\begin{array}{c}\begin{array}{cccccc}{\Phi }_n^{\prime }\left(\theta \right)-1& =\frac{{\alpha }_n}{\pi n}({\int }_{-1+\varepsilon }^{1-\varepsilon }\frac{Q^{\prime }\left({\alpha }_{n}u\right)\sqrt{1-u^2}-Q^{\prime }\left({\alpha }_{n}x\right)\sqrt{1-x^2}}{u-x}du& & +Q^{\prime }\left({\alpha }_{n}x\right)\sqrt{1-x^2}PV{\int }_{-1+\varepsilon }^{1-\varepsilon }\frac{du}{u-x}+{\int }_{-1}^{-1+\varepsilon }\frac{Q^{\prime }\left({\alpha }_{n}u\right)\sqrt{1-u^2}}{u-x}du& & +{\int }_{1-\varepsilon }^1\frac{Q^{\prime }\left({\alpha }_{n}u\right)\sqrt{1-u^2}}{u-x}du).\end{array}\end{array}$$

(41)

The first two terms within parentheses in the right hand side of  41 are uniformly bounded. Let us estimate $$\left|{\int }_{1-\varepsilon }^1\frac{Q^{\prime }\left({\alpha }_{n}u\right)\sqrt{1-u^2}}{u-x}du\right|\le \frac{1}{\varepsilon }{\int }_{1-\varepsilon }^1{Q}^{\prime }\left({\alpha }_{n}u\right)\sqrt{1-u^2}du$$  (the remaining term is analyzed in a similar fashion). Integrating by parts, $${\int }_{1-\varepsilon }^1{Q}^{\prime }\left({\alpha }_{n}u\right)\sqrt{1-u^2}du=-\frac{1}{{\alpha }_n}Q\left({\alpha }_n\right(1-\varepsilon \left)\right)\sqrt{\varepsilon (2-\varepsilon )}+\frac{1}{{\alpha }_n}{\int }_{1-\varepsilon }^1\frac{uQ\left({\alpha }_{n}u\right)}{\sqrt{1-u^2}}du,$$  and using  35 we get that $$\left|{\Phi }_n^{\prime }\left(\theta \right)-1\right|=\mathcal{O}\left(\frac{log\left(n\right)}{n}\right),$$  which proves i).

On the other hand, by [14,lemma4.2a)] , the sequence $\left|{\Phi }_n^{\prime }\right(\theta \left)\right|$  is uniformly bounded on $[0,\pi ]$  . Thus, by the dominated convergence theorem,

$$\begin{array}{cc}{\int }_0^{\pi }\left|\frac{1}{{\Phi }_n^{\prime }\left({\Phi }_n^{[-1]}\right(\eta \left)\right)}-1\right|d\eta & ={\int }_0^{\pi }\left|\frac{1}{{\Phi }_n^{\prime }\left(\theta \right)}-1\right|{\Phi }_n^{\prime }\left(\theta \right)d\theta \end{array}$$

$$\begin{array}{cc}& ={\int }_0^{\pi }\left|1-{\Phi }_n^{\prime }\left(\theta \right)\right|d\theta \to 0,n\to \infty .\end{array}$$

Furthermore, given $\theta \in [0,\pi ]$  , $$\left|{\Phi }_n\left(\theta \right)-\theta \right|=\left|{\int }_0^{\theta }\left({\Phi }_n^{\prime }\right(\eta )-1)d\eta \right|\le {\int }_0^{\pi }\left|{\Phi }_n^{\prime }\left(\eta \right)-1\right|d\eta =o\left(1\right),$$  and so the uniform convergence of ${\Phi }_n$  on $[0,\pi ]$  follows again by the dominated convergence theorem. Finally, if $\theta ={\Phi }_n\left(\eta \right)$  ,

$$\begin{array}{cc}& \left|{\Phi }_n^{[-1]}\left(\theta \right)-\theta \right|=\left|{\Phi }_n^{[-1]}\left({\Phi }_n\right(\eta \left)\right)-{\Phi }_n\left(\eta \right)\right|=\left|\eta -{\Phi }_n\left(\eta \right)\right|,\end{array}$$

showing that ${\Phi }_n^{[-1]}\left(\theta \right)$  converges uniformly to $\theta$  on $[0,\pi ]$  . This concludes the proof.

4 Asymptotics of $F_n$  : proof of Theorem  1.1 

We follow the scheme of proof of [3] . We have established already that the weight $w\in \mathcal{ℱ}(C^2+)$  ; furthermore, $w\left(x\right)>0$  on $(-1,1)$  , and in consequence, it is an Erdos-Turan weight. One of the most relevant facts about these weights is that the sequence $p_n^2\left(x\right)w\left(x\right)dx$  converges in the weak- $*$  topology to the equilibrium (Robin) measure $\mu$  of $[-1,1]$  (see [16, 18] ). In other words, for any $f\in C[-1,1]$  ,

$$\begin{array}{c}{lim}_{n\to \infty }{\int }_{-1}^{1}f\left(x\right)p_n^2\left(x\right)w\left(x\right)dx=\frac{1}{\pi }{\int }_{-1}^{1}f\left(x\right)\frac{dx}{\sqrt{1-x^2}}.\end{array}$$

(42)

We make use also of the following technical lemmas. In the sequel we write that $x_n\sim y_n$  if the ratios $x_n/y_n$  and $y_n/x_n$  are uniformly bounded in $n$  .

Lemma 4.1 $${\int }_{\Delta \backslash [-{\alpha }_n,{\alpha }_n]}log\left(p_n^2\right(x)w(x\left)\right)p_n^2\left(x\right)w\left(x\right)dx=o\left(1\right),n\to \infty .$$ 

Proof. ] Since $w\in \mathcal{ℱ}(C^2+)$  , by [13,Theorem1.18] , $${sup}_{x\in \Delta }\left|p_n\left(x\right)\sqrt{w\left(x\right)}\right|\sim n^{1/6}({\alpha }_n{)}^{-1/3}{\left(\frac{T\left({\alpha }_n\right)}{{\alpha }_n}\right)}^{1/6}\sim n^{1/3},$$  where we have taken into account Proposition  3.2 . Thus, for $\varepsilon >0$  , there exists $C_1>0$  such that $${\left(p_n^2\left(x\right)w\left(x\right)\right)}^{1+\varepsilon }\le C_1{n}^{(2+2\varepsilon )/3}$$  and so, if $\varepsilon <1/2$  , then

$$\begin{array}{cc}& {\int }_{\Delta \backslash [-{\alpha }_n,{\alpha }_n]}log\left(p_n^2\right(x)w(x\left)\right)p_n^2\left(x\right)w\left(x\right)dx\le {\int }_{\Delta \backslash [-{\alpha }_n,{\alpha }_n]}\left(p_n^2\right(x)w(x){)}^{1+\varepsilon }dx\end{array}$$

$$\begin{array}{cc}\le & C_1{n}^{(2+2\varepsilon )/3}{\int }_{\Delta \backslash [-{\alpha }_n,{\alpha }_n]}dx\le C_2{n}^{\frac{2+2\varepsilon }{3}-1}=C_2{n}^{\frac{-1+2\varepsilon }{3}}=o\left(1\right).\end{array}$$

For a lower bound it is sufficient to take into account that function

$$\begin{array}{c}\mathcal{ℛ}\left(y\right)=y^{2}log\left(y^2\right)\end{array}$$

(43)

is bounded from below on $[0,+\infty )$  , and thus $${\int }_{\Delta \backslash [-{\alpha }_n,{\alpha }_n]}log\left(p_n^2\right(x)w(x\left)\right)p_n^2\left(x\right)w\left(x\right)dx\ge C{\int }_{\Delta \backslash [-{\alpha }_n,{\alpha }_n]}dx=o\left(1\right),n\to \infty ,$$  which concludes the proof.

Lemma 4.2

$$\begin{array}{c}{lim}_{n\to \infty }{\int }_{-{\alpha }_n}^{{\alpha }_n}log\left(\sqrt{{\alpha }_n^2-x^2}\right)p_n^2\left(x\right)w\left(x\right)dx=-log\left(2\right).\end{array}$$

(44)

Proof. ] Let $${\ell }_n\left(x\right)\stackrel{\text{def}}{=}\{\begin{array}{cc}log\left(\sqrt{{\alpha }_n^2-x^2}\right),& x\in (-{\alpha }_n,{\alpha }_n),\\ 0,& x\in [-1,1]\backslash (-{\alpha }_n,{\alpha }_n).\end{array}$$  Obviously, ${\ell }_n\left(x\right)\to log\left(\sqrt{1-x^2}\right)$  pointwise for $x\in (-1,1)$  . Furthermore, there exists $C_1>0$  such that for $t>1$  it holds that $(log(t){)}^3<C_{1}t$  . Then, by Lebesgue dominated convergence theorem, $${\Vert {\ell }_n\left(x\right)-log\left(\sqrt{1-x^2}\right)\Vert }_{L^3}=o\left(1\right),n\to \infty$$  (here and in the sequel $\parallel \cdot {\parallel }_{L^p}$  denotes the $p$  -norm with respect to the Lebesgue measure on $\Delta$  ).

Furthermore, from [13,Theorem13.6] it follows that the sequence $\parallel p_n\sqrt{w}{\parallel }_{L^p}$  is uniformly bounded as long as $p<4$  (and in particular, for $p=3$  ). Thus, by Holder inequality,

$$\begin{array}{cc}\left|{\int }_{-1}^1\left({\ell }_n\left(x\right)-log\left(\sqrt{1-x^2}\right)\right)p_n^2\left(x\right)w\left(x\right)dx\right|& \le {\Vert {\ell }_n\left(x\right)-log\left(\sqrt{1-x^2}\right)\Vert }_{L^3}{\Vert p_n^{2}w\Vert }_{L^{3/2}},\end{array}$$

which implies that

$$\begin{array}{c}\begin{array}{cc}{\int }_{-{\alpha }_n}^{{\alpha }_n}log\left(\sqrt{{\alpha }_n^2-x^2}\right)p_n^2\left(x\right)w\left(x\right)dx& ={\int }_{-1}^{1}log\left(\sqrt{1-x^2}\right)p_n^2\left(x\right)w\left(x\right)dx+o\left(1\right),n\to \infty .\end{array}\end{array}$$

(45)

On the other hand, if for $\varepsilon \in (0,1)$  we denote $${log}_{\varepsilon }\left(x\right)\stackrel{\text{def}}{=}max\{log(\varepsilon ),log(\sqrt{1-x^2}\left)\right\}\in C[-1,1],$$  then

$$\begin{array}{cc}& \left|{\int }_{-1}^{1}log\left(\sqrt{1-x^2}\right)p_n^2\left(x\right)w\left(x\right)dx-\frac{1}{\pi }{\int }_{-1}^{1}log\left(\sqrt{1-x^2}\right)\frac{dx}{\sqrt{1-x^2}}\right|\end{array}$$

$$\begin{array}{cc}\le & \left|{\int }_{-1}^1\left(log\left(\sqrt{1-x^2}\right)-{log}_{\varepsilon }\left(x\right)\right)p_n^2\left(x\right)w\left(x\right)dx\right|\end{array}$$

$$\begin{array}{cc}+& \left|{\int }_{-1}^1{log}_{\varepsilon }\left(x\right)\left(p_n^2\left(x\right)w\left(x\right)-\frac{1}{\pi \sqrt{1-x^2}}\right)dx\right|\end{array}$$

$$\begin{array}{cc}+& \left|{\int }_{-1}^1\left({log}_{\varepsilon }\left(x\right)-log\left(\sqrt{1-x^2}\right)\right)\frac{1}{\pi \sqrt{1-x^2}}dx\right|=I_1+I_2+I_3.\end{array}$$

By  42 , $I_2=o\left(1\right)$  , as $n\to \infty$  , while $$I_3=-{\int }_{\sqrt{1-x^2}<\varepsilon }log\left(\sqrt{1-x^2}\right)\frac{1}{\pi \sqrt{1-x^2}}dx=\mathcal{O}\left(\varepsilon \right).$$  Using the same arguments as for  45 we find also that $I_1=\mathcal{O}\left(\varepsilon \right)$  . Taking into account that $\varepsilon >0$  is arbitrary and using  45 we obtain that

$$\begin{array}{cc}{\int }_{-{\alpha }_n}^{{\alpha }_n}log\left(\sqrt{{\alpha }_n^2-x^2}\right)p_n^2\left(x\right)w\left(x\right)dx& =\frac{1}{\pi }{\int }_{-1}^{1}log\left(\sqrt{1-x^2}\right)\frac{dx}{\sqrt{1-x^2}}+o\left(1\right),n\to \infty .\end{array}$$

The identity $$\frac{1}{\pi }{\int }_{-1}^{1}log\left(\sqrt{1-x^2}\right)\frac{dx}{\sqrt{1-x^2}}=-log\left(2\right)$$  is straightforward, which concludes the proof.

Remark 4.3This result is not surprising: if we denote by ${\widehat{\nu }}_n$  the absolutely continuous measure on $[-{\alpha }_n,{\alpha }_n]$  with ${\widehat{\nu }}_n^{\prime }\left(x\right)=p_n^2\left(x\right)w\left(x\right)$  , then ${\widehat{\nu }}_n\to \mu$  in the weak- $*$  topology, where $\mu$  is the Robin measure of $\Delta$  . The integral in the left hand side of  44 can be rewritten as $$-\frac{1}{2}\left(V^{{\widehat{\nu }}_n}(-{\alpha }_n)+V^{{\widehat{\nu }}_n}\left({\alpha }_n\right)\right)⟶-\frac{1}{2}\left(V^{\mu }(-1)+V^{\mu }\left(1\right)\right)=-log\left(2\right).$$ 

Now we turn to the proof of Theorem  1.1 . Using function ${\Phi }_n$  introduced in  40 , let us denote

$$\begin{array}{c}\begin{array}{cccc}{f}_n({\alpha }_{n}cos\theta )& \stackrel{\text{def}}{=}\sqrt{{\alpha }_n}{p}_n({\alpha }_{n}cos\theta )\sqrt{w({\alpha }_{n}cos\theta )}\sqrt{sin\theta },& g_n({\alpha }_{n}cos\theta )& \stackrel{\text{def}}{=}\sqrt{\frac{2}{\pi }}cos\left(\frac{\theta }{2}-\frac{\pi }{4}+n{\Phi }_n\left(\theta \right)\right).\end{array}\end{array}$$

(46)

Since $w\in \mathcal{ℱ}(C^2+)$  , by [13,Theorem15.1andLemma15.4] ,

$$\begin{array}{c}{\int }_0^{\pi }\left|f_n\right({\alpha }_{n}cos\theta )-g_n({\alpha }_{n}cos\theta \left)\right|d\theta =o\left(1\right),n\to \infty ,\end{array}$$

(47)

(where we have used that on the bounded interval convergence in $L^2$  is stronger than in $L^1$  ). Let us rewrite the definition of $F_n$  as

$$\begin{array}{cc}{F}_n=& -{\int }_{-{\alpha }_n}^{{\alpha }_n}log\left(p_n^2\right(x)w(x\left)\right)p_n^2\left(x\right)w\left(x\right)dx\end{array}$$

$$\begin{array}{cc}& -{\int }_{\Delta \backslash [-{\alpha }_n,{\alpha }_n]}log\left(p_n^2\right(x)w(x\left)\right)p_n^2\left(x\right)w\left(x\right)dx\end{array}$$

$$\begin{array}{cc}=& -{\int }_{-{\alpha }_n}^{{\alpha }_n}log\left(p_n^2\right(x)w(x\left)\right)p_n^2\left(x\right)w\left(x\right)dx+o\left(1\right),\end{array}$$

where we have used Lemma  4.1 . With the notation  43 and ( 46 ) it is equivalent to

$$\begin{array}{c}\begin{array}{cccc}{F}_n& =-{\int }_{-{\alpha }_n}^{{\alpha }_n}\mathcal{ℛ}\left(f_n\right(x\left)\right)\frac{dx}{\sqrt{{\alpha }_n^2-x^2}}+{\int }_{-{\alpha }_n}^{{\alpha }_n}log\left(\sqrt{{\alpha }_n^2-x^2}\right)p_n^2\left(x\right)w\left(x\right)dx+o\left(1\right)& & =-{\int }_0^{\pi }\mathcal{ℛ}\left(f_n\right({\alpha }_{n}cos\theta \left)\right)d\theta -log2+o\left(1\right),n\to \infty ,\end{array}\end{array}$$

(48)

(see Lemma  4.2 ). Since by [13,Theorem1.17] there exists a constant $M>\sqrt{2/\pi }$  such that for all $n\in \mathbb{N}$  , $\left|f_n\right({\alpha }_{n}cos\theta \left)\right|\le M$  , for $\theta \in [0,\pi ]$  , we get by ( 47 )

$$\begin{array}{cc}& {\int }_0^{\pi }\left|\mathcal{ℛ}\left(f_n\right({\alpha }_{n}cos\theta \left)\right)-\mathcal{ℛ}\left(g_n\right({\alpha }_{n}cos\theta \left)\right)\right|d\theta \end{array}$$

$$\begin{array}{cc}\le & {max}_{y\in [0,M]}{\mathcal{ℛ}}^{\prime }\left(y\right){\int }_0^{\pi }\left|f_n({\alpha }_{n}cos\theta )-g_n({\alpha }_{n}cos\theta )\right|d\theta =o\left(1\right),\end{array}$$

which yields $$F_n=-{\int }_0^{\pi }\mathcal{ℛ}\left(g_n({\alpha }_{n}cos\theta )\right)d\theta -log2+o\left(1\right),n\to \infty .$$  With the change of variable $\eta ={\Phi }_n\left(\theta \right)$  we rewrite

$$\begin{array}{cc}{F}_n& =-{\int }_0^{\pi }\mathcal{ℛ}\left(\sqrt{\frac{2}{\pi }}cos\left({\Phi }_n^{[-1]}\left(\eta \right)-\frac{\pi }{4}+n\eta \right)\right)\frac{d\eta }{{\Phi }_n^{\prime }\left({\Phi }_n^{[-1]}\right(\eta \left)\right)}-log2+o\left(1\right)\end{array}$$

$$\begin{array}{cc}& =-{\int }_0^{\pi }\mathcal{ℛ}\left(\sqrt{\frac{2}{\pi }}cos\left({\Phi }_n^{[-1]}\left(\eta \right)-\frac{\pi }{4}+n\eta \right)\right)d\eta -log2+o\left(1\right)\end{array}$$

$$\begin{array}{cc}& =-{\int }_0^{\pi }\mathcal{ℛ}\left(\sqrt{\frac{2}{\pi }}cos\left(\eta -\frac{\pi }{4}+n\eta \right)\right)d\eta -log2+o\left(1\right),\end{array}$$

(49)

where we have used Lemma  3.3 . It remains to use the following analogue of the Lebesgue lemma, proved in [2] under weaker conditions:

Lemma 4.4Let $g$  be a $\pi$  -periodic continuous function on $[0,+\infty )$  , and $h\in C[0,\pi ]$  .

Then $${\int }_0^{\pi }f(n\theta +h(\theta \left)\right)d\theta ={\int }_0^{\pi }f\left(\theta \right)d\theta +o\left(1\right).$$ 

Applying this lemma to  49 , we obtain finally $$F_n=-1+log\left(2\right)+log\left(\pi \right)-log\left(2\right)+o\left(1\right)=log\left(\pi \right)-1+o\left(1\right),n\to \infty .$$ 

5 Asymptotics of $G_n$  : proof of Theorem  1.3 

We start again with some technical results:

Lemma 5.1When $n\to \infty$  , $${\int }_0^{1}x(1-x^2)s^{\prime }\left(x\right)p_n^2\left(x\right)w\left(x\right)dx=B_2+o\left(1\right),$$  where

$$\begin{array}{c}{B}_2={\int }_0^{1}x(1-x^2)s^{\prime }\left(x\right)\frac{1}{\pi \sqrt{1-x^2}}dx.\end{array}$$

(50)

Proof. ] By  42 , it is sufficient to show that $(1-x^2)s^{\prime }\left(x\right)$  can be extended as a continuous function to the whole interval $\Delta$  ; for this purpose we only need to show that the limit $${lim}_{x\to 1^{-}}(1-x^2)s^{\prime }\left(x\right)$$  exists. From the explicit expression for $s$  it is easy to find that for $x\in (0,1)$  ,

$$\begin{array}{cc}(1-x^2)s^{\prime }\left(x\right)=& (1-x^2)\frac{w^{\prime }\left(x\right)}{w\left(x\right)}-(1-x^2)\frac{w_0^{\prime }\left(x\right)}{w_0\left(x\right)}=(1-x^2)\frac{w^{\prime }\left(x\right)}{w\left(x\right)}+\frac{2\pi a}{\sqrt{1-x^2}}\end{array}$$

$$\begin{array}{cc}=& -2x(\lambda -1/2)-2ax+2arccosx\frac{a}{\sqrt{1-x^2}}\end{array}$$

$$\begin{array}{cc}& -2\frac{t}{x}\left(Im\psi (\lambda +it)-\frac{\pi }{2}\right).\end{array}$$

It remains to use  26 , and the statement follows.

Let us denote

$$\begin{array}{c}{p}_n\left(x\right)={\gamma }_n{x}^n+{\beta }_n{x}^{n-2}+\text{lower degree terms};\end{array}$$

(51)

the explicit expression for ${\gamma }_n$  was given in  17 .

Lemma 5.2The followings identities hold:

$$\begin{array}{c}{\int }_{-1}^{1}x{p}_n\left(x\right)p_n^{\prime }\left(x\right)w\left(x\right)dx=n,\end{array}$$

(52)

and

$$\begin{array}{c}{\int }_{-1}^1{x}^3{p}_n\left(x\right)p_n^{\prime }\left(x\right)w\left(x\right)dx=n(a_{n+1}^2+a_n^2)-2a_n{a}_{n-1}\frac{{\beta }_n}{{\gamma }_{n-2}},\end{array}$$

(53)

where $a_n$  are the coefficients of the recurrence relation ( 6 ), and ${\gamma }_n$  , ${\beta }_n$  are the coefficients of $p_n$  defined in ( 51 ).

Proof. ] By the recurrence relation ( 6 ), $$x{p}_n\left(x\right)p_n^{\prime }\left(x\right)=a_{n+1}{p}_{n+1}\left(x\right)p_n^{\prime }\left(x\right)+a_n{p}_{n-1}\left(x\right)p_n^{\prime }\left(x\right),$$  so that

$$\begin{array}{cc}{\int }_{-1}^{1}x{p}_n\left(x\right)p_n^{\prime }\left(x\right)w\left(x\right)dx& =a_n{\int }_{-1}^1{p}_{n-1}\left(x\right)p_n^{\prime }\left(x\right)w\left(x\right)dx\end{array}$$

$$\begin{array}{cc}& =a_n\frac{n{\gamma }_n}{{\gamma }_{n-1}}{\int }_{-1}^1{p}_{n-1}^2\left(x\right)w\left(x\right)dx=n,\end{array}$$

where we have used the well known fact that $a_n={\gamma }_{n-1}/{\gamma }_n$  . This proves ( 52 ). Again, from ( 6 ) it is easy to find that $$x^2{p}_n\left(x\right)=a_{n+2}{a}_{n+1}{p}_{n+2}\left(x\right)+(a_{n+1}^2+a_n^2)p_n\left(x\right)+a_n{a}_{n-1}{p}_{n-2},$$  and we get

$$\begin{array}{cc}{\int }_{-1}^1{x}^3& p_n\left(x\right)p_n^{\prime }\left(x\right)w\left(x\right)dx=(a_{n+1}^2+a_n^2){\int }_{-1}^{1}x{p}_n\left(x\right)p_n^{\prime }\left(x\right)w\left(x\right)dx\end{array}$$

$$\begin{array}{cc}& +a_n{a}_{n-1}{\int }_{-1}^{1}x{p}_{n-2}{p}_n^{\prime }\left(x\right)w\left(x\right)dx.\end{array}$$

(54)

First integral in the right hand side was computed in ( 52 ), and it remains to concentrate our attention on the second one. Since

$$\begin{array}{cc}x{p}_n^{\prime }\left(x\right)& =n{\gamma }_n{x}^n+(n-2){\beta }_n{x}^{n-2}+\text{lower degree terms}\end{array}$$

$$\begin{array}{cc}& =n{p}_n\left(x\right)-2\frac{{\beta }_n}{{\gamma }_{n-2}}{p}_{n-2}\left(x\right)+\text{lower degree terms,}\end{array}$$

we have $${\int }_{-1}^{1}x{p}_{n-2}\left(x\right)p_n^{\prime }\left(x\right)w\left(x\right)dx=-2\frac{{\beta }_n}{{\gamma }_{n-2}}.$$  Substituting it in ( 54 ), we obtain ( 53 ).

We find next an expression for the ratio ${\beta }_n/{\gamma }_n$  in terms of the coefficients of the recurrence relation:

Lemma 5.3With the notations introduced in ( 6 ) and ( 51 ),

$$\begin{array}{c}\frac{{\beta }_{n+1}}{{\gamma }_{n+1}}=\frac{{\beta }_n}{{\gamma }_n}-a_n^2,n\in \mathbb{N};\end{array}$$

(55)

in particular

$$\begin{array}{c}\frac{{\beta }_{n+1}}{{\gamma }_{n+1}}=-{\sum }_{k=1}^n{a}_k^2.\end{array}$$

(56)

Proof. ] Comparing the coefficients of $x^{n-1}$  in both sides of ( 6 ) we obtain that $${\beta }_n=a_{n+1}{\beta }_{n+1}+a_n{\gamma }_{n-1},$$  so that $$\frac{{\beta }_{n+1}}{{\gamma }_{n+1}}=\frac{1}{a_{n+1}}\frac{{\gamma }_n}{{\gamma }_{n+1}}\frac{{\beta }_n}{{\gamma }_n}-\frac{a_n}{a_{n+1}}\frac{{\gamma }_{n-1}}{{\gamma }_n}\frac{{\gamma }_n}{{\gamma }_{n+1}},$$  and so, using the identity $$a_n=\frac{{\gamma }_{n-1}}{{\gamma }_n},$$  ( 55 ) holds. Formula ( 56 ) follows from ( 55 ) and the fact that ${\beta }_1=0$  .

Remark 5.4Observe that we have used only the symmetry of the recurrence relation, so these formulas are valid for any even weight function on $[-1,1]$  .

Corollary 5.5For the symmetric Pollaczek polynomials the following asymptotic formula is valid:

$$\begin{array}{c}\frac{{\beta }_n}{{\gamma }_n}=-\frac{n}{4}+\frac{a}{2}logn-\frac{a-\lambda }{4}-\frac{a}{2}\psi (a+\lambda )+\mathcal{O}\left(\frac{1}{n}\right),n\to \infty .\end{array}$$

(57)

Proof. ] From ( 7 ) it is easy to obtain that $$a_k^2=1+\frac{a^2-{\lambda }^2-a+\lambda }{k+\lambda +a-1}+\frac{{\lambda }^2-a^2-\lambda -a}{k+\lambda +a},$$  and ( 56 ) lead us to

$$\begin{array}{cc}\frac{{\beta }_n}{{\gamma }_n}=& -\frac{n-1}{4}+\frac{1}{4}{\sum }_{k=1}^{n-1}\frac{2a}{k+\lambda +a-1}+\frac{1}{4}\frac{-{\lambda }^2+a^2+\lambda +a}{n+\lambda +a-1}\end{array}$$

$$\begin{array}{cc}& -\frac{1}{4}\frac{-{\lambda }^2+a^2+\lambda +a}{\lambda +a}.\end{array}$$

Using that $${\sum }_{k=1}^{n-1}\frac{1}{k+\lambda +a-1}=\psi (n+\lambda +a-1)-\psi (\lambda +a)=log\left(n\right)-\psi (\lambda +a)+\mathcal{O}\left(\frac{1}{n}\right),$$  when $n\to \infty$  , we finally get $$\frac{{\beta }_n}{{\gamma }_n}=-\frac{n-1}{4}+\frac{2a}{4}(log(n)-\psi (\lambda +a\left)\right)+\frac{1}{4}(\lambda -a-1)+\mathcal{O}\left(\frac{1}{n}\right),$$  which is equivalent to the statement of the Lemma.

Now we can prove Theorem  1.3 . Remember that the technique of [14] is not valid here because an additional assumption on $w$  from [14] is not satisfied.

The central idea in our proof is to take advantage of the fact that the main contribution to the asymptotics of $G_n$  comes from the behavior of the weight $w$  at the endpoints of $\Delta$  (see Section  2 ). Using functions $w_0$  and $s$  introduced in ( 22 ), we write $G_n$  in the form

$$\begin{array}{cc}{G}_n& ={\int }_{-1}^{1}log\left(w_0\right(x\left)\right)p_n^2\left(x\right)e^{s\left(x\right)}{w}_0\left(x\right)dx+{\int }_{-1}^{1}s\left(x\right)p_n^2\left(x\right)w\left(x\right)dx.\end{array}$$

In particular, since $s\in C[-1,1]$  (see Lemma  2.1 ), applying ( 42 ) in the second integral we have

$$\begin{array}{c}{G}_n={\int }_{-1}^{1}log\left(w_0\right(x\left)\right)p_n^2\left(x\right)e^{s\left(x\right)}{w}_0\left(x\right)dx+B_1+o\left(1\right),\end{array}$$

(58)

where

$$\begin{array}{cc}{B}_1& ={lim}_{n\to \infty }{\int }_{-1}^{1}s\left(x\right)p_n^2\left(x\right)w\left(x\right)dx={\int }_{-1}^1\frac{s\left(x\right)}{\pi \sqrt{1-x^2}}dx.\end{array}$$

If we denote $$g\left(x\right)=log\left(w_0\right(x\left)\right)=-\frac{2\pi a\left|x\right|}{\sqrt{1-x^2}},$$  taking into account the symmetry, we can rewrite the integral in the right hand side of ( 58 ) as

$$\begin{array}{cc}& {\int }_{-1}^{1}g\left(x\right)p_n^2\left(x\right)e^{s\left(x\right)}{e}^{g\left(x\right)}dx=2{\int }_0^1\frac{g\left(x\right)}{g^{\prime }\left(x\right)}{p}_n^2\left(x\right)e^{s\left(x\right)}{e}^{g\left(x\right)}{g}^{\prime }\left(x\right)dx\end{array}$$

$$\begin{array}{cc}& =2{\int }_0^1\left(\frac{g\left(x\right)}{g^{\prime }\left(x\right)}{p}_n^2\left(x\right)e^{s\left(x\right)}\right)d{e}^{g\left(x\right)}.\end{array}$$

Observe that for $x\in [0,1]$  , $g\left(x\right)/g^{\prime }\left(x\right)=x(1-x^2)$  , so integrating by parts,

$$\begin{array}{cc}& {\int }_{-1}^{1}g\left(x\right)p_n^2\left(x\right)e^{s\left(x\right)}{e}^{g\left(x\right)}dx\end{array}$$

$$\begin{array}{cc}& =2{\left[x(1-x^2)p_n^2\left(x\right)w\left(x\right)\right]}_{x=0}^{x=1}-2{\int }_0^1{\left(x(1-x^2)p_n^2\left(x\right)e^{s\left(x\right)}\right)}^{\prime }{w}_0\left(x\right)dx\end{array}$$

$$\begin{array}{cc}& =-{\int }_{-1}^1(1-3x^2)p_n^2\left(x\right)w\left(x\right)dx-{\int }_{-1}^{1}x(1-x^2){\left(p_n^2\left(x\right)\right)}^{\prime }w\left(x\right)dx\end{array}$$

$$\begin{array}{cc}& -2{\int }_0^{1}x(1-x^2)p_n^2\left(x\right)s^{\prime }\left(x\right)w\left(x\right)dx.\end{array}$$

The asymptotics of each of these three integrals can be computed by means of ( 42 ), Lemma  5.2 , and Lemma  5.1 , respectively, obtaining that

$$\begin{array}{cc}{G}_n& =2n(a_{n+1}^2+a_n^2-1)-4a_n{a}_{n-1}\frac{{\beta }_n}{{\gamma }_{n-2}}+B_1-2B_2+\frac{1}{2}+o\left(1\right)\end{array}$$

$$\begin{array}{cc}& =2n(a_{n+1}^2+a_n^2-1)-4\frac{{\beta }_n}{{\gamma }_n}+B_1-2B_2+\frac{1}{2}+o\left(1\right),n\to \infty .\end{array}$$

Using that $a_n\to 1/2$  (see  7 ) and  57 , we get that

$$\begin{array}{c}{G}_n=-2alog\left(n\right)+B_1-2B_2+\frac{1}{2}+B_3+o\left(1\right),\end{array}$$

(59)

where

$$\begin{array}{c}{B}_3=-2a+(a-\lambda )+2a\psi (a+\lambda )=-a-\lambda +2a\psi (a+\lambda ).\end{array}$$

(60)

Let us simplify the expression of the constant term of this asymptotics.

First, $$B_1-2B_2={\int }_{-1}^1\frac{s\left(x\right)}{\pi \sqrt{1-x^2}}dx-2{\int }_0^{1}x\sqrt{1-x^2}{s}^{\prime }\left(x\right)\frac{1}{\pi }dx,$$  and integrating by parts the second integral,

$$\begin{array}{cc}{B}_1-2B_2& =2{\int }_0^1\frac{s\left(x\right)}{\pi \sqrt{1-x^2}}-\frac{2}{\pi }{\left[x\sqrt{1-x^2}s\left(x\right)\right]}_{x=0}^{x=1}+2{\int }_0^1\frac{1-2x^2}{\pi \sqrt{1-x^2}}s\left(x\right)dx\end{array}$$

$$\begin{array}{cc}& =\frac{4}{\pi }{\int }_0^{1}s\left(x\right)\sqrt{1-x^2}dx.\end{array}$$

(61)

Using the explicit expression for $s$  on $[0,1]$  , the right hand side in ( 61 ) is reduced to

$$\begin{array}{cc}& log\left(\frac{2^{2\lambda }(\lambda +a)}{2\pi \Gamma \left(2\lambda \right)}\right)\frac{4}{\pi }{\int }_0^1\sqrt{1-x^2}dx+\frac{4(\lambda -1/2)}{\pi }{\int }_0^{1}log(1-x^2)\sqrt{1-x^2}dx\end{array}$$

$$\begin{array}{cc}& +\frac{4}{\pi }{\int }_0^{1}2axarccosxdx+\frac{4}{\pi }{\int }_0^{1}log{\left|\Gamma (\lambda +it)\right|}^2\sqrt{1-x^2}dx+4{\int }_0^{1}axdx\end{array}$$

$$\begin{array}{cc}& =log\left(\frac{2^{2\lambda }(\lambda +a)}{2\pi \Gamma \left(2\lambda \right)}\right)+(1-2log(2\left)\right)(\lambda -1/2)\end{array}$$

$$\begin{array}{cc}& +a+\frac{4}{\pi }{\int }_0^{1}log{\left|\Gamma (\lambda +it)\right|}^2\sqrt{1-x^2}dx+2a\end{array}$$

$$\begin{array}{cc}& =log\left(\frac{(\lambda +a)}{\pi \Gamma \left(2\lambda \right)}\right)+\lambda +3a+\frac{4}{\pi }{\int }_0^{1}log{\left|\Gamma \left(\lambda +i\frac{ax}{\sqrt{1-x^2}}\right)\right|}^2\sqrt{1-x^2}dx.\end{array}$$

(62)

Let us compute now the value of this last integral. With the change of variables $u=x/\sqrt{1-x^2}$  we obtain that

$$\begin{array}{cc}{\int }_0^{1}log{\left|\Gamma \left(\lambda +i\frac{ax}{\sqrt{1-x^2}}\right)\right|}^2\sqrt{1-x^2}dx& ={\int }_0^{+\infty }log{\left|\Gamma \left(\lambda +iau\right)\right|}^2\frac{du}{(1+u^2{)}^2}\end{array}$$

$$\begin{array}{cc}& ={\int }_{-\infty }^{+\infty }log\left|\Gamma \left(\lambda +iau\right)\right|\frac{du}{(1+u^2{)}^2}.\end{array}$$

With $\lambda >0$  , $a\ge 0$  , function $$f\left(u\right)\stackrel{\text{def}}{=}\frac{log\left(\Gamma \left(\lambda +iau\right)\right)}{(1+u^2{)}^2}$$  is meromorphic and single valued in the lower half plane $\{Im(u)<0\}$  , with a double pole at $u=-i$  . Taking into account that by Stirling formula, $$log\left(\Gamma \left(\lambda +iau\right)\right)\sim (\lambda +iau)\left|u\right|log(\lambda +iau)\text{as}u\to \infty ,Im\left(u\right)<0,$$  we may apply the residue calculus to establish that $${\int }_{-\infty }^{+\infty }log\left(\Gamma \left(\lambda +iau\right)\right)\frac{du}{(1+u^2{)}^2}=-2\pi ire{s}_{u=-i}f\left(u\right)=\frac{\pi }{2}\left(a\psi (\lambda +a)-log(\Gamma (\lambda +a\left)\right)\right).$$  Taking the real part, we get that

$$\begin{array}{c}\frac{4}{\pi }{\int }_0^{1}log{\left|\Gamma \left(\lambda +i\frac{ax}{\sqrt{1-x^2}}\right)\right|}^2\sqrt{1-x^2}dx=2a\psi (\lambda +a)-2log(\Gamma (\lambda +a\left)\right).\end{array}$$

(63)

Gathering  60 – 63 in  59 we conclude the proof of Theorem  1.3 .

Remark 5.6The idea of this proof can be applied also to the case of a non-symmetric weight of the form $$w\left(x\right)=exp\left\{-\frac{4c}{(1-x{)}^{\alpha }}-\frac{4d}{(1+x{)}^{\alpha }}+s\left(x\right)\right\},$$  where $s\in C^1[-1,1]$  , and $\alpha \in [1/2,1]$  .

6 Proof of corollaries  1.6 and  1.8 

Proof of Corollary  1.6 : Consider the integral $$I_n=-\frac{2}{\pi }{\int }_{{\alpha }_{-n}}^{{\alpha }_n}\frac{Q\left(x\right)}{\sqrt{({\alpha }_n-x)(x-{\alpha }_{-n})}}dx=-\frac{4}{\pi }{\int }_0^{{\alpha }_n}\frac{Q\left(x\right)}{\sqrt{({\alpha }_n-x)(x-{\alpha }_{-n})}}dx.$$  By ( 22 ),

$$\begin{array}{cc}{I}_n=& -\frac{4}{\pi }{\int }_0^{{\alpha }_n}\left(\frac{\pi a\left|x\right|}{\sqrt{1-x^2}}+\frac{1}{2}log\left(w\right(0\left)\right)-\frac{s\left(x\right)}{2}\right)\frac{1}{\sqrt{{\alpha }_n^2-x^2}}dx\end{array}$$

$$\begin{array}{cc}=& -\frac{4}{\pi }{\int }_0^1\left(\frac{\pi a\left|{\alpha }_{n}x\right|}{\sqrt{1-{\alpha }_n^2{x}^2}}+\frac{1}{2}log\left(w\right(0\left)\right)-\frac{s\left({\alpha }_{n}x\right)}{2}\right)\frac{1}{\sqrt{1-x^2}}dx\end{array}$$

$$\begin{array}{cc}=& -\frac{4}{\pi }{\int }_0^1\frac{\pi a\left|{\alpha }_{n}x\right|}{\sqrt{1-{\alpha }_n^2{x}^2}}\frac{1}{\sqrt{1-x^2}}dx\end{array}$$

$$\begin{array}{cc}& -\frac{4}{\pi }{\int }_0^1\left(\frac{1}{2}log\left(w\right(0\left)\right)-\frac{s\left({\alpha }_{n}x\right)}{2}\right)\frac{1}{\sqrt{1-x^2}}dx.\end{array}$$

Since $s\in C[-1,1]$  , the second integral is bounded; hence, by  35 , $$I_n=-4a{\alpha }_n\frac{-1}{2}log(1-{\alpha }_n)+\mathcal{O}\left(1\right)=2alog(1-{\alpha }_n)+\mathcal{O}\left(1\right),n\to \infty .$$  Finally, from the asymptotics of ${\alpha }_n$  found in Proposition  3.2 we obtain $$I_n=2alog(a/n)+\mathcal{O}\left(1\right)=-2alog\left(n\right)+\mathcal{O}\left(1\right),n\to \infty ,$$  and comparing this expression with the result of Corollary  1.4 , the statement follows.

Proof of Corollary  1.8 : Taking into account the relation between the entropy and the mutual energy ( 16 ) we obtain that $$I[{\rho }_n,{\nu }_n]=\frac{E_n+2log\left({\gamma }_n\right)}{2n}.$$  By  17 , $$2log\left({\gamma }_n\right)=2nlog\left(2\right)+2alog\left(n\right)-log\left(\frac{\Gamma (\lambda +a+1)\Gamma (\lambda +a)}{\Gamma \left(2\lambda \right)}\right)+o\left(1\right),$$  so that by  9 , $$I[{\nu }_n,{\lambda }_n]=log\left(2\right)+\frac{1}{2n}\left(\tau (\lambda ,a)-log\left(\frac{\Gamma (\lambda +a+1)\Gamma (\lambda +a)}{\Gamma \left(2\lambda \right)}\right)\right)+o\left(\frac{1}{n}\right).$$  The use of the explicit expression for $\tau (\lambda ,a)$  in  10 concludes the proof of the Corollary. Acknowledgement The authors were supported, in part, by a research grant from the Ministry of Science and Technology (MCYT) of Spain, project code BFM2001-3878-C02, and by Junta de Andalucıa, Grupo de Investigacion FQM229. The research of A.M.F. was supported also by Research Network “Network on Constructive Complex Approximation (NeCCA)”, INTAS 03-51-6637, and by NATO Collaborative Linkage Grant “Orthogonal Polynomials: Theory, Applications and Generalizations”, ref. PST.CLG.979738. We are indebted also to Prof. D. S. Lubinsky for very interesting discussions.

References

1. M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, Dover Publ., New York, 1972.
2. A. I. Aptekarev, V. S. Buyarov, and J. S. Dehesa, Asymptotic behavior of the $L^p$  -norms and the entropy for general orthogonal polynomials, Russian Acad. Sci. Sb. Math., 82 (1995), pp. 373–395.
3. B. Beckermann, A. Martßnez-Finkelshtein, E. A. Rakhmanov, and F. Wielonsky, Asymptotic upper bounds for the entropy of orthogonal polynomials in the Szego class, J. Math. Physics, 45 (2004), pp. 4239–4254.
4. I. Bialynicki-Birula and J. Mycielsky, Uncertainty relations for information entropy in wave mechanics, Commun. Math. Phys., 44 (1975), pp. 129–132.
5. V. S. Buyarov, On information entropy of Gengenbauer polynomials, Vesntik Moskow Univ., Ser. 1 6 (1997), pp. 8–11 (in russian).
6. V. S. Buyarov, J. S. Dehesa, A. Martßnez-Finkelshtein, and E. B. Saœ, Asymptotics of the information entropy for Jacobi and Laguerre polynomials with varying weights, J. Approx. Theory, 99 (1999), pp. 153–166.
7. V. S. Buyarov, P. Lopez-Artes, A. Martßnez-Finkelshtein, and W. Van Assche, Information entropy of Gegenbauer polynomials, J. Phys. A, 33 (2000), pp. 6549–6560.
8. V. Buyarov, J. S. Dehesa, A. Martßnez-Finkelshtein, and J. Sanchez-Lara, Computation of the entropy of polynomials orthogonal on an interval, SIAM J. Sci. Comp., 26 (2004), pp. 488–509.
9. J. S. Dehesa, W. Van Assche, R. J. Yanez, Information entropy of classical orthogonal polynomials and their application to the harmonic oscillator and Coulomb potentials, Methods Appl. Anal. 4 (1997), pp. 91–110.
10. J. S. Dehesa, A. Martßnez-Finkelshtein, and J. Sanchez-Ruiz, Quantum information entropies and orthogonal polynomials, J. Comput. Appl. Math., 133 (2001), pp. 23–46.
11. J. S. Dehesa, R. J. Yanez, A. I. Aptekarev, and V. S. Buyarov, Strong asymptotics of Laguerre polynomials and information entropies of 2D harmonic oscillator and 1D Coulomb potentials, J. Math. Physics, 39 (1998), pp. 3050–3060.
12. I. S. Gradshtein and I. M. Ryzhik, Table of Integrals, Series and Products, Academic Press, San Diego, CA, fifth ed., 1995.
13. E. Levin and D. S. Lubinsky, Orthogonal Polynomials for Exponential Weights, vol. 4 of CMS Books in Mathematics, Springer Verlag, 2001.
14. E. Levin and D. S. Lubinsky, Asymptotics for entropy integrals associated with exponential weights, J. Comput. Appl. Math., 156 (2003), pp. 265–283.
15. Y. L. Luke, Mathematical Functions and their Approximations, Academic Press, New York, 1975.
16. A. Mate, P. Nevai, V. Totik, Strong and weak convergence of orthogonal polynomials, Am. J. Math. 109 (1987), pp. 239–282.
17. R. G. Parr and W. Yang, Density Functional Theory of Atoms and Molecules, Oxford University Press, New York, 1989.
18. E. A. Rakhmanov, On the asymptotics of the ratio of orthogonal polynomials, Math. USSR Sb., 32 (1977), pp. 199–213.
19. E. B. Saœ and V. Totik, Logarithmic Potentials with External Fields, vol. 316 of Grundlehren der Mathematischen Wissenschaften, Springer-Verlag, Berlin, 1997.
20. J. F. Sanchez Lara, On the asymptotic expansion of the entropy of Gegenbauer polynomials, J. Comput. Appl. Math., 142 (2002), pp. 401–409.
21. G. Szego, Orthogonal Polynomials, vol. 23 of Amer. Math. Soc. Colloq. Publ., Amer. Math. Soc., Providence, RI, fourth ed., 1975.
22. R. J. Yanez, W. Van Assche, and J. S. Dehesa, Position and momentum information entropies of the D-dimensional harmonic oscillator and hydrogen atom, Physical Rev. A, 50 (1994), pp. 3065–3079.
23. R. J. Yanez, W. Van Assche, R. Gonzalez-Ferez, J. S. Dehesa, Entropic integrals of hyperspherical harmonics and spatial entropy of $D$  -dimensional central potentials, J. Math. Phys. 40 (1999), pp. 5675–5686.