Lee–Yang Property and Gaussian Multiplicative Chaos

The Lee–Yang property of certain moment generating functions having only pure imaginary zeros is valid for Ising type models with one-component spins and XY models with two-component spins. Villain models and complex Gaussian multiplicative chaos are two-component systems analogous to XY models and related to Gaussian free fields. Although the Lee–Yang property is known to be valid generally in the first case, we show that is not so in the second. Our proof is based on two theorems of general interest relating the Lee–Yang property to distribution tail behavior.


Overview
The original theorem of Lee and Yang [17] on zeros of the partition function of a spin-1 2 ((±1) −valued) Ising model {σ v } with pair ferromagnetic interactions implies that the moment generating function E e z X of X = λ v σ v with λ v ≥ 0 has only pure imaginary zeros (PIZ) in the complex z−plane. As we explain in more detail below, this has been extended [2,18] for certain twocomponent XY models where S v takes values in the unit circle of R 2 . It also extends to D−component models with D = 3 [2], known as classical Heisenberg models, but, as far as we know, has not been extended to D ≥ 4-see Remark 4 in [2].
In this paper, we consider two interesting two-component systems, namely the Villain model [27] and complex Gaussian multiplicative chaos, related to the Gaussian Free Field (GFF) [15] and whether they have a Lee-Yang PIZ property. For Villain models, like for XY models, such a property is expressed in terms of the cosines of the angular variables; for Gaussian multiplicative chaos, it is expressed in terms of the real part of the complexvalued field and hence in terms of a sine-Gordon field-since the partition function of subcritical imaginary Gaussian multiplicative chaos in an external field equals that of a sine-Gordon field. A PIZ property seems plausible for both systems because of the connection between GFF and the XY model via the spin-wave conjecture [4,21] on the one hand and the connection between GFF and the Villain Gaussian-like distribution on the other. Although Villain models are known to have the PIZ property (see p. 636 of [8] and Theorem 3 below) we prove that complex Gaussian multiplicative chaos in general does not (see Propositions 18 and 20 below).
Our analysis relies on two results relating for (symmetric) random variables X n , the PIZ property, tail behavior and convergence in distribution of X n to X , which may be of independent interest. The first of these (Theorem 7) is a new result which states that the PIZ property plus a sub-Gaussian tail bound for each symmetrically distributed X n (with no uniformity in n) implies not only that X has the PIZ property but also, surprisingly, sub-Gaussian tail behavior. The second (Theorem 11) states that any X with tail behavior strictly between Gaussian and exp −c |x| 1+ε cannot have the PIZ property; it follows directly from a result of Goldberg and Ostrovskii [10] (see Theorem 14.4.2 of [19]).

