                                                DENSITY OF EIGENVALUES AND ITS PERTURBATION
                                                INVARIANCE IN UNITARY ENSEMBLES OF RANDOM
                                                                 MATRICES




arXiv:0810.5425v1 [math-ph] 30 Oct 2008
                                                          DANG-ZHENG LIU, ZHENG-DONG WANG AND KUI-HUA YAN


                                                    Abstract. We generally study the density of eigenvalues in unitary ensembles
                                                    of random matrices from the recurrence coefficients with regularly varying
                                                    conditions for the orthogonal polynomials. First we calculate directly the
                                                    moments of the density. Then, by studying some deformation of the moments,
                                                    we get a family of differential equations of first order which the densities satisfy
                                                    (see Theorem 1.2), and give the densities by solving them. Further, we prove
                                                    that the density is invariant after the polynomial perturbation of the weight
                                                    function (see Theorem 1.5).




                                                                        1. Introduction and main results
                                             Statistical behavior of eigenvalues of random matrices was first studied by E.P.
                                          Wigner in 1950s to get an information about spectra of heavy nuclei, and the famous
                                          semicircle-law was found then. In 1960s random matrices were intensively developed
                                          by E.P. Wigner, F.J. Dyson, M.L. Mehta and others for a better understanding of
                                          statistical behavior of energy levels in nuclear Physics, see [Po] as a collection of
                                          early papers. Later more and more importance was gained in other areas of Physics
                                          and Mathematics (see [Me] for a comprehensive introduction, or J. Phys. A 36,
                                          2003).
                                             It is of special importance to observe that the eigenvalue distribution of unitary
                                          invariant ensembles which can be described in [FK] by weight functions from the
                                          families of classical orthogonal polynomials. Such ensembles can be defined by the
                                          probability distribution on a space HN of N -order Hermitian matrices as

                                                                                            N
                                                                                            Y           Y
                                                                        P (dA) = F (A)           dajj         dℜajk dℑajk                  (1.1)
                                                                                           j=1          j<k

                                          where F is an integrable nonnegative class function on HN .
                                            Using Proposition 2 in [HZ] we obtain

                                          Z
                                               F (A)dA =
                                          HN
                                                   Z +∞         Z +∞                                  Y
                                              ZN          ···          F (diag(x1 , . . . , xN ))              (xj − xk )2 dx1 . . . dxN   (1.2)
                                                   −∞           −∞                                  1≤j<k≤N



                                              2000 Mathematics Subject Classification. Primary 15A52; Secondary 60F99, 30E05.
                                                                                                 1
2                DANG-ZHENG LIU, ZHENG-DONG WANG AND KUI-HUA YAN


where
                                                   π k(k−1)/2
                                     ZN =                           .                                     (1.3)
                                                k!(k − 1)! · · · 1!
Especially when
                                                                    N
                                                                    Y
                                                     −1
                       F (diag(x1 , . . . , xN )) = ZN                    ω(xj )χI (xj ),                 (1.4)
                                                                    j=1

one introduces the joint distribution function for the eigenvalues

                                                  N
                                                  Y                     Y
                  PN (x1 , x2 , . . . , xN ) =             ω(xj )                (xj − xk )2              (1.5)
                                                  j=1               1≤j<k≤N

where I = (α, β), −∞ ≤ α < β ≤ +∞ and ω(x) is a weight function on the interval
I all finite moments of which exist.

Remark 1.1 (1.5) was imposed on a definition of a matrix ensemble in [Me], and
obtained in [TW] for
                                                  
                         F (A) = exp − Tr(V (A))                          (1.6)
where V (x) is a real-valued function such that ω(x) = exp(−V (x)) defines a weight
function. In addition, it is easy to draw the conclusion similar to (1.5) for the
orthogonal and symplectic ensembles.

Remark 1.2 The appearance of the characteristic function χI of the set I means
that we just consider sub-ensemble of Hermitian matrices space whose eigenvalues
are all in I. Moreover, I can be given by the union of some disjoint intervals and
then ω(x) is the associated weight function.
   From (1.5) the n−point correlation function is defined in [Me] by
                                        Z         Z
                               N!
     Rn (x1 , . . . , xn ) =                ···           PN (x1 , x2 , . . . , xN ) dxn+1 . . . dxN .    (1.7)
                             (N − n)!   I             I
    If we introduce the orthogonal polynomials
                     Z
                        pj (x)pk (x)ω(x)dx = δjk , j, k = 0, 1, . . .                                     (1.8)
                          I
and the associated functions
                                                    p
                                     ϕk (x) = pk (x) ω(x),                                                (1.9)
then by the property of Vandermonde determinant the joint distribution function
(1.5) reads:

                                              1                         2
                 PN (x1 , x2 , . . . , xN ) =     det[ϕj−1 (xk )]|N
                                                                  j,k=1
                                              N!
                                              1
                                            =    det[KN (xj , xk )]|N
                                                                    j,k=1                                (1.10)
                                              N!
where
                                                          N
                                                          X −1
                                 KN (x, y) =                     ϕj (x)ϕj (y).                           (1.11)
                                                          j=0
                 DENSITY OF EIGENVALUES IN UNITARY ENSEMBLES                         3


Because of the orthonormality of the ϕj (x)’s one can show just as in [Me, Chap.6]
that
                     Rn (x1 , . . . , xn ) = det[KN (xj , xk )]|nj,k=1 .    (1.12)
   It is of considerable interest to obtain the behavior of n−point correlation of
(1.12) after some appropriate scaling in the limit of large N . In particular, putting
n = 1, we get the density of eigenvalues (also called level density )
                                                 N
                                                 X −1
                          R1 (x) = KN (x, x) =          ϕ2j (x)                (1.13)
                                                 j=0

and denote the normalized density by
                                           N −1
                              1          1 X 2
                                R1 (x) =        ϕ (x).                         (1.14)
                              N          N j=0 j

   Remark 1.3 For the eigenvalues x1 , x2 , . . . , xN of a random Hermitian matrix
A, the integral of the normalized counting function over an interval ∆ = (c, d) is
given by
                                     Z X N
                                   1
                         νN (∆) =              δ(x − xj )dx.                 (1.15)
                                   N ∆ j=1
A calculation in [PS] shows that
                                          Z
                                             1
                            E(νN )(∆) =        R1 (x)dx                        (1.16)
                                           ∆ N
where E(νN )(∆) denotes the expectation of νN (∆) with respect to the probability
measure of (1.5).
   Our motivation is to obtain the density in the scaling limit of large N from
(1.14) and further study its polynomial-perturbation invariance. Before our results
are stated we review some known results about the density.
   As we know, Wigner in [Wig1, 2] not only got his famous semicircle law (here
                                         2
corresponding to the weight ω(x) = e−x called Gauss unitary ensemble, denoting
GUE)
                                   2p
                            σ(x) =     1 − x2 χ[−1,1] (x),                     (1.17)
                                   π
