arxiv: v2 [math.pr] 2 Nov 2009

Transcription

1 arxiv: v2 [math.pr] 2 Nov 2009 Limit Distribution of Eigenvalues for Random Hankel and Toeplitz Band Matrices Dang-Zheng Liu and Zheng-Dong Wang School of Mathematical Sciences Peking University Beijing, 0087, P. R. China September, 2009 Abstract Consider real symmetric, complex Hermitian Toeplitz and real symmetric Hankel band matrix models, where the bandwidth b N but b N /N b, b [0, ] as N. We prove that the distributions of eigenvalues converge weakly to universal, symmetric distributions γ T (b) and γ H (b). In the case b > 0 or b = 0 but with the addition of b N C N 2 +ǫ0 for some positive constants ǫ 0 and C, we prove almost sure convergence. The even moments of these distributions are the sum of some integrals related to certain pair partitions. In particular, when the bandwidth grows slowly, i.e. b = 0, γ T (0) is the standard Gaussian distribution and γ H (0) is the distribution x exp( x 2 ). In addition, from the fourth moments we know that the γ T (b) s are different for different b s, the γ H (b) s different for different b [0, 2 ] and the γ H (b) s different for different b [ 2, ]. Introduction In Random Matrix Theory, the most important information is contained in the eigenvalues of matrices and the most prominent analytical object is the distribution of eigenvalues. That is, for a real symmetric or complex Hermitian N N matrix A with eigenvalues λ,,λ N, the distribution of 0 Key words: Random matrix theory; distribution of eigenvalues; Toeplitz band matrices; Hankel band matrices.

2 its eigenvalues is the normalized probability measure µ A := N N δ λj. (.) j= If A is a random matrix, then µ A is a random measure. The limit distribution of (.) is of much interest in Random Matrix Theory. In [8], Wigner found his famous semicircular law and then in [9] pointed out that the semicircular law was valid for a wide class of real symmetric random matrices. Since then much more has been done on various random matrix models, a standard reference is Metha s book [6]. Recently in his review paper [], Bai proposes the study of random matrix models with certain additional linear structure. In particular, the properties of the distributions of eigenvalues for random Hankel and Toeplitz matrices with independent entries are listed among the unsolved random matrix problems posed in [], Section 6. Bryc, Dembo and Jiang [7] proved the existence of limit distributions γ H and γ T for real symmetric Hankel and Toeplitz matrices. The moments of γ H and γ T are the sum of volumes of solids, from which we can see that γ H and γ T are symmetric and of unbounded support. At the same time Hammond and Miller [2] also independently proved the existence of limit distribution for symmetric Toeplitz matrices. In the present paper we shall prove the existence of limit distribution of eigenvalues for real symmetric and complex Hermitian Toeplitz band matrices, and also for real symmetric Hankel band matrices. Random band matrices have arisen in connection with the theory of quantum chaos [8, ], and the limit distribution of eigenvalues of band matrices has been studied in [3, 7]. Note that Toeplitz matrices emerge in many aspects of mathematics and physics and also in plenty of applications, e.g. [6], especially [2, 9] for connections with random matrices. Hankel matrices arise naturally in problems involving power moments, and are closely related to Toeplitz matrices. Explicitly, for a Toeplitz matrix of the form T N = (a i j ) N i,j= and a Hankel matrix of the form H N = (a i+j 2 ) N i,j=, if P N = (δ i,n j ) N i,j= the backward identity permutation, then P N T N is a Hankel matrix and P N H N is a Toeplitz matrix. In this paper we always write a Hankel matrix H N = P N T N where we assume the matrix entries a N+,,a 0,,a N are real-valued, thus H N is a real symmetric matrix. In addition, if we introduce the Toeplitz or Jordan matrices B = (δ i+,j ) N i,j= and F = (δ i,j+) N i,j=, called the backward shift and forward shift because of their effect on the elements of the standard basis {e,,e N } of R N, then an N N matrix 2