Introduction
The Lee-Yang theorem, first obtained by Lee and Yang when studying phase transitions in the classical Ising model, states that all the zeroes of the partition function of the Ising model as a function of the external magnetic field lie on the imaginary axis [17]. Lee and Yang proved this result only for the spin-1 2 (i.e., (±1) −valued) Ising model, but later it was extended by Griffiths to spin-n 2 models [11], and then by Newman to ferromagnetic Ising models with quite general single spin distributions [22] (see also [18]). Since then, it has been applied to prove the properties of the infinite volume limit and existence of a mass gap under an external magnetic field [7,12,24], and to prove correlation inequalities [3,23]-see also [6] for a general review of the Lee-Yang type theorem and their applications.
For two-component ferromagnets, analogous Lee-Yang type theorems were first proved by Dunlop and Newman [2] for the classical XY model (see also [26]) and then by Lieb and Sokal [18] for generic two-component ferromagnets with quite general single spin distributions. We now recall the Lee-Yang theorem for the classical XY model. Given a finite graph G = (V (G) , E (G)), we denote the spin variable at each v ∈ V (G) by S v = S 1 v , S 2 v in the unit circle. Given real numbers {J e } e∈E , the classical XY model on G is defined by a Gibbs measure with the Hamiltonian H (S) given by where S i · S j = S 1 i S 1 j + S 2 i S 2 j , and with Z G = Z G (T, {J e }) the normalization constant that makes (1) a probability measure. This is an XY model on G with free boundary conditions; for a discussion of other boundary conditions where the conclusions of the next theorem remain valid, see Remark 5 below.
Theorem 1 ([2,18]). Suppose that J e ≥ 0 for all e ∈ E and λ v ≥ 0 for all v ∈ V(G), then has only pure imaginary zeros (namely, it is not zero when z > 0 or z < 0). Here E dentoes the expectation with respect to the probability measure (1).
The classical Villain model is another two-component spin model which is closely related to the XY model [27]. Given any finite graph G = (V (G) , E (G)), and positive numbers {J e } e∈E(G) , a Villain model on G is defined by the Gibbs measure where is a periodized Gaussian, and Z G is the normalizing constant. Let ( v : v ∈ V (G)) be jointly distributed by the Gibbs measure (3). Given non-negative real numbers The validity of the Lee-Yang type property as well as correlation inequalities for Villain models, due to Fröhlich and Spencer and to Bellissard ( [8] and ref. 34 there), is discussed and a proof is sketched on pp. 635-636 of [8]. Since that Lee-Yang property does not seem to be widely known, we present it as the next theorem and provide a detailed proof in Sections 4 and 5. The Lee-Yang property does not seem to follow from the Lieb-Sokal approach since unlike ( 2), the Villain model (3) is not written as a Gibbs measure with ferromagnetic pair interactions.
has only pure imaginary zeros as a function of complex z. Here E dentoes the expectation with respect to the probability measure (3). Remark 4. By extending Theorem 7 below to a multivariable version, the same proof leads to a more general result, that is, for any {z v } v∈V(G) such that z v > 0 ∀v, Remark 5. Theorem 3 is stated for the Villain model (3) with free boundary condition. The same result holds for some other boundary conditions, including periodic (i.e., the Villain model on a torus), or for Dirichlet-type boundary condition where v for all v in the boundary are equal with a value that is uniformly distributed on the unit circle. Periodic boundary condition on a torus simply correspond to a different choice of G than do free boundary conditions, while Dirichlet type boundary conditions can be handled as a corollary of Theorem 3 by adding couplings between boundary vertices and letting the coupling magnitudes tend to infinity. We note that the Lee-Yang property of Theorem 1 for XY models will also be valid for such boundary conditions.
Obtaining new Lee-Yang type theorems for lattice models can be useful to derive new results for continuum field theories. For example, Simon and Griffiths [25] proved a Lee-Yang result for the continuum φ 4 2 Euclidean field theory by first obtaining a new Lee-Yang result for φ 4 lattice models. Indeed, part of our motivation comes from the so-called spin-wave conjecture [4,21], which states that for both the XY and the Villain model at temperature T less than some critical value, on large scales the angular variables v behave like a Gaussian Free Field (GFF) modulo 2π . This suggests that the spin field (S · in the XY model and exp (i · ) in the Villain model) may behave like a version of complex Gaussian multiplicative chaos (see e.g., [15]); or, after a duality transformation, a version of the Sine-Gordon field [9]. We have not obtained a Lee-Yang property for complex Gaussian multiplicative chaos and even ruled it out (as we discuss below) in a certain parameter range. However that parameter range does not correspond to very low temperature T , so it may still be that there is a low T Lee-Yang property.
Our approach is based on showing that a certain set of properties of moment generating functions (see Definition 6) is preserved under convergence in distribution (Theorem 7). This is more than the requirement that all zeros are pure imaginary-it further requires a sub-Gaussian tail for the distribution. In Sections 4 and 5 we prove Theorem 3 by approximating the Villain model on any finite graph by one dimensional XY spin chains. Since the XY models do satisfy the Lee-Yang property, we obtain the result by applying the convergence Theorem 7.
Based on the spin-wave conjecture, it is natural to ask whether there is a Lee-Yang property for complex Gaussian multiplicative chaos, namely exp (iβh) where h is a two dimensional GFF. In Section 6 we give a negative answer to this question when β ∈ 1, √ 2 -the Lee-Yang property does not hold for complex Gaussian multiplicative chaos with such values of β. The proof is based on a general theorem that shows that the pure imaginary zeros property does not hold for random variables with tail slower than Gaussian (Theorem 11), and an explicit computation of the tail probability for the integral of complex Gaussian multiplicative chaos. We also study the so-called discrete complex Gaussian multiplicative chaos on finite graphs, which roughly speaking, is exp (iβh) where h is a discrete Gaussian free field (DGFF) with certain boundary condition where the Lee-Yang property might be expected to hold. We show (see Proposition 20) that for discrete complex Gaussian multiplicative chaos with β ∈ 1, √ 2 , the Lee-Yang property cannot hold on all finite graphs. Interestingly, we use there a corollary of the weak convergence result of Theorem 7 to rule out the Lee-Yang property.
Complex Gaussian multiplicative chaos corresponds to a type of a Sine-Gordon field (see, e.g., [9] and [1]). It has been pointed out to us by T. Spencer [personal communication] that because of the relation of the Lee-Yang property to exponential decay or existence of a mass gap for Sine-Gordon fields, how the validity of the Lee-Yang property depends on the parameter β is of some interest. The Sine-Gordon betaparameter used in [9], which we denote here byβ, is related to the β we use for complex Gaussian multiplicative chaos byβ = 2πβ 2 . Thus our region β ∈ 1, √ 2 of nonvalidity of the Lee-Yang property corresponds toβ ∈ (2π, 4π ).