but also invented the calculation method which had some independent interest as
he thought. Afterwards, Bronk got the associated density for ω(x) = xρ e−x , ρ >
−1, 0 < x < +∞ for Laguerre ensembles in [Br] and Leff for ω(x) = (1 − x)a (1 +
x)b , a, b > −1, −1 < x < 1 for Jacobian ensembles in [Lef]. Recently, an irradiative
idea in [HT] was introduced by Haagerup and Thorbjornsen to give a short proof of
Wigner’s semicircle law, using the Laplace transform of (1.14) and the property of
Hermitian polynomials. In [Lef], Ledoux pushed forward the investigation in [HT]
and in a slightly different way obtained Wigner’s semicircle law, also the densities
for Laguerre and Jacobian unitary ensembles but whose expressions appeared com-
plicated. In addition, based on the spirit of statistical mechanics the method of
equilibrium measure is used to obtain the density, mainly for the weight function
ω(x) = exp(−V (x)),
                       V (x) = γ2m x2m + · · · + γ0 , γ2m > 0.                 (1.18)
4             DANG-ZHENG LIU, ZHENG-DONG WANG AND KUI-HUA YAN


See [Jo], [DMK] or [De] for the equilibrium measure method. In a recent survey on
orthogonal polynomial ensembles ([Kö]), König obtained semicircle law for GUE
respectively by moment method and equilibrium measure method.
   In this paper, we will generally deal with the weight function ω(x) on the inter-
val I and obtain the density of eigenvalues σ(x). In fact, Professor W. Van Assche
told us that P. Nevai and W. Van Assche had shown that the density and zero
distribution of the orthogonal polynomials were the same ([Va1], Theorem 5.3), see
Remark 1.8 below. However, our method is to calculate directly the moments with-
out other knowledge than three-term recurrence formula , and further the density
was obtained using a different method. To our knowledge, the family of densities
satisfying differential equations of first order is first obtained in the present paper.
   It is a well known result that three-term recurrence formula holds for the orthog-
onal polynomials defined by (1.8)
         xpn (x) = an+1 pn+1 (x) + bn pn (x) + an pn−1 (x), n = 0, 1, · · ·      (1.19)
where an > 0, p−1 (x) = 0. We will assume that there is a positive and non-
decreasing sequence cn such that
                                an              bn
                            lim    = a > 0, lim    = b.                   (1.20)
                        n→∞ cn              n→∞ cn

An extra condition on the contraction sequence is that cn is a regularly varying
sequence with index λ ≥ 0, i.e.
                                         cn = nλ L(n)                            (1.21)
where L : (0, +∞) −→ (0, +∞) is slowly varying , that is,
                                   L(xt)
                               lim       = 1,       ∀t > 0.                      (1.22)
                               x→∞ L(x)

The condition of (1.21) was first introduced by W. Van Assche to study the asymp-
totics for orthogonal polynomials, see [Va2] as a general survey for the condition of
(1.20) and [Fe, Chap.VIII] for regular functions (1.22).

  Remark 1.4 Assuming c(x) is a positive, non-decreasing differentiable function
on the interval (0, +∞) with
                                            ′
                                      xc (x)
                                   lim       = λ ≥ 0,                            (1.23)
                                  x→+∞ c(x)

then one easily knows that c(x) can be represented as
                                     c(x) = xλ L(x).                             (1.24)
A discrete version of this result is that (1.21) holds if
                                    c         
                                       n+1
                              lim n         − 1 = λ ≥ 0.                         (1.25)
                             n→∞       cn

   Remark 1.5 The conditions of (1.20) and (1.21) are connected close with the
asymptotic problems of the orthogonal polynomials. The class of Freud weights
plays a most close role on (1.21), for examples, ω(x) = exp(−Q(x)) where Q(x)
grows like a power at infinity, in particular,
                         Q(x) = γ2m x2m + · · · + γ0 , γ2m > 0                   (1.26)
                 DENSITY OF EIGENVALUES IN UNITARY ENSEMBLES                         5


with λ = 1/(2m), see [DKMVZ] ;
                                               α
                         ω(x) = |x|β e−|x| , β > −1, α > 0                      (1.27)
with λ = 1/α, see [LMS] or [Va1]. In addition, λ = 1 for Laguerre weights
                                ω(x) = xα e−x , α > −1                          (1.28)
and λ = 0 for Jacobi weights
                        ω(x) = (1 − x)α (1 + x)β , α, β > −1.                   (1.29)
A classic result of E.A. Rakhmanov in [Si] asserts that if the weight ω on [−1, 1]
satisfying ω > 0 a.e. on [−1, 1], then ω(x) belongs to the Nevai-Blumenthal class,
that is
                            lim an = 1/2, lim bn = 0.                       (1.30)
                              n→∞                   n→∞
Obviously, λ = 0 for the Nevai-Blumenthal class (see [Va2]). We strongly refer the
reader to [Lub] for a recent survey for a wide variety of weights on finite or infinite
intervals.
   Now we can state our main results. In Section 2 we will introduce ascending,
equilibrating and descending operators which describe the transforming of polyno-
mials, and explicitly calculate the moments of the density. Then in Section 3 we
consider a simple deformation of the moments for any given density and obtain a
corresponding density determined by a differential equation with respect to the new
moments. In section 4 using the results in Section 3 we give the proofs of Theorems
1.2 and 1.5 below.
   First let us rescale the density of (1.14) by
                                           cN
                                σN (x) =      R1 (cN x).                        (1.31)
                                           N
Note that σN (x) is our main object and we will study its limit behavior.

                                                                                  (N )
Theorem 1.1 Denote the kth moment of the scaling density σN (x) by Mk                    .
Under the contraction conditions of (1.20) and (1.21), we have

                                    (N )
                             lim Mk        = Mk ,     k = 0, 1, . . .           (1.32)
                         N →∞
where
                                       [k/2]                  
                                   1 X j j
                       Mk =                  Ck Ck−j a2j bk−2j .                (1.33)
                                1 + λk j=0
  Remark 1.6 One sets a = 1/2 in the following since cn can be chosen freely from
some constant. Observe
                     [k/2]
                      X                              a    k
                          Ckj Ck−j
                               j
                                   a2j bk−2j = L0 az + + b                      (1.34)
                      j=0
                                                      z

where the operator L0 (f ) represents the constant term of Laurent series, using
Cauchy contour integral, that is
                                           I
                                        1    f (z)
                             L0 (f ) =             dz.                    (1.35)
                                       2πi     z
6              DANG-ZHENG LIU, ZHENG-DONG WANG AND KUI-HUA YAN


Theorem 1.2 A probability density σ(x) with its kth moment Mk exists, and is
uniquely determined by the following differential equation of first order
                                 (1)    1        1
                   σ(x) − λ xσ(x)      = p                 χI ,           (1.36)
                                          π 1 − (x − b)2 b
with the following conditions
                                         Z +∞
                             σ(x) ≥ 0,          σ(x)dx = 1                    (1.37)
                                         −∞
where Ib = (−1 + b, 1 + b) and χIb is a characteristic function of Ib . Exactly, the
support of σ(x) can be restricted to a finite interval, that is, for λ = 0
                                     1         1
                             σ(x) = p                                         (1.38)
                                     π 1 − (x − b)2
while for λ > 0
                                 supp(σ) = [B1 , B2 ]                         (1.39)
where
                      B1 = min{b − 1, 0}, B2 = max{b + 1, 0}.                 (1.40)

    One directly solves the equation of (1.36) and easily obtains