3 A can be written in the form A = N j=0 N a j B j + j= a j F j (.2) if and only if A is a Toeplitz matrix where a N+,,a 0,,a N are complex numbers [3]. It is worth emphasizing that this representation of a Toeplitz matrix is of vital importance and the starting point of our method. The shift matrices B and F exactly present the information of the traces. Consider a Toeplitz band matrix as follows. Given a band width b N < N, let {, i j b N ; η ij = (.3) 0, otherwise. Then a Toeplitz band matrix is T N = (η ij a i j ) N i,j=. (.4) Moreover, the Toeplitz band matrix T N can be also written in the form T N = b N j=0 a j B j + b N j= a j F j. (.5) Obviously, a Toeplitz matrix can be considered as a band matrix with the bandwidth b N = N. Note that when referring to a Hankel band matrix H N, we always mean H N = P N T N where T N is a Toeplitz band matrix. In this paper the three models under consideration are : (i) Hermitian Toeplitz band matrices. The model consists of N-dimensional random Hermitian matrices T N = (η ij a i j ) N i,j= in Eq.(.4). We assume that Re a j = Re a j and Im a j = Im a j for j =,2,, and {a 0 } {Re a j,im a j } j N is a sequence of independent real random variables such that E[a j ] = 0, E[ a j 2 ] = for j Z (.6) and further Notice that sup j Z sup E[ a j k ] = C k < for k N. (.7) j Z E[ a j k ] < sup{e[ Re a j k ], E[ Im a j k ]} <, j Z 3

4 and the assumption (.7) shows that when fixing k, these moments of the independent random variables whose orders are not larger than k can be controlled by some constant only depending on k. The case of real symmetric Toeplitz band matrices is very similar to the Hermitian case, except that we now consider N-dimensional real symmetric matrices T N = (η ij a i j ) N i,j=. (ii) Symmetric Toeplitz band matrices. We assume that a j = a j for j = 0,,, and {a j } j=0 is a sequence of independent real random variables such that and further E[a j ] = 0, E[ a j 2 ] = for j = 0,,, (.8) sup E[ a j k ] = C k < for k N. (.9) j 0 (iii) Symmetric Hankel band matrices. Let T N = (η ij a i j ) N i,j= be a Toeplitz band matrix and H N = P N T N be the corresponding Hankel band matrix. We assume that {a j } j Z is a sequence of independent real random variables such that and further E[a j ] = 0, E[ a j 2 ] = for j Z (.0) sup E[ a j k ] = C k < for k N. (.) j Z We assume that b N grows either as ) b N /N b as N, b (0,] (proportional growth), or as 2) b N as N, b N = o(n) (slow growth). Under these assumptions we will establish the limit distributions γ T (b) for Toeplitz band matrices and γ H (b) for Hankel band matrices. In the case b > 0, we prove that the distribution of eigenvalues for Toeplitz (Hankel) band matrices converges almost surely to γ T (b) (γ H (b)), the moments of which are the sum of some integrals depending on b. For b = 0, we obtain that γ T (0) is the standard Gaussian distribution and γ H (0) is the distribution x exp( x 2 ). In the case b = 0, with the addition of b N C N 2 +ǫ 0 for some positive constants ǫ 0 and C, we also prove almost sure convergence. The plan of the remaining part of our paper is the following. Sections 2 and 3 are devoted to the statements and proofs of the main theorems for band matrices with proportionally growing and slowly growing bandwidths b N, respectively. In Section 4, the fourth moments of the limit distributions γ T (b) and γ H (b) are calculated, which shows their difference for different b s. 4

