arxiv: v1 [math.pr] 9 Sep 2014

Transcription

1 MOD-GAUSSIAN CONVEGENCE AND ITS APPLICATIONS FO MODELS OF STATISTICAL MECHANICS PIEE-LOÏC MÉLIOT AND ASHKAN NIKEGHBALI arxiv: v1 [math.p] 9 Sep 014 ABSTACT. In this paper we complete our understanding of the role played by the limiting or residue function in the context of mod-gaussian convergence. The question about the probabilistic interpretation of such functions was initially raised by Marc Yor. After recalling our recent result which interprets the limiting function as a measure of "breaking of symmetry" in the Gaussian approximation in the framework of general central limit theorems type results, we introduce the framework of L 1 -mod-gaussian convergence in which the residue function is obtained as up to a normalizing factor the probability density of some sequences of random variables converging in law after a change of probability measure. In particular we recover some celebrated results due to Ellis and Newman on the convergence in law of dependent random variables arisisng in statistical mechanics. We complete our results by giving an alternative approach to the Stein method to obtain the rate of convergence in the Ellis-Newman convergence theorem and by proving a new local limit theorem. More generally we illustrate our results with simple models from statistical mechanics. In memoriam, Marc Yor. 1. INTODUCTION Let X n n N be a sequence of real-valued random variables. In the series of papers [JKN11, DKN11, KN10, KN1, FMN13], we introduced the notion of mod-gaussian convergence and more generally of mod-convergence with respect to an infinitely divisible law φ: Definition 1. The sequence X n n N is said to converge in the mod-gaussian sense with parameters t n + and limiting or residue function θ if, locally uniformly in, E[e itx n ] e t n t = θt 1 + o1, where θ is a continuous function on with θ0 = 1. A trivial situation of mod-gaussian convergence is when X n = G n + Y n is the sum of a Gaussian variable of variance t n and of an independent random variable Y n that converges in law to a variable Y with characteristic function θ. More generally X n can be thought of as a Gaussian variable of variance t n, plus a noise which is encoded by the multiplicative residue θ in the characteristic function. In this setting, θ is not necessarily the characteristic function of a random variable the residual noise. For instance, consider X n = 1 n 1/3 n Y i, i=1 Date: September 10,

2 PIEE-LOÏC MÉLIOT AND ASHKAN NIKEGHBALI where the Y i are centred, independent and identically distributed random variables with convergent moment generating function. Then a Taylor expansion of E[e ity ] shows that X n n N converges in the mod-gaussian sense with parameters n 1/3 VarY and limiting function E[Y θt = 3 ] it 3 exp, 6 which is not the characteristic function of a random variable, since it does not go to zero as t goes to infinity. In 008, during the workshop "andom matrices, L-functions and primes" held in Zürich, Marc Yor asked the second author A. N. about the role of the limiting function θ. In [KNN13] it is proved that the set of possible limiting functions is the set of continuous functions θ from to C such that θ0 = 1 and θ t = θt for t. But this characterization does not say anything on the probabilistic information encoded in θ. We now wish to develop more on probabilistic interpretations of the limiting function and the implications of mod-gaussian convergence in terms of classical limit theorems of probability theory. We first note that by looking at E[e itx n/ t n ], one immediately sees that mod-gaussian convergence implies a central limit theorem for the sequence X n tn : X n tn n N 0, 1, 1 where the convergence above holds in law see [JKN11, -3] for more details on this. On the other hand, with somewhat stronger hypotheses on the remainder o1 that appears in Definition 1, a local limit theorem also holds, see [KN1, Theorem 4] and [DKN11, Theorem 5]: P[X n B] = P[G n B] 1 + o1 = mb πtn 1 + o1 for relatively compact sets B with m B = 0, m denoting the Lebesgue measure. In [FMN13], it is then explained that by looking at Laplace transforms instead of characteristic functions, and by assuming the convergence holds on a whole band of the complex plane, one can obtain in the setting of mod-gaussian convergence precise estimates of moderate or large deviations. In fact these results provide a new probabilistic interpretation of the limiting function as a measure of the "breaking of symmetry" in the Gaussian approximation of the tails of X n see 1.1 for more details. The goal of this paper is threefold: to propose a new interpretation of the limiting function in the framework of mod-gaussian convergence with Laplace transforms; these results allow us in particular to recover some well known exotic limit theorems from statistical mechanics due to Ellis and Newman [EN78] and similar one for other models or in higher dimensions.

3 MOD-GAUSSIAN CONVEGENCE AND APPLICATIONS FO MODELS OF STATISTICAL MECHANICS 3 to show that once one is able to prove mod-gaussian convergence, then one can expect to obtain finer results than merely convergence in law, such as speed of convergence and local limit theorems. esults on the rate of convergence in the Curie-Weiss model were recently obtained using Stein s method see e.g. [EL10] while the local limit theorem, to the best of our knowledge, is new. to explore the applications of the results obtained in [FMN13] on the "breaking of symmetry" in the central limit theorem to some classical models of statistical mechanics. In particular our approach determines the scale up to which the Gaussian approximation for the tails is valid and its breaking at this critical scale. Our results are best illustrated with some classical one-dimensional models from statistical mechanics, such as the Curie-Weiss model or the Ising model. To illustrate the flexibility of our approach, we shall also prove similar results for weighted symmetric random walks in dimensions and 3. The statistics of interest to us will be the total magnetization, which can be written as a sum of dependent random variables. These examples add to the already large class of examples of sums of dependent random variables we have already been able to deal with in the context of mod-φ convergence. In the remaining of the introduction we recall the results obtained in [FMN13] which led us to the "breaking of symmetry" interpretation, as well as an underlying method of cumulants that enabled us to establish the mod-gaussian convergence for a large family of sums of dependent random variables. The important aspect of the cumulant method is that it provides a tool to prove mod-gaussian convergence in situations where one cannot explicitly compute the characteristic function. We eventually give an outline of the paper Complex convergence and interpretation of the residue. We consider again a sequence of real-valued random variables X n n N, but this time we assume that their Laplace transforms E[e zx n] are convergent in an open disk of radius c > 0. In this case, they are automatically well-defined and holomorphic in a band of the complex plane B c = {z C, ez < c} see [LS5, Theorem 6], and [Ess45] for a general survey of the properties of Laplace and Fourier transforms of probability measures. Definition. The sequence X n n N is said to converge in the complex mod-gaussian sense with parameters t n and limiting function ψ if, locally uniformly on B c, E[e zx n ] e t n z = ψz 1 + o1, where ψ is a continuous function on B c with ψ0 = 1. Then, one has in particular convergence in the sense of Definition 1, with θt = ψit. In this setting which is more restrictive than before, the residue ψ has a natural interpretation as a measure of "breaking of symmetry" when one tries to push the estimates