Corollary 1.3 For b = 0 and m ∈ N, we have
  (1) for λ = 1/2m
                                  m
                        4p       X              j−1 
                 σ(x) =   1 − x2     (2x)2m−2j C2j−2   m
                                                     /C2m                     (1.41)
                        π        j=1

and
  (2) for λ = 1/(2m − 1)
             m 
                            √        √
                                       1 − x2 X 1 x 2m−2−2j 
                                              m−1
         mC2m    x 2m−2 1 + 1 − x2
 σ(x) =         ( )      ln        +                j
                                                       ( )    . (1.42)
           2π     2         |x|          2    j=1 jC2j
                                                        2

    For b = −1 and q ∈ N ∪ {0}, we have
    (3) for λ = 1/(q + 1)
                                                     r
                           q + 1 x q X Cqj             2 − x 1+2j
                                          q
                    σ(x) =                                                    (1.43)
                             π   2 j=0 1 + 2j              x
and
  (4) for λ = 1/(q + 12 )
                             r       r      r
         2q + 1  x q− 12                     2 − x X 1 x q−j  q
                                                    q
                               2−x     2
 σ(x) =          ( )      ln       +      +             j (8)  C2q . (1.44)
           4π     8             x      x        x j=1 jC2j

   Remark 1.7 For the weight ω(x) = exp(−V (x)) from (1.18) which   p is correspond-
ing with λ = 1/2m, the density was given in [Jo] in the form of r(x) (x − x1 )(x2 − x).
Here x1 < x2 and r(x) is a polynomial of degree 2m − 2 depending on V (x); when
λ = 1 and b = 0, it was obtained in [CI] using Mathematica for Meixner-Pollaczek
polynomials. It is just the case where λ = 1 and b = −1 for Laguerre weights of
(1.28), see [Sze]. The authors believe that the densities in Corollary 1.3, especially
for small integrals m and q, may appear in some other random matrix models.
                 DENSITY OF EIGENVALUES IN UNITARY ENSEMBLES                                       7


  From Theorems 1.1 and 1.2, it is obvious that

Theorem 1.4 Denote σN (x) and σ(x) as above , then
                                                   W
                                         σN (x) −→ σ(x)                                       (1.45)
where W means in the weak sense.

  Remark 1.8 Let x1,n , x2,n , · · · , xn,n be zeros of pn . P. G. Nevai and J. S. De-
hesa in [ND] also got the same moments of zero distribution under the contraction
condition of (1.20) (this fact was pointed out by W. Van Assche in [Va2]), namely,
                                            n
                                       1 X xj,n k
                                   lim       (    ) = Mk .                                    (1.46)
                                   n→∞ n       cn
                                         j=1

A beautiful probabilistic interpretation in [Va1] shows Mk is the kth moment of
Nevai-Ullman measure. In addition, Nevai and Van Assche proved that weak con-
vergence of zero distribution implied some weak convergence to the same probability
                                                         NX−1
measure for Christoffel functions, which are defined by (     p2j (x))−1 (see Theorem
                                                                       j=0
5.3 in [Va1]). It is obvious that the density of eigenvalues and zero distribution are
the same according to (1.46) and Theorem 1.1, thus we give a new proof.

  Let p(x) be a fixed lth order polynomial, and we consider a new weight function
                                      ω̂(x) = p2 (x)ω(x)                                      (1.47)
with                Z
                         p̂j (x)p̂k (x)ω̂(x)dx = δjk , j, k = 0, 1, . . . .                   (1.48)
                     I
Associated 1-point correlation function is given by
                                           N
                                           X −1
                              R̂1 (x) =           p̂2 (x)p2 (x)ω(x).                          (1.49)
                                           j=0

Under the same scaling we write
                                                  cN
                                    σ̂N (x) =        R̂1 (cN x).                              (1.50)
                                                  N
  Calculate the moments of (1.50) for the new weight function, we find

                                                                                       (N )     (N )
Theorem 1.5 Denote the kth moments of σ̂N (x) and σN (x) by Mk and M̂k                                 ,
respectively. Under the contraction conditions of (1.20) and (1.21), we have

                            (N )                (N )
                  lim M̂k          = lim Mk            = Mk ,      k = 0, 1, . . . .          (1.51)
                 n→∞                 n→∞

Namely,
                                                   W
                                         σ̂N (x) −→ σ(x)                                      (1.52)
where σ(x) is the density determined by (1.36) and (1.37).
8              DANG-ZHENG LIU, ZHENG-DONG WANG AND KUI-HUA YAN


                           2. Calculation of the moments
  By the recursion formula of (1.19), we regard the multiplication by x as an
operator Ax , and it can be represented as
                                         Ax = A+ + A0 + A−                                   (2.1)
where A+ , A0 and A− are called ascending, equilibrating and descending operators
respectively, defined by
    A+ pn (x) = an+1 pn+1 (x),           A0 pn (x) = bn pn (x),   A− pn (x) = an pn−1 (x).   (2.2)
Thus we calculate the kth moment of σN (x) using (1.31), (2.1) and (2.2) as follows:
                           Z
                (N )
              Mk       =       xk σN (x)dx
                                  N
                                  X −1 Z
                            1
                       =                 xk p2j ω(x)dx
                         N (cN )k j=0 I

                              1
                                    N
                                    X −1 D             E
                                           k
                       =                  x  p j , p j
                           N (cN )k j=0                  L2 (ω)


                              1
                                    N
                                    X −1 D                           E
                                                         k
                       =                  (A+ + A0 + A− )  p j , p j          .              (2.3)
                           N (cN )k j=0                                L2 (ω)


   Let Λqk be a set composed of those terms in the expansion of (A+ + A0 + A− )k ,
in which the operators A+ S   and A− exactly appear q times respectively. Note that
 T pj , pj L2 (ω) = 0 for T 6∈ Λqk , then one obtains the following lemma:
                                    q


Lemma 2.1
                                            N
                                            X −1 [k/2]
                                                  X X D            E
                        (N )          1
                   Mk          =          k
                                                         T pj , pj         .                 (2.4)
                                   N (cN ) j=0 q=0     q            L2 (ω)
                                                        T ∈Λk

   For convenience of the following calculation, we introduce one property of the
regular varying functions.

Lemma 2.2 Let c(x) be a positive and non-decreasing function on (0, +∞), and

                                              c(x) = xλ L(x)                                 (2.5)
where L(x) satisfies
                                         L(xt)
                                        lim    = 1,         ∀t > 0.                          (2.6)
                                    x→+∞ L(x)
Then
                                          1
                                              Rn
                                          n   1
                                                (c(x))k dx      1
                                   lim                     =        .                        (2.7)
                                   n→∞         (c(n))k       1 + λk

Proof By Lebesgue’s dominated theorem and exchanging limits and integrals , it
is easy to prove using (2.5) and (2.6). 
                  DENSITY OF EIGENVALUES IN UNITARY ENSEMBLES                                        9


  Now we give a proof of Theorem 1.1 by two cases of b.

Case 1 b 6= 0.

  Assume that b > 0 for convenience. Write
                                an               bn
                                   = a(1 + ξn ),    = b(1 + ηn ).                                (2.8)
                                cn               cn
