arxiv: v1 [math.pr] 15 May 2018

Transcription

1 Ferromagnetic to paramagnetic transition in spherical spin glass Jinho Baik Ji Oon Lee Hao Wu October 8, 08 arxiv: v [math.pr] 5 May 08 Abstract We consider the spherical spin glass model defined by a combination of the pure -spin spherical Sherrington-Kirkpatrick Hamiltonian and the ferromagnetic Curie-Weiss Hamiltonian. In the large system limit, there is a two-dimensional phase diagram with respect to the temperature and the coupling strength. The phase diagram is divided into three regimes; ferromagnetic, paramagnetic, and spin glass regimes. The fluctuations of the free energy are known in each regime. In this paper, we study the transition between the ferromagnetic regime and the paramagnetic regime in a critical scale. Introduction We consider a disordered system defined by random Gibbs measures whose Hamiltonian is the sum of a spin glass Hamiltonian and a ferromagnetic Hamiltonian. Depending on the strength of the coupling constant and the temperature, the system may exhibit several phases in the large system limit. The paper is concerned with the fluctuations of the free energy near the boundary between two phases known as ferromagnetic and paramagnetic regimes. Consider the sum of the pure -spin spherical Sherrington-Kirkpatrick SSK Hamiltonian and the Curie-Weiss CW Hamiltonian. We call this sum the SSK+CW Hamiltonian. We denote the coupling constant by J and the inverse temperature by β. We consider the random Gibbs measure with the SSK+CW Hamiltonian. The focus of this paper is on the free energy. The limiting free energy was obtained non-rigorously by Kosterlitz, Thouless, and Jones [] in 976. When J = 0, this formula is the explicit evaluation of the Crisanti Sommers formula [5] which was proved rigorously by Talagrand [6] in the case of the pure -spin SSK. The Crisanti Sommers formula is the spherical version of the Parisi formula [4, 7]. The formula of Kosterlitz, Thouless, and Jones shows a two-dimensional phase transition: see Figure. The three regimes are determined by the condition that max{, β, J} is equal to spin glass regime, β paramagnetic Department of Mathematics, University of Michigan, Ann Arbor, MI, 4809, USA baik@umich.edu Department of Mathematical Sciences, KAIST, Daejeon, 344, Korea jioon.lee@kaist.edu Department of Mathematics, University of Michigan, Ann Arbor, MI, 4809, USA lingluan@umich.edu

2 β paramagnetic spin glass ferromagnetic 0 J Figure : Phase diagram for SSK+CW model. Here, β is the inverse temperature and J is the coupling constant. regime or J ferromagnetic regime. The limiting free energy is analytic with respect to both β and J in each regime, but not on the boundary. Recently, the authors of [6] showed that the result of Kosterlitz, Thouless, and Jones is rigorous. Furthermore, the authors also evaluated the distribution of the fluctuations of the free energy in each regime. The case when J = 0 was obtained earlier in [4]. The order of the fluctuations are /3,, / and the limiting distributions are Tracy-Widom, Gaussian, and Gaussian in the spin glass, paramagnetic regime, ferromagnetic regime, respectively. In the same paper, the transition between the spin glass regime and the ferromagnetic regime was also studied. However, the other two transitions and the triple point were left open. The goal of this paper is to describe the transition between the paramagnetic regime and and the ferromagnetic regime. Another system which combines a spin glass and a ferromagnetic model is the SSK with an external field. The difference between the CW Hamiltonian and an external field is that one is a quadratic function and the other is a linear function of the spin variables. These two models are related; see [] for a one-sided inequality. For the spin glass with external field, the fluctuations of the free energy were computed recently in [3, 4] when the coupling constant is positive for both SSK and SK Sherrington-Kirkpatrick cases with general spin interactions. However, the transitions are not obtained except for certain large deviation results [8, 6]. One of the interests of the SSK+CW model is that it is an easier model which can be analyzed in detail in the transitional regimes.. Model Let S = {σ = σ,, σ R : σ + + σ = }. be a sphere in R of radius. Define the SSK+CW Hamiltonian by H σ = H SSK σ + H CW σ, σ S.

3 where H SSK σ = i,j= A ij σ i σ j, H CW σ = J σ i σ j = J σ i..3 Here J is the coupling constant. The random coefficients A ij satisfy A ij = A ji and A ij, i j, are independent centered random variables. We call A ij disorder variables. The precise conditions are given in Definition. below. ote that as a function of σ, H CW σ is large when the coordinates of σ have same sign. On the other hand, the maximizers σ of H SSKσ depend highly on {A ij}. With β > 0 representing the inverse temperature, the free energy and the partition function are defined by F = log Z, Z = e βh σ dω σ.4 S i,j= i= where ω is the normalized uniform measure on S. ote that F and Z are random variables since they depend on the disorder variables A ij. The free energy and the partition function depend on the parameters β and J, F = F β, J, Z = Z β, J..5 Since the Curie-Weiss Hamiltonian is a quadratic function of the spin variable, we can write the SSK+CW Hamiltonian as H σ = i,j= M ijσ i σ j where M ij = A ij + J are non-centered random variables. In terms of matrix notations, H σ = σ T Mσ, M = A + J T.6 with A = A ij i,j, =,, T, M = M ij i,j, and σ = σ,, σ T. The noncentered random symmetric matrix M is an example of a real Wigner matrix perturbed by a deterministic finite rank matrix. Such matrices are often called spiked random matrices. We will use the eigenvalues of spiked random matrices in our analysis of the free energy. We assume the following conditions on the disorder variables. Definition. Assumptions on disorder variables. Let A ij, i j, be independent real random variables satisfying the following conditions: All moments of A ij are finite and E[A ij ] = 0 for all i j. For all i < j, E[A ij ] =, E[A3 ij ] = W 3, and E[A 4 ij ] = W 4 for some constants W 3 R and W 4 0. For all i, E[A ii ] = w for a constant w 0. Set A ij = A ji for i > j. Let A = A ij i,j= and we call it a Wigner matrix of zero mean. 3