5 2 Proportional Growth of b N In order to calculate the moments of the limit distribution, we review some basic combinatorical concepts. Definition 2.. Let the set [n] = {,2,,n}. () We call π = {V,,V r } a partition of [n] if the blocks V j ( j r) are pairwise disjoint, non-empty subsets of [n] such that [n] = V V r. The number of blocks of π is denoted by π, and the number of the elements of V j is denoted by V j. (2) Without loss of generalization, we assume that V,,V r have been arranged such that s < s 2 < < s r, where s j is the smallest number of V j. Therefore we can define the projection π(i) = j if i belongs to the block V j ; furthermore for two elements p,q of [n] we write p π q if π(p) = π(q). (3) The set of all partitions of [n] is denoted by P(n), and the subset consisting of all pair partitions, i.e. all V j = 2, j r, is denoted by P 2 (n). The subset of P 2 (n) consisting of such pair partitions that each contains exactly one even number and one odd number is denoted by P 2 (n). Note that P 2 (n) is an empty set if n is odd. We can formulate our results for Hankel and Toeplitz band matrices with the proportional growth of b N as follows. By convention, I B represents the characteristic function of a set B. Theorem 2.2. Let T N be either a Hermitian ((.6) (.7)) or real symmetric ((.8) (.9)) Toeplitz random band matrix, where b N /N b as N, b (0,]. Take the normalization X N = T N / (2 b)bn, then µ XN converges almost surely to a symmetric probability distribution γ T (b) which is determined by its even moments m 2k (γ T (b)) = (2 b) k π P 2 (2k) 2 [0,] [,] k j= I [0,] (x 0 +b j ǫ π (i)x π(i) ) d x l (2.) where ǫ π (i) = if i is the smallest number of π (π(i)); otherwise, ǫ π (i) =. Theorem 2.3. Let H N be a real symmetric ((.0) (.)) Hankel random band matrix, where b N /N b as N, b (0,]. Take the normalization Y N = H N / (2 b)bn, then µ YN converges almost surely to a symmetric i= l=0 5

6 probability distribution γ H (b) which is determined by its even moments m 2k (γ H (b)) = (2 b) k π P 2 (2k) 2 [0,] [,] k j= I [0,] (x 0 b j ( ) i x π(i) ) d x l. i= l=0 (2.2) Let us first give two basic lemmas about traces of Toeplitz and Hankel band matrices. Although their proofs are simple, they are very useful in treating random matrix models closely related to Toeplitz matrices. Lemma 2.4. For a Toeplitz band matrix T = (η ij a i j ) N i,j= with the bandwidth b N where a N+,,a N are complex numbers, we have the trace formula tr(t k ) = N b N i= j,,j k = b N l= a jl l I [,N] (i + j q ) P, k N. δ0, k j q l= (2.3) Proof. For the standard basis {e,,e N } of the Euclidean space R N, we have T e i = b N j=0 a j B j e i + b N j= a j F j e i = Repeating T s effect on the basis, we have T k e i = b N j,,j k = b N l= a jl b N I [,N] (i + l= By tr(t k ) = k e t i T k e i, we complete the proof. i= j= b N a j I [,N] (i + j)e i+j. l j q P. )ei+ k j q Lemma 2.5. Given a Toeplitz band matrix T = (η ij a i j ) N i,j= with the bandwidth b N where a N+,,a N are real numbers, let H = PT be the corresponding Hankel band matrix where P = (δ i,n j ) N i,j=. We have the trace formula 6

7 tr(t k )= N b N a jl i= j,,j k = b N l= l= N b N a jl i= j,,j k = b N l= l= I [,N] (i I [,N] (i l l ( ) q j q ) δ0, k P, ( ) q j q ( ) q j q ) δ2i N, k P, ( ) q j q k even, k odd. (2.4) Proof. Follow a similar procedure as in the proof of Lemma 2.4. We are now ready to prove the main results of this section. Proof of Theorem 2.2. We prove the theorem by the following steps: Step. Calculation of the Moments Let m k,n = E[ x k µ XN (dx)] = N E[tr(Xk N)]. Using Lemma 2.4, we have m k,n = N k 2 (2 b) k 2 b k 2 N b N l I [,N] (i + j q ) P E[ a δ0, k jl ]. j q i= j,,j k = b N l= l= For j=(j,,j k ), we construct a set of numbers S j = { j,..., j k } with multiplicities. Note that the random variables whose subscripts have different absolute values are independent. If S j has one number with multiplicity, by independence of the random variables, we have E[ k l= a j l ] = 0. Thus the only contribution to the k-th moment comes when each of S j has multiplicity at least 2, which implies that S j has at most [ k 2 ] distinct numbers. Once we have specified the distinct numbers of S j, the subscripts j,,j k are determined in at most 2 k ([ k 2 ])k ways. By independence and the assumptions (.7) and (.9), we find Therefore, for odd k m k,n = O((b N ) k 2 +[ k 2 ] ) = O(N k 2 +[ k 2 ] ). lim m k,n = 0. N It suffices to deal with m 2k,N. If each of S j = { j,, j 2k } (here j=(j,,j 2k )) has multiplicity at least 2 but one of which at least 3, then S j has at most k distinct numbers. Thus, the contribution of such terms to m 2k,N is 7