If one writes for j > k
                                  uj =        max       {|ξm |, |ηm |},                          (2.9)
                                          j−k≤m≤j+k

then
                                               lim uj = 0.                                      (2.10)
                                              j→∞

Thus we can suppose uj < 1, ∀j = 0, 1, . . . and by the definition of uj one obtains
       max      {am } ≤ a cj+k (1 + uj ),               max         {bm } ≤ b cj+k (1 + uj ).   (2.11)
   j−k≤m≤j+k                                        j−k≤m≤j+k

Furthermore, for T ∈ Λqk we have
   D          E                    2q                                            k−2q
    T pj , pj     ≤       max {am }                             max         {bm }
             L2 (ω)            j−k≤m≤j+k                    j−k≤m≤j+k
                              2q k−2q
                        ≤ a b            (cj+k ) (1 + uj )k .
                                                k
                                                                                                (2.12)
Summing by q and j, we get

                                                             1 P
                                                               N −1
                                                                          k         k
                             [k/2]
                               X                            N j=0 (cj+k ) (1 + uj )
                 (N )
             Mk         ≤             Ckq Ck−q
                                           q
                                               a2q bk−2q                                    .   (2.13)
                              q=0
                                                                          (cN )k

By Cauchy-Maclaurin summation formula and Lemma 2.2 one obtains

                     NP
                      −1                                            NP
                                                                     −1
                 1                                              1
                 N          (cj+k )k (1 + uj )k                 N         (cj+k )k
                     j=0                                            j=0
          lim                                       =   lim
         N →∞                 (cN )k                    N →∞       (cN )k
                                                                1 P
                                                                  N
                                                                N    (cj )k
                                                                    j=1
                                                    =   lim
                                                        N →∞   (cN )k
                                                            1
                                                               RN        k
                                                            N 1 (c(x)) dx
                                                    = lim
                                                      N →∞      (c(N ))k
                                                         1
                                                    =        .                                  (2.14)
                                                      1 + λk

  Analogously, for j > k, one obtains
                                                                NP
                                                                 −1
                                                            1
                             [k/2]                        N         (cj−k )k (1 − uj )k
                 (N )
                               X                                j=k
                Mk ≥                  Ckq Ck−q
                                           q
                                               a2q bk−2q                                        (2.15)
                               q=0
                                                                          (cN )k
10               DANG-ZHENG LIU, ZHENG-DONG WANG AND KUI-HUA YAN


and
                                         NP
                                          −1
                                    1
                                    N          (cj−k )k (1 − uj )k
                                         j=k            1
                             lim                             .          =                          (2.16)
                      N →∞          (cN                  )k
                                                     1 + λk
     Combining (2.13) — (2.16) we complete the proof of this case. 

Case 2 b = 0.

     Similarly, write
                                     an               bn
                                        = a(1 + ξn ),    = ηn                                      (2.17)
                                     cn               cn
and
                                    vj =         max          {|ξm |, |ηm |},                      (2.18)
                                           j−k≤m≤j+k
then
                                                 lim vj = 0.                                       (2.19)
                                                 j→∞
Thus we can suppose vj < 1, ∀j = 0, 1, . . .. For 0 ≤ j − k ≤ m ≤ j + k, one obtains
                 a cj−k (1 − vj ) ≤ am ≤ a cj+k (1 + vj ), |bm | ≤ cj+k vj .                       (2.20)

                                                  (N )
  Now we give an estimation of Mk when k is odd. For T ∈ Λqk , by (2.20) we
have
            D           E
              T pj , pj     ≤ a2q vjk−2q (cj+k )k (1 + vj )k−2q
                               L2 (ω)

                                           ≤ a2q vjk−2q (cj+k )k (1 + vj )k .                      (2.21)
If denoting
                                              [k/2]
                                               X
                                     vj =             Ckq Ck−q
                                                           q
                                                               a2q vjk−2q ,                        (2.22)
                                               q=0
then
                                                 lim v j = 0.                                      (2.23)
                                                 j→∞
Thus by the Stoltz formula one obtains
                                              NP
                                               −1
                                          1
                                          N           (cj+k )k (1 + vj )k v j
                             (N )              j=0
                            Mk       ≤                                          −→ 0.              (2.24)
                                                         (cN )k

     When k is even. For T ∈ Λqk , q 6= k/2,
                     D          E
                       T pj , pj 2     ≤ a2q vjk−2q (cj+k )k (1 + vj )k                            (2.25)
                                         L (ω)
                 k/2
and for T ∈ Λk , j > k,
                                       D         E
               ak (cj−k )k (1 − vj )k ≤ T pj , pj 2                   ≤ ak (cj+k )k (1 + vj )k .   (2.26)
                                                              L (ω)

If one writes
                        (N )          1
                                            N
                                            X −1 X X D             E
                       Mk      =                       T p j , p j                                 (2.27)
                                   N (cN )k j=0     q               L2 (ω)
                                                         q<k/2 T ∈Λk
                   DENSITY OF EIGENVALUES IN UNITARY ENSEMBLES                            11


and
                                  1
                                        N
                                        X −1 X D              E
                     f(N ) =
                     M                            T p j , p j         ,              (2.28)
                       k
                               N (cN )k j=0   k/2
                                                               L2 (ω)
                                                T ∈Λk
then
                                (N )    (N )
                             Mk = M k + M     f(N ) .                                (2.29)
                                                 k
  Similar to the case where k is odd, by (2.25) one obtains
                                                (N )
                                        lim M k        =0                            (2.30)
                                       N →∞
and similar to the case where b 6= 0, by (2.26) one obtains
                                 f(N ) =     1     k/2
                           lim M   k             C ak .                   (2.31)
                          N →∞            1 + λk k
  Combining (2.24), (2.29), (2.30) and (2.31) we complete the proof of this case
where b = 0. 
  Therefore, Theorem 1.1 has been proved.

                         3. Deformation of the moments
   Following Wigner’s original papers [Wig1] and [Wig2], to obtain the distribution
of eigenvalues in random matrix models, a standard procedure is first to calculate
the moments of all orders after some appropriate scaling of eigenvalues and then
to determine an explicit distribution function of the moments. In this section we
discuss some deformations of the moments of a given density function and determine
the corresponding density function with respect to the new moments, and then
give an application in the following section. A natural question: which kind of
deformations of the moments has a corresponding distribution function? We will
make an interesting try.
   Let f (x) be a continuous density function on the interval I = (α, β), −∞ ≤
α < β ≤ +∞, and its moments
                              Z
                        mk = xk f (x)dx, k = 0, 1, 2, . . .                   (3.1)
                                  I
exist. Now consider the following deformation of moments with a parameter λ > 0,
                                  1
                       Mk =           mk , k = 0, 1, 2, . . . .            (3.2)
                               1 + λk
   Question: When is there a unique density function σ(x) which is a solution of
the moment problem of (3.2), i.e.
                     Z +∞
                            xk σ(x)dx = Mk , k = 0, 1, 2, . . . ,          (3.3)
                         −∞
and further how to determine σ(x) from the given density f (x)?
  Note that twoZdensity functions
                                Z f (x) = g(x) are said to be equal if their distri-
                     x                   x