4 Definition. Eigenvalues of non-zero mean Wigner matrices. Let M be the symmetric matrix defined in.6. We call it a Wigner matrix of non-zero mean. Its eigenvalues are denoted by λ λ λ..7 We introduce the following terminology. Definition.3 High probability event. We say that an -dependent event Ω holds with high probability if, for any given D > 0, there exists 0 > 0 such that for any 0.. Previous results in each regime PΩ c D We review the results on the fluctuations in each regime obtained in [6]. We state two types of results: one in terms of the eigenvalues of M and the other in terms of limiting distributions. Set J := max{j, }..8 It was shown in [6] that the following holds with high probability. In both ferromagnetic and the spin glass regimes given by J > β, with any ɛ > 0, F = F + β J λ J J + o +ɛ..9 In the paramagnetic regime given by J < β, F = F i= log β + β λ i + o..0 Here, F is a deterministic function of, β, J. The above results show that the fluctuations of F are determined, to the leading order, by the top eigenvalue λ in the ferromagnetic and spin glass regimes, while they are determined by all eigenvalues in the paramagnetic regime. A limit theorem for F follows if we use limit theorems for the eigenvalues of random matrices. The relevant random matrices are Wigner matrices of non-zero mean in.6. For such random matrices, the following is known [5, ] see [3] for complex matrices: { /3 λ TW if J <, / λ J J W3 J J 4, J if J >,. In [6], we consider the case when the diagonal entries of M have mean J and the off-diagonal entries have mean J where J and J are allowed to be different. However, in this case, M = + J T + J J I where I is the identity matrix. This only shifts all eigenvalues by a deterministic small number. As we will see in Remark., it is not more general than the case with J = J. 4

5 where the convergences are in distribution. Here TW denotes the GOE Tracy-Widom distribution and a, b denotes the Gaussian distribution of mean a and variance b. The dichotomy is due to the effect of the non-zero mean; if J is not large enough i.e. J <, then the influence of the non-zero mean is negligible to contribute to the fluctuations of the top eigenvalue. For J <, the top eigenvalue is close to the second eigenvalue with order O /3+ɛ. But for J >, the difference of the top eigenvalue and the second eigenvalue is of order O. On the other hand, the following is also known see Theorem.6 of [6]: if a function ϕ is smooth in an open interval containing the interval [, J + J ], then ϕλ i i= ϕxdσ scl x f, a, dσ scl x := 4 x dx,. π for some explicit constants f, a. This result is applicable to the paramagnetic regime. Together, we have the following asymptotic results obtained in Theorem.4 of [6] with a small correction in [5]: i Spin glass regime If β > and J <, then ii Paramagnetic regime If β < and β < J, then β /3 F F TW..3 F F f, α..4 iii Ferromagnetic regime If J > and β > J, then F F f, α..5 for some deterministic function F = F β, J and some explicit constants f, α, f and α depending on β and J..3 Results We state the results on the transition between the paramagnetic regime and the ferromagnetic regime. The boundary between these two regimes is given by the equation β = J with J >. In the transitional regime, the correct scaling turns out to be the following: let J > be fixed and let β = β be given by β = J + B.6 with fixed B R. The following is the first main result of this paper. This relates the free energy with the eigenvalues of M. 5

6 Theorem.4. Let β be given by.6. Then, for every 0 < ɛ < 8, F = F gλ i + Qχ + O 3/+4ɛ, χ := λ J + J,.7 i= with high probability as, where F = βj + J logβ + Also, with and 4 log + log β, gz := logj + J z..8 π sx Qx = sx x sx logsx x sx x 4J + + log I J sx = x BJ + x + BJ + 4J Iα = where the square root denotes the principal branch..9.0 e α 4 t + it + it dt,. The formula.7 shows a combined contribution from λ,, λ and a distinguished contribution from λ. Compare the formula with.9 and.0. ow we state a result analogous to.4 and.5. This follows if we have limit theorems for Qχ and i= gλ i. From the second part of., Qχ converges to an explicit function of a Gaussian random variable. On the other hand, i= gλ i is different from i= gλ i by one term. It is not difficult to show that removing one term does not affect the fluctuations much and the fluctuations are still given by a Gaussian random variable similar to.; see Theorem. in the next section. In random matrix theory, these sums are known as partial linear statistic and linear statistic, respectively. The main technical part of this paper is to evaluate the joint distribution of Qχ and i= gλ i. We show that jointly they converge in distribution to a bivariate Gaussian variable with an explicit covariance. See the next section for the precise statement. These results are interesting on their own in random matrix theory. Putting together, we obtain the following result. Theorem.5. We have F 4J B J log 4 B J G + QG. 4 in distribution as where G and G are bivariate Gaussian random variables with E[G ] = 4 logj + w 4J + W 4 3 8J 4 + log πj,.3 Var[G ] = log J + w 4J + W 4 3 8J 4,.4 6

7 a pdf of QG b pdf of normalized QG Figure : a Probability density function of QG for B =, 0,, b Probability density function of normalized QG resembles a Gaussian density as B +. and E[G ] = W 3 J J 4, Var[G ] = J,.5 CovG, G = W 3J J 4..6 ote that G and G do not depend on B. The function Q is defined in.9. ote that if the third moment W 3 of A ij with i j is zero, then G and G are independent Gaussians. The above result is consistent with the results on ferromagnetic and paramgnatic regimes if we let formally B + and B, respectively. One can show that when B +, QG dominates G. Furthermore, while QG is not Gaussian, upon proper normalization, it converges to a Gaussian as B +. See Figure. On the other hand, when B, the leading two terms of QG are constants and the random part is smaller than G. See Section 6 for details. Let us comment on the other transitions in the phase digram in Figure. As mentioned before, the transition between the spin glass and ferromagnetic regimes was discussed in [6]. ote that.9 is valid in both regimes. It was shown that if we let β > / be fixed and consider -dependent J = + w /3, then for each w R,.9 still holds. ow, for such J, it was shown in [9] that /3 λ TW,w where TW,w is a one-parameter family of random variables interpolating TW and Gaussian distributions. Hence, we obtain the fluctuations for the transitional regime. On the other hand, the transition between the spin glass and paramagnetic regimes is an open question. By matching the fluctuation scales in both regimes, we expect that the critical scale is β = + O log /3..4 Organization The rest of the paper is organized as follows. In Section, we first state new results on random matrices. They are given in Theorem. partial linear statistics and Theorem.3 joint convergence. Using them, we derive Theorem.5 from Theorem.4. In Section 3, we prove Theorem 7

8 .4. In the next two sections, we prove the random matrix results stated in Section ; Theorem. in Section 4 and Theorem.3 in Section 5. In Section 6, we show that Theorem.5 is consistent with the previous results on ferromagnetic and paramagnetic regimes. Acknowledgments The work of Jinho Baik was supported in part by SF grant DMS and the Simons Fellows program. The work of Ji Oon Lee was supported in part by Samsung Science and Technology Foundation project number SSTF-BA Results on Wigner matrices with non-zero mean In order to prove Theorem.5 from Theorem.4, we need some new results on random matrices. We need i a limit theorem for partial linear statistics i= gλ i and ii a joint convergence of the large eigenvalue and partial linear statistics. These results are interesting on their own in random matrix theory. We state them here and prove them in Section 4 and Section 5 below. Using these results, we prove Theorem.5 in Subsection.3. Recall that the symmetric matrix M is given by M = A+ J T where A = A ij is a symmetric matrix with independent entries for i j satisfying the conditions given in Definition. and =,, T. The matrix M is called a Wigner matrix with a non-zero mean J. Recall that we assume J >.. The eigenvalues of M are denoted by λ λ. It is known that λ is close to J+J with high probability and λ,, λ are in a neighborhood of [, ] with high probability. See Lemma 3. below for the precise statement.. Partial linear statistics A linear statistic is the sum of a function of the eigenvalues. The fluctuations of linear statistics for Wigner matrices and other random matrix ensembles are of central interest in the random matrix theory; see, for example, [9,, ]. For Wigner matrices with non-zero mean, the following result was obtained in Theorem.6 and Remark.7 of [6]. Set Ĵ = J + J.. Let ϕ : R R be a function which is analytic in an open neighborhood of [, Ĵ] and has compact support. Then, as, the random variable ϕ := ϕλ i i= ϕxdσ scl x Mϕ, V ϕ.3 8

9 where Mϕ = 4 ϕ + ϕ 3 τ 0ϕ J τ ϕ + w τ ϕ + W 4 3τ 4 ϕ + ϕĵ J l τ l ϕ, l= V ϕ =w τ ϕ + W 4 3τ ϕ + Here, W 4 = E[A 4 ], w = E[A ], and lτ l ϕ. l=.4 τ l ϕ = π where T l t are the Chebyshev polynomials of the first kind. ϕx T lx/ dx = π ϕ cosθ coslθdθ,.5 4 x π π We are interested in a partial linear statistic, i= ϕλ i. See [7, 3] for other types of partial linear statistics. The partial linear static i= ϕλ i is the linear statistic minus one term ϕλ. Since λ Ĵ in probability see the second part of., by.3, Slutsky s theorem implies that ϕλ i ϕxdσ scl x Mϕ ϕĵ, V ϕ. i= Since this follows from.3, this is true assuming that ϕ is analytic in an open neighborhood of [, Ĵ]. However, we are interested in the test function ϕx = gx = logĵ x see.7. Since this function is not analytic at x = Ĵ, the above simple argument does not apply. onetheless, if we adapt the proof of.3, one can show that it is enough to assume that the test function is analytic in a neighborhood of the interval [, ], not of [, Ĵ]. Theorem.. Let J >. Then for every test function ϕ which is analytic in a neighborhood of [, ], ϕ := ϕλ i ϕxdσ scl x M ϕ, V ϕ.6 as with i= M ϕ = 4 ϕ + ϕ 3 τ 0ϕ J τ ϕ + w τ ϕ + W 4 3τ 4 ϕ J l τ l ϕ, l=.7 and V ϕ = V ϕ where V ϕ is defined in.4. ote that M ϕ = Mϕ ϕĵ.8 for ϕ analytic in a neighborhood of [, Ĵ]. 9

10 Remark.. We comment on a case when the test function depends on. Consider the function ϕ defined by ϕ x = ϕx + φx + O uniformly for x in a neighborhood of [, ] for analytic functions ϕ and φ. Define the corresponding linear statistic ϕ = i= ϕ λ i ϕ xdσ scl x, then ϕ = ϕ λ i i= = ϕ + ϕ xdσ scl x φλ i i= φxdσ scl x + O..9 By Theorem., the second order term converges to zero in probability. Thus, ϕ and ϕ converge to the same Gaussian distribution. The same argument also applies to full linear statistics; this is used in Remark 5.4 below. ow, the claim in footnote is verified by noting that ϕx + J J = ϕx + ϕ xj J + O.. Joint convergence of the largest eigenvalue and linear statistics By Theorem. and the second part of., the partial linear statistic and the largest eigenvalue each converge to Gaussian distributions individually. The following theorem shows that they converge jointly to a bivariate Gaussian with an explicit covariance. Theorem.3. Let J >. Then for ϕx which is analytic in a neighborhood of [, ], ϕ := i= ϕλ i ϕxdσ sclx and χ := λ Ĵ converges jointly in distribution to a bivariate Gaussian variable with mean M ϕ, W 3 J J 4.0 and covariance V ϕ W 3 τ ϕ J W 3 τ ϕ J J.. The proof of this theorem, given in Section 5, is the main technical part of this paper. We prove the theorem first for the Gaussian case, and then use an interpolation argument..3 Proof of Theorem.5 We now derive Theorem.5 from Theorem.4 using the results on the eigenvalues stated in the previous two subsections. The term Qχ converges to QG in distribution from Theorem.3. Consider the rest. It was shown in A.5 of [4] that for gz = logj + J z, gzdσ scl x = + log J.. J 0

11 Inserting β = J +B / and using the Taylor expansion log+ BJ = BJ B J +O 3/, F gzdσ scl x = 4J + B J + log 4 + [ B J 4 + log ] + O 3/..3 πj We can evaluate M g using.7 of [6] which evaluated the Mh with hx = logβ + β x: note that J = J here M g = lim β J Mh logβ + β J J = logj w J W 4 3 4J 4..4 The variance V g = V g, which is independent of J, is given by 4 times 3.3 of [4] if we replace β by J : V g = log J + J w + J 4 W For the covariance term, we have τ g = J.3, we obtain the result. from A.7 of [4]. Hence, from Theorem. and 3 Proof of Theorem.4 The proof follows the steps for the proof of the Theorem.5 of [6] for paramagnetic and ferromagnetic regimes with necessary adjustments. The analysis is based on applying a method of steepest-descent to a random integral. The location of the critical point is important. In the transitional regime, the critical point is close to the largest eigenvalue but not as close as the ferromagnetic case. On the other hand, the critical point is away from the largest eigenvalue in the paramagnetic case. See Subsection 3. below for details. 3. Preliminaries The following formula is a simple result in []. Lemma 3. []; also Lemma.3 of [4]. Let M be a real symmetric matrix with eigenvalue λ λ λ. Then for fixed β > 0, S e βσt Mσ dw σ = C γ+i γ i e Gz dz, Gz = βz logz λ i, 3. where γ is any constant satisfying γ > λ, the integration contour is the vertical line from γ i to γ + i, the log function is defined in the principal branch, and Here Γz denotes the Gamma function. C = i= Γ/. 3. πiβ /

12 Let M be a Wigner matrix with non-zero mean as in.6. Then its eigenvalues λ i are random variables, and hence the above result gives a random integral representation of the partition function. In [6, 4], the above random integral was evaluated using the method of steepest-descent for different choices of random matrices. The key ingredient in controlling the error term is a precise estimate for the eigenvalues which are obtained in the random matrix theory. Lemma 3. Rigidity of eigenvalues: Theorem.3 of [7] and Theorem 6.3 of [0]. For each positive integer k [, ], set ˆk := min{k, + k}. Let γ k be the classical location defined by dσ scl x = k. 3.3 γ k Then, for every 0 < ɛ <, λ k γ k ˆk /3 /3+ɛ 3.4 for all k =, 3,, with high probability. Furthermore, for fixed J >, recall Ĵ = J + J, holds with high probability. λ Ĵ /+ɛ 3.5 From the rigidity, it is easy to obtain the following law of large numbers for eigenvalues. Corollary 3.3 c.f. Lemma 5. of [4]. Fix δ > 0, let {f α } α I C [ δ, + δ] be a family of monotonic increasing functions satisfying sup α I max x f α x C 0 and sup α I max x f αx C. Then, for every 0 < ɛ <, with high probability. sup α I i= f α λ i f α xdσ scl x = O +ɛ 3.6 Proof. Let f = f α for some α I. The absolute value on the left hand-side is bounded above by fλ i fγ i + fγ i fxdσ scl x. 3.7 i= i= i= By Lemma 3., fλ i fγ i max f x i= i= λ i γ i C 0 ɛ 3.8 with high probability. On the other hand, set ˆγ j by dσ scl x = j, j =,,,, 3.9 ˆγ j

13 and by convention ˆγ 0 =. As fx is a monotonic increasing function, for i =, 3,,, Thus, ˆγi ˆγ i+ fxdσ scl x fγ i i= fγ i fxdσ scl x ˆγi ˆγ i fxdσ scl x max fx Since the upper bounds are independent of f, we obtain the result. 3C Steepest-descent analysis We now apply steepest descent analysis to the integral in Lemma 3.. We deform the contour to pass a critical point and show that the main contribution to the integral comes from a small neighborhood of the critical point. For Gz given in 3., it is easy to check that all solutions of G z = 0 are real-valued, and there is a unique critical point γ which lies in the interval λ, see Lemma 4. of [6]. ote that since G is random, the critical point is also random. For the paramagnetic regime, it was shown in [6] that γ λ = O with high probability. In the same paper, it was also shown that in the ferromagnetic regime, γ λ = O +ɛ with high probability. The following lemma establishes a corresponding result for the transitional regime; it shows that γ λ = O +ɛ with high probability. Lemma 3.4 Critical point. Recall that see.6 J > is fixed and β = β = J + B with fixed B R. Then, for every 0 < ɛ < 4, γ = λ + χ BJ + χ + J B + 4J + O +ɛ 3. with high probability, where we set χ := λ Ĵ. ote that γ given above is larger than λ with high probability since the term in the big parenthesis is positive. Proof. Set θ := χ BJ + χ + J B + 4J. 3.3 ote that θ > 0. By the rigidity of λ, we have χ ɛ 4 and hence, θ ɛ 3 with high probability. On the other hand, using a + a + b b =, a +b +a θ = J J B + χ + 4J + J B + χ, and hence θ C ɛ 4 for some constant C > 0 with high probability. Hence, ɛ 3 θ ɛ

14 with high probability. Set γ ± := λ + θ ± +ɛ. 3.5 By the above properties of θ, we have γ ± > λ with high probability. We will show that G γ < 0 and G γ + > 0 with high probability. Since G z is a monotone increasing function for real z in the interval λ,, this shows that γ < γ < γ + with high probability, proving the lemma. Recall that λ Ĵ in probability. Let us write γ ± = J + J + φ ± +ɛ, φ := θ + χ 3.6 where χ = λ Ĵ. ote that φ = O ɛ 3 with high probability. ow, notice that G z = β i= z λ i z λ. 3.7 We apply Corollary 3.3 to the family of the function { z x } z>+c for some constant c > 0 and obtain γ ± γ G ± 4 γ ± = β + O + ɛ 3 γ ± λ with high probability. By 3.6, γ ± γ± 4 By 3.5, = J φ J ± +ɛ ɛ + + O 3. γ ± λ = θ +ɛ + O θ +ɛ. Using the formula of β and the estimate 3.4 for θ, we find that G γ ± = B + φ J ± θ J + θ +ɛ ɛ + + O with high probability since 0 < ɛ < 4. By the definition of θ, the leading term is zero. The coefficient of the second term is positive. Hence we find that G γ < 0 and G γ + > 0, and we obtain the lemma. Then we have the following lemma. Lemma 3.5. Set Then, for every ɛ > 0, s = s := γ J J and = := γ λ = s χ. 3.9 s = χ BJ + χ + J B + 4J 4 θ + O +ɛ 3.0

15 with high probability. We also have with high probability. s ɛ and ɛ ɛ 3. Proof. The previous lemma implies 3.0. The first part of 3. follows from the fact that χ = O ɛ with high probability. The second part is the estimate 3.4 in the proof of the previous lemma. We also need the following lemma. Lemma 3.6. For every 0 < ɛ <, with high probability. Proof. This follows from Corollary 3.3 applied to fx = γ x. i= γ λ i = J + O +ɛ 3. The following auxiliary lemma is used to estimate an error in the steepest descent analysis. Lemma 3.7. Define I m α := t m + it e α 4 t + it dt 3.3 for non-negative integers m and α > 0, where the square root is the defined on the principal branch. We set Iα := I 0 α; see.. Then, 4π Iα = α + Oα as α +, 3.4 and for every m 0, Iα = 8π e + Oα as α 0 +, 3.5 I m α is uniformly bounded for α 0,. 3.6 A particular consequence is that the derivative I α = 4 I α is uniformly bounded for α > 0. Furthermore, Iα > 0 for all α > 0. Proof. Consider 3.4. Applying the method of steepest-descent to Iα = gteαht dt with hz = z 4 and gz = +iz e iz, we find that [ ] Iα = eαhzc π α h z c gz c + Oα = 4π α + Oα 3.7 5

16 as α +. For I m α, using ym e αy dy = Oα m+/, we find that I m α = Oα m+ as α Consider the limit α 0 +. After the change of the variables t = z/α, Iα = e 4α α e z i 4α dz. 3.9 α + iz The integrand is analytic in the complex plane minus the vertical line from iα to i. ote that the saddle point is i and it is on the branch cut. We show that the main contribution to the integral comes from the branch point z = iα. We deform the contour so that it consists of the following four line segments: L from i to i on the left half-plane, L from i to iα lying on the left of the branch cut, L 3 from iα to i lying on the right of the branch cut, and L 4 from i to i + lying on the right-half plane. On L 4, setting z = i + αx, e L4 z i 4α dz = e x 4 α dx = O α 3.30 α + iz 0 α + i αx as α 0. Similarly, the integral over L is also of the same order. On the other hand, setting z = iα + iy, e L L3 z i 4α α dz = α + iz 0 e α+y 4α y dy = e α 4α α 0 e y + y y 4α dy. 3.3 y The function y y decreases as y increases from y = 0 to y =. Hence the main contribution to the integral comes near the point y = 0. Using Watson s lemma, α 0 e y + y y 4α dy = Γ/ α + Oα. 3.3 y Combining together and using Γ/ = π, we obtain 3.5. For I m α, the analysis is same except that we use α iα + iy m e y + y y 4α dy = Oα m+/ y 0 Hence, we find that for m 0, I m α = O as α 0 +. Together with 3.8, this implies the uniform boundness of I m α. For the positiveness of Iα, we first write it as e α 4 t + i t arctan t e α 4 t Iα = + t /4 dt = 0 + t cos t arctan t dt /4 The function θt = t arctan t is monotone increasing. We use the inverse function, t = tθ, to change the variables and find that Iα = e α + 4 t t 3/4 θ t cos dθ, t = tθ

17 Since e α 4 tθ is positive and monotone decreasing in θ, we obtain Iα > 0 for every α > 0 if we show that i π + t 3/4 θ 3π + t 3/4 θ t cos dθ π t cos dθ, 3.36 and ii 0 k+π k + t 3/4 θ t cos dθ, k =,, 3,, 3.37 k π is decreasing in k. i can be verified numerically. On the other hand, ii follows immediately from the fact + t 3/4 /t is a decreasing function of t. This completes the proof. We now evaluate the integral in 3. using the steepest descent analysis. Lemma 3.8. Fix J > and let β = J + B /. Consider Gz defined in 3.. Then, for every 0 < ɛ < 8, γ+i γ i e Gz dz = i e Gγ IF γ + O +4ɛ 3.38 with high probability, where F z = βz logz λ i logγ λ z γ γ λ i= 3.39 and Iα is defined in.. Recall that = γ λ see Lemma 3.5. Proof. We choose the γ, which defines the contour, as the critical point of Gz. The path of steepest-descent is locally a vertical line near the critical point. It turns out that, instead of using the path of steepest-descent, it is enough to proceed the analysis using the straight line γ + ir globally. This choice was also made for the analysis in the paramagnetic regime in [6]. We first write, using the function F z, γ+i γ i γ+i e Gz dz = e Gγ e Gz F z+ F z Gγ dz γ i From the definitions of Gz and F z, e γ Gz F z λ = e z γ γ λ. 3.4 z λ Changing the variables z = γ + it / and using the notation = γ λ, γ+i γ i e Gz dz = ie Gγ e it e F γ+it / Gγ dt it 7

18 It is easy to check that the part of the integral with t ɛ is small. To show this, we first note that R F γ + it γ λ i + it / Gγ = R log log + c t γ λ i i= with high probability for some constant c > 0, since there is a constant c > 0 such that c γ λ i c for all i =,,, with high probability. Hence, it e e F γ+it / Gγ dt ɛ + it e log + c t dt ɛ 3.43 c 8 ɛ dt + c t /4 dt = Oe ɛ + O /8 ɛ e with high probability. Consider the part t ɛ. ote that F z satisfies F γ = Gγ, F γ = G γ = 0, and for each m, F m z = O uniformly for z in a small neighborhood of γ by Corollary 3.3. For m =, by Lemma 3.6, c F γ c 3.44 for some constants 0 < c < c, uniformly in. By Taylor expansion, for t ɛ, and hence, Therefore, F γ + it / Gγ = F γt if γt 3 e F γ+it / Gγ = e F γt 4 6 3/ + O +4ɛ 3.45 if γt 3 / + O +6ɛ ɛ = ɛ it e e F γ+it / Gγ dt + it e it e F γ t 4 dt if γ + it / = IF γ if γ 4 t 3 e it + it / I 3 F γ + O +6ɛ. e F γ 4 t dt + O +6ɛ By 3.44 and Lemma 3.5, c ɛ F γ c ɛ. Hence, Lemma 3.7 implies that 3.47 IF γ c ɛ 3.48 for some constant c > 0. Hence, using Lemma 3.5, Lemma 3.7, and the uniform boundedness of F γ, we find that 3.47 is equal to IF γ + O +4ɛ

19 if 0 < ɛ < 8. Thus, using 3.48 and Lemma 3.5 again, we conclude that γ+i γ i e Gz dz = i e Gγ IF γ + O /+4ɛ Proof of Theorem.4 Proof of Theorem.4. From Lemma 3. and Lemma 3.8, for every 0 < ɛ < 8, with high probability. Using Stirling s formula, Z = C i e Gγ IF γ + O +4ɛ 3.5 C = Γ/ πiβ / = thus we find that F = log Z satisfies F = Gγ logβ + β i πβe / + O, 3.5 β log π + log IF γ + O 3 +4ɛ 3.53 with high probability. Let us consider Gγ. Since γ and Ĵ = J +J are away from λ,, λ with high probability, logγ λ i = logĵ λ i log γ Ĵ γ λ i = logĵ λ i + γ Ĵ γ Ĵ + γ λ i γ λ i + O γ Ĵ for i =,,, where we also use that γ Ĵ = O +ɛ with high probability see Lemma 3.5. Then, using Lemma 3.6 and the fact that G γ = β i= γ λ i = 0, logγ λ i = i= logĵ λ i + βγ Ĵ γ Ĵ γ Ĵ + γ λ J + O 3 +3ɛ i= with high probability. Hence, from the formula of Gz in 3., Gγ = βĵ = βĵ logĵ λ i logγ λ + γ Ĵ γ Ĵ γ λ i= J + O 3 +3ɛ logĵ λ i log + s s J + O 3 +3ɛ i= 9

20 using the notations s = γ Ĵ and = γ λ in Lemma 3.5. Thus, F =βĵ logβ + s i= logĵ λ i + 4 log s 4J + log + log β π + log IF γ + O 3 +4ɛ To conclude Theorem.4, we use i the fact that = s χ, ii the asymptotic 3.0 of s in terms of χ, iii the fact that F γ = J + O +ɛ which follows from Lemma 3.6, and iv the fact that I α is uniformly bounded for α > 0 see Lemma Partial linear statistics This section is devoted to a proof of Theorem. on partial linear statistics. The proof is a simple modification of [6] for the linear statistics of all eigenvalues, which, in turn, follows the proof of [, ] for the case when the random matrix has zero mean. 4. Proof of Theorem. Recall Ĵ := J + J denotes the classical location of the largest eigenvalue of a Wigner matrix of non-zero mean. Fix -independent constants a < and < a + < Ĵ. Let Γ be the rectangular contour whose vertices are a ± iv 0 and a + ± iv 0 for some v 0 0, ]. The contour is oriented counter-clockwise. For a test function ϕx which is analytic in a neighborhood of [, ], we consider where ϕ := ϕλ i = i= i= R ϕxdσ scl x ϕz dz πi Γ z λ i πi ξ z := i= Decompose Γ into Γ u Γ d Γ l Γ r Γ 0, where R Γ λ i z R ϕz z x dzdσ sclx = ϕzξ πi dz, Γ 4. x z dσ sclx. 4. Γ u ={z = x + iv 0 : a x a + }, 4.3 Γ d ={z = x iv 0 : a x a + }, 4.4 Γ l ={z = a + iy : δ y v 0 }, 4.5 Γ r ={z = a + + iy : δ y v 0 }, 4.6 Γ 0 ={z = a + iy : y < δ } {z = a + + iy : y < δ } 4.7 0

