The early study of large random matrices was stimulated by analysis of high-dimensional data. One example is Wishart’s (1928) investigation on large covariance matrices whose statistical properties are mainly determined by eigenvalues and eigenvectors from the point view of a principal components analysis. Since then, the random matrix theory has been developed very rapidly and found many applications in areas such as heavy-nuclei atoms (Wigner, 1955), number theory (Mezzadri and Snaith, 2005), quantum mechanics (Mehta, 2005), condensed matter physics (Forrester, 2010), wireless communications (Couillet and Debbah, 2011).
The study of random matrices has greatly been motivated by Tracy and Widom’s (1994, 1996) work. They show that the largest eigenvalues of the three Hermitian matrices (Gaussian orthogonal ensemble, Gaussian unitary ensemble and Gaussian symplectic ensemble) converge to some special distributions that are now known as the Tracy-Widom laws. Subsequently, the Tracy-Widom laws have found their applications in the study of problems such as the longest increasing subsequence (Baik et al., 1999), combinatorics, growth processes, random tilings and the determinantal point processes (see, e.g., Tracy and Widom (2002), Johansson (2007) and references therein) and the largest eigenvalues in the high-dimensional statistics (see, e.g., Johnstone (2001, 2008) and Jiang (2009)). Some recent research focuses on the universality of the largest eigenvalues of matrices with non-Gaussian entries; see, for example, Tao and Vu (2011), Erdős et al. (2012) and the references therein.
A very recent paper by Jiang and Qi (2017) studies the largest radii of three rotation-invariant and non-Hermitian random matrices: the spherical ensemble, the truncation of circular unitary ensemble and the product of independent complex Ginibre ensembles. It is proved in the paper that the spectral radii converge to the Gumbel distribution and some new distributions.
where C is a normalizing constant. See, e.g., Zyczkowski and Sommers (2000).
Assume p=pn depends on n and set c=limn→∞npn. Życzkowski and Sommers (2000) show that the empirical distribution of zi’s converges to the distribution with density proportional to (1−∣z∣2)21 for ∣z∣≤c if c∈(0,1). Dong et al. (2012) prove that the empirical distribution goes to the circular law and the arc law as c=0 and c=1, respectively. See also Diaconis and Evans (2001) and Jiang (2009, 2010) and references therein for more results.
Jiang and Qi (2017) have proved that the spectral radius max1≤j≤p∣zj∣ for the truncated circular unitary ensemble converges to the Gumbel distribution when the dimension of the truncated truncated circular unitary matrix is of the same order as the dimension of the original circular unitary matrix, see Theorem 1 in section 2.
In this paper we consider heavily truncated and lightly truncated circular unitary matrices and investigate the limiting distribution of the spectral radii for those truncated circular unitary matrices. Our results complement that in Jiang and Qi (2017).
The rest of the paper is organized as follows. The main results in this paper are given in section 2 and their proofs are provided in section 3.
Main Results
Consider the pn×pn submatrix A, truncated from a n×n circular unitary matrix U in section 1. Denote the pn eigenvalues as z1,⋯,zpn with the joint density function given by (1.1).
For completeness, we first quote a theorem in Jiang and Qi (2017) on the limiting distribution of the spectral radii max1≤j≤pn∣zj∣ before we give our results in the paper.
The main results of the paper are the following theorems:
Under condition (2.4), ((k+1)!)1/k2n1+1/k(max1≤j≤pn∣zj∣−1) converges weakly to the reversed Weibull distribution Wk(x) defined as
We notice that the limiting distribution of the spectral radii depends on the dimension of truncated matrices. Our results in Theorems 2, 3 and 4 indicate that the limiting distribution of the spectral radii of the truncated circular unitary matrices is Gumbel distribution Λ if the parameter kn=n−pn, the number of truncated columns and rows diverges. When the truncation is very light, that is, kn=n−pn=k≥1 is a fixed integer, the limiting distribution of the spectral radii of the truncated matrices is the reversed Weibull distribution Wk.
It is obvious that the case when kn=n−pn is of order between logn and (logn)3 has not been covered in Theorems 1 to 4. We conjecture that max1≤j≤pn∣zj∣, after properly normalized, converges in distribution to the Gumbel distribution in this case.
Proofs
For random variables {Xn;n≥1} and constants {an;n≥1}, we write Xn=Op(an) if limx→+∞limsupn→∞P(∣anXn∣≥x)=0. It is well known that anbnXn→0 in probability as n→∞ if Xn=Op(an) and {bn;n≥1} is a sequence of constants with limn→∞bn=∞.
Let Ui, i≥1 be a sequence of i.i.d. random variables uniformly distributed over (0,1), and U1:n≤U2:n≤⋯≤Un:n be the order statistics of U1,U2,⋯,Un for each n≥1. Then from page 14 on the book by Balakrishnan and Cohen (1991), we know that the cumulative distribution function of Ui:n is given by
for each 1≤i≤n, and the probability density function (pdf) of Ui:n is given by
This is the so-called Beta distribution, denoted by Beta(i,n−i+1).
From (3.2), Upn−j+1:n−j has a Beta(pn−j+1,kn) distribution with pdf given by
For each n≥2, let {Ynj;1≤j≤pn} be independent random variables such that Ynj and (Upn−j+1:n−j)1/2 have the same distribution. Jiang and Qi (2017) have shown that max1≤j≤pn∣zj∣ and max1≤j≤pnYnj have the same distribution, that is
for x∈(0,1). See the proof of Theorem 2 in Jiang and Qi (2017).
We will present some useful lemmas before we prove our main results.
Suppose {ln;n≥1} is sequence of positive integers. Let znj∈[0,1) be real numbers for 1≤j≤ln such that max1≤j≤lnznj→0 as n→∞. Then n→∞limj=1∏ln(1−znj)∈(0,1) exists if and only if the limit n→∞limj=1∑lnznj=:z∈(0,∞) exists and the relationship of the two limits is given by
which implies log(1−znj)=−znj+O(znj2) uniformly over 1≤j≤ln since max1≤j≤lnznj→0 as n→∞. Therefore,
The lemma can be easily concluded from the above expression. ■
Let {ln} be a sequence of positive integers such that ln→∞ and for each n, {znj,1≤j≤ln} are non-negative numbers such that znj is non-increasing in j with zn1>0. Then for any sequence of positive integers {rn}satisfying that rn<ln for all large n and rn/ln→1 as n→∞, we have
as n→∞. In fact, from the monotonicity of znj, we have znj≤rn1∑j=1rnznj for rn+1≤j≤ln. Hence
where tn∈(0,1) will be specified later in the proof of each theorem. From (3.4),
Obviously, we have for 1≤j≤pn.
Assume that 1≤pn<n and pn→∞ as n→∞. Let {rn} satisfy the condition in Lemma 3.2 with ln=pn. Assume αn>0 and βn are real numbers such that limn→∞P(Yn12>βn+αnx)=0 for any x∈R. If (max1≤j≤rnYnj2−βn)/αn converges in distribution to a cdf G, then (max1≤j≤pnYnj2−βn)/αn converges in distribution to the same distribution G.
Proof. Note that (max1≤j≤rnYnj2−βn)/αn converges in distribution to the cdf G if and only if
for every continuity point x of G with G(x)∈(0,1). We need to prove the above expression is still true when rn is replaced by pn. Now fix x, a continuity point of G with G(x)∈(0,1). Set tn=tn(x)=βn+αnx and define znj as in (3.7). Note that (3.8) holds, zn1→0 as n→∞,
By using Lemma 3.1 and (3.10) we have ∑j=1rnznj→z=−logG(x) which, together with Lemma 3.2, implies ∑j=1pnznj→z=−logG(x). Once again we have from Lemma 3.1 that
This completes the proof of the lemma. ■
Let Zn be nonnegative random variables such that (Zn2−βn)/αn converges weakly to a cdf G(x), where αn>0 and βn>0 are constants satisfying that limn→∞αn/βn=0. Then
Proof. Set Wn=(Zn2−βn)/αn. We have Zn2=βn+αnWn=βn(1+βnαnWn). Then by Taylor’s expansion
where for i=1,2, li(t) is a polynomial in t of degree ≤3i, depending on r and k, and all of its coefficients are of order O(\big{(}\frac{r}{(r-k)k}\big{)}^{i/2}).
Define Vpn−j+1:n−j as in (3.25). Assume that kn=n−pn→∞ and kn/n→0 as n→∞. Then for any δn>0 such that δn→∞ and δn=o(kn1/6)
uniformly over 0≤x≤δn, 1≤j≤pn−kn as n→∞.
Proof. Set βnj(x)=n−jpn−j+(n−j)3/2((pn−j)kn)1/2x=n−jpn−j(1+(n−j)1/2(pn−j)1/2kn1/2x). Then 1−βnj(x)=n−jn−pn−(n−j)3/2((pn−j)kn)1/2x=n−jkn(1−(n−j)1/2kn1/2(pn−j)1/2x), and the density function of Vpn−j+1:n−j is given by
To estimate hj(x), we need Stirling’s formula:
and Taylor’s expansion: 1−t=exp(log(1−t))=exp(−t−21t2+O(t3)) as t→0. By applying Stirling’s formula to kn!, (pn−j)! and (n−j)!, the product in (3.16) is equal to 2π1+o(1)(1+o(1)) for ∣x∣=o(kn1/6) uniformly over 1≤j≤pn−kn as n→∞. By applying Taylor’s expansion to (1+(n−j)1/2(pn−j)1/2kn1/2x)pn−j and (1−(n−j)1/2kn1/2(pn−j)1/2x)kn, the product in (3.16) is equal to
for ∣x∣=o(kn1/6) uniformly over 1≤j≤pn−kn as n→∞. Therefore, we have
uniformly over 1≤j≤pn−kn as n→∞. Next we will give an estimate of the upper bound of hj(x) for large x. Note that βnj(x)<1 if and only if x<(n−j)1/2kn1/2(pn−j)1/2=O(kn1/2) uniformly over 1≤j≤pn−kn and thus (n−j)1/2(pn−j)1/2kn1/2x≤O(n1/2kn1/2)→0 uniformly over 0<x<(n−j)1/2kn1/2(pn−j)1/2, 1≤j≤pn−kn as n→∞. Now by applying Taylor’s expansion to (1+(n−j)1/2(pn−j)1/2kn1/2x)pn−j and inequality
to (1−(n−j)1/2kn1/2(pn−j)1/2x)kn we get
uniformly over 1≤j≤pn−kn as n→∞.
Assume that 0≤x≤δn. From (3.19) we have
uniformly over 0≤x≤δn, 1≤j≤pn−kn as n→∞. Therefore, to complete the proof of the lemma, it suffices to show that
uniformly over 0≤x≤δn, 1≤j≤pn−kn as n→∞.
uniformly if y→∞. Since (δnkn1/6)1/2=o(kn1/6), by using (3.18) we have
proving (3.20). This completes the proof of the lemma. ■
2 Proofs of the Theorems
where An=cn+21(1−cn2)1/2(n−1)−1/2an, Bn=21(1−cn2)1/2(n−1)−1/2bn,
with a(x)=(logx)1/2−(logx)−1/2log(2πlogx) and b(x)=(logx)−1/2 for x>3.
Since Ynj2 and Upn−j+1:n−j are identically distributed, we have
Part 1. First we show (3.21) under condition (2.1). We will prove that
Let jn=[pn5/8], the integer part of pn5/8. For 1≤j≤jn, define
Then we see that uniformly over 1≤j≤jn,
uniformly for all 1≤j≤jn as n→∞.
In Lemma 3.6, take r=n−j and k=n−pn to have
uniformly over 1≤j≤jn as n→∞, where
and where, for i=1,2, li(t) is a polynomial in t of degree ≤3i, depending on n, and all of its coefficients are of order O((1/pn)i/2) uniformly over 1≤j≤jn as n→∞. Now, by taking B=(unj,∞) we obtain
uniformly for 1≤j≤jn as n→∞. From L’Hospital’s rule, we have that for any r≥0
Since min1≤j≤jnunj→∞ as n→∞ by (3.24), it follows from (3.27) that
holds uniformly over 1≤j≤jn. Furthermore, since the coefficients of li(t) are uniformly bounded by O((1/pn)i/2) for i=1,2, we have
In Lemma 3.5, by taking xn=ncn2/(1−cn2) and and define cn,j such that unj=(j−1)cn,j+a(xn)+b(xn)x for 1≤j≤jn with cn,1=xn−1/2. It follows from (3.24) that cn,j=xn−1/2(1+o(1)) uniformly over 1≤j≤jn as n→∞, which implies limn→∞max1≤j≤jn∣cn,jxn1/2−1∣=0. Then we get
We will show that the second term and the third term on the line below (3.2) converge to zero as n→∞. By noting that n1/2cn(1−cn2)1/2∼pn−1/2 we have
uniformly over 1≤j≤jn. Thus, it follows from (3.13) that
as n→∞. Therefore, by combining (3.2), (3.26) and (3.29) we get
for all large n. Then we estimate znjn by using (3.26) with j=jn and (3.30)
as n→∞. From (3.8) we have ∑j=jn+1pnznj≤pnznjn=O(pn−1/2)→0 as n→∞, which together with (3.31) yields
We can also prove from (3.26) that zn1→0 as n→∞. In view of (3.22) and Lemma 3.1 we conclude (3.23), ie.,
where αn=cn(1−cn2)1/2(n−1)−1/2bn and βn=cn2+cn(1−cn2)1/2(n−1)−1/2an. Since
as n→∞, we can apply Lemma 3.4 and get that
Recall that An=cn+21(1−cn2)1/2(n−1)−1/2an=cn(1+o(1)) and Bn=21(1−cn2)1/2(n−1)−1/2bn∼21(n−1)−1/2(logpn)−1/2. Then
ie., (3.21) holds. The proof of Part 1 is completed.
Part 2. We will show (3.21) under condition (2.2). First, it follows from condition (2.2) that
as n→∞. Noting that n≤n2/kn≤n2, we get that
is of order (logn)1/2 as n→∞.
Use the same notation as in Part 1. Recall that pn/n→1, kn=n−pn=o(n) and logn=o(kn1/3) as n→∞. In order to use both Lemmas 3.5 and 3.7, we take xn=1−cn2ncn2. Define jn=[5(logn)1/2xn]+1. Then jn∼kn5n(logn)1/2=o(n), which implies 1≤jn≤pn−kn for all large n. Define cn,j for 1≤j≤jn in the same way as in Part 1. Similar to the proof of (3.24) we can show that
uniformly for 1≤j≤jn as n→∞. Then unj=O((logn)1/2)=o(kn1/6) uniformly for 1≤j≤jn as n→∞. We can also verify that all conditions in Lemma 3.5 are satisfied. Thus, from (3.14) and (3.12) we have
Note that (3.14) holds uniformly over 1≤j≤pn−kn and unjn≥4(logn)1/2 for all large n. By employing (3.14) with j=jn and x=4(logn)1/2 and using equation (3.8) and Lemma 3.7 we have
Thus, we obtain that ∑j=1pnznj→e−x for any x. Then equation (3.23) follows from equation (3.22) and Lemma 3.1. The rest of the proof will follow from the same lines in the proof of the first part. Again Lemma 3.4 will be used. The details are omitted. ■
Proof of Theorem 3. Recall that an is given by
By using Stirling’s formula (3.17), we have under condition (2.3) that
Since te−t is strictly increasing in (0,1), for all large n such that yn<1 define εn as the unique solution to te−t=yn in (0,1), that is, εne−εn=yn, which implies that εn→0 and εn∼yn as n→∞ and
for all large n. Then it follows from the the first inequality in (3.37) that
for all large n, which together with (3.36) implies that an≤knεn for all large n,and thus an=o(kn) as n→∞. By plugging y=an in (3.37) and using an≤knεn for large n we conclude
Define znj as in (3.7) with tn=tn(x)=1−nan(1−knx) for any fixed x. Then 1−tn=nan(1−knx)=o(nkn)=o(nlogn) as n→∞, where we have used the fact that kn=n−pn→∞ and kn=o(logn) from (2.3). This implies n(1−tn)2→0 as n→∞.
It is easy to verify the following expression
where 0≤d(t)≤t2 for 0≤t≤1/2. Then
Furthermore, by using Stirling’s formula (3.17), we get
uniformly over 1≤j≤jn, where jn:=pn−kn3. Then from (3.9) we obtain that
Note that n(1−tn)=o(kn). Then it follows from (3.39) and (3.37) that
uniformly over 1≤j≤jn, and thus
The second integral above is dominated by the first one since tkne−t is increasing over (0,kn) and kn3(1−tn)/(n(1−tn))=kn3/n→0 as n→∞. Therefore, in view of (3.37) and (3.38) we have
as n→∞. Therefore, it follows from Lemma 3.2 that ∑j=1pnznj→e−x as n→∞. It is easy to conclude that zn1→0 as n→∞ from (3.39), (3.37) and the above estimates. Accordingly, by taking βn=1−nan and αn=nknan with G(x)=Λ(x) in Lemmas 3.3 and 3.4 we conclude that
as n→∞. This completes the proof of the theorem. ■
Fix x<0. Let tn=tn(x)=1+n1+1/k((k+1)!)1/kx. Then tn∈(0,1) for all large n. Since n(1−tn)→0 as n→∞, we have
Therefore, we have from (3.7) and (3.9) that
uniformly over 1≤j≤pn=n−k. Since
we have max1≤j≤pnznj=O(1/n)→0 as n→∞. To complete the proof of (3.40), by using (3.5) we need to show that
Let {jn} be a sequence of integers such that jn→∞ and jn/n→0 as n→∞. Then
uniformly over 1≤j≤n−jn, which implies that
as n→∞. This proves (3.41) and thus we obtain (3.40).
Finally, the theorem follows from Lemma 3.4 with αn=n1+1/k((k+1)!)1/k and βn=1. This completes the proof. ■
Acknowledgements We would like to thank an anonymous referee for his/her careful reading of the original version of the paper and pointing out some imperfections in the proofs. Gui’s work was partially supported by the program for the Fundamental Research Funds for the Central Universities (2014RC042).