bution functions         f (s)ds and         g(s)ds are equal at their points of continuity.
                    −∞                  −∞
To make sure that the density is unique we assume the moments mk , k = 0, 1, 2, . . .
satisfy Carleman’s condition
                                 X∞
                                      −1/2k
                                    m2k     = ∞.                             (3.4)
                                       k=0
12             DANG-ZHENG LIU, ZHENG-DONG WANG AND KUI-HUA YAN


It is obvious that the moments Mk also satisfy Carleman’s condition
                               X∞
                                     −1/2k
                                  M2k       = ∞.                                  (3.5)
                                  k=0
Due to Carleman’s famous theorem (see [ST] or [Fe]) (3.5) assures that the moment
problem of (3.2) is determined, i.e. numbers Mk determine a unique density whose
kth moment is Mk .
   To prove that the moment problem of (3.2) has a solution σ(x) whose spectrum
support suppσ(x) is to be contained in the interval J, given in advance, we first
introduce an important theorem. Let P (u) be any polynomial in u,
                                          X
                                  P (u) =   xk uk                            (3.6)
                                            k
where numbers xk are real constants. Introduce the functional µ(P ) defined by
                                       X
                               µ(P ) =    Mk xk .                            (3.7)
                                            k


Theorem 3.1 ([ST], Theorem 1.1) A necessary and sufficient condition that the
moment problem defined by the sequence of moments Mk shall have a solution on
J is that the functional µ(P ) be non-negative, that is
                     µ(P ) ≥ 0,    whenever P (u) ≥ 0 on J.                       (3.8)
   Note that Theorem 3.1 can be applied to derive explicit necessary and sufficient
conditions by a special choice of J, which depend on representations of non-negative
polynomials on J. In particular, it is a well-known fact (see [PóS]) any polynomial
P (u) ≥ 0 for all real u can be presented as
                               P (u) = P1 (u)2 + P2 (u)2                          (3.9)
where P1 (u) and P2 (u) are polynomials with real coefficients.
  If we take for P (u) the particular polynomial P (u) = (x0 + x1 u + · · · + xn un )2 ,
we have
                                       Xn
                             µ(P ) =       Mj+k xj xk .                        (3.10)
                                        j,k=0
From the theory of quadratic forms and (3.9) it is well known that the conditions
of (3.8) are equivalent to
                          |Mj+k |nj,k=0 > 0, n = 0, 1, 2, · · ·                  (3.11)
if the spectrum of the solution is not reducible to a finite set of points.

     Now we can state our theorem as follows:

Theorem 3.2 A probability density σ(x) on the real axis with its kth moment
Mk of (3.2) exists, and is uniquely determined by the differential equation of first
order
                            σ(x) − λ[xσ(x)](1) = f (x)χI                     (3.12)
with the following conditions
                                      Z +∞
                            σ(x) ≥ 0,       σ(x)dx = 1.                      (3.13)
                                         −∞
                   DENSITY OF EIGENVALUES IN UNITARY ENSEMBLES                         13


  Proof. It is sufficient to prove (3.11) if (Mj+k )nj,k=0 is a positive definite matrix.
Taking f (x)χI for the solution of the moments mk on the real axis, we have
                            |mj+k |nj,k=0 > 0, n = 0, 1, 2, · · · ,               (3.14)
that is,
        ∆n = (mj+k )nj,k=0 is a positive definite matrix, n = 0, 1, 2, · · · .    (3.15)
   The formula of (3.2) shows that
                                   (Mj+k )nj,k=0 = ∆n ∗ Λn                        (3.16)
where ∗ represents Schur product and
                                   1      n
                       Λn =                      , n = 0, 1, 2, · · · .           (3.17)
                              1 + λ(j + k) j,k=0
Thus by the property of Schur product it is sufficient to prove that Λn is a positive
definite matrix. Note that
                                          Z 1
                                  1
                                        =     tλk dt,                         (3.18)
                               1 + λk      0
then one can obtain
                             Z 1     Z 1 Y                       n
                        1                                    2 Y
            detΛn =              ···               tλj − tλk        dtj > 0.  (3.19)
                     (n + 1)! 0       0                         j=0
                                               0≤j<k≤n

   To derive the equation of (3.12) we make a Fourier transform and write
                        Z +∞                       Z
                 H(t) =        eitx σ(x)dx, F (t) = eitx f (x)dx.         (3.20)
                              −∞                                    I

Thus,
                   Z +∞                        ∞
                                               X (it)k+1
                          itx eitx σ(x)dx =                         Mk+1
                     −∞                                   k!
                                               k=0
                                               X∞
                                                     (it)k+1
                                           =                  (k + 1)Mk+1
                                                     (k + 1)!
                                               k=0
                                               X∞
                                                     (it)k
                                           =               kMk .                  (3.21)
                                                       k!
                                               k=0

Obviously,
               Z +∞                                 ∞
                                                    X (it)k
                      (1 + λitx) eitx σ(x)dx =                      (1 + λk)Mk
                −∞                                             k!
                                                    k=0
                                                    X∞
                                                          (it)k
                                               =                mk
                                                            k!
                                                 k=0
                                                 Z
                                               =   eitx f (x)dx.                  (3.22)
                                                     I

Namely,
                                                ′
                                   H(t) + λtH (t) = F (t).                        (3.23)
14              DANG-ZHENG LIU, ZHENG-DONG WANG AND KUI-HUA YAN


     We also make an inverse Fourier transform. Combining
                Z +∞                            Z
              1         −itx    ′           1 h +∞ −itx ′           i(1)
                      e      tH (t)dt = −              e    H (t)dt
             2π −∞                        2πi −∞
                                                   Z +∞
                                            1 h                       i(1)
                                      = −       ix       e−itx H(t)dt
                                          2πi       −∞
                                                 (1)
                                      = − xσ(x)                              (3.24)
and (3.23) one obtains the differential equation on real axis
                                         (1)
                          σ(x) − λ xσ(x)       = f (x)χI                     (3.25)
whose possible singular points are 0, α and β (if α and β are finite numbers). 

  Note that if α or β is finite we have an exact information about the spectrum
support J of σ(x).

     Case 1. −∞ < α ≤ 0, β = +∞ (or α = −∞, 0 ≤ β < +∞ similarly ). J = I.

     Any non-negative polynomial P (u) on I = (α, +∞) can be represented by
                 P (u) = P1 (u)2 + P2 (u)2 + (u − α)(P3 (u)2 + P4 (u)2 )     (3.26)
where P1 (u), P2 (u), P3 (u) and P4 (u) are polynomials with real coefficients.
   An analogous procedure to the above arguments of (3.9) — (3.11) one obtains
the conditions of (3.8) are equivalent to the matrices
                               (Mj+k )nj,k=0 , n = 0, 1, 2, · · ·            (3.27)
and
                       (Mj+k+1 − αMj+k )nj,k=0 , n = 0, 1, 2, · · ·          (3.28)
are positive definite .
   Note that
                          mk+1 − αmk                       mk
        Mk+1 − αMk =                   + (−λα)                               (3.29)
                          1 + λ(k + 1)          (1 + λk)(1 + λ(k + 1))