21 for some δ > 0. In the proof of Theorem.6 in [6], the authors showed that ξ z := i= λ i z R x z dσ sclx = ξ z + λ z 4.8 converges weakly to a Gaussian process with mean bz = b z + and covariance Γz Ĵ z i, z j = Γ z i, z j where b z and Γ z i, z j are given in the proposition below. Since for each fixed z C +, λ z in probability by Lemma 3., it is natural to expect the following result for Ĵ z a partial sum. Proposition 4.. Let sz = x z dσ sclx = z + z 4 be the Stieltjes transform of the semicircle measure. Fix a constant c > 0 and a path K C + such that Iz > c for z K. Then the process {ξ z : z K} converges weakly to a Gaussian process with the mean b z = sz J sz + Jsz + w sz + s zsz + W 4 3sz Ĵ z and the covariance matrix Γ z i, z j = s z i s z j Remark 4.. ote that as z Ĵ, 4.9 w + W 4 3sz i sz j + sz i sz j. 4. sz J sz + Jsz = s z J + sz = s Ĵ + + Oz Ĵ. 4. z Ĵ s Ĵ Hence, b z is analytic near Ĵ and thus analytic for z C \ [, ]. In order to complete the proof of Theorem., we will prove the following lemma. Lemma 4.3. Define the events Ω := {λ Ĵ /3, λ + /3 } 4.3 which satisfies PΩ c < D for any fixed large D > 0. Then for some δ > 0, lim lim sup E ξ v z Ω dz = 0, 4.4 Γ # where Γ # can be Γ r, Γ l or Γ 0. From the explicit formulas 4.0 and 4., it is easy to check that lim E ξ z dz = v Γ #