Lee-Yang Property and Weak Convergence
Let μ be a probability measure on R and X be a random variable on some probability space ( , F, P) with distribution μ. Definition 6. We say μ (or X ) is of Lee-Yang type (and write μ ∈ L) if (1) X has the same distribution as −X only has zeros on the pure imaginary axis.
The next theorem states that the Lee-Yang type property is preserved under weak convergence and helps explain why the sub-Gaussian property (2) is built into Definition 6.
Theorem 7. Suppose for each n ∈ N, μ n ∈ L, and μ n converges weakly to the probability measure μ. Then μ ∈ L.
Remark 8. We note in particular that the limiting measure must satisfy μ (r, ∞) ≤ exp −cr 2 for some c > 0. We will use this fact later. This property of μ seems a-priori surprising, because if one does not assume μ n has the pure imaginary zero property (Definition 6,(3)), the conclusion is not true since the constant b in Definition 6,(2) may depend on n.
Remark 9. Approximation schemes play a major role in classical work (e.g., [22,25]) on Lee-Yang type theorems for Ising-like systems. It seems that Theorem 7 fits together with that work primarily as a helpful tool-i.e., it suffices to show weak convergence, without any extra moment or moment generating function estimates, to guarantee that the limit system will have the desired Lee-Yang property. (1) X n has the same distribution as −X n , and If X n converges in distribution to X and E exp bX 2 = ∞ for all b > 0, it follows that for all but finitely many n, E exp (z X n ) has some zeros that are not purely imaginary.
Proof. The proof is by contradiction. If the conclusion were not valid, then there would be a subsequence X n k with distributions μ k ∈ L which would converge weakly to the distribution μ of X with μ / ∈ L. That would contradict Theorem 7, which completes the proof.
The next theorem relates the Lee-Yang property to the distribution tail behavior and explains further why the sub-Gaussian property (2) is natural in Definition 6. As we explain below, it follows directly from a theorem of Goldberg and Ostrovskii [10] (see Theorem 14.4.2 of [19]). Theorem 11. Suppose the random variable X satistifes the following two properties: (1) Ee b|X | a < ∞ for some b > 0 and a > 1, Then Ee z X has some zeroes that are not purely imaginary.

