Spectral Radii of Truncated Circular Unitary Matrices

Wenhao Gui, Yongcheng Qi

Introduction

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 CC is a normalizing constant. See, e.g., Zyczkowski and Sommers (2000).

Assume p=pnp=p_{n} depends on nn and set c=lim⁡n→∞pnnc=\lim_{n\to\infty}\frac{p_{n}}{n}. Życzkowski and Sommers (2000) show that the empirical distribution of ziz_{i}’s converges to the distribution with density proportional to 1(1−∣z∣2)2\frac{1}{(1-|z|^{2})^{2}} for ∣z∣≤c|z|\leq c if c∈(0,1)c\in(0,1). Dong et al. (2012) prove that the empirical distribution goes to the circular law and the arc law as c=0c=0 and c=1c=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 max⁡1≤j≤p∣zj∣\max_{1\leq j\leq p}|z_{j}| 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×pnp_{n}\times p_{n} submatrix A\mathbf{A}, truncated from a n×nn\times n circular unitary matrix U\mathbf{U} in section 1. Denote the pnp_{n} eigenvalues as z1,⋯ ,zpnz_{1},\cdots,z_{p_{n}} 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 max⁡1≤j≤pn∣zj∣\max_{1\leq j\leq p_{n}}|z_{j}| before we give our results in the paper.

The main results of the paper are the following theorems:

Under condition (2.4), 2n1+1/k((k+1)!)1/k(max⁡1≤j≤pn∣zj∣−1)\frac{2n^{1+1/k}}{((k+1)!)^{1/k}}(\max_{1\leq j\leq p_{n}}|z_{j}|-1) converges weakly to the reversed Weibull distribution Wk(x)W_{k}(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 Λ\Lambda if the parameter kn=n−pnk_{n}=n-p_{n}, the number of truncated columns and rows diverges. When the truncation is very light, that is, kn=n−pn=k≥1k_{n}=n-p_{n}=k\geq 1 is a fixed integer, the limiting distribution of the spectral radii of the truncated matrices is the reversed Weibull distribution WkW_{k}.

It is obvious that the case when kn=n−pnk_{n}=n-p_{n} is of order between log⁡n\log n and (log⁡n)3(\log n)^{3} has not been covered in Theorems 1 to 4. We conjecture that max⁡1≤j≤pn∣zj∣\max_{1\leq j\leq p_{n}}|z_{j}|, after properly normalized, converges in distribution to the Gumbel distribution in this case.

Proofs

For random variables {Xn; n≥1}\{X_{n};\,n\geq 1\} and constants {an; n≥1}\{a_{n};\,n\geq 1\}, we write Xn=Op(an)X_{n}=O_{p}(a_{n}) if lim⁡x→+∞lim sup⁡n→∞P(∣Xnan∣≥x)=0\lim_{x\to+\infty}\limsup_{n\to\infty}P(|\frac{X_{n}}{a_{n}}|\geq x)=0. It is well known that Xnanbn→0\frac{X_{n}}{a_{n}b_{n}}\to 0 in probability as n→∞n\to\infty if Xn=Op(an)X_{n}=O_{p}(a_{n}) and {bn; n≥1}\{b_{n};\,n\geq 1\} is a sequence of constants with lim⁡n→∞bn=∞\lim_{n\to\infty}b_{n}=\infty.

Let UiU_{i}, i≥1i\geq 1 be a sequence of i.i.d. random variables uniformly distributed over (0,1)(0,1), and U1:n≤U2:n≤⋯≤Un:nU_{1:n}\leq U_{2:n}\leq\cdots\leq U_{n:n} be the order statistics of U1,U2,⋯ ,UnU_{1},U_{2},\cdots,U_{n} for each n≥1n\geq 1. Then from page 14 on the book by Balakrishnan and Cohen (1991), we know that the cumulative distribution function of Ui:nU_{i:n} is given by

for each 1≤i≤n1\leq i\leq n, and the probability density function (pdf) of Ui:nU_{i:n} is given by

This is the so-called Beta distribution, denoted by Beta(i,n−i+1i,n-i+1).

From (3.2), Upn−j+1:n−jU_{p_{n}-j+1:n-j} has a Beta(pn−j+1,knp_{n}-j+1,k_{n}) distribution with pdf given by

For each n≥2n\geq 2, let {Ynj; 1≤j≤pn}\{Y_{nj};\,1\leq j\leq p_{n}\} be independent random variables such that YnjY_{nj} and (Upn−j+1:n−j)1/2(U_{p_{n}-j+1:n-j})^{1/2} have the same distribution. Jiang and Qi (2017) have shown that max⁡1≤j≤pn∣zj∣\max_{1\leq j\leq p_{n}}|z_{j}| and max⁡1≤j≤pnYnj\max_{1\leq j\leq p_{n}}Y_{nj} have the same distribution, that is

for x∈(0,1)x\in(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}\{l_{n};\,n\geq 1\} is sequence of positive integers. Let znj∈[0,1)z_{nj}\in[0,1) be real numbers for 1≤j≤ln1\leq j\leq l_{n} such that max⁡1≤j≤lnznj→0\max_{1\leq j\leq l_{n}}z_{nj}\to 0 as n→∞n\to\infty. Then lim⁡n→∞∏j=1ln(1−znj)∈(0,1)\displaystyle\lim_{n\to\infty}\prod^{l_{n}}_{j=1}(1-z_{nj})\in(0,1) exists if and only if the limit lim⁡n→∞∑j=1lnznj=:z∈(0,∞)\displaystyle\lim_{n\to\infty}\sum^{l_{n}}_{j=1}z_{nj}=:z\in(0,\infty) exists and the relationship of the two limits is given by

which implies log⁡(1−znj)=−znj+O(znj2)\log(1-z_{nj})=-z_{nj}+O(z_{nj}^{2}) uniformly over 1≤j≤ln1\leq j\leq l_{n} since max⁡1≤j≤lnznj→0\max_{1\leq j\leq l_{n}}z_{nj}\to 0 as n→∞n\to\infty. Therefore,

The lemma can be easily concluded from the above expression. ■\blacksquare

Let {ln}\{l_{n}\} be a sequence of positive integers such that ln→∞l_{n}\to\infty and for each nn, {znj, 1≤j≤ln}\{z_{nj},~{}1\leq j\leq l_{n}\} are non-negative numbers such that znjz_{nj} is non-increasing in jj with zn1>0z_{n1}>0. Then for any sequence of positive integers {rn}\{r_{n}\}satisfying that rn<lnr_{n}<l_{n} for all large nn and rn/ln→1r_{n}/l_{n}\to 1 as n→∞n\to\infty, we have

as n→∞n\to\infty. In fact, from the monotonicity of znjz_{nj}, we have znj≤1rn∑j=1rnznjz_{nj}\leq\frac{1}{r_{n}}\sum^{r_{n}}_{j=1}z_{nj} for rn+1≤j≤lnr_{n}+1\leq j\leq l_{n}. Hence

where tn∈(0,1)t_{n}\in(0,1) will be specified later in the proof of each theorem. From (3.4),

Obviously, we have for 1≤j≤pn1\leq j\leq p_{n}.

Assume that 1≤pn<n1\leq p_{n}<n and pn→∞p_{n}\to\infty as n→∞n\to\infty. Let {rn}\{r_{n}\} satisfy the condition in Lemma 3.2 with ln=pnl_{n}=p_{n}. Assume αn>0\alpha_{n}>0 and βn\beta_{n} are real numbers such that lim⁡n→∞P(Yn12>βn+αnx)=0\lim_{n\to\infty}P(Y_{n1}^{2}>\beta_{n}+\alpha_{n}x)=0 for any x∈Rx\in R. If (max⁡1≤j≤rnYnj2−βn)/αn(\max_{1\leq j\leq r_{n}}Y_{nj}^{2}-\beta_{n})/\alpha_{n} converges in distribution to a cdf GG, then (max⁡1≤j≤pnYnj2−βn)/αn(\max_{1\leq j\leq p_{n}}Y_{nj}^{2}-\beta_{n})/\alpha_{n} converges in distribution to the same distribution GG.

Proof. Note that (max⁡1≤j≤rnYnj2−βn)/αn(\max_{1\leq j\leq r_{n}}Y_{nj}^{2}-\beta_{n})/\alpha_{n} converges in distribution to the cdf GG if and only if

for every continuity point xx of GG with G(x)∈(0,1)G(x)\in(0,1). We need to prove the above expression is still true when rnr_{n} is replaced by pnp_{n}. Now fix xx, a continuity point of GG with G(x)∈(0,1)G(x)\in(0,1). Set tn=tn(x)=βn+αnxt_{n}=t_{n}(x)=\beta_{n}+\alpha_{n}x and define znjz_{nj} as in (3.7). Note that (3.8) holds, zn1→0z_{n1}\to 0 as n→∞n\to\infty,

By using Lemma 3.1 and (3.10) we have ∑j=1rnznj→z=−log⁡G(x)\sum^{r_{n}}_{j=1}z_{nj}\to z=-\log G(x) which, together with Lemma 3.2, implies ∑j=1pnznj→z=−log⁡G(x)\sum^{p_{n}}_{j=1}z_{nj}\to z=-\log G(x). Once again we have from Lemma 3.1 that

This completes the proof of the lemma. ■\blacksquare

Let ZnZ_{n} be nonnegative random variables such that (Zn2−βn)/αn(Z_{n}^{2}-\beta_{n})/\alpha_{n} converges weakly to a cdf G(x)G(x), where αn>0\alpha_{n}>0 and βn>0\beta_{n}>0 are constants satisfying that lim⁡n→∞αn/βn=0\lim_{n\to\infty}\alpha_{n}/\beta_{n}=0. Then

Proof. Set Wn=(Zn2−βn)/αnW_{n}=(Z_{n}^{2}-\beta_{n})/\alpha_{n}. We have Zn2=βn+αnWn=βn(1+αnβnWn)Z_{n}^{2}=\beta_{n}+\alpha_{n}W_{n}=\beta_{n}(1+\frac{\alpha_{n}}{\beta_{n}}W_{n}). Then by Taylor’s expansion

where for i=1,2i=1,2, li(t)l_{i}(t) is a polynomial in tt of degree ≤3i\leq 3i, depending on rr and kk, and all of its coefficients are of order O(\big{(}\frac{r}{(r-k)k}\big{)}^{i/2}).

Define Vpn−j+1:n−jV_{p_{n}-j+1:n-j} as in (3.25). Assume that kn=n−pn→∞k_{n}=n-p_{n}\to\infty and kn/n→0k_{n}/n\to 0 as n→∞n\to\infty. Then for any δn>0\delta_{n}>0 such that δn→∞\delta_{n}\to\infty and δn=o(kn1/6)\delta_{n}=o(k_{n}^{1/6})

uniformly over 0≤x≤δn0\leq x\leq\delta_{n}, 1≤j≤pn−kn1\leq j\leq p_{n}-k_{n} as n→∞n\to\infty.

Proof. Set βnj(x)=pn−jn−j+((pn−j)kn)1/2(n−j)3/2x=pn−jn−j(1+kn1/2(n−j)1/2(pn−j)1/2x)\beta_{nj}(x)=\frac{p_{n}-j}{n-j}+\frac{((p_{n}-j)k_{n})^{1/2}}{(n-j)^{3/2}}x=\frac{p_{n}-j}{n-j}(1+\frac{k_{n}^{1/2}}{(n-j)^{1/2}(p_{n}-j)^{1/2}}x). Then 1−βnj(x)=n−pnn−j−((pn−j)kn)1/2(n−j)3/2x=knn−j(1−(pn−j)1/2(n−j)1/2kn1/2x)1-\beta_{nj}(x)=\frac{n-p_{n}}{n-j}-\frac{((p_{n}-j)k_{n})^{1/2}}{(n-j)^{3/2}}x=\frac{k_{n}}{n-j}(1-\frac{(p_{n}-j)^{1/2}}{(n-j)^{1/2}k_{n}^{1/2}}x), and the density function of Vpn−j+1:n−jV_{p_{n}-j+1:n-j} is given by

To estimate hj(x)h_{j}(x), we need Stirling’s formula:

and Taylor’s expansion: 1−t=exp⁡(log⁡(1−t))=exp⁡(−t−12t2+O(t3))1-t=\exp(\log(1-t))=\exp(-t-\frac{1}{2}t^{2}+O(t^{3})) as t→0t\to 0. By applying Stirling’s formula to kn!k_{n}!, (pn−j)!(p_{n}-j)! and (n−j)!(n-j)!, the product in (3.16) is equal to 1+o(1)2π(1+o(1))\frac{1+o(1)}{\sqrt{2\pi}}(1+o(1)) for ∣x∣=o(kn1/6)|x|=o(k_{n}^{1/6}) uniformly over 1≤j≤pn−kn1\leq j\leq p_{n}-k_{n} as n→∞n\to\infty. By applying Taylor’s expansion to (1+kn1/2(n−j)1/2(pn−j)1/2x)pn−j(1+\frac{k_{n}^{1/2}}{(n-j)^{1/2}(p_{n}-j)^{1/2}}x)^{p_{n}-j} and (1−(pn−j)1/2(n−j)1/2kn1/2x)kn(1-\frac{(p_{n}-j)^{1/2}}{(n-j)^{1/2}k_{n}^{1/2}}x)^{k_{n}}, the product in (3.16) is equal to

for ∣x∣=o(kn1/6)|x|=o(k_{n}^{1/6}) uniformly over 1≤j≤pn−kn1\leq j\leq p_{n}-k_{n} as n→∞n\to\infty. Therefore, we have

uniformly over 1≤j≤pn−kn1\leq j\leq p_{n}-k_{n} as n→∞n\to\infty. Next we will give an estimate of the upper bound of hj(x)h_{j}(x) for large xx. Note that βnj(x)<1\beta_{nj}(x)<1 if and only if x<kn1/2(pn−j)1/2(n−j)1/2=O(kn1/2)x<\frac{k_{n}^{1/2}(p_{n}-j)^{1/2}}{(n-j)^{1/2}}=O(k_{n}^{1/2}) uniformly over 1≤j≤pn−kn1\leq j\leq p_{n}-k_{n} and thus kn1/2x(n−j)1/2(pn−j)1/2≤O(kn1/2n1/2)→0\frac{k_{n}^{1/2}x}{(n-j)^{1/2}(p_{n}-j)^{1/2}}\leq O(\frac{k_{n}^{1/2}}{n^{1/2}})\to 0 uniformly over 0<x<kn1/2(pn−j)1/2(n−j)1/20<x<\frac{k_{n}^{1/2}(p_{n}-j)^{1/2}}{(n-j)^{1/2}}, 1≤j≤pn−kn1\leq j\leq p_{n}-k_{n} as n→∞n\to\infty. Now by applying Taylor’s expansion to (1+kn1/2(n−j)1/2(pn−j)1/2x)pn−j(1+\frac{k_{n}^{1/2}}{(n-j)^{1/2}(p_{n}-j)^{1/2}}x)^{p_{n}-j} and inequality

to (1−(pn−j)1/2(n−j)1/2kn1/2x)kn(1-\frac{(p_{n}-j)^{1/2}}{(n-j)^{1/2}k_{n}^{1/2}}x)^{k_{n}} we get

uniformly over 1≤j≤pn−kn1\leq j\leq p_{n}-k_{n} as n→∞n\to\infty.

Assume that 0≤x≤δn0\leq x\leq\delta_{n}. From (3.19) we have

uniformly over 0≤x≤δn0\leq x\leq\delta_{n}, 1≤j≤pn−kn1\leq j\leq p_{n}-k_{n} as n→∞n\to\infty. Therefore, to complete the proof of the lemma, it suffices to show that

uniformly over 0≤x≤δn0\leq x\leq\delta_{n}, 1≤j≤pn−kn1\leq j\leq p_{n}-k_{n} as n→∞n\to\infty.

uniformly if y→∞y\to\infty. Since (δnkn1/6)1/2=o(kn1/6)(\delta_{n}k_{n}^{1/6})^{1/2}=o(k_{n}^{1/6}), by using (3.18) we have

proving (3.20). This completes the proof of the lemma. ■\blacksquare

2 Proofs of the Theorems

where An=cn+12(1−cn2)1/2(n−1)−1/2anA_{n}=c_{n}+\frac{1}{2}(1-c_{n}^{2})^{1/2}(n-1)^{-1/2}a_{n}, Bn=12(1−cn2)1/2(n−1)−1/2bnB_{n}=\frac{1}{2}(1-c_{n}^{2})^{1/2}(n-1)^{-1/2}b_{n},

with a(x)=(log⁡x)1/2−(log⁡x)−1/2log⁡(2πlog⁡x)a(x)=(\log x)^{1/2}-(\log x)^{-1/2}\log(\sqrt{2\pi}\log x) and b(x)=(log⁡x)−1/2b(x)=(\log x)^{-1/2} for x>3x>3.

Since Ynj2Y_{nj}^{2} and Upn−j+1:n−jU_{p_{n}-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]j_{n}=[p_{n}^{5/8}], the integer part of pn5/8p_{n}^{5/8}. For 1≤j≤jn1\leq j\leq j_{n}, define

Then we see that uniformly over 1≤j≤jn1\leq j\leq j_{n},

uniformly for all 1≤j≤jn1\leq j\leq j_{n} as n→∞n\to\infty.

In Lemma 3.6, take r=n−jr=n-j and k=n−pnk=n-p_{n} to have

uniformly over 1≤j≤jn1\leq j\leq j_{n} as n→∞n\to\infty, where

and where, for i=1,2i=1,2, li(t)l_{i}(t) is a polynomial in tt of degree ≤3i\leq 3i, depending on nn, and all of its coefficients are of order O((1/pn)i/2)O((1/p_{n})^{i/2}) uniformly over 1≤j≤jn1\leq j\leq j_{n} as n→∞n\to\infty. Now, by taking B=(unj,∞)B=(u_{nj},\infty) we obtain

uniformly for 1≤j≤jn1\leq j\leq j_{n} as n→∞n\to\infty. From L’Hospital’s rule, we have that for any r≥0r\geq 0

Since min⁡1≤j≤jnunj→∞\min_{1\leq j\leq j_{n}}u_{nj}\to\infty as n→∞n\to\infty by (3.24), it follows from (3.27) that

holds uniformly over 1≤j≤jn1\leq j\leq j_{n}. Furthermore, since the coefficients of li(t)l_{i}(t) are uniformly bounded by O((1/pn)i/2)O((1/p_{n})^{i/2}) for i=1,2i=1,2, we have

In Lemma 3.5, by taking xn=ncn2/(1−cn2)x_{n}=nc_{n}^{2}/(1-c_{n}^{2}) and and define cn,jc_{n,j} such that unj=(j−1)cn,j+a(xn)+b(xn)xu_{nj}=(j-1)c_{n,j}+a(x_{n})+b(x_{n})x for 1≤j≤jn1\leq j\leq j_{n} with cn,1=xn−1/2c_{n,1}=x_{n}^{-1/2}. It follows from (3.24) that cn,j=xn−1/2(1+o(1))c_{n,j}=x_{n}^{-1/2}(1+o(1)) uniformly over 1≤j≤jn1\leq j\leq j_{n} as n→∞n\to\infty, which implies lim⁡n→∞max⁡1≤j≤jn∣cn,jxn1/2−1∣=0\lim_{n\to\infty}\max_{1\leq j\leq j_{n}}|c_{n,j}x_{n}^{1/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→∞n\to\infty. By noting that (1−cn2)1/2n1/2cn∼pn−1/2\frac{(1-c_{n}^{2})^{1/2}}{n^{1/2}c_{n}}\sim p_{n}^{-1/2} we have

uniformly over 1≤j≤jn1\leq j\leq j_{n}. Thus, it follows from (3.13) that

as n→∞n\to\infty. Therefore, by combining (3.2), (3.26) and (3.29) we get

for all large nn. Then we estimate znjnz_{nj_{n}} by using (3.26) with j=jnj=j_{n} and (3.30)

as n→∞n\to\infty. From (3.8) we have ∑j=jn+1pnznj≤pnznjn=O(pn−1/2)→0\sum^{p_{n}}_{j=j_{n}+1}z_{nj}\leq p_{n}z_{nj_{n}}=O(p_{n}^{-1/2})\to 0 as n→∞n\to\infty, which together with (3.31) yields

We can also prove from (3.26) that zn1→0z_{n1}\to 0 as n→∞n\to\infty. 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\alpha_{n}=c_{n}(1-c_{n}^{2})^{1/2}(n-1)^{-1/2}b_{n} and βn=cn2+cn(1−cn2)1/2(n−1)−1/2an\beta_{n}=c_{n}^{2}+c_{n}(1-c_{n}^{2})^{1/2}(n-1)^{-1/2}a_{n}. Since

as n→∞n\to\infty, we can apply Lemma 3.4 and get that

Recall that An=cn+12(1−cn2)1/2(n−1)−1/2an=cn(1+o(1))A_{n}=c_{n}+\frac{1}{2}(1-c_{n}^{2})^{1/2}(n-1)^{-1/2}a_{n}=c_{n}(1+o(1)) and Bn=12(1−cn2)1/2(n−1)−1/2bn∼12(n−1)−1/2(log⁡pn)−1/2B_{n}=\frac{1}{2}(1-c_{n}^{2})^{1/2}(n-1)^{-1/2}b_{n}\sim\frac{1}{2}(n-1)^{-1/2}(\log p_{n})^{-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→∞n\to\infty. Noting that n≤n2/kn≤n2n\leq n^{2}/k_{n}\leq n^{2}, we get that

is of order (log⁡n)1/2(\log n)^{1/2} as n→∞n\to\infty.

Use the same notation as in Part 1. Recall that pn/n→1p_{n}/n\to 1, kn=n−pn=o(n)k_{n}=n-p_{n}=o(n) and log⁡n=o(kn1/3)\log n=o(k_{n}^{1/3}) as n→∞n\to\infty. In order to use both Lemmas 3.5 and 3.7, we take xn=ncn21−cn2x_{n}=\frac{nc_{n}^{2}}{1-c_{n}^{2}}. Define jn=[5(log⁡n)1/2xn]+1j_{n}=[5(\log n)^{1/2}\sqrt{x_{n}}]+1. Then jn∼5n(log⁡n)1/2kn=o(n)j_{n}\sim\frac{5n(\log n)^{1/2}}{\sqrt{k_{n}}}=o(n), which implies 1≤jn≤pn−kn1\leq j_{n}\leq p_{n}-k_{n} for all large nn. Define cn,jc_{n,j} for 1≤j≤jn1\leq j\leq j_{n} in the same way as in Part 1. Similar to the proof of (3.24) we can show that

uniformly for 1≤j≤jn1\leq j\leq j_{n} as n→∞n\to\infty. Then unj=O((log⁡n)1/2)=o(kn1/6)u_{nj}=O((\log n)^{1/2})=o(k_{n}^{1/6}) uniformly for 1≤j≤jn1\leq j\leq j_{n} as n→∞n\to\infty. 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−kn1\leq j\leq p_{n}-k_{n} and unjn≥4(log⁡n)1/2u_{nj_{n}}\geq 4(\log n)^{1/2} for all large nn. By employing (3.14) with j=jnj=j_{n} and x=4(log⁡n)1/2x=4(\log n)^{1/2} and using equation (3.8) and Lemma 3.7 we have

Thus, we obtain that ∑j=1pnznj→e−x\sum^{p_{n}}_{j=1}z_{nj}\to e^{-x} for any xx. 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. ■\blacksquare

Proof of Theorem 3. Recall that ana_{n} is given by

By using Stirling’s formula (3.17), we have under condition (2.3) that

Since te−tte^{-t} is strictly increasing in (0,1)(0,1), for all large nn such that yn<1y_{n}<1 define εn\varepsilon_{n} as the unique solution to te−t=ynte^{-t}=y_{n} in (0,1)(0,1), that is, εne−εn=yn\varepsilon_{n}e^{-\varepsilon_{n}}=y_{n}, which implies that εn→0\varepsilon_{n}\to 0 and εn∼yn\varepsilon_{n}\sim y_{n} as n→∞n\to\infty and

for all large nn. Then it follows from the the first inequality in (3.37) that

for all large nn, which together with (3.36) implies that an≤knεna_{n}\leq k_{n}\varepsilon_{n} for all large nn,and thus an=o(kn)a_{n}=o(k_{n}) as n→∞n\to\infty. By plugging y=any=a_{n} in (3.37) and using an≤knεna_{n}\leq k_{n}\varepsilon_{n} for large nn we conclude

Define znjz_{nj} as in (3.7) with tn=tn(x)=1−ann(1−xkn)t_{n}=t_{n}(x)=1-\frac{a_{n}}{n}(1-\frac{x}{k_{n}}) for any fixed xx. Then 1−tn=ann(1−xkn)=o(knn)=o(log⁡nn)1-t_{n}=\frac{a_{n}}{n}(1-\frac{x}{k_{n}})=o(\frac{k_{n}}{n})=o(\frac{\log n}{n}) as n→∞n\to\infty, where we have used the fact that kn=n−pn→∞k_{n}=n-p_{n}\to\infty and kn=o(log⁡n)k_{n}=o(\log n) from (2.3). This implies n(1−tn)2→0n(1-t_{n})^{2}\to 0 as n→∞n\to\infty.

It is easy to verify the following expression

where 0≤d(t)≤t20\leq d(t)\leq t^{2} for 0≤t≤1/20\leq t\leq 1/2. Then

Furthermore, by using Stirling’s formula (3.17), we get

uniformly over 1≤j≤jn1\leq j\leq j_{n}, where jn:=pn−kn3j_{n}:=p_{n}-k_{n}^{3}. Then from (3.9) we obtain that

Note that n(1−tn)=o(kn)n(1-t_{n})=o(k_{n}). Then it follows from (3.39) and (3.37) that

uniformly over 1≤j≤jn1\leq j\leq j_{n}, and thus

The second integral above is dominated by the first one since tkne−tt^{k_{n}}e^{-t} is increasing over (0,kn)(0,k_{n}) and kn3(1−tn)/(n(1−tn))=kn3/n→0k_{n}^{3}(1-t_{n})/(n(1-t_{n}))=k_{n}^{3}/n\to 0 as n→∞n\to\infty. Therefore, in view of (3.37) and (3.38) we have

as n→∞n\to\infty. Therefore, it follows from Lemma 3.2 that ∑j=1pnznj→e−x\sum^{p_{n}}_{j=1}z_{nj}\to e^{-x} as n→∞n\to\infty. It is easy to conclude that zn1→0z_{n1}\to 0 as n→∞n\to\infty from (3.39), (3.37) and the above estimates. Accordingly, by taking βn=1−ann\beta_{n}=1-\frac{a_{n}}{n} and αn=annkn\alpha_{n}=\frac{a_{n}}{nk_{n}} with G(x)=Λ(x)G(x)=\Lambda(x) in Lemmas 3.3 and 3.4 we conclude that

as n→∞n\to\infty. This completes the proof of the theorem. ■\blacksquare

Fix x<0x<0. Let tn=tn(x)=1+((k+1)!)1/kn1+1/kxt_{n}=t_{n}(x)=1+\frac{((k+1)!)^{1/k}}{n^{1+1/k}}x. Then tn∈(0,1)t_{n}\in(0,1) for all large nn. Since n(1−tn)→0n(1-t_{n})\to 0 as n→∞n\to\infty, we have

Therefore, we have from (3.7) and (3.9) that

uniformly over 1≤j≤pn=n−k1\leq j\leq p_{n}=n-k. Since

we have max⁡1≤j≤pnznj=O(1/n)→0\max_{1\leq j\leq p_{n}}z_{nj}=O(1/n)\to 0 as n→∞n\to\infty. To complete the proof of (3.40), by using (3.5) we need to show that

Let {jn}\{j_{n}\} be a sequence of integers such that jn→∞j_{n}\to\infty and jn/n→0j_{n}/n\to 0 as n→∞n\to\infty. Then

uniformly over 1≤j≤n−jn1\leq j\leq n-j_{n}, which implies that

as n→∞n\to\infty. This proves (3.41) and thus we obtain (3.40).

Finally, the theorem follows from Lemma 3.4 with αn=((k+1)!)1/kn1+1/k\alpha_{n}=\frac{((k+1)!)^{1/k}}{n^{1+1/k}} and βn=1\beta_{n}=1. This completes the proof. ■\blacksquare

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).

References