22 Proposition 4., Lemma 4.3 and 4.5 imply that ϕ converges in distribution to a Gaussian random variable with the following mean and variance: ϕzb zdz, πi πi ϕz ϕz Γz, z dz dz. 4.6 Γ Γ Γ It is direct to check that these are equal to M ϕ and V ϕ see Section 4. in [6]. We thus obtain Theorem.. 4. Proof of Proposition 4. From Theorem 7. of [8], we need to show i the finite-dimensional convergence of ξ z to a Gaussian vector with desired mean and variance, and ii the tightness of ξ z. We will base our proof on the corresponding properties of ξ z obtained in [6]. Let us first recall the limit theorem for ξ z. Lemma 4.4 Proposition 4. in [6]. Let sz and K defined in the same way as in Proposition 4.. Then, the process {ξ z : z K} converges weakly to a Gaussian process {ξz : z K} with the mean bz = sz J sz + Jsz + w sz + s zsz + W 4 3sz and the covariance matrix Γz i, z j = s z i s z j w + W 4 3sz i sz j + sz i sz j. 4.8 Let z, z,, z p are p distinct points in K. The above lemma implies that the random vector ξ z i p i= converges weakly to a p-dimensional Gaussian distribution with the mean bz i p i= and the covariance matrix Γz i, z j. Since the distance between K and λ is bounded below, λ z i in probability for i =,, p. Hence, by Slutsky s theorem, ξ Ĵ z i z i p i= converges weakly to a p-dimensional Gaussian distribution vector with the mean b z i p i= and the covariance matrix Γ z i, z j, where b z = bz Ĵ z, 4.9 and Γ z i, z j = Γz i, z j. From Theorem.3 of [8], in order to show the tightness of a random process ζ z z K, it is sufficient to show that i ζ z is tight for a fixed z, and ii the following Hölder condition holds: for some -independent constant K > 0, E ζ z ζ z K z z, z, z K. 4.0 In [6], the authors considered the random process ζ z := ξ z E[ξ z], and proved that it satisfies conditions i and ii. ow, we consider ξ z := ζ + E[ξ z], where ζ z := ζ z λ z. Since E[ξ z] converges, it is enough to check that ζ z satisfies conditions