Remark 12.
A natural question is how much Property (1) of this theorem can be weakened without changing the conclusions. The answer is not very much. This can be seen by constructing examples of random variables X where E e b|X | log|X | < ∞ for some b > 0 but E e z X has no zeros that are not purely imaginary. Probably the simplest example is a Poisson random variable X , where E e z X = exp (λ (e z − 1)) has no zeros at all. One can also extend this example to obtain symmetric random variables with Poissontype tail behavior whose moment generating functions do have many zeros, all purely imaginary.
We now turn to the proof of Theorems 7 and 11. The proof of Theorem 7 is based on the uniform convergence of entire functions and an application of Hurwitz Theorem. The proof will be given in two steps. In the first step we prove Theorem 7 under the additional assumption that sup n E X 2 n < ∞. In the second step we show sup n E X 2 n < ∞ automatically holds.
The key ingredient for the proof is the following Proposition (see Proposition 2 of [23]). We include a proof for completeness.
Proof. We first note that f is entire of (exponential) order 2 (and finite type), since .
By the Hadamard factorization theorem, where P 2 is a quadratic polynomial, m 0 is the degree of zero of f at the origin, z j are the other zeros and j z j −2 < ∞. Since μ is a symmetric probability measure, Proof of Theorem 7. Let f n (z) = E exp (z X n ) . We first prove Theorem 7 assuming sup n E X 2 n < ∞. We claim that it suffices to prove sup n | f n (z)| < ∞ uniformly on compact sets of z (6) and that We now explain why (6) and (7) suffice to imply the conclusion of Theorem 7. First note that the validity of (1) of Definition 6 for each μ n implies it for μ. As X n converges in distribution to X , f n → f on the pure imaginary axis, and (6) implies that f extends to an entire function with f n → f uniformly on compact sets. Moreover, by Hurwitz' Theorem, open zero-free regions for all f n (e.g., C\iR) are zero-free for f . This verifies (3) of Definition 6 for X . Finally, (7) implies (2) of Definition 6 for X (e.g., by taking any b ∈ 0, b ) and thus X ∈ L. We next claim that (6) and (7) are direct consequences of Proposition 13. Apply Proposition 13 and use the fact that 1 + z 2 /y 2 ≤ exp |z| 2 /y 2 to see that where B (n) and y (n) k for X n correspond to B and {y k } for X in Proposition 13. Since we assumed sup n E X 2 n < ∞, we conclude (6). To prove (7), note that (8) implies that the tail of X n is dominated by the tail of Y n ∼ N (0, Var [X n ]). Therefore we conclude (7) with any b < 2 sup n E X 2 n −1 .
Finally we prove that convergence of X n to some X in distribution implies that sup n E X 2 n < ∞. We will argue by contradiction. Suppose that X n ∈ L, sup n E X 2 n = ∞ and X n converges to some X in distribution. By taking a subsequence and applying Proposition 13 we may assume that We also know from the convergence in distribution of X n that f n (it) → f (it) uniformly on compact subsets of t ∈ R. Both f n (it) and f (it) are real and continuous in t so that there exists ε > 0 such that f (it) > 0 for t ∈ [0, ε] . Since f n iy However, by Proposition 13, for t ∈ (0, y 1 ] , By (9) this goes to zero as n → ∞. This contradicts f (it) > 0 for t ∈ [0, ε] and completes the proof.
Finally we note that Theorem 11 is an immediate consequence of the following proposition, which follows from the Goldberg-Ostrovskii result [10] stated as Theorem 14.4.2 in [19]. We also note that there is a typographical error in [19] and k a 2 k < ∞ should be replaced there by k 1/a 2 k < ∞. Proposition 14. Suppose that the random variable Y satisfies Property 1 of Theorem 11, and Ee zY has only pure imaginary zeroes. Then Y − EY belongs to the class L.
Proof. The proof follows directly from a result of Goldberg and Ostrovskii (see Theorem 14.4.2 of [19]). Arguing as in the proof of Proposition 13, let f (z) = Ee zY ; then by Young's inequality, By Property 1 of Theorem 11 we have that Therefore f is an entire function of finite (exponential) order a . By the Goldberg-Ostrovskii result it follows that with α ∈ R, B ≥ 0 and k |y k | −2 < ∞. But this implies that EY = α and Y − EY satisfies all the properties to be in L.