8 O(N ). So it suffices to consider all pair partitions of [2k] = {,2,,2k}. That is, for π P 2 (2k), if p π q, then it is always the case that j p = j q or j p = j q. Under the condition 2k j q = 0 according to (2.3), considering the main contribution to the trace, we should take j p = j q ; otherwise, there exists p 0,q 0 [2k] such that j p0 = j q0 = 2 (j p 0 + j q0 We can choose other k distinct numbers, which determine j p0 = j q0. This shows that there is a loss of at least one degree of freedom and the contribution of such terms is O(N ). Thus a jp and a jq are conjugate. So we can write m 2k,N = o()+ N k (2 b) k b k i.e. N b N 2k j q ). 2 π P 2 (2k) i= j,,j k = b N l= I [,N] (i+ Follow a similar argument from [7, 3]: for fixed π P 2 (2k), N k+ N k+ N i= N i= b N 2 j,,j k = b N l= b N j,,j k = b N l= 2 I [ I [,N] (i + l ǫ π (q)j π(q) ). l ǫ π (q)j π(q) ), N,]( i N + l ǫ π (q) j π(q) N ), can be considered as a Riemann sum of the definite integral 2 [0,] [ b,b] k l= I [0,] (x 0 + l ǫ π (q)x π(q) ) d x l. l=0 (2.5) Given π P 2 (2k), as N, each term in (2.5) can be treated as an integral. Thus we have m 2k (γ T (b)) = lim m 2k,N is the representation of N the form (2.). Step 2. Carleman s Condition 8

9 Obviously, from (2.) we have m 2k (γ T (b)) (2 b) k π P 2 (2k) [0,] [,] k d x l = l=0 2 k (2k )!!. (2 b) k By Carleman s theorem (see [0]), the limit distribution γ T (b) is uniquely determined by the moments. Step 3. Almost Sure Convergence It suffices to show that N= N 4 E[(tr(Xk N ) E[tr(Xk N )])4 ] < (2.6) for each fixed k. Indeed, by Lemma 2.4, we have where A[p 0 ;p] = tr(x k N) = l= a pl N k 2 A[p 0 ;p] p 0,p I [,N] (p 0 + l= l p q P, )δ0, k p q p = (p,...,p k ), and the summation p 0,p runs over all possibilities that p { b N,...,b N } k and p 0 {,...,N}. Thus, N 4 E[(tr(Xk N) E[tr(X k N)]) 4 ] = (2 b) 2k b 2k N 4+2k i 0,j 0,s 0,t 0 i,j,s,t E[ p {i,j,s,t} (A[p 0 ;p] E[A[p 0 ;p]])] (2.7) where the summation i 0,j 0,s 0,t 0 runs over all possibilities that i,j,s,t i,j,s,t { b N,...,b N } k and i 0,j 0,s 0,t 0 {,...,N}. For i = (i,...,i k ), as in Step we still construct a set of numbers S i = { i,..., i k } with multiplicities. Obviously, if one of S i,...,s t, for example, S i does not have any number coincident with numbers of the other three sets, then the term in (2.7) equals 0 by independence. Also, if the union S = S i S t has one number with multiplicity, the term in (2.7) equals 0. Now, let us estimate the non-zero term in (2.7). Assume that S has p distinct numbers with multiplicities ν,...,ν p, subject to the constraint 9