23 i and ii. ow for a fixed z, the tightness of ζ z and the boundedness of λ z imply that ζ z is tight. On the other hand, since ζ z satisfies the Hölder condition and Iz c for z K, E ζ z ζ z E ζ z ζ z + E λ z λ z K z z + z z c 4 = K + c 4 z z. Thus {ξ z, z K} is tight. This completes the proof of Proposition Proof of Lemma For z Γ 0, we notice that ξ Ω C and then Γ 0 E ξ Ω C δ. 4. Thus 4.4 holds for Γ 0 with δ >. For Γ r and Γ l, it is sufficient to show E ξ < K for some -independent constant K > 0. The authors in [6] showed that E ξ z < K. Hence, for z Γ r, ξ z Ω ξ z Ω + λ z Ω. 4.3 The lemma then follows from the fact that λ z Ω is bounded. 5 Joint Distribution of χ and ϕ As before, let A be a random symmetric matrix of size whose entries are up to the symmetry condition independent centered random variables satisfying Definition.. Let M = A+ J T where J >. Let λ λ be the eigenvalues of M. Let χ = λ Ĵ denoting the rescaled largest eigenvalue. Given an analytic function ϕx, recall the partial linear statistics ϕ = i= ϕλ i ϕxdσ sclx. We saw in the previous sections that χ and ϕ converge individually to Gaussian random variables. In this section, we consider the joint distribution and prove Theorem.3. In Subsection 5., we first prove Theorem.3 assuming that the disorder variables are Gaussian random variables. In Subsection 5., the general disorder variables are considered using an interpolation trick. 5. Asymptotic Independence for the GOE case Let the off-diagonal entries of A be Gaussian random variables of variance and the diagonal entries be Gaussian random variables of variance. In random matrix theory, the random symmetric Even though it is stated in Lemma 4. of [6] that the lemma holds for sufficiently small δ > 0, the proof of it is valid for any δ > 0, and we use δ > for our purpose. 3