and
               (mj+k )nj,k=0 , (mj+k+1 − αmj+k )nj,k=0 , n = 0, 1, 2, · · ·  (3.30)
are positive definite, by the property of Schur product one obtains the matrices of
(3.27) and (3.28) are positive definite.

  Case 2. −∞ < α ≤ 0, 0 < β < +∞ (or −∞ < α < 0, 0 ≤ β < +∞ similarly ).
J = I.

     Any non-negative polynomial P (u) on I = (α, β) can be represented by
                          P (u) = P1 (u)2 + (u − α)(β − u)P2 (u)2            (3.31)
where P1 (u) and P2 (u) are polynomials with real coefficients.
  In this case the conditions of (3.8) are equivalent to the matrices
                               (Mj+k )nj,k=0 , n = 0, 1, 2, · · ·            (3.32)
and
         (−Mj+k+2 + (β − α)Mj+k+1 − αβMj+k )nj,k=0 , n = 0, 1, 2, · · ·      (3.33)
are positive definite .
                  DENSITY OF EIGENVALUES IN UNITARY ENSEMBLES                        15


  Note that −Mk+2 + (β − α)Mk+1 − αβMk can be rewritten by
−mk+2 + (β − α)mk+1 − αβmk                    mk
                           + (−αβ)
       1 + λ(k + 1)                 (1 + λk)(1 + λ(k + 1))
                                          mk+2
                           + λ                                                   (3.34)
                               (1 + λ(k + 1))(1 + λ(k + 2))
and

      (mj+k )nj,k=0 , (−mk+2 + (β − α)mk+1 − αβmk )nj,k=0 , n = 0, 1, 2, · · ·   (3.35)

are positive definite, again by the property of Schur product one obtains the ma-
trices of (3.32) and (3.33) are positive definite.

  Remark 3.1 If αβ > 0, we cannot make sure that the matrices of (3.33) are
positive, e.g., α = 1, β = 2, n = 0
                                               m2    m1
          −M2 + (β − α)M1 − αβM0 = −               +      − 2m0 < 0              (3.36)
                                             1 + 2λ 1 + λ
for the sufficiently large λ.

  Case 3. 0 < α < β (or α < β < 0 similarly ). J = (0, β) (or J = (α, 0)).

  We take f (x) for the density function on (0, β)
                           
                              f (x)     α<x<β
                   fe(x) =
                              0          0<x≤α                                   (3.37)

and by using Case 1 and Case 2 it is obvious.

  Anyway, we get

Theorem 3.3 A probability density σ(x) with its kth moment Mk of (3.2) exists,
and is uniquely determined by the following differential equation
  (1) α ≤ 0 ≤ β,
                          σ(x) − λ[xσ(x)](1) = f (x), x ∈ I;                     (3.38)
  (2) α < β < 0,

                       σ(x) − λ[xσ(x)](1) = f (x)χI , x ∈ (α, 0);                (3.39)

  (3) 0 < α < β,

                       σ(x) − λ[xσ(x)](1) = f (x)χI , x ∈ (0, β).                (3.40)


                       4. Proofs of Theorems 1.2 and 1.5
   Proof of Theorem 1.2.
   In [Wig1] Wigner made use of the integral representation of Bessel function of
order 1 (pointed out by W.Feller to him) to get his semicircle law. However, in the
original analysis he got the semicircle law by leading a differential equation. In the
following we will combine these two kinds of method to derive the density.
16               DANG-ZHENG LIU, ZHENG-DONG WANG AND KUI-HUA YAN


                                               P
                                               k
                                                      j   j
     Uniqueness: Note that using (1.33)              C2k C2k−j ≤ 32k holds. Writing B =
                                               j=0
a + |b|, then we have
                                k
                                X      j   j
                        M2k ≤         C2k C2k−j a2j b2k−2j ≤ (3B)2k .                  (4.1)
                                j=0

Thus, the Carleman’s condition is satisfied and the density function is determined
by the moments.
   Derivation of the differential equation: Putting 2a = 1, using the integral repre-
sentation of Bessel function of order zero (see [Sze] ) and calculating directly

          Z +∞                                 ∞
                                               X (it)k
                 (1 + λitx) eitx σ(x)dx =                   (1 + λk)Mk
           −∞                                         k!
                                               k=0
                                               ∞       [k/2]
                                               X (it)k X
                                           =                      Ckj Ck−j
                                                                       j
                                                                           a2j bk−2j
                                                      k!    j=0
                                               k=0

                                                     (it)k           k
                                               X∞
                                                                 a
                                           =              L0 az + + b
                                                       k!        z
                                               k=0

                                                       (it)k            k
                                                     ∞
                                                     X             a
                                           = L0                   az ++b
                                                        k!         z
                                                  k=0
                                                                a     
                                           = L0 exp it az + + b
                                                                 z
                                              itb
                                                                   a 
                                           = e L0 exp it az +
                                                                     z
                                                   ∞
                                                   X         2k
                                                       (ita)
                                           = eitb               Ck
                                                        (2k)! 2k
                                                      k=0
                                           = eitb J0 (2at)
                                                   Z
                                                  1 1         eitx
                                           = eitb         √         dx
                                                  π −1       1 − x2
                                             Z b+1
                                           =         eitx fb (x)dx                     (4.2)
                                                b−1

where
                                           1      1
                                fb (x) =     p             .                           (4.3)
                                           π  1 − (x − b)2
     By using Theorems 3.2 and 3.3 in Section 3 one obtains
                                         (1)
                            σ(x) − λ xσ(x)     = fb (x)χIb                             (4.4)

and supp(σ) = [B1 , B2 ] for λ > 0 while for λ = 0, obviously
                                           1      1
                                σ(x) =       p             .                           (4.5)
                                           π  1 − (x − b)2
                   DENSITY OF EIGENVALUES IN UNITARY ENSEMBLES                         17


   Now we give an exact solution of the equation of (4.4). Note that when λ > 0
the equation of (4.4) is known as a Cauchy-Euler equation of order 1 (see [GN],
P99). Thus we have
   (1) −1 < b < 1. 0 is a singular point of the equation of (4.4), and we have to
determine how the solutions for x < 0 and x > 0 can be pieced together to give
solutions valid on the whole interval Ib (see [GN], P22 ). So one obtains
                                Z
                     
                        1 1 −1 b+1 − 1
                          xλ           s λ fb (s)ds             x > 0,
             σ(x) =     λ         x Z x                                         (4.6)
                     
                        1      1                 1
                          (−x) λ −1      (−s)− λ fb (s)ds     x < 0.
                        λ            b−1
                                                           Z b+1
   Note that xσ(x) is absolutely continuous on Ib . Thus         σ(x)dx = 1 from the
                                                               b−1
equation of (4.4) by using σ(b − 1) = σ(b + 1) = 0. Besides, σ(x) is continuous on
Ib when 0 < λ < 1, while 0 is a singular point of σ(x) when λ ≥ 1.
  (2) b = ±1.
  If b = 1, then
                             Z 2
                   1 1                1    1
             σ(x) = x λ −1         s− λ p         ds,            x ∈ (0, 2).        (4.7)
                   λ          x          s(2 − s)