1 D XY Spin Chain
In this and the next sections we will show that the distribution of (5) from the Villain model is of Lee-Yang type as defined in Definition 6, which implies Theorem 3. The proof uses the following scaling limit result for a one dimensional XY model. Let G n denote the graph whose vertex set V (G n ) is {0, 1/n, 2/n, ..., 1} and edge set is {{( j − 1) /n, j/n} : j = 1, ..., n}. We assign to each i ∈ V (G n ) a spin variable S n i in the unit circle, with the corresponding angle θ n i ∈ (−π, π], and consider the XY model on G n defined by the Gibbs measure with Hamiltonian For the remainder of this section, we write B or B n for 1/T or 1/T n . We note that in the next proposition, the parameter b > 0 will eventually be proportional to one of the J e 's in (4) of the Villain model-see Remark 16 below.
Proof. One can view θ n with S n = (cos θ n , sin θ n ) as a one dimensional Markov process with transition density given by and initial distribution uniform on (−π, π]. We provide a sketch of a proof of the scaling limit result. Based on a standard convergence result for discrete time Markov chains (see, e.g., Theorem 17.28 of [14]), it suffices to prove that for any twice differentiable function f : S 1 → R, the transition operator K B defined by Laplace's method (see e.g., [5]) yields π −π e B n cos φ dφ = e B n 2π Also, the value of θ that minimizes K B θ, θ f θ is given by Let g θ = B −1 n log exp B n cos θ − θ f θ , then Laplace's method yields In particular, for any f ∈ H 1 S 1 , We will also need the following version of Proposition 15 with Dirichlet boundary conditions. The proof follows from essentially the same arguments as for Proposition 15-see also Remark 16 for the last statement of the following proposition.

Proof of Theorem 3
We now apply Theorem 7, Proposition 15 and Proposition 17 to prove Theorem 3. We first prove Theorem 3 when G is a single edge, i.e., G = ({x, y} , {e}), and without loss of generality we can identify x and y with 0 and 1 in the unit interval [0, 1]. Consider an XY model on G n of Section 4, defined by the Gibbs measure (10) with B = B n = n J −1 e . By Proposition 15, S n 0 , S n n converges in distribution to (exp (i (0)) , exp (i (1))), where (0) is uniform in [−π, π) and exp (i ( (1) − (0))) is distributed as exp i B (1) √ J e . In other words, (exp (i (0)) , exp (i (1))) has the probability density which has the same distribution as for the Villain model on G. By Theorem 1, for all λ 0 , λ 1 ≥ 0, the distribution of  (1) is also of Lee-Yang type. This proves Theorem 3 in the special case when G is a single edge.
In the general case (see Figure 1 for an illustration), given G = (V (G) , E (G)) we first replace each vertex v ∈ V (G) by several new vertices-one denoted v * and then one more denoted (v, e, 0) for each e ∈ E (G) incident on v, which we write as e ∼ v. (We will also denote (v, e, 0) for e = {v, w} by (w, e, n).) We then create one new edge between {v * , (v, e, 0)} v∈V and each (v, e, 0). Each e = {v, w} ∈ E (G) is replaced by a collection V (e, n) of n − 1 new vertices which will be labelled (v, e, 1) , (v, e, 2) , ..., (v, e, n − 1) (and also in opposite order (w, e, 1) , (w, e, 2) , ..., (w, e, n − 1)) and a collection E (e, n) of n new edges: , 0) , (w, e, 1)} , ...). This defines a new graph G * n = V G * n , E G * n (Fig. 1). As indicated by Figure 1, to prove Theorem 3, it suffices to prove that for all non- has only pure imaginary zeroes.
For any fixed J > 0 and {J e } e∈E(G) , let B e = n J −1 e , and consider an (inhomogeneous) XY model on G * n , defined by the Gibbs measure with Hamiltonian Applying the Lee-Yang theorem for the XY model (Theorem 1), we see that for all has only pure imaginary zeros. Given e = (u, v) ∈ E (G), and any θ u , θ v ∈ (−π, π], we also define Using the definition of H J XY , we see that has only pure imaginary zeros. Applying Proposition 17, by taking n → ∞ we have Therefore, omitting all the superscripts n in θ n · , the integral (13) converges to By Theorem 7, we see that (14), as a function of z, only has pure imaginary zeroes. Applying Laplace's method, we can multiply (14) by the right J −dependent factor and let J → ∞ to recover the Villain model on G, finishing the proof of (12).