24 matrix H = A is said to belong to the Gaussian orthogonal ensemble GOE. A special property of GOE, compared with general random symmetric matrices, is that the probability measure of GOE is invariant under orthogonal conjugations. The following result is basically in [0]. Lemma 5.. Let A ii, A ij, y i i<j be i.i.d. standard Gaussian random variables. Let H = A with A = A ij i,j and let Y = y,, y T. Define Gz = H zi for z C\[ δ, + δ], which is well defined with high probability for fixed δ > 0. Then, for z R\[ δ, + δ], where n = nz := n z := Y GzY TrGz nz 5. 0, dσ scl x is a Gaussian random variable. x z Proof of Lemma 5.. We follow the idea presented in [0]. By Theorem 5. of [0], it is enough to check the following three conditions for G: i There exists an -independent constant a such that G a with high probability, ii Tr G converges to a constant in probability, and iii i= G ii converges to a constant in probability. They follow from rigidity of eigenvalue Lemma 3., law of large numbers Corollary 3.3, and local law Theorem.9 of [7], respectively. We are now ready to prove the following property of GOE matrices. Proposition 5.. For H defined in Lemma 5., denote its eigenvalues by ρ ρ ρ. For fixed k, consider a random vector X, X,, Xk whose entries are real measurable functions of those eigenvalues, i.e., X i = Xi ρ, ρ,, ρ for i =,,, k. Suppose there is a random vector X i k i= such that Xi k i= Xi k i= as. Then for n and n defined as in 5., X, X,, Xk, n X, X, X k, n, where n is independent from X, X,, X k. Proof. For the convergence, it is enough to show i X, X, Xk, n is tight, and ii convergence of characteristic function. The tightness follows from the tightness of individual random vector variable, which is a consequence of individual convergence. For ii, consider the eigenvalue decomposition H = OP O T, where P = diagρ, ρ,, ρ and O is an orthogonal matrix. Since the H is orthogonal invariant, P and O are independent. Set X = O T Y. Then X = x,, x where x,, x are i.i.d standard Gaussian X is also independent with P. ow, n = Y GzY TrGz = x i i= ρ i z. Since E[etx ] = t, we find that for any t ir, the conditional expectation over X given P satisfies E X [ e tn P ] = EX [ e t i= x i ρ i z P ] = i= e log t t ρi z ρi z. ote that X, X,, Xk only depends on the eigenvalues, and hence it is independent of X. Thus, for any u, u,, u k, t ir, E [e ] k j= u jx j +tn = E [ e k j= u jx j 4 i= e log t t ρi z ρi z ]. 5.