While b = −1, we have
                       Z x
            1     1
                    −1           1    1
     σ(x) = (−x)  λ        (−s)− λ p          ds,                x ∈ (−2, 0).       (4.8)
            λ           −2          −s(s + 2)

  (3) b > 1. The equation of (4.4) can be rewritten by one nonhomogeneous linear
equation on Ib = (b − 1, b + 1)
                                     (1)         1
                       σ(x) − λ xσ(x)      =p                              (4.9)
                                               1 − (x − b)2
and the other homogeneous linear equation on (0, b − 1)
                                          (1)
                            σ(x) − λ xσ(x)      = 0.                               (4.10)
  Next we solve the related nonhomogeneous linear equation of (4.9)
                                      Z b+1
                               1 1             1
                        σ(x) = x λ −1       s− λ fb (s)ds                          (4.11)
                               λ       x
and the related homogeneous linear equation of (4.10)
                                                 1
                                   σ(x) = C+ x λ −1                                (4.12)
where                                 Z b+1
                                  1              1
                             C+ =             s− λ fb (s)ds.                       (4.13)
                                  λ     b−1
                                                                     R b+1
Note that C+ assures that σ(x) is continuous on (0, b + 1) and         0
                                                                             σ(x)dx = 1.
  (4) b < −1. Similar to the case where b > 1. One obtains the solutions of the
equation of (4.4) on Ib = (b − 1, b + 1)
                                         Z x
                              1      1              1
                      σ(x) = (−x) λ −1        (−s)− λ fb (s)ds           (4.14)
                              λ           b−1
18                 DANG-ZHENG LIU, ZHENG-DONG WANG AND KUI-HUA YAN


and on (b + 1, 0)
                                                                 1
                                          σ(x) = C− (−x) λ −1                                                        (4.15)
where                                             Z b+1
                                              1                  1
                                   C− =                   (−s)− λ fb (s)ds.                                          (4.16)
                                              λ    b−1

     In the end we will complete the proof of Theorem 1.5. First, setting
                               Hn = span{p0 (x), p1 (x), · · · , pn−1 (x)},                                          (4.17)
                                                            2
then Hn is an n-dimensional subspace of L (ω). It is obvious that p̂0 (x)p(x), · · · ,
p̂n−l−1 (x)p(x) is a family of normalized orthogonal vectors in Hn , extended by
                        (n)             (n)
                       e0 (x), · · · , el−1 (x), p̂0 (x)p(x), · · · , p̂n−l−1 (x)p(x)                                (4.18)
to a normalized orthogonal base of Hn .
   Let Pn be a projective operator from L2 (ω) to Hn . We construct an operator
from Hn to itself as follows,
                                    Tn(k) = Pn ◦ Akx : Hn −→ Hn                                                      (4.19)
where Ax is the multiplication by x.

                                                                                                     (N )               (N )
Lemma 4.1 Denote the kth moments of σ̂N (x) and σN (x) by Mk                                                and M̂k
respectively, then
                           (N )   (N )
                         M̂k = Mk + ΘN                                                                               (4.20)
where
             1      NX
                      −1                                         l−1
                                                                 X                                              
                                                                         (k)    (N )          (N )
 ΘN =                              Akx (p̂j p), p̂j p L2 (ω) −          TN (ej         ), ej         L2 (ω)
                                                                                                                    . (4.21)
          N (cN )k                                               j=0
                         j=N −l



Proof. We first point out that
                             Z
                     (N )
                   Mk     =    xk σN (x)dx
                                               N
                                               X −1 Z
                                         1
                                    =                 xk p2j (x)ω(x)dx
                                      N (cN )k j=0 I

                                           1
                                                 N
                                                 X −1 D             E
                                                        k
                                    =                  x  p j , p j
                                        N (cN )k j=0                  L2 (ω)


                                         1
                                                 NX−1 D                  E
                                                         (k)
                                    =                  T N   (p j ), p j
                                      N (cN )k j=0                         L2 (ω)

                                           (k) 
                                      Tr TN
                                    =            .                                                                    (4.22)
                                      N (cN )k
On the other hand, by the normalized orthogonal base of (4.18) one obtains
                       NX
                        −l−1                                    l−1
                                                                X
           (k)                   (k)                                   (k)    (N )      (N )
      Tr TN        =            TN (p̂j p), p̂j p L2 (ω) +             TN (ej         ), ej      L2 (ω)
                                                                                                            .        (4.23)
                        j=0                                     j=0
                     DENSITY OF EIGENVALUES IN UNITARY ENSEMBLES                           19


    Thus,
                Z
        (N )
      M̂k =         xk σ̂N (x)dx
              N
              X −1 Z
        1
 =                   xk p̂2j (x)p2 (x)ω(x)dx
     N (cN )k j=0 I
               −l−1 Z
             NX                                 l−1 Z
                                                 X                        
      1                                                     (N )
 =        k
                       xk 2
                         p̂ j (x)p 2
                                     (x)ω(x)dx +       xk (ej (x))2 ω(x)dx + ΘN
   N (cN )     j=0   I                           j=0 I

      1      NX
               −l−1                                l−1
                                                   X                              
                      (k)                                (k) (N )      (N )
 =                  T N   (p̂ j p), p̂ j p  2
                                           L (ω)
                                                 +     T N  (e j  ), e j     2
                                                                            L (ω)
                                                                                    + ΘN
   N (cN )k    j=0                                 j=0
        (N )
 = Mk          + ΘN .                                                                (4.24)


                                                 1/2
Lemma 4.2 Write k · k = h·, ·iL2 (ω) , and for k = 1, 2, · · · ,
                                                   k−1
                                                    Y            
                           kAkx f k ≤ 3k                     Dn+j kf k, ∀f ∈ Hn      (4.25)
                                                       j=0

where
                                            Dn = max {aj , |bj |}.                   (4.26)
                                                       0≤j≤n


                         NP
                          −1
Proof. Set f (x) =             lj pj (x) and note that
                         j=0

                                   N
                                   X −1
                 kAx f k2 = k                lj xpj k2
                                   j=0
                                   N
                                   X −1
                            = k              lj (aj+1 pj+1 + bj pj + aj pj−1 )k2
                                   j=0
                                   X
                            = k             (lj−1 aj + lj bj + lj+1 aj+1 )pj k2
                                       j
                                X
                            =              |lj−1 aj + lj bj + lj+1 aj+1 |2
                                   j
                                            X
                                                    2
                            ≤ 3Dn2                (lj−1 + lj2 + lj+1
                                                                 2
                                                                     )
                                             j

                            ≤ 9Dn2 kf k2 .                                           (4.27)
Thus, kAx f k ≤ 3Dn kf k and (4.25) is easily proved. 

  Proof of Theorem 1.5: By (1.21) and Lemma 4.2, there exist constants C1 and
C2 which only depend on k such that
            Akx (p̂j p), p̂j p L2 (ω) ≤ kAkx (p̂j p)k ≤ 3k DN
                                                            k
                                                              +l+k−2 ≤ C1 (cN )
                                                                               k
                                                                                     (4.28)
20                DANG-ZHENG LIU, ZHENG-DONG WANG AND KUI-HUA YAN