4 4 PIEE-LOÏC MÉLIOT AND ASHKAN NIKEGHBALI of the central limit theorem from the scale t n to the scale t n. The previously mentioned central limit theorem 1 tells us that: P [ X n a t n ] = 1 π a e x dx 1 + o1 for any a. In the setting of complex mod-gaussian convergence, this estimate remains true with a = o t n, so that if ε = o1, then 1 P [X n ε t n ] = π ε e x dx 1 + o1, t n = tnε e o1 if 1 ε, πtn ε tn where the notation a n b n stands for b n = oa n. Then, at scale t n, the limiting residue ψ comes into play, with the following estimate that holds without additional hypotheses than those in Definition : x 0, c, P[X n xt n ] = tnx e ψx 1 + o1, πtn x the remainder o1 being uniform when x stays in a compact set of + 0, c. This estimate of positive large deviations has the following counterpart on the negative side: x 0, c, P[X n xt n ] = tnx e ψ x 1 + o1. πtn x So for instance, if Y n n N is a sequence of i.i.d. random variables with convergent moment generating function, mean 0, variance 1 and third moment E[Y] > 0, then X n = 1 n 1/3 i=1 n Y i converges in the complex mod-gaussian sense with parameters n 1/3 and limiting function ψz = expe[y 3 ] z 3 /6, and therefore for x > 0, P [ n i=1 Y i xn /3 ] [ = P N 0, 1 xn 1/6] exp E[Y 3 ] x o1. 6 Thus, at scale n /3, the fluctuations of the sum of i.i.d. random variables are no more Gaussian, and the residue ψx measures this "breaking of symmetry": in the previous example, it makes moderate deviations on the positive side more likely than moderate deviations on the negative side, since ψx > 1 > ψ x for x > 0. emark. The problem of finding the normality zone, i.e. the scale up to which the central limit theorem is valid, is a known problem in the case of i.i.d. random variables see e.g. [IL71]. The description of the "symmetry breaking" is new and moreover the mod-gaussian framework covers many examples with dependent random variables see also [FMN13] for more examples. Thus, the observation of large deviations of the random variables X n provides a first probabilistic interpretation of the residue ψ in the deconvolution of a sequence of characteristic functions of random variables by a sequence of large Gaussian variables.

5 MOD-GAUSSIAN CONVEGENCE AND APPLICATIONS FO MODELS OF STATISTICAL MECHANICS 5 In Section 3, we shall provide another interpretation of ψ, which is inspired by some classical results from statistical mechanics cf. [EN78, EN80]. 1.. The method of joint cumulants. The appearance of an exponential of a monomial Kx r 3 as the limiting residue in mod-gaussian convergence is a phenomenon that occurs not only for sums of i.i.d. random variables, but more generally for sums of possibly non identically distributed and/or dependent random variables. For instance, 1 the number of zeroes of a random Gaussian analytic function k=0 N C k z k in the disk of radius 1 1 n ; the number of triangles in a random Erdös-ényi graph Gn, p; are both mod-gaussian convergent after proper rescaling, and with limiting function of the form explz 3, with the constant L depending on the model see again [FMN13]. The reason behind these universal asymptotics lies in the following method of cumulants. If X is a random variable with convergent Laplace transform E[e zx ] on a disk, we recall that its cumulant generating function is log E[e zx κ ] = r X z r, 3 r! r 1 which is also well-defined and holomorphic on a disk around the origin. Its coefficients κ r X are the cumulants of the variable X, and they are homogenenous polynomials in the moments of X; for instance, κ 1 X = E[X], κ X = E[X ] E[X], and κ 3 X = E[X 3 ] 3 E[X ] E[X] + E[X] 3. Consider now a sequence of random variables W n n N with κ 1 W n = 0, and for r, κ r W n = K r α n 1 + o1, 4 with α n +. This assumption is inspired by the case of a sum W n = i=1 n Y i of centred i.i.d. random variables for which κ r W n = n κ r Y. If it is satisfied, then one can formally write [ log E e z ] Wn αn 1/3 = α n /3 κ W n z α n 1/3 K z α n 1/3 K z + K 3 z K 3 z α n 1 κ3 W n z r 4 K r z r α n 1 r/3 r! + r 4 κ r W n α n 1/3 z r r! whence the mod-gaussian convergence of X n = α n 1/3 W n, with parameters t n = K α n 1/3 and limiting function expk 3 z 3 /6. The approximation is valid if the o1