25 Since log z z = z + Oz 3 as z 0, using Corollary 3.3, i= e log t t ρi z ρi z = e i= t ρ i z +O = e t x z dσ sclx+o [ = E e tnz] e O 5.3 with high probability. Denote this high probability event by Ω. Then, [ lim E e ] k j= u jx j +tn = lim E [e k j= u jx j +tn ] Ω PΩ + E = E [e k j= u jx j] [ E e tnz], [e k j= u jx j +tn Ω c ] PΩ c 5.4 since t, u, u,, u k ir and hence all exponents are pure imaginary. ote that the characteristic function of X,, X k, n is equal to the product of the characteristic functions of individual random vector variable. Thus nz is independent from X,, X k. This completes the proof. Corollary 5.3. Fix δ > 0, consider z C\R and z R\[ δ, + δ]. Recall sz defined in 4.9. Then TrGz sz, n z converges in distribution to independent Gaussian random variables. Proof. ote that TrGz sz is complex, we consider the random vector RTrGz sz, ITrGz sz. By Proposition 5., it is enough to show that RTrGz s, ITrGz sz converges to a Gaussian random vector. Consider the expression z = E + iη for ɛ, η R and η 0. Recalling the definition of linear statistics ϕ defined in.3, we have and RTrGz sz = ϕ r, ϕ r x = ITrGz sz = ϕ i, ϕ i x = x E x E + η, η x E + η. That is, they are both linear statistics. Then Corollary then follows from Theorem. of []. Remark 5.4. When we prove Theorem.3 for GOE, we use Proposition 5. and Corollary 5.3 with + -dependent z i. First, for a fixed z R\[ δ, +δ] for some δ > 0, let z = z := z. Using the exactly same argument in the proof of Lemma 5., one can show n z nz. Since the 5.3 still holds for z and nz, the asymptotic independence in Proposition 5. is still valid, i.e. X, X,, X, k n z X, X, X k, nz, where nz is independent from X, X,, X k. Second, for z C\R, consider z = z := + z. otice that = + x z x z z x z + O. 5

26 Then, by the discussion in Remark., TrG z s z = x z converges to a Gaussian random variable. ow, putting together, for z and z defined as above, TrG z s z, n z converge jointly to independent Gaussian random variables. We now prove Theorem.3 for the case where the disorder belongs to GOE. Proof of Theorem.3 when A belongs to GOE. Recall that λ i are the eigenvalues of M = A + J T with A from the GOE. Since the means and variances follow from [0] and Theorem., it is enough to prove the asymptotic independence of χ and ϕ. otice that that W 3 = 0 for Gaussian A ij. ow, for any analytic test function ϕ, the partial linear statistics can be expressed as see 4. an integral of ξ z = i= λ i z R x z dσ sclx, z C\R. 5.5 Then according to Lemma 4.3 and what follows, it is enough to prove that χ and ξ z are asymptotically independent for fixed z C \ R. Let ξ z = ξ z + λ z = TrM zi sz. Since λ z Ĵ z independent. in probability, it is enough to prove that χ and ξ z are asymptotically Since the GOE is orthogonal invariant, for every deterministic matrix U, the eigenvalues of A + U have the same distribution as A + OUO T for any orthogonal matrix O. Thus, we may consider the following equivalent model: Following the proof of Theorem. in [0], we write [ A + J M = Y Since detm zi = det ˆM zi A + J z Y ĜzY the largest eigenvalue of M satisfies M = A + diagj, 0,, Y ˆM ]. 5.7 with Ĝz := ˆM zi = Â zi, 5.8 λ = J + A Y Ĝλ Y 5.9 if λ is not an eigenvalue of twice, we write ˆM, which holds with high probability. Using the resolvent formula Ĝλ = ĜĴ + Ĝλ ĜĴ = ĜĴ + λ ĴĜλĜĴ = ĜĴ + λ ĴĜĴ + λ Ĵ Ĝλ ĜĴ. 6

27 Hence, λ Ĵ = A J Y Ĝλ Y = A J Y ĜĴY + λ ĴY ĜĴ Y + λ Ĵ Y Ĝλ ĜĴ Y with high probability. Moving all terms with factor λ Ĵ to the left and taking it out as a common factor, we arrive at χ = λ Ĵ = A J + Y ĜĴY + Y ĜĴ Y + λ ĴY Ĝλ ĜĴ Y 5.0 with high probability. ote that ˆM and Y satisfy the setting of Corollary 5.3 up to the scaling factor. Set Ỹ = Y, Gz = ˆM zi 5. Then, Ỹ and G satisfy the setting of Corollary 5.3, and ow, by Corollary 3.3, Y ĜĴY = Ỹ G J Ỹ, J := 5. Ĵ. Tr G J = sĵ + O +ɛ = J + O +ɛ 5.3 with high probability. By Lemma 5., Corollary 3.3 and Lemma 3., Y ĜĴ Y J, λ ĴY Ĝλ ĜĴ Y in probability. Using 5.4, 5.3 and denoting the denominator in 5.0 by D, we write χ = D A ñ J + O +ɛ, 5.5 where n J = Ỹ G J Ỹ Tr G J see 5. and D J in probability. J ote that A and n J are independent, the distribution of χ is governed by their convolution. We now turn to the linear statistic ξ z. Using Schur complement of M with block structure in 5.7, for any z C\R, TrM zi =J + A z Y ĜzY + Y Ĝz Y + TrĜz 5.6 Using Lemma 5. and Lemma 3.3, D = D := + Y Ĝz Y J + A z Y ĜzY + s z J z sz 7