10 ν ν p = 4k. If p 2k 2, by (.6) (.9), the contribution of the terms corresponding to S is bounded by (2 b) 2k b 2k C k N 4+2k N 4 b 2k 2 N C k,b N 2 for some constants C k and C k,b. For the case of p = 2k, we have two numbers shared by three of four sets S i,...,s t each or one number shared by the four sets. It is not hard to see that there always exists one of S i,...,s t, for example, S i, which has one number with multiplicity, denoted by i q0 for some q 0 {,...,k}. Consequently, under the constraint k i q = 0 according to (2.3), we have i q0 = i q0 k i q. Hence we can choose the other 2k 2 distinct numbers of S, which determine the subscript i q0. This shows that there is a loss of at least one degree of freedom and the contribution of the terms corresponding to S is bounded by C k,b N 2. The case of p = 2k implies each number in S occurs exactly two times. Thus there exist two of S i,...,s t, for example, S i and S j, which share one number with one of S s and S t respectively. Consequently, under the constraints k k i q = 0 and j q = 0 according to (2.3), there is a loss of at least two degrees of freedom and the contribution of the terms corresponding to S is bounded by C k,b N 2. Therefore, N 4 E[(tr(Xk N ) E[tr(Xk N )])4 ] C k,b N 2 and almost sure convergence is proved. Consequently, the proof of Theorem 2.2 is complete. Proof of Theorem 2.3. Using Lemma 2.5 and our assumptions, the proof is very similar to that of of Theorem 2.2. Here we don t repeat the process. 0

11 Note that in the representation of even moments in Eq.(2.2) the pair partition π P2 (2k) comes from the constraint k ( ) q j q = 0 (2.8) in the form of the trace for Hankel matrices in Lemma Slow Growth of b N We give our results for Hankel and Toeplitz band matrices with slow growth of the bandwidth as follows. Theorem 3.. Let T N be either a Hermitian ((.6) (.7)) or real symmetric ((.8) (.9)) Toeplitz random band matrix, where b N but b N /N 0 as N. Take the normalization X N = T N / 2b N, then µ XN converges weakly to the standard Gaussian distribution. In addition, if there exist some positive constants ǫ 0 and C such that b N C N 2 +ǫ 0, then µ XN converges almost surely to the standard Gaussian distribution. Theorem 3.2. Let H N be a real symmetric ((.0) (.)) Hankel random band matrix, where b N but b N /N 0 as N. Take the normalization Y N = H N / 2b N, then µ YN converges weakly to the distribution f(x) = x exp( x 2 ). In addition, if there exist positive constants ǫ 0 and C such that b N C N 2 +ǫ 0, then µ YN converges almost surely to the distribution f(x) = x exp( x 2 ). Remark 3.3. For real symmetric palindromic Toeplitz Matrices, i.e. real symmetric Toeplitz matrices under extra conditions: a j = a N j,0 < j < N, Massey, Miller and Sinsheimer [5] have obtained the Gaussian normal distribution for eigenvalues. And for random reverse circulant matrices, i.e., Hankel matrices under extra conditions: a j = a N j, < j < N, Bose and Mitra [5] obtained the same distribution of eigenvalues f(x) = x exp( x 2 ) as in Theorem 3.2. Since proofs of both theorems are very similar, we prove only the Toeplitz case. Proof of Theorem 3.. The proof is quite similar to that of Theorem 2.2 and some details will be omitted. We will lay a strong emphasis on the derivation of the Gaussian distribution. Here the slow growth of b N leads to an easy