and
         (k)   (N )      (N )                           (N )
       TN (ej         ), ej     L2 (ω)
                                         ≤ kPN k kAkx ej       k ≤ 3 k DN
                                                                        k                k
                                                                          +k−1 ≤ C2 (cN ) .   (4.29)

Thus, using (4.21), for the large N ,

                     1      NX
                              −1           l−1
                                           X             (C + C )l
                                                            1   2
        |ΘN | ≤          k
                                 C (c
                                  1 N )k
                                         +     C2 (cN )k =          .                         (4.30)
                  N (cN )                  j=0
                                                              N
                                  j=N −l

     Again by Lemma 4.1 we get
                                             (N )               (N )
                                    lim M̂k         = lim Mk           = Mk .                 (4.31)
                                   n→∞               n→∞



                                          Acknowledgements
  The authors would like to thank Professors Y. Chen, D. S. Lubinsky and W. Van
Assche for helpful discussions on the relationship between the density of eigenvalues
and zero distribution of the orthogonal polynomials. The first author is grateful to
Professor Y. Chen for his encouragement.

                                              References
[Br]    B. V. Bronk, Exponential ensemble for random matrices, J. Math. Phys. 6 (1965),
        228-237.
[CI]    Y. Chen and M.E.H. Ismail, Asymptotics of extreme zeros of the Meixner-Pollaczek
        polynomials, J. Comput. Appl. Math. 82 (1997), 59-78.
[De]    P. Deift, Orthogonal polynomials and random matrices: a Riemann-Hilbert approach,
        Courant Lecture Notes in Mathematics. 3, Courant Institute, New York, 1999.
[DKM]   P. Deift, T. Kriecherbauer and K.T-R. McLaughlin, New results for the asymptotics
        of polynomials and related problems via the Lax-Levermore method, in Recent ad-
        vances in partial differential equations, Venice 1996, 87- 104. Proc.Sympos.Appl.Math.,
        54.Amer.Math.Soc., Providence, R.I., 1998.
[DKMVZ] P. Deift, T. Kriecherbauer, K.T.-R. McLaughlin, S. Venakides and X. Zhou, Uniform
        asymptotics for polynomials orthogonal with respect to varying exponential weights
        and applications to universality questions in random matrix theory, Comm. Pure Appl.
        Math. 52 (1999), 1335-1425.
[Fe]    W. Feller, An Introduction to Probability Theory and its Applications, Vol.II, 2nd Ed,
        John Wiley & Sons, Inc., New York, 1971.
[FK]    D. Fox and P. B. Kahn, Higher order spacing distributions for a class of unitary en-
        sembles, Phys. Rev. 134 (1964), B1151-1155.
[GN]    M. M. Guterman and Z. H. Nitecki, Differential Equations: A First Course, Saunders
        College Pub., New York, 1984.
[HT]    U. Haagerup and S. Thorbjornsen, Random matrices with complex Gaussian entries,
        Expo. Math. 21 (2003), 293-337.
[HZ]    J. Harer and D. Zagier, The Euler characteristic of the moduli space of curves, Invent.
        math. 85 (1986) , 457-485.
[Jo]    K. Johansson, On fluctuations of eigenvalues of random Hermitian matrices, Duke
        Math. J. 91 (1998), no.1,151-204.
[Kö]   W. König, Orthogonal polynomial ensembles in probability theory, Probab. Surveys 2
        (2005), 385-447.
[Led]   M. Ledoux, Differential operators and spectral distributions of invariant ensembles
        from the classical orthogonal polynomials, the continuous case, Electron. J. Probab. 9
        (2004), 177-208.
[Lef]   H. S. Leff, Class of ensembles in the statistical theory of energy level spectra, J. Math.
        Phys. 5 (1964), 763-768.
                 DENSITY OF EIGENVALUES IN UNITARY ENSEMBLES                                   21


[Lub]    D. S. Lubinsky, Asymptotics of orthogonal polynomials: some old, some new, some
         identities. Acta Appl. Math. 61 (2000), 207- 256.
[LMS]    D. S. Lubinsky, H. N. Mhaskar and E. B. Saff, A proof of Freud’s conjecture for
         exponential weights, Constr. Approx. 4 (1988), 65-83.
[Me]     M. L. Mehta, Random Matrices, Elsevier Academic Press, 3rd ed, New York, 2004.
[ND]     P. G. Nevai and J. S. Dehesa, On asymptotic average properties of zeros of orthogonal
         polynomials, SIAM J. Math. Anal. 10 (1979), 1184-1192.
[PS]     L. Pastur and M. Shcherbina , Universality of the Local Eigenvalue statistics for a class
         of unitary invariant random matrix ensembles , J. Stat. Phys. 86 (1997), 109-147.
[PóS]   G. Pólya and G. Szeg ö, Problems and Theorems in Analysis, Vol.II, Springer-Verlag,
         New York,1976.
[Po]     C. E. Porter (Ed.), Statistical Theories of Spectra, Academic Press, New York, 1965.
[ST]     J. A. Shohat and J. D. Tamarkin, Mathematical Surveys, the Problem of Moments,
         Amer. Math. Soc., New York, 1943.
[Si]     B. Simon, Orthogonal Polynomials on the Unit Circle, Parts 1 and 2, American Math-
         ematical Society, Providence, 2005.
[Sze]    G.Szegö, Orthogonal Polynomials, American Mathematical Society Colloquium Pub-
         lications, 23. Amer. Math. Soc., Providence, R.I., 1939.
[TW]     C. A. Tracy and H. Widom, Introduction to random matrices, In: Geometric and
         Quantum Aspects of Integrable Systems, Lecture Notes in Physics, Vol.424 (1993),
         103-130.
[Va1]    W. Van Assche, Asymptotics for Orthogonal Polynomials, Lecture Notes in Mathe-
         matics Vol.1265, Springer, Berlin, 1987.
[Va2]    W. Van Assche, Asymptotics for orthogonal polynomials and three-term recurrences,
         in Orthogonal Polynomials: Theory and Practice (P. Nevai, ed.), Vol.294 of NATO-ASI
         Series C: Mathematical and Physical Sciences, Kluwer Academic Publishers, Dordrecht,
         1990, 435-462.
[Wig1]   E. P. Wigner, Characteristic vectors of bordered matrices with infinite dimensions,
         Ann. Math. 62 (1955), 548-564.
[Wig2]   E. P. Wigner, Characteristic vectors of bordered matrices with infinite dimensions II,
         Ann. Math. 65 (1957), 203-207.

   SCHOOL OF MATHEMATICAL SCIENCES, PEKING UNIVERSITY, BEIJING, 100871,
P. R. CHINA
   E-mail address: DZLIUMATH@GMAIL.COM

   SCHOOL OF MATHEMATICAL SCIENCES, PEKING UNIVERSITY, BEIJING, 100871,
P. R. CHINA
   E-mail address: ZDWANG@PKU.EDU.CN

  SCHOOL OF MATHEMATICS AND PHYSICS, ZHEJIANG NORMAL UNIVERSITY,
ZHEJIANG JINHUA, 321004, P. R. CHINA
  E-mail address: YANKH@ZJNU.CN