6 6 PIEE-LOÏC MÉLIOT AND ASHKAN NIKEGHBALI in the asymptotics of κ W n is small enough namely oα n 1/3, and if the series r 4 can be controlled, which is the case if r, κ r W n Cr r α n 5 for some constant C. The method of cumulants in the setting of mod-gaussian convergence amounts to prove 4 for the first cumulants of the sequence X n n N, and 5 for all the other cumulants. From such estimates one then obtains mod-gaussian convergence for an appropriate renormalisation of W n r 3, with limiting function expk r z r /r!, where r is the smallest integer greater or equal than 3 such that K r = 0. This method of cumulants works well with sequences W n n N that write as sums of weakly dependent random variables. Indeed, cumulants admit the following generalization to families of random variables, see [LS59]. Denote Q r the set of partitions of [[1, r]] = {1,, 3,..., r}, and µ the Möbius function of this poset see [ot64] for basic facts about Möbius functions of posets: µπ = 1 lπ 1 lπ 1! where lπ = s if Π = π 1 π π s has s parts. The joint cumulant of a family of r random variables with well defined moments of all order is [ ] lπ κx 1,..., X r = µπ E X j. Π Q r i=1 j π i It is multilinear and generalizes Equation 3, since κx 1,..., X r = r z 1 z z r κ X,..., X = κ }{{} r X. r times z1 = =z r =0 log E[e z 1 X 1 + +z r X r ] Suppose now that W n = W = i=1 n Y i is a sum of dependent random variables. By multilinearity, κ r W = κy i1,..., Y ir, 6 i 1,...,i r so in order to obtain the bound 5, it suffices to bound each "elementary" joint cumulant κy i1,..., Y ir. To this purpose, it is convenient to introduce the dependency graph of the family of random variables Y 1,..., Y n, which is the smallest subgraph G of the complete graph on n vertices such that the following property holds: if Y i i I and Y j j J are disjoint subsets of random variables with no edge of G between a variable Y i and a variable Y j, then Y i i I and Y j j J are independent. Then, in many situations, one can write a bound on the elementary cumulant κy i1,..., Y ir that only depends on the induced subgraph G[i 1,..., i r ] obtained from the dependency graph by keeping only the vertices i 1,..., i r and the edges between them. In particular: 1 κy i1,..., Y ir = 0 if the induced graph G[i 1,..., i r ] is not connected. if Y i 1 for all i, then κy i1,..., Y ir r 1 STG[i 1,..., i r ], where STH is the number of spanning trees on a connected graph H.

7 MOD-GAUSSIAN CONVEGENCE AND APPLICATIONS FO MODELS OF STATISTICAL MECHANICS 7 By gathering the contributions to the sum of Formula 6 according to the nature and position of the induced subgraph G[i 1,..., i r ] in G, one is able to prove efficient bounds on cumulants of sums of dependent variables, and to apply the method of cumulants to get their mod-gaussian convergence. We refer to [FMN13] for precise statements, in particular in the case where each vertex in G has less than D neighbors, with D independent of the vertex and of n. In Section 5, we shall apply this method to a case where G is the complete graph on n vertices, but where one can still find correct bounds and in fact exact formulas for the joint cumulants κy i1,..., Y ir : the one-dimensional Ising model Basic models. As mentioned above, the goal of the paper is to study the phenomenon of mod-gaussian convergence for probabilistic models stemming from statistical mechanics; this extends the already long list of models for which we were able to establish this asymptotic behavior of the Fourier or Laplace transforms [JKN11, KN1, FMN13]. More precisely, we shall focus on one-dimensional spin configurations, which already yield an interesting illustration of the theory and technics of mod- Gaussian convergence. Given two parameters α and β +, we recall that the Curie-Weiss model and the one-dimensional Ising model are the probability laws on spin configurations σ : [[1, n]] {±1} given by CW α,β σ = I α,β σ = 1 Z n CW, α, β exp α 1 Z n I, α, β exp α n i=1 n i=1 σi + β σi + β n n i=1 σi n 1 σi σi + 1 i=1 ; 7. 8 The coefficient α measures the strength and direction of the exterior magnetic field, whereas β measures the strength of the interaction between spins, which tend to align in the same direction. This interaction is local for the Ising model, and global for the Curie-Weiss model. Set M n = i=1 n σi: this is the total magnetization of the system, and a random variable under the probabilities CW α,β and I α,β. In Section, we quickly establish the mod-gaussian convergence of the magnetization for the Ising model, using the explicit form of the Laplace transform of the magnetization, which is given by the transfer matrix method. Alternatively, when α = 0, in the appendix, we apply the cumulant method and give an explicit formula for each elementary cumulant of spins see Section 5. This allows us to prove the analogue for joint cumulants of the well-known fact that covariances between spins decrease exponentially with distance in the 1D-Ising model. This second method is much less direct than the transfer matrix method, but we consider the Ising model to be a very good illustration of the method of joint cumulants. Moreover it illustrates the fact that one does not necessarily need to be able to compute precisely the moment generating function of the random variables. In Section 3, we focus on the Curie-Weiss model, and we interpret the magnetization as a change of measure on a sum of i.i.d. random variables. Since these sums converge in the mod-gaussian sense, it leads us to study the effect of a change of measure on a mod-gaussian convergent sequence. We prove that in the setting of

8 8 PIEE-LOÏC MÉLIOT AND ASHKAN NIKEGHBALI L 1 -mod-gaussian convergence, such changes of measures either conserve the mod- Gaussian convergence with different parameters, or lead to a convergence in law, with a limiting distribution that involves the residue ψ. We thus recover the results of [EN78, EN80], and extend them to the setting of L 1 -mod-gaussian convergence. In Section 4, using Fourier analytic arguments, we quickly recover the optimal rate of convergence of the Ellis-Newman limit theorem for the Curie-Weiss model which was recently obtained in [EL10] using Stein s method, and then we establish a local limit theorem, thus completing the existing limit theorems for the Curie-Weiss model at critical temperature CW 0,1.. MOD-GAUSSIAN CONVEGENCE FO THE ISING MODEL: THE TANSFE MATIX METHOD In this section, σi i [[1,n]] is a random configuration of spins under the Ising measure 8, and M n = i=1 n σi is its magnetization. The mod-gaussian convergence of M n after appropriate rescaling can be obtained by two different methods: the transfer matrix method, which yields an explicit formula for E[e zm n]; and the cumulant method, which gives an explicit combinatorial formula for the coefficients of the power series log E[e zm n]. We use here the transfer matrix method, and refer to the appendix Section 5 for the cumulant method. FIGUE 1. Two configurations of spins under the Ising measures of parameters α = 0, β = 0.3 and α = 0, β = 1. The Laplace transform E[e zm n] of the magnetization of the one-dimensional Ising model is well-known to be computable by the following transfert matrix method, see [Bax8, Chapter ]. Introduce the matrix e α+β e T = α β, e α β e α+β and the two vectors V = e α, e α and W = 1 1. If the rows and columns of T correspond to the two signs +1 and 1, then any configuration of spins σ = σi i [[1,n]] has under the Ising measure I α,β a probability proportional to V σ1 T σ1,σ T σ,σ3 T σn 1,σn. Therefore, the partition function Z n I, α, β is given by V σ1 T σ1,σ T σ,σ3 T σn 1,σn = VT n 1 W σ1,...,σn = a + λ + n 1 + a λ n 1,

9 MOD-GAUSSIAN CONVEGENCE AND APPLICATIONS FO MODELS OF STATISTICAL MECHANICS 9 where a + = cosh α + λ + = e β cosh α + eβ sinh α + e β ; a = cosh α eβ sinh e β sinh α + e β e β sinh α + e β ; λ = e β cosh α α + e β e β sinh α + e β e β sinh α + e β. Indeed, λ ± are the two eigenvalues of T, and a + and a are obtained by identification of coefficients in the two formulas Z 1 I, α, β = e α + e α Z I, α, β = e α+β + e α+β + e β. Then, the Laplace transform of M n is given by In particular, E α,β [M n ] = E[ezM n] z E α,β [e zm n ] = Z ni, α + z, β. Z n I, α, β = z=0 α log Z ni, α, β = n whence a formula for the asymptotic mean magnetization by spin: m = e β sinh α e β sinh α + e β. e β sinh α + O1. e β sinh α + e β A more precise Taylor expansion of Z n I, α + z, β leads to the following: Theorem 3. Under the Ising measure I α,β, M n nm n 1/3 sense with parameters converges in the complex mod-gaussian and limiting function ψz = exp t n = n 1/3 e β cosh α e β sinh α + e β 3/ eβ sinh 3 α + 3e β e 3β sinh α 6e β sinh z 3 α + e β 5/ Proof. In the following, we are dealing with square roots and logarithms of complex numbers, but each time in a neighborhood of +, so there is no ambiguity in the choice of the branches of these functions. That said, it is easier to work with log-laplace.

10 10 PIEE-LOÏC MÉLIOT AND ASHKAN NIKEGHBALI transforms: [ ] log E e z M n nm n 1/3 = log Z n I, α + z n 1/3, β log Z n I, α, β zn /3 m log Z n I, α, β = log a + α, β + n 1 log λ + α, β + o1 log Z n I, α + z n 1/3, β = log a + α + z n 1/3, β + n 1 log λ + α + z n 1/3, β + o1 = log a + α, β + n 1 log λ + α, β + zn /3 α log λ +α, β + z n 1/3 α log λ +α, β + z3 3 6 α 3 log λ +α, β + o1. Thus, it suffices to compute the first derivatives of log λ + α, β with respect to α: log λ + α, β = log e β cosh α + e β sinh α + e β α log λ +α, β = α log λ +α, β = e β sinh α e β sinh α + e β = m e β cosh α e β sinh = α + e β σ 3/ 3 α 3 log λ +α, β = eβ sinh 3 α + 3e β e 3β sinh α e β sinh α + e β 5/ = K 3. We therefore get [ ] log E e z M n nm n 1/3 = n 1/3 σ z + K 3 z 3 + o1. 6 By using Formula, this result leads to new estimates of moderate deviations for the probability P α,β [M n nm + n 1/3 x]. In the special case when α = 0, the limiting function ψz of Theorem 3 is equal to 1, and one has to push the expansion of log Z n I, 0, β to order 4 to get a meaningful mod-gaussian convergence the same phenomenon will occur in the case of the Curie-Weiss model: M Theorem 4. Under the Ising measure I 0,β, n converges in the complex mod-gaussian sense n 1/4 with parameters t n = n 1/ e β and limiting function ψz = exp 3e6β e β z 4. 4

11 MOD-GAUSSIAN CONVEGENCE AND APPLICATIONS FO MODELS OF STATISTICAL MECHANICS11 3. MOD-GAUSSIAN CONVEGENCE IN L 1 AND THE CUIE-WEISS MODEL In this Section, X n n N is a sequence of random variables with entire moment generating series E[e zx n], and we assume the following: A One has mod-gaussian convergence of the Laplace transforms, i.e., there is a sequence t n + and a function ψ continuous on such that ψ n t = E[e tx n ] e t n t converges locally uniformly on to ψt. B Each function ψ n, and their limit ψ are in L 1. We denote P n the law of X n, e tx n t n t Q n [dx] = E e tn x [e X n tn ] P n [dx], 9 and Y n a random variable under the new law Q n. Note that hypothesis B implies that Z n = E[e X n /t n ] is finite for all n N. Indeed, [ ] ] ] ψ n t dt = E dt = E [e X n tn e X n tnt π tn dt = E [e X n tn. Therefore the new probability measures Q n are well defined. The goal of this section is to study the asymptotics of the new sequence Y n n N. As we shall see in 3.3, the Curie-Weiss model defined by Equation 7 is one of the main examples in this framework. However, it is more convenient to look at the general problem, and we shall introduce later other models concerned by our general results Ellis-Newman lemma and deconvolution of a large Gaussian noise. Suppose for a moment that hypothesis A is replaced by the stronger hypotheses of Definition, with c = + and therefore B c = C. Fix then 0 < a < b, and consider the partial integral E[e X n /t n 1 tn a X n t n b]. By integration by parts of iemann-stieltjes integrals, one has: tn b t n a [ ] tn e tn x P n [dx] = e tn x b P n [X n x] = = + t n a ] b [ e t n x P n [X n t n x] [ ψx ] b + πtn x a = a tn π tn b ψx dx π a tn b + t n a b b a a t n x t n e x tn P n [X n x] dx t n x e t n x P n [X n t n x] dx ψx dx 1 + o a,b o a,b 1 because of the estimates of precise deviations. In this computation, o a,b 1 is uniform for a, b in compact sets of 0, +. In fact this estimate remains true for a, b in a

12 1 PIEE-LOÏC MÉLIOT AND ASHKAN NIKEGHBALI compact set of ; hence, a and b can be possibly negative. If the estimate is also true with a = and b = +, then Q n [t n a Y n t n b] = E[eX n /t n 1 tn a X n t n b] E[e X n /t n] t n b π a ψx dx = t n 1 + o1 + π ψx dx = b a + ψx dx 1 + o1, ψx dx so Y n t n n N converges in law to the density ψx/ ψx dx. We now wish to identify the most general conditions under which this convergence in law happens. To this purpose, it is useful to produce random variables with density ψ n x/ ψ nx dx. They are given by the following Proposition, which appeared in [EN78] as Lemma 3.3: Proposition 5. Let G n be a centred Gaussian variable with variance t 1 n, and independent from Y n. The law of W n = G n + Y n t n has density ψ n x/ ψ nx dx. Proof. Denote Z n = E[e X n /t n ], and f X x dx respectively, P X the density respectively, the law of a random variable X. One has w P[W n w] = f Gn x u P Yn [du] dx t n tn w = e t y n x tn P Yn [dy] dx π tn w y yx = tn Q n [dy] e t n x dx π = 1 Z n tn π w w e e yx P n [dy] e t n x dx = 1 tn ψ n x dx. Z n π t Making w go to + gives an equation for Z n = n π ψ nx dx. One concludes that: w P[W n w] = ψ nx dx ψ nx dx. This important property was not used in our previous works: to get the residue of deconvolution ψ n of a random variable X n by a large Gaussian variable of variance t n that is to say that one wants to remove a Gaussian variable of variance t n from X n, one can make the exponential change of measure 9, and add an independent Gaussian

13 MOD-GAUSSIAN CONVEGENCE AND APPLICATIONS FO MODELS OF STATISTICAL MECHANICS13 variable of variance t n : the random variable thus obtained, which is t n W n with the previous notation, has density proportional to ψ n w/t n dw. 3.. The residue of mod-gaussian convergence as a limiting law. We can now state and prove the main result of this Section. We assume the hypotheses A and B, and keep the same notation as before. Theorem 6. The following assertions are equivalent: i The sequence Y n t n n N is tight. ii The sequence Y n t n n N converges in law to a variable with density ψx/ ψx dx. iii The convergence ψ n ψ, which is supposed locally uniform on, also occurs in L 1. We shall then say that X n n N converges in the L 1 -mod-gaussian sense with parameters t n and limiting function ψ. In this setting, the residue ψ can be interpreted as the limiting law of X n n N after an appropriate change of measure. Proof. Since the Gaussian variable G n of variance t 1 n converges in probability to 0, Y n t n n N converges to a law µ if and only if W n n N converges to the law µ. If iii is satisfied, then by Proposition 5, lim P[W n w] = n w ψx dx ψx dx, so the cumulative distribution functions of the variables W n converge to the cumulative distribution function of the law ψx/ ψx dx, and ii is established. Obviously, one also has ii i. Finally, if iii is not satisfied, then by Scheffe s lemma one also has ψ n x dx ψx dx. However, by Fatou s lemma, ψx dx lim inf n ψ nx dx. Therefore, the nonconvergence in L 1 is only possible if ψx dx < lim sup n ψ nx dx. Thus, there is an ε > 0 and a subsequence n k k N such that k N, ψ nk x dx ε + ψx dx. Then, for all a, b, lim sup k b P[a W nk b] = lim sup a ψ n k x dx k ψ n k x dx ψx dx ε + < 1 ψx dx = b a ψx dx lim inf k ψ n k x dx which amounts to saying that W n n N and therefore Y n t n n N is not tight; hence, i implies iii. To complete this result, it is important to compare the two notions of complex mod- Gaussian convergence and of integral L 1 -mod-gaussian convergence. Though there

14 14 PIEE-LOÏC MÉLIOT AND ASHKAN NIKEGHBALI are no direct implication between these two assumptions, the following Proposition shows that the latter notion is a stronger type of convergence: Proposition 7. Let X n n N be a sequence that converges in the L 1 -mod-gaussian sense with parameters t n and limiting function ψ L 1. The estimate of precise large deviations is then satisfied. Proof. ecall than Z n = E[e X n /t n ] = P[X n t n x] = t n x P n [dy] = Z n t n x t n π ψ nx dx. We want to compute e y tn Q n [dy] = Z n x e t n u P Yn t n [du]. Suppose for a moment that we can replace the law of Y n t n by the one of W n = G n + Y n t n the previous computation. Then, one obtains from Proposition 5 in Z n x e t n u P Wn [du] = tn π x e t nu ψ n u du. Fix ε > 0. Since ψ n converges locally uniformly to the continuous function ψ, there is an interval [x, x + η] such that for n large enough and u [x, x + η], Therefore, for n large enough, ψx ε x+η x ψx ε ψx ε < ψ n u < ψx + ε. e t nu x+η du e t nu x+η ψ n u du ψx + ε e t nu du x x tnx e t n x ψx + ε tnx e t n x. Indeed, by integration by parts, x+η x e t n u t du is asymptotic to e n x t n x. On the other hand, since ψ n L 1 ψ, for the remaining part of the integral, x+η e t nu ψ n u du e t nx+η x+η ψ n u du e t nx+η x+η ψu du which is much smaller than the previous quantities. Therefore, assuming that one can replace Y n t n by W n, we obtain the asymptotics P[X n t n x] = tn x e ψx 1 + o1 πtn x

15 MOD-GAUSSIAN CONVEGENCE AND APPLICATIONS FO MODELS OF STATISTICAL MECHANICS15 for all x > 0; this is what we wanted to prove. Finally, the replacement Y n t n W n is indeed valid, because [ ] e t n u P Wn [du] = e t n u P[W n u] + t n u e t n u P[W n u] du x [ e t n u P[Y n /t n u] e t n u P Yn x t n [du] x ] x x + x t n u e t n u P[Y n /t n u] du by using on the second line the fact that both Y n t n and W n converge in law to the same limit, and therefore have equivalent cumulative distribution function on +. In the same setting of L 1 -mod-gaussian convergence, one has similarly the estimates on the negative part of the real line, and around 0, as described on page 4 in the setting of complex mod-gaussian convergence Application to the Curie-Weiss model. Consider i.i.d. Bernoulli random variables σi i 1 with P[σi = 1] = 1 P[σi = 1] = eα cosh α for some α. We set U n = i=1 n σi, so that E[e zu n coshz + α n ] = = cosh z + sinh z tanh α n cosh α [ ] n E e z U n n tanh α coshzn n 1/3 = 1/3 + sinhzn 1/3 tanh α e zn 1/3 tanh α [ ] log E e z U n n tanh α n 1/3 = n1/3 cosh α z sinh α 3 cosh 3 α z3 + o1 so one has complex mod-gaussian convergence of U n n tanh α n 1/3 with parameters n 1/3 cosh α and limiting function exp sinh α 3 cosh 3 α z3. If α = 0, then the term of order 3 disappears in the Taylor expansion of the characteristic function, and one obtains instead [ ] log E e z U n n 1/4 = n1/ z z4 1 + o1, hence a complex mod-gaussian convergence of X n = U n with parameters n 1/ and n 1/4 limiting function exp z 4 /1. Since this function restricted to is integrable, this leads us to the following result, which originally appeared in [EN78] without the mod- Gaussian interpretation: Theorem 8. Let X n = n 1/4 i=1 n σi be a rescaled sum of centred ±1 independent Bernoulli random variables. It converges in the L 1 -mod-gaussian sense, with parameters n 1/ and limiting function exp 1 z4. As a consequence, if Y n = n 1/4 M n is the rescaled magnetization of

16 16 PIEE-LOÏC MÉLIOT AND ASHKAN NIKEGHBALI a Curie-Weiss model CW 0,1 of parameters α = 0 and β = 1, then Y n /n 1/ converges in law to the distribution exp 1 x4 dx x4 exp 1 dx. Proof. The function ψ n t is in our case ψ n t = e t n 1/ cosh t n n 1/4, and we have seen that it converges locally uniformly to ψt = exp 1 t4. By Scheffe s lemma, to obtain the L 1 -mod-convergence, it is sufficient to prove that ψ nt dt converges to t4 exp 1 dt. This is a simple application of Laplace s method: ψ n t dt = e t n 1/ cosh t n n n 1/4 dt = n 1/4 e u cosh u du and the function u e u cosh u attains its global maximum at u = 0, with a Taylor expansion 1 u4 1 + ou4, see Figure FIGUE. The function f u = e u cosh u. Then, the exponential change of measure 9 gives a probability measure on spin configurations proportional to Yn 1 exp = exp n M n, n 1/ so it is indeed the Curie-Weiss model CW 0,1.

17 MOD-GAUSSIAN CONVEGENCE AND APPLICATIONS FO MODELS OF STATISTICAL MECHANICS17 It is easily seen that the proof adapts readily to the case where Bernoulli variables are replaced by so-called pure measures, so we recover all the limit theorems stated in [EN78, EN80]. However, by choosing the setting of mod-gaussian convergence, we also obtain new limit theorems for models that do not fall in the Curie-Weiss setting. The following result explains how it would work to replace the Bernoulli distribution by more general ones; cf. [KNN13, Proposition.]. Proposition 9. Let k be an integer, and let B n n 1 be a sequence of i.i.d random variables in L r for some r > k + 1, such that the first k moments of B 1 are the same as the corresponding moments of the Standard Gaussian distribution. Then the sequence of random variables 1 n 1/k+1 n B k k=1 n 1 converges in the mod-gaussian sense with parameters and limiting function t n = n k 1/k+1, θt = e itk+1 c k+1 k+1!, where c k+1 denotes the k + 1-th cumulant of B 1. When the random variables B n have an entire moment generating function, then one can replace t with it to obtain mod-gaussian convergence with the Laplace transforms. If B 1 is symmetric, then k is necessarily an odd number of the form s 1 and hence ψt = e 1s t s c s s!. In the case of the Bernoulli random variables, s = and c 4 = 1/1. In order to have our theorem of L 1 -mod-gaussian convergence to hold, we need to find conditions on the distribution of B 1 such that c s is negative and that ψ n converges to ψ. The conditions in [EN78] and [EN80] precisely imply these. But within our more general framework, following the discussion in Section 1., we could well imagine a situation which fulfils the assumptions of Theorem 6 but where the initial symmetric random variables are not necessarily i.i.d but simply independent or even weakly dependent. The following paragraph yields an example of such a setting Mixed Curie-Weiss-Ising model. Consider the one-dimensional Ising model of parameter α = 0, and β arbitrary. We have shown in Section the complex mod- Gaussian convergence of n 1/4 M n n N with parameters n 1/ e β and limiting function ψz = exp 3e 6β e β z 4 /4. estricted to, this limiting function is integrable, and again one has L 1 -mod-convergence. Indeed, recall that E[e tm n ] = Z ni, t, β Z n I, 0, β = 1 a + t, β λ+ t, β n 1 + a t, β cosh β λ t, β n 1. cosh β

18 18 PIEE-LOÏC MÉLIOT AND ASHKAN NIKEGHBALI It will be convenient to work with n 1/4 M n+1 instead of n 1/4 M n in order to work with n-th powers. Then, [ ] ψ n t = E e t M n+1 n 1/4 e n1/ e β z ψ n tdt = n1/4 λ+ u, β n eβu λ u, β n eβu a + u, β e + a u, β e du cosh β cosh β and for every parameter β 0, the functions u λ +u, β cosh β eβu e and u λ u, β eβu e cosh β attain their unique maximum at u = 0, see Figure 3 for the graph of the first function. FIGUE 3. The function f u, β = λ +u,β cosh β e eβ u using MATHEMATICA. Their Taylor expansions at u = 0 are respectively 1 3e6β e β 4 u 4 + ou 4 and tanh β + o1, so again by the Laplace method we get lim n ψ nt dt = ψt dt and the L1 - mod-convergence. As a consequence, consider the random configuration of spins σ on [[1, n]] with probability proportional to exp β n 1 σiσi + 1 i=1 n + 1 ne β σi. i=1 This model has a local interaction with coefficient β and a global interaction with coefficient 1 e β, so it is a mix of the Ising model and of the Curie-Weiss model. The previous

19 MOD-GAUSSIAN CONVEGENCE AND APPLICATIONS FO MODELS OF STATISTICAL MECHANICS19 discussion and Theorem 6 show that its magnetization satisfies the non standard limit theorem M n n 3/4 ψx dx n ψx dx with ψx = exp 3e6β e β x Sub-critical changes of measures. In the mixed Curie-Weiss-Ising model, one may ask what happens if instead of β and 1 one puts arbitrary coefficients for the e β local and the global interaction. More generally, given a sequence X n n N that converges in the L 1 -mod-gaussian sense with parameters t n and limiting function ψ, one can look at the change of measure Q γ n [dx] = E e γx tn [e γx n tn ] P n [dx] with γ 0, 1 for γ > 1, the change of measure is not necessarily well-defined, since the hypotheses A and B do not ensure that E[e γx n /t n ] < +. These subcritical changes of measures do not modify the order of magnitude of the fluctuations of X n, and more precisely: Theorem 10. Suppose that X n n N converges in the L 1 -mod-gaussian sense with parameters t n and limiting function ψ. Then, if X n γ n N is a sequence of random variables under the new probability measures Q γ n, it converges in the L 1 -mod-gaussian sense with parameters 1 γt n and limit ψ. Example. Consider a random configuration of spins σ on [[1, n]] with probability proportional to exp β n 1 n σiσi γ σi, n i=1 i=1 with γ < e β. The total magnetization of the system has order of magnitude n 1/, and more precisely, one has the central limit theorem M n n 1/ n N 0, 1 γe β e β, and in fact a L 1 -mod-gaussian convergence of M n n 1/4, with a limiting function that does not depend on γ. ψx = exp 3e 6β e β x 4 /4 Proof of Theorem 10. We denote as before Y n n N a sequence of random variables under the laws Q n = Q 1 n. We first compute the asymptotics of Z n γ = E[e γx n /t n ]: Z γ n = Z n E[e 1 γy n /t n ]

20 0 PIEE-LOÏC MÉLIOT AND ASHKAN NIKEGHBALI tn = π tn = π = o1 1 γ [ ] ψx dx E e t n1 γ Y n t n 1 + o1 [ ] ψx dx E e t n1 γ W n 1 + o1 by using on the third line the same argument as in the proof of Proposition 7 to replace Y n t n by W n ; and by using the Laplace method on the fourth line in order to compute e t n1 γx / ψ n x dx. The same computations give the asymptotics of E[e tx n+γx n /t n ] = Z n E[e ty n 1 γy n /t n ] [ tn = ψx dx E π tn = π = e t nt 1 γ e t nt Y n t n t n1 γ Y n t n ] 1 + o1 [ ] ψx dx E e t nt W n t n1 γ W n 1 + o1 1 ψt 1 + o1 1 γ with again a Laplace method on the fourth line. Since E[e txγ n ] = E[etX n+γx n /t n ], this shows the hypotheses A and B for the sequence X n γ n N. Then, since Y n t n n N converges in law, by using the implication ii iii in Theorem 6 for the sequence X n γ n N, we see that the mod-gaussian convergence of Laplace transforms necessarily happens in L 1. Z γ n 3.6. andom walks changed in measure. In this section, we shall make a brief excursion in the higher dimensions. Since we do not want to enter details on mod- Gaussian convergence for random vectors for which we refer the reader to [KN1] and [FMN13], we shall only consider the simple case X = X 1,..., X d is a random vector with values in d such that E[expz 1 X z d X d ] is entire in C d. We shall say that the sequence X n of random vectors converges in the complex mod- Gaussian sense with parameter t n and limiting function ψz 1,, z d if the following convergence holds locally uniformly on compact subsets of C d : ψ n t = E[expz 1 X n z d X n d z ] exp t z d n ψz 1,..., z d. In this vector setting, the assumptions A and B of Section 3 now simply amount to the fact that the convergence above holds locally uniformly for t = t 1,, t d d and that ψ n and ψ are both in L 1 d.

21 MOD-GAUSSIAN CONVEGENCE AND APPLICATIONS FO MODELS OF STATISTICAL MECHANICS1 Following the case d = 1 we denote P n the law of X n on d, Q n [dx] = [ E e x tn e X n tn ] P n [dx], and Y n a random variable under the new law Q n. Note that here again hypothesis B implies that Z n = E[e X n /t n ] is finite for all n N. Indeed, with the notation u, v = u 1 v u d v d, we have [ ] ψ tn t nt dt = E e t,x n dt d d [e X n π d/ E[e X n = E tn d e Xn tnt tn ] dt = Therefore, the new probabilities Q n are well-defined and d/ Z n = E[e X n /t tn n ] = π ψ nt dt. d Then it is clear that Proposition 5 holds with G n being a Gaussian vector with covariance matrix 1/t n I d where I d is the identity matrix of size d. Similarly one can establish an analogue of Theorem 6 in d. Let W n be a simple random walk on the lattice Z d : at each step, each of the d neighbors of the state that is occupied has the same probability of transition d 1. The d-dimensional characteristic function of W n = W n 1,..., W n d is Therefore, one has the asymptotics [ log E e ] z 1 W n 1 + +z d Wd n n 1/4 E[e z 1W n 1 + +z d W d n cosh z1 + + cosh z ] = d d = n log 1 + z z d dn 1/ + z z d 4 + o 4dn = n 1/ z z d d t n n. 1 n 3z z d dz z d 4 4d tn ]. + o1. One obtains a d-dimensional complex mod-gaussian convergence of X n = n 1/4 W n with parameters n1/ d and limiting function ψz 1,..., z d = exp 3 z z d d z z d 4 4d. In [FMN13], we used this mod-convergence to prove quantitative estimates regarding the breaking of the radial symmetry when one considers random walks conditioned to be of large size of order n 3/4 instead of the expected order n 1/. With the notion of L 1 -mod-gaussian convergence, one can give another interpretation, but only for d =

22 PIEE-LOÏC MÉLIOT AND ASHKAN NIKEGHBALI or d = 3. estricted to d, the limiting function is indeed not integrable for d 4: if t,..., t d [ 1, 1], then 3 t t d d t t d 4 3 t 1 + d 1 d t dt d 1t 1 + 3d 1. So, restricted to the domain [ 1, 1] d 1, ψt 1,..., t d K expat 1 4 bt 1 for some positive constants a, b and K; therefore, this function is not integrable. On the other hand, if d = or d = 3, then ψ is integrable on d, and one has L 1 - mod-gaussian convergence. Indeed, when d =, the limiting function is ψt 1, t = exp t t 4 + 6t 1 t, which is clearly integrable; and the residues [ ψ n t 1, t = E ψ nt 1, t dt 1 dt = e ] t 1 W n 1 +t W n n 1/4 e n1/ t 1 +t 4 converge locally uniformly on to ψt 1, t, but also in L 1. Indeed, cosh t 1 + cosh t n 1/4 n 1/4 = n 1/ 4 reaches its unique global maxi- and the function u 1, u cosh u 1+cosh u mum at u 1 = u = 0, with Taylor expansion cosh u1 + cosh u e u 1 +u 1 u u 4 + 6u 1 u + o u 4 96 around this point see Figure 4. n e n1/ t 1 +t 4 dt 1 dt e u 1 +u n 4 du 1 du, Thus, by using the multi-dimensional Laplace method, one sees that the limit of the integral ψ n t 1, t dt 1 dt is ψt 1, t dt 1 dt, and the L 1 convergence is shown. Similarly, when d = 3, the limiting function is ψt 1, t, t 3 = exp t 1t + t 1 t 3 + t t 3, and the following computation shows that it is integrable: 3 ψx, y, z dx dy dz = yz x dx dy dz e = 6 π = 3 π e y +z e yz 36 dy dz y + z r=0 π θ=0 e r4 sin θ 144 dr dθ

23 MOD-GAUSSIAN CONVEGENCE AND APPLICATIONS FO MODELS OF STATISTICAL MECHANICS3 FIGUE 4. The function f u 1, u = cosh u 1+cosh u e u 1 +u 4. = 1 3π r=0 π e r4 dθ dr < + θ=0 sin θ 1 since is integrable at 0. On the other hand, the residues sin θ [ ψ n t 1, t, t 3 = E e t 1 W 1 n ] +t W n +t 3 W3 n n 1/4 e n1/ t 1 +t +t 3 6 converge to ψt 1, t, t 3 locally uniformly on 3 and in L 1 3. Indeed, one has again cosh ψ nt 1, t, t 3 dt = n 1/ u1 + cosh u + cosh u 3 e u 1 +u +u 3 n 6 du and the function in the brackets reaches its unique maximum at u 1 = u = u 3 = 0, with Taylor expansion corresponding to the limiting function ψ after application of the Laplace method. The multidimensional analogue of Theorem 6 thus yields the following multidimensional extension of the limit theorem for the Curie-Weiss model: Theorem 11. Let W n be a simple random walk in dimension d 3. If V n is obtained from W n by a change of measure by the factor expd W n /n, then V n n 3/4 n ψx dx 3 ψx dx, where ψx = exp x 4 /1 in dimension 1, and ψ is given by Formulas 10 and 11 in dimension and 3. emark. Suppose d =. Then, there is a limit in law not only for V n n 3/4, but in fact for the whole random walk V k n 3/4 k n, viewed as a random element of C +, or of the Skorohod space D +,, see Figure 5.

24 4 PIEE-LOÏC MÉLIOT AND ASHKAN NIKEGHBALI FIGUE 5. A -dimensional random walk changed in measure by e W n /n, here with n = LOCAL LIMIT THEOEM AND ATE OF CONVEGENCE IN THE ELLIS-NEWMAN LIMIT THEOEM We keep the same notation as before and note I n = ψ nx dx and I = ψx dx. In this section we wish to provide a quick approach based on Fourier analysis, 1 to compute the Kolmogorov distance between the rescaled magnetization Y n /n 1/ = M n /n 3/4 in the Curie-Weiss model and the random variable W with density ψx/i, where ψx = exp x 4 /1. This problem was recently solved in [EL10] using Stein s method. As in [EL10], our method would cover many more general models as well: it is just a matter of specializing Lemma 1 and Lemma 13 below which are stated in all generality. to prove a new local limit theorem for the rescaled magnetization n 1/4 M n in the Curie-Weiss model. Here again we shall indicate how one can establish local limit theorems in more general situations Speed of convergence. Getting back to our special case of the Curie-Weiss model, we denote X n = 1 n 1/4 i=1 n B i a scaled sum of ±1 independent Bernoulli random variables; Y n the random variable with modified law ψ n x I n Q n [dy] = e y n 1/ P n [dy] ] ; E [e X n n 1/ G n an independent Gaussian random variable of variance 1 ; and W n 1/ n = Y n + G n 1/ n. It follows from the previous results that the law of W n has density cosh = 1 e n1/x I n x n 1/4 n,

25 MOD-GAUSSIAN CONVEGENCE AND APPLICATIONS FO MODELS OF STATISTICAL MECHANICS5 which converges in L 1 towards the law ψx I = I 1 e x4 1. We hence wish for an upper bound for the Kolmogorov distance between Y n and W n 1/. For this we shall need the following general lemmas. Lemma 1. Consider the two distributions W n = ψ nx dx I n and W = distance between them is smaller than ψ ψ n L o1. I ψx dx I. The Kolmogorov Proof. Fix a, and suppose for instance that ψx dx ψ nx dx. We have a F Wn a F W a = ψ a a a nx dx ψx dx ψx dx ψx dx + I n I n a a = ψ nx ψx dx + ψx dx ψx ψ nx dx I n I I n a ψx ψ nx dx + ψx ψ nx dx I n I n a ψx ψ nx dx ψ ψ n L 1. I n I n Writing F W a F Wn a = 1 F Wn a 1 F W a, one sees that the inequality is in fact valid with an absolute value on the left-hand side. Since I n = I 1 + o1, this shows the claim. If ψ nx dx ψx dx, it suffices to exchange the roles played by ψ n and ψ to get the inequality. I n I The asymptotics of the L 1 -norm ψ ψ n L 1 in the Curie-Weiss model are computed as follows. Noting that one always has ψ n x ψx, it suffices to compute ψ n x dx = e n1/ x n n coshx n 1/4 dx = n 1/4 e u coshu du. By the Laplace method see [Zor04, Formula 19.17, p ], the asymptotics of the integral is n /4 Γ n /4 Γ smaller terms. 10 The first term corresponds to I = ψx dx = e x4 /1 dx. As a consequence, ψ ψ n L 1 = 1 1 Γ 3 4 I n 1/ 5 Γ o1. The main work now consists in computing d Kol Y n n 1/, W n. We start by a Lemma which is a variation of arguments used for i.i.d. random variables in [Tao1, p. 87].

26 6 PIEE-LOÏC MÉLIOT AND ASHKAN NIKEGHBALI In the following, given a function f L 1, we write its Fourier transform f ξ = f x eiξx dx. ecall that the function υξ = {e 1 1 4ξ if ξ < 1, 0 otherwise. is even, of class C and with compact support [ 1, 1 ]. We set ρ = υ, so that ρ x = 1 1 υξ e ixξ dξ π 1 by the Fourier inversion theorem. By construction, the Fourier transform of ρ has support equal to [ 1, 1 ]. Set now ρx = ρ x ρ y dy. By construction, ρ is smooth, even, non-negative and with integral equal to 1. Moreover, ρ is up to a constant equal to ν νξ, so it has support included into [ 1, 1]. The convolution of ρ with characteristic functions of intervals will allow us to transform estimates on test functions into estimates on cumulative distribution functions. More precisely, for a and ε > 0, set ρ ε x = 1 ε ρ x ε, and φ a,εx = φ ε x a, where φ ε is the function 1,0] ρ ε. One sees φ a,ε as a smooth approximation of the characteristic function 1,a]. [ ] For all a, ε, φ a,ε has Fourier transform compactly supported on 1 ε, 1 ε. Moreover, it has negative derivative, and decreases from 1 to 0. Later, we will use the identity φ ε εx = φ 1 x = φx. On the other hand, we have the following estimates for K > 0 we used SAGE for numerical computations: ρ K = 1 1 πk υ ξ e ikξ dξ 1 1 πk υ ξ dξ = K ; ρ y dy = 1 1 υξ dξ = π 0 Therefore, for any K > 0, 1 ρk = ρ K = ρ K ρ y dy 99 K 4 ; φk = 1 φ K = 0 ρk + y dy 33 K 3. 0 Lemma 13. Let V and W be two random variables with cumulative distribution functions F V and F W. Assume that for some ε > 0 E[φ a,ε V] E[φ a,ε W] Bε,

27 MOD-GAUSSIAN CONVEGENCE AND APPLICATIONS FO MODELS OF STATISTICAL MECHANICS7 where the positive constant B is independent of a. We also suppose that W has a density w.r.t. Lebesgue measure that is bounded by m. Then, sup F V a F W a B + 10m ε. a Proof. Fix a positive constant K, and denote δ = sup a F V a F W a the Kolmogorov distance between V and W. One has F V a = E[1 V a ] E[φ a+kε,ε V] + E[1 φ a+kε,ε V 1 V a ] E[φ a+kε,ε W] + E[1 φ a+kε,ε V 1 V a ] + Bε. The second expectation writes as E[1 φ a+kε,ε V 1 V a ] = 1 φ a+kε,ε x 1,a] x f V x dx = 1 φ a+kε,ε x 1,a] x F V x dx = φ a+kε,ε x 1,a]x F V x dx + 1 φ a+kε,ε x 1 a x F V x dx. For the first integral, since F V x F W x δ and the derivative of φ a+kε,ε is negative, an upper bound on I 1 is φ a+kε,ε x 1,a]x F W x dx δ φ a+kε,ε x 1 V ax = φ a+kε,ε x 1,a]x F W x dx + 1 φ a+kε,ε a δ = φ a+kε,ε x 1,a]x F W x dx + 1 φ K δ. As for the second integral, it is simply 1 φ a+kε,εaf V a, and by writing F V a F W a + δ, one gets the upper bound on I 1 φ a+kε,ε x 1 a x F W x dx + 1 φ a+kε,ε a δ = 1 φ a+kε,ε x 1 a x F W x dx + 1 φ K δ. One concludes that E[1 φ a+kε,ε V 1 V a ] E[1 φ a+kε,ε W 1 W a ] + 1 φ Kδ. On the other hand, if m is a bound on the density f W of W, then E[φ a+kε,ε W 1 W a ] = m mε a φ a+kε,ε y f W y dy a 0 φ ε y a Kε dy = m φu K du mε K + 4.8, 0 φ ε y Kε dy