12 calculation of the complicated integrals. We now complete the proof of the Gaussian law by showing () m k,n = N E[tr(Xk N )] converges to the k-th moment of the standard Gaussian distribution. (2) Assume that b N C N 2 +ǫ 0 for some positive constants ǫ 0 and C. For each fixed k, N 4 E[(tr(Xk N ) E[tr(Xk N )])4 ] <. (3.) N= By Lemma 2.4 and our assumptions (.6) (.9), it follows that Since b N as N, for odd k m k,n = O((b N ) k 2 +[ k 2 ] ). lim m k,n = 0. N It suffices to deal with m 2k,N. However, again by Lemma 2.4 and our assumptions, the contribution with the exception of all pair partitions to m 2k,N is O(b N ). So it suffices to consider all pair partitions of [2k] = {,2,,2k}. That is, for π P 2 (2k), if p π q, then it is always the case that j p = j q or j p = j q. Under the condition 2k j q = 0 according to (2.3), considering the main contribution to the trace, we should take j p = j q (otherwise, there is a loss of at least one degree of freedom). Thus a jp and a jq are conjugate. So we can write m 2k,N = o() + N i.e. 2 k b k N N b N 2 π P 2 (2k) i= j,,j k = b N l= I [,N] (i + Note that the constraints l i + ǫ π (q)j π(q) N, l 2k, N i l N + ǫ π (q) j π(q) N 2, l 2k. l ǫ π (q)j π(q) ). (3.2)

13 Since as N, we have l ǫ π (q) j π(q) N lb N N 0 Hence I [,N] (i + l ǫ π (q)j π(q) ) = + o(). m 2k,N = o() + (2k )!!( 2b N + 2b N ) k ( + o()) 2k (3.3) converges to (2k )!! as N. Assertion () is then proved. We now prove (2). The same argument as in the proof of almost sure convergence in Theorem 2.2 shows that Since b N C N 2 +ǫ 0, we have N 4 E[(tr(Xk N ) E[tr(Xk N )])4 ] = O( b 2 ). N N 4 E[(tr(Xk N ) E[tr(Xk N )])4 ] C k N +2ǫ, 0 where C k is a constant depending on k only. This completes the proof of assertion (2). The proof of Theorem 3. is then complete. 4 First Four Moments of γ T (b) and γ H (b) We will compute the second and fourth moments of γ T (b) and γ H (b). From the fourth moment we can read that the γ T (b) s are different for different b s, the γ H (b) s different for different b [0, 2 ] and the γ (b) s different for H different b [ 2,]. However, for the 2k-th moments (k 3), the integrals on the right-hand sides of (2.) and (2.2) are in general different for different pair partitions, which makes difficult the explicit calculations of the higher moments. Observe that for π P2 (2k), the corresponding integrals on the righthand sides of (2.) and (2.2) are in fact the same. Therefore, it is sufficient 3