Continuum complex Gaussian multiplicative chaos.
In this subsection we apply Theorem 11 to complex Gaussian multiplicative chaos. Let D r ⊂ C, r > 0 be the disk of radius r centered at the origin. For β > 0, complex Gaussian multiplicative chaos in D r is defined as exp (iβh r ), where h r is a Dirichlet zero boundary condition GFF in D r .
As was discussed in [15], for β ∈ 0, √ 2 the (complex-valued) measure e iβh r (x) dx is well-defined and is absolutely continuous with respect to Lebesgue measure.
We will obtain results for complex Gaussian multiplicative chaos in the whole plane, which is defined informally as exp (iβh), where h is a GFF in all of R 2 . Mathematically, the GFF in R 2 is not well-defined as a random generalized function-it is only defined up to an additive constant. Therefore, a priori, exp (iβh) is only defined up to a multiplicative constant (on the unit circle in C). In this section, we obtain the measure exp (iβh) dx by defining for any bounded simply connected domain U ⊂ R 2 , the random variable as r → ∞. The existence of the limit of the modulus of U exp (iβh r (x)) dx as r → ∞ is proved in the appendix of [16] (and an alternative proof may be obtained by using a construction similar to that in [15]). To give a precise definition of the limiting field without analyzing multiplicative factors, we may simply proceed as follows. Let be a random variable uniformly distributed on (−π, π] that is independent of the entire collection of random fields {h r : r > 0}. Then define U exp (iβh (x)) dx as the r → ∞ limit in distribution of Based on the spin-wave picture of the XY and Villain models, it a priori seems reasonable to conjecture that complex Gaussian multiplicative chaos with any of the boundary conditions mentioned in Remark 5 also satisfies a Lee-Yang property. Here we restrict attention first to the whole-plane field (without boundary). Let U ⊂ R 2 be a bounded simply connected domain with smooth boundary, and for any non-negative bounded continuous function λ : D r → R, define the random variable If the Lee-Yang property were valid, then E exp (z X ∞ ) as a function of z would only have pure imaginary zeroes. The next theorem disproves the Lee-Yang property for complex Gaussian multiplicative chaos, when β ∈ 1, √ 2 .
Proposition 18. Let X ∞ be as in (15) with λ (x) ≡ 1 for any bounded simply connected domain U with smooth boundary. Then for any β ∈ 1, √ 2 , E exp (z X ∞ ) has some zeroes that are not purely imaginary.
To prove Proposition 18, we start with the following tail estimate for
For a lower bound, we use the fact that W U is equidistributed with e iφ W U for any real φ, which implies that X ∞ = W U and Y ∞ = W U are equidistributed. Thus Then by Proposition 19, E e b |W U | 2 /2 = ∞ for any b > 0 and now by Theorem 11 we conclude that E exp (z X ∞ ) must have some zeroes that are not purely imaginary.