28 in probability. Then, by setting z := z = z, we write ξ z =TrM zi sz = D + TrĜz sz + O +ɛ =D sz + zs z + Tr G z s z + O +ɛ. 5.7 That is, the fluctuation of ξ z is govern by Tr G z s z. ow using Corollary 5.3 and Remark 5.4, one can conclude that Tr G z s z, n J converge to independent Gaussian random variables. Furthermore, A is independent of both Y and ˆM. Thus by 5.5 and 5.7, ξ z, χ converge to independent random variables. Theorem.3 then follows. 5. Proof of Theorem.3 for general case We prove Theorem.3 for general disorders, where the disorder matrix A is a Wigner matrix and satisfies Definition.. Unlike the GOE, Wigner matrices are not orthogonal invariant, hence we cannot apply 5.6 where we replaced the rank- perturbation in M by a diagonal matrix. To overcome the difficulty, we use an interpolation method. It has been successfully applied in many works in random matrix theory, where a given matrix and a reference matrix such as GOE are interpolated. We refer to [] for its application in the analysis of linear eigenvalue statistics. Let V = A be a normalized Wigner matrix and V G be a normalized GOE matrix independent from V. Define Ht = V cos t + V G sin t 5.8 so that H0 = V and H π = V G. ote that E[H ij ] = and for i j. Let e = T =,,..., T R 5.9 Mt = Ht + Jee T, 5.0 whose eigenvalues are denoted by λ λ λ. Define the resolvents Gz = M zi, Ĝz = H zi. 5. Here, we omit the dependence on t for the ease of notation. We note that G and Ĝ are symmetric not Hermitian. For any small fixed δ > 0, Ĝz is well-defined for z C \ [ δ, + δ] with high probability. For χ = λ Ĵ, we notice that Ĝ ee λ := e, Ĝλ e = J 5. with high probability. The claim holds since 0 = detm λ I = deth λ I deti + JĜλ ee T = deth λ I deti + Je T Ĝλ e = deth λ I JĜeeλ 8

29 and λ is not an eigenvalue of H with high probability See Lemma 6. of [0]. Furthermore, by Taylor expansion, J = Ĝeeλ = ĜeeĴ + Ĝ eeĵλ Ĵ + O +ɛ 5.4 with high probability, since λ Ĵ = O +ɛ and Ĝ z = O with high probability. From the isotropic local law, Theorem. of [0], we find that Ĝ ee Ĵ = sĵ + O +ɛ, Ĝ eeĵ = s Ĵ + O +ɛ 5.5 with high probability. Thus, using Lemma 5., χ = λ Ĵ = J + ĜeeĴ s Ĵ + O +ɛ 5.6 with high probability. That is, the behavior of χ is governed by the fluctuation of ĜeeĴ. To prove the Theorem.3, as in the Gaussian disorder case, it is enough to show the convergence of the joint distribution of χ and the full linear statistics ξ z = TrGz sz for fixed z C \ R. Under the light of 5.6, we set out to calculate the following characteristic function involving ξ z and ĜeeĴ. Explicitly, for t, t, t 3 ir and z = E + iη with E R and η > 0, we define E [ e P t] := E [e t Rξ +t Iξ +t 3 n ], P t := t Rξ z + t Iξ z + t 3 n, 5.7 where n = Ĝ ee Ĵ J ote that n is real, the exponent P t is pure imaginary and thus e P t. For our purpose, it is desired to estimate E[e P 0 ]. At t = π, the disorder H π reduces to the GOE case. From Subsection 5., χ and ξ are asymptotically independent in the GOE case, then [ lim E e P π ] [ ] = E e t Rξ+t Iξ E [ e t 3n ] 5.9 for some Gaussian random variables ξ, n with known mean and variance. Thus, it only remains to estimate the t-derivative of E[e P t ]. Here, we recall the following identity for the derivative of the resolvent G. For i, j, a, b =,,,, M ij G ab = β jk G aj G kb + G ak G jb 5.30 with β jk = { j k, / j = k

30 We note that the above identity also holds if one replace G by Ĝ. Thus for any fixed event Ω, d [ ] dt E e P t Ω = E dm ij e P t dt M ij Ω i j = [ Vij E sin t Vij G cos t t R G ij + t I ] G ij + t 3 Ĝ pi Ĝ jq e P t Ω. i,j p,q 5.3 The reason for the introduction of Ω will be revealed in a minute. The right hand side of 5.3 motivates us to apply the generalized Stein s lemma. More precisely, we will use Proposition 3. of [] with a small modification as follows: Proposition 5.5. Given an event Ω, let X be a random variable such that E[ X p+ Ω] < for a certain non-negative integer p. Denote the conditional cumulants of X by κ l := κ l Ω, l =,..., p +. Then for any function Φ : R C of the class C p+ with bounded derivatives Φ l, l =,..., p +, we have E[XΦX Ω] = p l=0 where the remainder term ɛ p admits the bound ɛ p C p E X p+ + max for some constant C p that depends only on p. κ l+ l! j p+ E[Φ l X Ω] + ɛ p, 5.33 Φ p+ vx dv 0 p+ j Ω 5.34 Proof. We basically follow the proof of Proposition 3. of []. Let π p be the degree p Taylor polynomial of Φ and let r p = Φ π p. Then, as in the proof of Proposition 3. of [], p κ j+ E[Xπ p X Ω] = E[π p j X Ω] j! j=0 Thus p E[XΦX Ω] κ l+ E[Φ l p X Ω] l! E[Xr κ l+ [ px Ω] + E r l p X Ω] l! l=0 l=0 Since r p X = Xp+ Φ p+ vx v p dv, 5.37 p! l=0 0 by the estimate κ j j j E[ X j Ω] and Hölder s inequality, p κ l+ [ ] E r l p κ ] l+ p X Ω [ X l! l!p l! E p+ l Φ p+ vx dv Ω l=0 l=0 0 p l + l+ p+ E X p+ + Φ p+ p+ l vx dv l!p l! Ω

31 As E[Xr p Ω] can also be bounded by the right hand side of 5.38, the proof is complete. In order to apply Proposition 5.5 to 5.3, we need prior bounds of P t and its derivatives to bound ɛ p in As we will see later, it is enough to bound G ij, G ij, Ĝ ij and p Ĝip. In the following, we are going to introduce a high probability event Ω, on which we have the desired bounds. With the trivial bound G η recall that z = E + iη, we have that G ij η and G ij G Ĵ/}. It is easy to. For η Ĝij, we introduce the high probability event Ω = {λ + and thus check that Ĝ Ω Ĵ Ĝij Ω Ĵ, 5.39 For p Ĝip, we recall the following concentration theorem for the quadratic function of Ĝ: Proposition 5.6 Theorem.3 and Remark.4 of [0]. Fix Σ 3. Set ϕ = log log log. Then there exist constants C and C such that for any E [Σ, ϕ C 3 ] [ + ϕ C 3, Σ], and any η 0, Σ], and any deterministic v, w C, v, Ĝzw sz v, w ϕc with high probability, uniformly on z = E + iη. Isz v w 5.40 η Let e i := 0,,,, 0. oting that p= Ĝpi = e, Ĝe i, we can derive a prior bound for p= Ĝpi, which is summarized in the following Corollary. Corollary 5.7. For any fixed E R\[, ], the tail bound p ĜE pi ɛ 5.4 holds simultaneously for i =,, with high probability. We also have that with high probability. v, ĜEw se v, w v w +ɛ 5.4 Proof. We first prove 5.4. Consider z = E + i /. Using Proposition 5.6, we find there exists some C > 0 such that v, ĜEw se v, w v, Ĝz ĜEw + v, Ĝzw sz v, w + sz se v, w C / v w + Cϕ C v w + C / Here we also use the fact that Ω holds with high probability. Since ϕ ɛ, 5.4 then follows. The tail bound 5.4 can be obtained from 5.4 by setting v = e and w = e i. 3