14 to calculate the integral in the Toeplitz case for every pair partition. For the pair partition π P 2 (2k), we introduce the symbol p π (b) = 2 [0,] [,] k j= I [0,] (x 0 + b j ǫ π (i)x π(i) ) d x l. i= For k= it is easy to obtain the second moments m 2 (γ T (b)) =, m 2 (γ H (b)) =. Thus the fourth moment is the first free moment in seeing the shape of the distribution. When k = 2, a direct calculation (this can be checked using Mathematica) for all pair partitions π = {{,2}, {3,4}}, π 2 = {{,4}, {2,3}},π 3 = {{,3}, {2,4}} of P 2 (4) shows l=0 and { 2 p π (b) = p π2 (b) = 3 (6 5b), b [0, 2 ]; +6b 2b 3, b ( 3b 2 2,] (4.) p π3 (b) = { 4( b), b [0, 2 ]; 2( +6b 6b 2 +2b 3 ) 3b 2, b ( 2,]. (4.2) Thus, for b [0, 2 ], 4(9 8b) m4(γt (b)) = 3(2 b) 2 which strictly decreases on [0, 2 ]. When b ( 2,], m 4 (γ T (b)) = 4( + 6b 3b2 ) 3b 2 (2 b) 2 which also strictly decreases on ( 2,]. Therefore, one knows that m 4 (γ T (b)) strictly decreases on [0,]. Further, the distributions γ T (b) s are different. In particular, γ T (b) (0 < b ) is indeed different from the normal distribution. Note that for b = the fact that the limit distribution is not Gaussian has been observed in [4, 7, 2]. We turn to the Hankel-type distribution γ H (b). From (4.), one obtains that for b [0, 2 ] 4(6 5b) m4(γh (b)) = 3(2 b) 2 4

15 which strictly decreases on [0, 2 ] while for b ( 2,], m 4 (γ H (b)) = 2( + 6b 2b2 ) 3b 2 (2 b) 2 which strictly increases on ( 2,]. Thus, according to the fourth moments, we only know that the γ H (b) s (0 b 2 ) are different and the γ (b) s H ( 2 b ) are different. An added note. After the paper as submitted we learned from the Associate Editor and the referee about a related preprint Limiting Spectral Distribution of Some Band Matrices by Basak and Bose at the site statmath/html/publication/techreport.html. Although their paper and ours contain the same results for Hankel and Toeplitz band matrices, the former assumes less integrability on the entries of a matrix, allows more general rates for the bandwidth, and also covers more ensembles of structured matrices related to Toeplitz matrices. On the other hand, our paper covers Hermitian Toeplitz matrices and gives low order moment calculations. Besides, our paper gives a different method for analyzing Toeplitz matrices by treating them as a linear combination of deterministic matrices with independent coefficients. This method can be used to derive other results, including those that deal with semicircle law. Acknowledgements The authors thank the Associate Editor for valuable comments, the referees for very helpful comments and careful writing instruction on the draft of this paper, and Steven Miller for information on his paper. References [] Z. D. Bai, Methodologies in spectral analysis of large-dimensional random matrices, a review, Statist. Sinica. 9 (999), [2] E. L. Basor, Toeplitz determinants, Fisher Hartwig symbols and random matrices, Recent Perspectives in Random Matrix Theory and Number Theory, London Math. Soc. Lecture Note Ser. 322 (2005), [3] L. V. Bogachev, S. A. Molchanov and L. A. Pastur, On the level density of random band matrices, Math. Notes. 50 (99), [4] A. Bose, S. Chatterjee, and S. Gangopadhyay, Limiting spectral distribution of large dimensional random matrices, J. Indian Statist. Assoc. 4 (2003),

16 [5] A. Bose and J. Mitra, Limiting spectral distribution of a special circulant, Stat. Probab. Letters 60 (2002), no., 20. [6] A. Bottcher and S. M. Grudsky, Spectral Properties of Banded Toeplitz Matrices, Society for Industrial Mathematics, Philadelphia, [7] W. Bryc, A. Dembo and T. Jiang, Spectral measure of large random Hankel, Markov and Toeplitz matrices, Ann. Probab. 34 (2006), no., 38. [8] G. Casati, L. Molinari and F. Izrailev, Scaling properties of band random matrices, Phys. Rev. Lett. 64 (990), [9] P. Diaconis, Patterns in eigenvalues: The 70th Josiah Willard Gibbs lecture, Bull. Amer. Math. Soc. 40 (2003), [0] W. Feller, An Introduction to Probability Theory and Its Applications, Volum 2, Wiley, New York, 97. [] Y. V. Fyodorov and A. D. Mirlin, Scaling properties of localization in random band matrices: A σ-model approach, Phys. Rev. Lett. 67 (99), [2] C. Hammond and S. J. Miller, Distribution of eigenvalues for the ensemble of real symmetric Toeplitz matrices, J. Theoret. Probab. 8 (2005), [3] R. A. Horn and C. R. Johnson, Matrix Analysis, Cambridge University Press, New York, 985. [4] A. Khorunzhy and W. Kirsch, On asymptotic expansions and scales of spectral universality in band random matrix ensembles, Commun. Math. Phys. 23 (2002), [5] A. Massey, S. J. Miller and J. Sinsheimer, Distribution of Eigenvalues of Real Symmetric Palindromic Toeplitz Matrices and Circulant Matrices, J. Theoret. Probab. 20 (2007), [6] M. L. Mehta, Random Matrices, 3nd ed., Academic Press, San Diego, [7] S. A. Molchanov, L. A. Pastur and A. M. Khorunzhy, Eigenvalue distribution for band random matrices in the limit of their infinite rank, Teor. Matem. Fizika. 90 (992), [8] E. P. Wigner, Characteristic vectors of bordered matrices with infinite dimensions, Ann. Math. 2 (955), [9] E. P. Wigner, On the distribution of the roots of certain symmetric matrices, Ann. Math. 67 (958),