Discrete complex Gaussian multiplicative chaos.
Let G = (V, E) be a planar domain, i.e., V ⊂ Z 2 and E is the set of nearest neighbor edges between vertices of V. We use ∂G to denote the boundary vertices of G, i.e., those that are connected to Z 2 \ V by one or more edges. Consider a random function h (i), i ∈ V, such that h (i) ∈ (−π, π], distributed according to the (conditional on φ) joint density where φ is the value of a uniformly distributed in (−π, π] random variable and B i j are some positive numbers. Heuristically, one can construct h by first choosing φ uniformly in (−π, π], then sampling a discrete Gaussian free field (DGFF) on G with Dirichlet boundary condition φ (denoted as h φ ); finally, we obtain h by taking the value of h φ modulo 2π . We call the random function h thus defined the discrete complex Gaussian multiplicative chaos on G, since it is a discrete analogue of the complex Gaussian multiplicative chaos considered in Section 6.1.
Motivated by the spin-wave picture of the XY and Villain models, it is natural to ask whether the discrete complex Gaussian multiplicative chaos with any boundary conditions mentioned in Remark 5 also satisfies a Lee-Yang theorem. Here we restrict attention to Dirichlet boundary conditions, with h sampled from the discrete complex Gaussian multiplicative chaos (19), and we define This is a discrete analogue of the random variable X ∞ defined in ( 15). In this section we show that the Lee-Yang property cannot be valid in general, by giving a family of counter-examples in the next proposition. We first introduce some notation. Recall that D r is the disk of radius r centered at the origin. Given x ∈ D r , we denote by C r (x) the conformal radius of D r from x. Proposition 20. Take B i j = β −2 , G n,r = D nr ∩ Z 2 and g = 2/π . Let h n,r be sampled from (19) on G n,r and M n,r = x∈D n ∩Z 2 1 n 2 exp ih n,r (x) + β 2 2 g log n − gC r x n .
Then for β ∈ 1, √ 2 , large enough r and then large enough n, E exp z M n,r has some zeros that are not purely imaginary.
We start by proving the following limit theorem for M n,r .

Lemma 21.
For any r > 0, as n → ∞, M n,r converges in distribution to D 1 e i(βh r (x)+ ) dx, (20) where h r is a continuum GFF on D r with zero boundary, and is uniformly distributed in (−π, π]. Proof. Recall that for any c > 0, both the discrete and continuum GFF have the Gibbs-Markov property: a GFF/DGFF with Dirichlet boundary condition c equals in law a zero boundary GFF/DGFF plus the constant c. Therefore it suffices to prove the lemma with = 0 in both (19) and (20). Next definê M n,r = x∈D n ∩Z 2 1 n 2 exp ih n,r (x) + β 2 2 g log n − gC r x n .
Since M n,r = M n,r is a bounded continuous function ofM n,r , it suffices to show thatM n,r converges in distribution to To see this, we compute the moments ofM n,r as E M k n,r = Let G nr be the Dirichlet Green's function on D nr ∩ Z 2 . We have We now use the standard estimates for the lattice Green's function (see, e.g., [20], Chapter 1.6), where g D r is the Dirichlet Green's function in D r with zero boundary condition. Substituting these into (21) and using the convergence of Riemann sums to integrals concludes the proof.
Proof of Proposition 20. Notice that when B i j = β −2 , Lemma 21 implies that as n → ∞, M n,r converges in distribution to the integral X r := D 1 e i(βh r (x)+ ) dx.
Moreover, by the construction and the rotational invariance of the whole plane complex Gaussian multiplicative chaos, as r → ∞, D 1 e i(βh r (x)+ ) dx converges in distribution to W D 1 e i , which is equidistributed with W D 1 . Thus X r converges in distribution to X ∞ . However, by Eq (18), when β > 1, E e b|X ∞ | 2 = ∞ for all b > 0. Applying Corollary 10 completes the proof.