arxiv: v1 [math.pr] 27 May 2013

Transcription

1 Gaussian multiplicative chaos and applications: a review arxiv: v [math.pr] 27 May 203 May 28, 203 Rémi Rhodes Université Paris-Dauphine, Ceremade, F-7506 Paris, France rhodes@ceremade.dauphine.fr Vincent Vargas CNRS, UMR 7534, F-7506 Paris, France Université Paris-Dauphine, Ceremade, F-7506 Paris, France vargas@ceremade.dauphine.fr Abstract In this article, we review the theory of Gaussian multiplicative chaos initially introduced by Kahane s seminal work in 985. Though this beautiful paper faded from memory until recently, it already contains ideas and results that are nowadays under active investigation, like the construction of the Liouville measure in 2d-Liouville quantum gravity or thick points of the Gaussian Free Field. Also, we mention important extensions and generalizations of this theory that have emerged ever since and discuss a whole family of applications, ranging from finance, through the Kolmogorov-Obukhov model of turbulence to 2d-Liouville quantum gravity. This review also includes new results like the convergence of discretized Liouville measures on isoradial graphs (thus including the triangle and square lattices) towards the continuous Liouville measures (in the subcritical and critical case) or multifractal analysis of the measures in all dimensions. Key words or phrases: Random measures, Gaussian processes, Multifractal processes. MSC 2000 subject classifications: 60G57, 60G5, 60G25, 28A80 Contents Introduction 2. A word on Quantum field theory and the Hoeg-Krohn model Notations State of the art since Kahane 6 2. The seminal work of Kahane in Examples of kernels of σ-positive type Exact kernels

2 Star scale invariant kernels Green functions Different notions of stochastic scale invariance Extensions of the theory 9 3. Limitations of Kahane s theory Generalized Gaussian multiplicative chaos Extension to open domains and non-stationary smooth fields Duplantier-Sheffield s approach in dimension Multifractal analysis of the measures The Peyrière probability measure and the thick points General multifractal formalism: a heuristic introduction Applications of Gaussian multiplicative chaos Volatility of a financial asset or boundary Liouville measure MRM with infinite correlation length or boundary Liouville measure Liouville Quantum Gravity and KPZ Convergence of discrete Liouville measures on isoradial graphs Kolmogorov-Obhukov model in turbulence Decaying Burgers turbulence Generalizations of the theory Gaussian multiplicative chaos at criticality Atomic Gaussian multiplicative chaos Conjectures in connection with lognormal star scale invariance About the maximum of log-correlated Gaussian fields and discrete GFF Matrix-valued Gaussian multiplicative chaos Beyond Gaussianity A Proof of lemma B Proof of convergence of Theorem C Multifractal formalism 6 C. Proof of Theorem Introduction Log-normal multiplicative martingales were introduced by Mandelbrot [] in order to build random measures describing energy dissipation and contribute explaining intermittency effects in Kolmogorov s theory of fully developed turbulence (see [34, 33, 37, 35, 65] and references therein). However, his model was difficult to define mathematically and this is why he proposed in [2] the simpler model of random multiplicative cascades whose detailed study started with Kahane s and Peyrière s notes [86, 22], improved and gathered in their joint paper [87]. 2

3 From that moment on, multiplicative cascades have been widely used as reference or toy models in many applications as they feature beautiful stochastic scaling relations, modeling a phenomenon that is commonly called intermittency. Let us roughly explain this point. Consider a multiplicative cascade M constructed on a dyadic tree. It is a random measure over the interval [0, ]. If you look at the measure at a dyadic scale 2 n, you observe the same object as M up to an independent stochastic factor: M[ k 2 n, k + 2 n ] law = e Ωn M[0, ], where Ω n is a random variable independent of M, the law of which depends on the scale 2 n. However, multiplicative cascades are constructed on a dyadic (or p-adic) tree and therefore possess many drawbacks: they do not possess stationary fluctuations and present discrete (p-adic) scaling relations. Gaussian multiplicative chaos, introduced by Kahane [85] in 985, is born from the need of making rigorous Mandelbrot s initial model of energy dissipation [], the so-called Kolmogorov-Obukhov model. It is about constructing a continuous parameter theory of suitable multifractal random measures. Kahane s efforts were followed by several authors [5, 2, 2, 63, 26, 27, 32, 33] coming up with various generalizations at different scales. These measures have found many applications in various fields of science, especially in mathematical finance (or equivalently boundary Liouville quantum gravity), 2d-Liouville quantum gravity and 3d-turbulence. In dimension d, a standard Gaussian multiplicative chaos is a random measure on a given domain D of R d that can be formally written, for any Borelian set A as: M γ (A) = A γ2 γx(x) e 2 E[X2 (x)] σ(dx) (.) where X is a centered Gaussian field and σ is a Radon measure on D. In the situations of interest, X is rather badly behaved and cannot be defined as a random function: it is a random distribution (in the sense of Schwartz), like the Gaussian Free Field (GFF for short, see [34] for an overview about the GFF) for instance. In his seminal work, Kahane focused on the case where X possesses a covariance kernel of the form: E[X(x)X(y)] = ln + + g(x, y), (.2) x y with ln + (u) = max(ln u, 0) and g a continuous bounded function over D D. Surprisingly, it turns out that this is the only situation of interest since this family of kernels can be thought of as a transition, separating the family of kernels for which (.) is trivially converging from the family of kernels for which (.) is trivially vanishing. The covariance kernel thus possesses a singularity and it is now clear that giving sense to (.) is not straightforward (how do you define the exponential of a distribution?). The standard approach consists in applying a cut-off to the distribution X, that is in regularizing the field X in order to get rid of the singularity of the covariance kernel and get a nicer field. The regularization usually depends on a small parameter that stands for the extent to which the field has been regularized. The measure (.) is naturally understood as the limit that you get when the regularization parameter goes to 0. Kahane s paper [85] is about making this sketch of construction rigorous. When σ is the Lebesgue measure for 3

4 instance, it turns out that it produces non trivial limiting objects when the parameter γ is less than some critical value γ c = 2d. These are the foundations of Kahane s theory. Several questions are then raised like: - Does the limiting measure M γ depend on the chosen cut-off procedure? - What are the geometrical and statistical properties of the measure M γ? - What are the regularity properties of the measure M γ? - How can we characterize the measure M γ? - What happens at γ c? and what about γ > γ c? In this review, we will discuss the above questions as well as possible generalizations and applications of the theory of Gaussian multiplicative chaos in the light of recent progresses. Among the main points that we will address are the convolution techniques used to produce a Gaussian multiplicative chaos. The situation may be summarized as follows: given a Gaussian distribution X with covariance kernel of the type (.2), what if we regularize X by convolution with a smoothing family of functions (θ ɛ ) ɛ, converging towards the Dirac mass at 0 when ɛ 0? If we formally define: X ɛ (x) = X(y)θ ɛ (x y) dy, it turns out that, under weak conditions, the family of random measures: γ2 γxɛ(x) M ɛ,γ (A) = e 2 E[X2 ɛ (x)] dx (.3) A weakly converges in law towards the same measure as that produced by Kahane s theory. Convolution techniques were first introduced in [3, 32], where convergence in law is established. In the particular case when X is a Gaussian Free Field (GFF), a particular convolution technique was also studied in [56], i.e. convolution of the field by a circle, where almost sure convergence along subsequences is established (see section 3). We will also discuss multifractality of Gaussian multiplicative chaos and related scaling relations. Roughly speaking, multifractal analysis is the study of objects, like measures or functions, possessing several levels of local regularity: for instance, the local Hölder exponent may vary spatially. In particular, we will explain why the (nowadays called) thick points of the GFF are very closely related to a general theory called multifractal analysis, which at least goes back to Kahane s paper when applied to Gaussian multiplicative chaos (see subsection 2.3). Nowadays, there is a huge amount of literature on multifractal analysis and it is far beyond the scope of this review to cite or discuss all the related mathematical achievements (in the case of Mandelbrot s multiplicative cascades see [4, 5, 22] and references therein). The applications that we will mention range from mathematical finance to fully developed turbulence through 2d-Liouville quantum gravity or decaying Burgers turbulence. On the one hand, some of them are rather well established so that we will just recall the basic framework and give references. On the other hand, we feel important to devote a considerable part of the paper to 2d-Liouville quantum gravity. Roughly, it can be seen as an attempt to construct a canonical Riemannian random metric on the sphere. Since Polyakov s original work [25], physicists have understood that such a metric takes on the form of the exponential of a GFF [47, 25]. Interestingly, they have also understood that such metrics can be discretized by randomly triangulated surfaces (see [6] for instance), a 4

5 recent field of mathematical research that culminated with the construction of the so-called Brownian map (see [05, 06, 07, 4]) in the special case of pure gravity. While constructing properly a metric is presently out of reach, Gaussian multiplicative chaos theory gives straightforwardly a way of defining the associated volume form, called the Liouville measure. Non specialists of Gaussian multiplicative chaos will find in subsection 5.2 the different constructions suggested in the literature, i.e. white noise decomposition/circle average/h -expansion of the GFF, as well as a proof that they all produce the same measure in law. Let us mention here the remarkable works of Duplantier-Sheffield [56] and Sheffield [35] where the authors state an impressive series of conjectures, which can be seen as a starting point to understand the physicist picture (see also the nice review for mathematicians [70]). In particular, this ambitious program could give a rigorous geometrical framework to the Knizhnik-Polyakov-Zamolodchikov formula (KPZ for short, see [94]). For instance, such a KPZ formula has already been used by physicists [54] to predict the exact values of the Brownian intersection exponents and has been rigorously checked to hold in some special cases (see [26]). We will explain the geometrical KPZ formulae rigorously proved in [56, 28] (see also [27, 7]). Furthermore, as pointed out in [28], the theory of Gaussian multiplicative chaos allows us to deal with much more general situations concerning constructions of measures or the KPZ formula. We stress here that the theory of Gaussian multiplicative chaos in [85] and the KPZ relation proved in [28] (or [7]) is valid in any dimension when applied to logcorrelated Gaussian fields. In particular, boundary Liouville measures are discussed since they are nothing but Gaussian multiplicative chaos along d Riemannian manifolds. Also, in dimension d one can consider situations where the field X has correlations given by the kernel (m 2 ) d/2 (with possibly m = 0) since such correlations are logarithmic (see [59]). We do not detail this situation here, first because we cannot explain in great details all the situations where the theory of Gaussian multiplicative chaos applies and second, because it could be instructive for the reader to check that the framework drawn in [85] for the construction of measures or in [28, 7] for the KPZ formula applies (the reader may also consult [38] on this topic). Finally, as a new result, we explain how to combine the results in [37, 32] to prove that the discrete Liouville measures on isoradial graphs converge towards the Liouville measure as the mesh of the graph converges to 0. The last part of this review (section 6) will be devoted to possible generalizations of Kahane s theory. We will discuss how to renormalize the vanishing measure (.) for γ 2 2d (when σ is the Lebesgue measure). This yields new qualitative behaviours of the limiting measure that may be classified in two categories: Critical Gaussian Multiplicative Chaos or Atomic Gaussian Multiplicative Chaos. We discuss the associated relations with duality in 2d-Liouville quantum gravity, the frozen phase of logarithmically correlated Gaussian potentials or the maximum of log-correlated Gaussian fields (including in particular the maximum of the GFF). Other possible generalizations are mentioned, like taking complex-valued fields X or matrix-valued fields X in (.). Finally, we mention that a preceeding review on multiplicative chaos has already appeared [62]. The reader may find in [62] some further fields of applications including Dvoretzy covering, percolation on trees, random cascades and Riesz products that we do not review here for the sake of non-overlapping. Concerning multiplicative cascades, the reader may consult the recent review [6]. 5

6 . A word on Quantum field theory and the Hoeg-Krohn model To our knowledge, the first mathematical occurence of measures of the form (.) appeared in Hoeg-Krohn s work [8]. In the context of Quantum field theory, Hoeg-Krohn focused on the case where X is the two dimensional massive free field and σ the Lebesgue measure (to be precise, their point of view is not exactly that of random measures). He showed that the measure M γ is non trivial for γ 2 < 2, thus working below the L 2 -threshold; under the assumption γ 2 < 2, one can perform L 2 -computations which considerably simplifies the study of Gaussian multiplicative chaos measures (see subsection 2.). This work led to other works [2, 3] in dimension 2 which generalized some of the initial results of [8]. Nonetheless, it seems that none of these works focused on building a general theory applicable to a wide range of random measures in all dimensions..2 Notations The truncated logarithm ln + is the function ln + (x) = max(0, ln x). The relation f g means that there exists a positive constant c > 0 such that c f(x) g(x) cf(x) for all x under consideration. We denote by (D, ρ) a metric space D equipped with its metric ρ. This metric space is endowed with its Borelian sigma algebra B(D). Given a domain D of R d, we denote by H0 (D) the classical Sobolev space defined as the Hilbert space closure with respect to the Dirichlet inner product of the set of smooth compactly supported functions on D. Finally, γ will denote the intermittency parameter. We will suppose that γ 0. Acknowledgements We would like to thank M. Bauer, F. David, B. Eynard, J.F. Le Gall for interesting discussions on these topics. Special thanks are addressed to all those who read and suggested improvements of a prior draft of this manuscript: J. Barral, N. Curien and C. Garban. 2 State of the art since Kahane 2. The seminal work of Kahane in 985 The theory of multiplicative chaos was first defined rigorously by Kahane in 985 in the article [85] to which the reader is referred for further details or definitions. More specifically, Kahane built a theory relying on the notion of σ-positive type kernel. Consider a locally compact metric space (D, ρ). A function K : D D R + { } is of σ-positive type if there exists a sequence (K k ) k of continuous nonnegative and positive definite kernels K k : D D R + such that: x, y D, K(x, y) = k K k (x, y). (2.) It is worth pointing out here that Kahane s theory uses nonnegativeness of the kernels (K k ) k for the only sake of an easy formulation of a uniqueness criterion (see Theorem 2.3 below). If the reader is not interested in the uniqueness part of Gaussian multiplicative 6

7 chaos theory, he may skip this assumption of nonnegativeness as we will discuss in section 3 a more elaborate uniqueness criterion. If K is a kernel of σ-positive type with decomposition (2.), one can consider a sequence of independent centered Gaussian processes (Y k ) k with covariance kernels (K k ) k. Then the Gaussian process X n = n k= has covariance kernel n k= K k. Given a Radon measure σ on D, it is proved in [85] that the sequence of random measures (M n ) n given by: A B(D), M n,γ (A) = A Y k γ2 γxn(x) e 2 E[Xn(x)2] σ(dx) (2.2) converges almost surely in the space of Radon measures (equipped with the topology of weak convergence) towards a random measure M, which is called Gaussian multiplicative chaos with kernel K acting on σ. Basically, this convergence relies on the fact that for each compact set A, the sequence (M n,γ (A)) n is a nonnegative martingale. This martingale structure ensuring the almost sure convergence of (2.2) at low cost is the main motivation for considering kernels of σ-positive type. Then Kahane established a whole set of properties of this chaos that he derived from the following comparison principle: Theorem 2.. Convexity inequalities. [Kahane, 985]. Let (A i ) i n and (B i ) i n be two centered Gaussian vectors such that: i, j, E[A i A j ] E[B i B j ]. Then for all combinations of nonnegative weights (p i ) i n and all convex (resp. concave) functions F : R + R with at most polynomial growth at infinity E [F ( n p i e A i 2 E[A2])] i (resp. ) E [F ( n p i e B i 2 E[B2])] i. (2.3) i= Remark 2.2. Note that theorem 2. is very general and can be useful for instance in the study of random Gibbs measures where the Hamiltonian is a Gaussian variable. This is for instance the case in the Sherrington-Kirkpatrick (SK) model of statistical physics. In an important work, the authors of [80] rederived (2.3) with F (x) = ln x in the context of the SK model and used this inequality to compare the model of size N with two independent subsystems of size N and N 2 with N = N + N 2. As a consequence, they noticed that the expected finite size free energy was subadditive with respect to its size hence obtaining the existence of the limiting free energy as the size of the system goes to infinity. Private communication with J.. Kahane: The terminology multiplicative chaos was adopted since this theory may be seen as a multiplicative counterpart of the additive Wiener chaos theory. Actually, this is Paul Lévy himself who suggested to J.P. Kahane in the seventies to construct a multiplicative theory of random variables, arguing that this should be as fundamental as the additive theory of random variables. It took almost ten years to Kahane to build his theory of Gaussian multiplicative chaos. i= 7

8 This ingenious inequality sheds some light on the mechanism of Gaussian multiplicative chaos. For instance, let us stress that a kernel K of σ-positive type admits infinitely many decompositions of the form (2.): you can obtain other decompositions by changing the order of the kernels K k, by gathering them, etc... so there are possibly quite different kernels (K k ) k whose sum is K. The important question thus is: does the law of the limiting measure M depend on the choice of the decomposition (K k ) k in (2.)? Theorem 2.3. Uniqueness. [Kahane, 985]. The law of the limiting measure M γ does not depend on the sequence of nonnegative and positive definite kernels (K k ) k used in the decomposition (2.) of K. Thus, the theory enables to give a unique and mathematically rigorous definition to a random measure M in D defined formally by: γ2 γx(x) A B(D), M γ (A) = e 2 E[X(x)2] σ(dx). (2.4) where (X(x)) x D is a centered Gaussian field whose covariance K is a σ-positive type kernel. To show the usefulness of Theorem 2., it is worth giving a few words about the proof of Theorem 2.3. It works roughly as follows. Assume that you have two decompositions (K k ) k and (K k ) k of K with associated Gaussian process sequences (X n ) n and (X n) n and associated measures (M n ) n and (M n) n. Both sequences ( n k= K k) n and ( n k= K k ) n converge pointwise towards K in a nondecreasing way. Therefore, if we choose a compact set T D then, for each fixed p and ɛ > 0, the Dini theorem entails that A p K k ɛ + k= for q large enough on T T. Since q k= K k (resp. p k= K k) is the covariance kernel of X q (resp. X p ), we can apply Kahane s convexity inequalities and get, for each bounded convex function F : R + R: q k= K k E[F (M p (A))] E[F (e ɛγz ɛγ2 2 M q (A))], where Z is a standard Gaussian random variable independent of M. By taking the limit as q tends to, and then p, we obtain E[F (M(A))] E[F (e ɛγz ɛγ2 2 M (A))]. Since ɛ > 0 can be chosen arbitrarily small, we deduce E[F (M(A))] E[F (M (A))]. The converse inequality is proved in the same way, showing E[F (M(A))] = E[F (M (A))] for each bounded convex function F. By choosing F (x) = e λx for λ > 0, we deduce that the measures M and M have the same law. However, the simplicity of the convergence does not solve the question of non-degeneracy of the limiting measure M: it is possible that M identically vanishes. A 0 law argument 8

9 straightforwardly shows that the event M is identically null has probability 0 or. It seems difficult to state a general decision rule to decide whether M is degenerate or not. It depends in an intricate way on the covariance structure, i.e. the kernel K, and on the measure σ. So Kahane focused on the situation when the kernel K and the measure σ are intertwined via the metric structure of D. More precisely, he assumed that K can be written as T x, y D, K(x, y) = ln + + g(x, y) (2.5) ρ(x, y) where T > 0, g : D D R is a bounded continuous function and σ is in the class R α + (denoted M α + in Kahane s paper): Definition 2.4. For α > 0, a Borel measure σ is said to be in the class R + α if for all ɛ > 0 there is δ > 0, C < and a compact set A ɛ D such that σ(d \ A ɛ ) ɛ and: O open set, σ(o A ɛ ) C ɛ diam ρ (O) α+δ, (2.6) where diam ρ (O) is the diameter of O with respect to ρ. For instance, the Lebesgue measure of R d, restricted to any bounded domain of R d equipped with the Euclidean distance ρ, is in the class R α + for all α < d. The above definition looks like a Hölder condition for measures. It is intimately related to the notion of measure with finite β-energy: a Borel measure σ is said to be of finite β-energy if I β (σ) = σ(dx)σ(dy) < +. (2.7) ρ(x, y) β Indeed, if σ has a finite β-energy then σ R α + for all α < β. Conversely, if σ is in the class R α +, then the measure σ Aɛ (dx) = Aɛ (x)σ(dx) has finite β-energy for all β < α + ɛ. To have a flavor of the forthcoming results, let us treat the following simple situation, which we call the below L 2 -threshold case. It is about formulating a criterion ensuring that the martingale (M n (A)) n is bounded in L 2 for some given bounded set A, and therefore uniformly integrable. A straightforward computation shows that: E[M n,γ (A) 2 ] = A A e γ2k(x,y) σ(dx)σ(dy) A A C σ(dx)σ(dy). ρ(x, y) γ2 A γ2 γxn(x) E[e 2 E[Xn(x)2] γ2 γxn(y) e 2 E[Xn(y)2] ] σ(dx)σ(dy) A Therefore, if the measure A (x)σ(dx) has finite γ 2 -energy, the martingale (M n (A)) n is bounded in L 2 and therefore converges towards a non trivial limit. In the case when σ is the Lebesgue measure of R d and ρ the Euclidean distance, this condition simply reads γ 2 < d. Kahane proved the following highly deeper result: Theorem 2.5. Non-degeneracy. [Kahane, 985]. Assume that the kernel K takes on the form (2.5) and that the measure σ is in the class R + α for some α > 0. Then, for each compact set A, the sequence (M n,γ (A)) n is a uniformly integrable martingale. Hence γ 2 < 2α M γ is non degenerate. 9

10 As a by-product of his proof, Kahane also shows the following result concerning the dimension of the carrier of the measure M γ : Theorem 2.6. Structure of the carrier. [Kahane, 985]. Assume that the kernel K takes on the form (2.5) and that the measure σ is in the class R α + for some α > 0. If γ 2 < 2α, the measure M γ is non degenerate and is almost surely in the class R +. α γ2 2 Note that Theorem 2.6 implies that the measure M cannot give positive mass to a set of Hausdorff dimension less or equal than α γ2 2. Therefore the measure M cannot possess atoms if σ is in some class R α +. We will see in subsection 4. that the Hausdorff dimension of the carrier is exactly d γ2 2 when σ is the Lebesgue measure on a domain of Rd. Remark 2.7. Stronger versions of Theorem 2.5 and Theorem 2.6 were proved recently in [29]. Remark 2.8. If σ is the Lebesgue measure on some domain D R d equipped with the Euclidean metric, then σ is in the class R + d ɛ for all ɛ > 0 hence, if γ2 < 2d, the measure M γ is non degenerate and in the class R + for all ɛ > 0. d γ2 ɛ 2 Partial converses of Theorem 2.5 are more intricate. Kahane first gave a general necessary condition: Theorem 2.9. Necessary condition of non-degeneracy. [Kahane, 985]. Assume (D, ρ) is a locally compact metric space and: -the function (t, s) ρ(t, s) 2 is of negative type, -σ has the doubling property, namely that there exists a constant C such that x T, r > 0, σ(b(x, 2r)) Cσ(B(x, r)), -the kernel K takes on the form (2.5). Denote by dim(d) the Hausdorff dimension of D. If γ 2 > 2dim(D) then M γ is degenerate. As pointed out by Kahane, for a squared distance of negative type, the assumptions of the above theorem are satisfied when the triple (D, ρ, σ) admits a Lipschitz immersion into a finite-dimensional space. Let us also stress that the critical situation γ 2 = 2dim(D) is not settled by this theorem. Nevertheless, he reinforced his assumptions to prove: Theorem 2.0. Necessary and sufficient condition of non-degeneracy. [Kahane, 985]. Assume (D, ρ) is a d-dimensional manifold of class C and let σ be its volume form (or any Radon measure absolutely continuous w.r.t the volume form with a bounded density). Assume that the kernel K takes on the form (2.5). Then More results in the Euclidean space M γ is non degenerate γ 2 < 2d. When the metric space (D, ρ) is an open subset of R d for some d equipped with the Euclidian distance and σ is the Lebesgue measure, the previous results can be strengthened. We first point out that Theorem 2.0 applies and the non-degeneracy necessary and sufficient condition reads γ 2 < 2d. Kahane s convexity inequalities (Theorem 2.) allow us to give a complete description of the moments of M γ : 0

11 (a) γ = 0.2 (b) γ = (c) γ =.8 Figure : Influence of the intermittency parameter γ: we plot the density profile of a 2-dimensional multiplicative chaos illustrating the clustering of the activity of the chaos when γ grows. Theorem 2.. Positive Moments. [Kahane, 985]. If the measure M γ is non degenerate, that is γ 2 < 2d, the measure M γ admits finite positive moments of order p for all p ]0, 2d [. More precisely, for all compact set A D and p ]0, 2d [, we have γ 2 γ 2 E[M γ (A) p ] < +. Basically, it suffices to prove the finiteness of the moments for your favorite kernel K of the type (2.5) and deduce that the conclusions remain valid for all the kernels of the type (2.5) via Theorem 2.. We know turn to the existence of negative moments which was not investigated by Kahane. We have: Theorem 2.2. Negative Moments. If the measure M γ is non degenerate, that is γ 2 < 2d, the measure M γ admits finite negative moments of order p for all p ], 0[. More precisely, for any compact nonempty Euclidean ball A D and p ], 0[, we have E[M γ (A) p ] < +.

12 The finiteness of moments of negative order is proved in [5] in the case of discrete cascades. By adapting the argument of [5] and using the convexity inequalities 2.3, theorem 2.2 is proved in [32]. Let us also point out an important result in [8] where the authors compute the tail distributions of the measure M γ in dimension with the kernel K(x, y) = ln + x y (which is of σ-positive type, see Proposition 2.5 below): Theorem 2.3. Distribution tails. [Barral, Jin, 202]. segment of R then there exists a constant c > 0 such that If A is some nonempty lim x 2 γ 2 P(M γ (A) > x) = c. x + We conclude this theoretical background by pointing out some further interesting properties that can be exhibited when the integrating measure σ in (2.4) is the Lebesgue measure. The first observation that can be made is that the measure M γ is stationary in space as soon as we consider a stationary Gaussian distribution X. This is due to the translation invariance of the Lebesgue measure and the stationarity of X. Furthermore, when M is non-degenerate, it can be shown that the support of M γ is almost surely the whole of R d. This results from the 0 law of Kolmogorov: if you consider a ball B, the 0 law tells you that the event {M γ (B) > 0} has probability 0 or. Indeed, we have inf γ 2 x B eγxn(x) 2 E[X2 n (x)] M γ,n (B) M γ (B) sup x B where M γ,n (B) is the Gaussian multiplicative chaos Since for k > n: M γ,n (dx) = lim k eγ(x k X n)(x) γ 2 2 E[(X k X n) 2 (x)] dx. X k X n = k p=n+ γ2 γxn(x) e 2 E[X2 n (x)] M γ,n (B), (2.8) we may say that we have just removed the dependency on the first n fields Y,..., Y n. Therefore (2.8) entails that for any n 0 Y p, {M γ (B) > 0} = {M γ,n (B) > 0}. Since the event {M γ,n (B) > 0} is independent of the fields Y,..., Y n, we deduce that the event {M γ (B) > 0} belongs to the asymptotic sigma-algebra generated by the fields (Y n ) n in such a way that it has probability 0 or. When M γ is non-degenerate, it has clearly probability. This argument can be reproduced for every ball chosen among a countable family of balls (B n ) n generating the open sets of R d. Therefore, the event n {M γ(b n ) > 0} has probability, proving that almost surely the support of M γ is the whole of R d. 2.2 Examples of kernels of σ-positive type In this section, we give a few important examples of σ-positive kernels. 2

13 Exact kernels We consider for T > 0 the kernel x, y R d, K(x, y) = ln + T x y. (2.9) It is of σ-positive type in dimension d =, 2 and is involved in exact scaling relations as explained in subsection 2.3. Star scale invariant kernels A simple way of constructing σ-positive kernels on R d is to consider x, y R d k ( (x y)u ), K(x, y) = du, (2.0) u where k is a continuous function of positive type with k(0) =. Such kernels are related to the notion of -scale invariance (see subsection 2.3). Whole plane massive Green functions are -scale invariant. Green functions If we consider a bounded domain D of R 2, the Green function G of the Laplacian with 0-boundary condition is of σ-positive type. The corresponding Gaussian distribution X with covariance G is the Gaussian Free Field (GFF for short). The associated Gaussian multiplicative chaos is called the Liouville measure. More details about these claims are given in subsection Different notions of stochastic scale invariance In this subsection, we consider the Euclidian framework. More precisely, we consider an open set D of R d and Gaussian multiplicative chaos of the type γ2 γx(x) A B(D), M γ (A) = e 2 E[X(x)2] dx (2.) where the Gaussian distribution X has covariance kernel of the form: T K(x, y) = ln + + g(x, y) (2.2) x y for some continuous and bounded function g over D 2. Then the power-law spectrum of such Gaussian multiplicative chaos presents some interesting features, such as non-linearity (in the parameter q in the following theorem): Theorem 2.4. Assume that the kernel K takes on the form (2.5) with γ 2 < 2d. Choose a point x D. For each q [0, 2d γ 2 [ we have E [ M γ (B(x, r)) q] r ξ(q) as r 0 where ξ is the structure exponent of the measure M: q [0, 2d γ 2 [, A ξ(q) = ( d + γ2 ) γ 2 q 2 2 q2. 3

14 Heuristic proof. For simplicity, assume that T = in (2.2). By making a change of variables, we get: γ2 γx(y) M γ (B(x, r)) = e 2 E[X(y)2] dy B(x,r) =r d γ2 γx(x+ry) e 2 E[X(x+ry)2] dy. B(0,) Then we observe that the field (X(x + ry)) y has a covariance structure approximatively given for r by: E[X(x + ry)x(x + rz)] ln + r y z = ln r + ln + y z ln r + E[X(x + y)x(x + z)]. The above relation gives us the following (good) approximation in law (X(x + ry)) y B(0,) Ω r + (X(x + y)) y B(0,) where Ω r is a centered Gaussian random variable independent of the field (X(x+y)) y B(0,) and with variance ln r. Therefore M γ (B(x, r)) =r d B(0,) r d γ2 γωr e 2 E[Ω2 r] γ2 γx(x+ry) e 2 E[X(x+ry)2] dy B(0,) =r d γ2 γωr e 2 E[Ω2 r ] M γ (B(x, )). By taking the q-th power and integrating, we get as r 0 γ2 γx(x+y) e 2 E[X(x+y)2] dy E [ M γ (B(x, r)) q] r ξ(q) E[M γ (B(x, ) q ], thus explaining the theorem. Actually, the rigorous proof of this result is very close to the heuristic developed here. Notice that the quadratic structure of the structure exponent is intimately related to the Gaussian nature of the random distribution X. Random measures with a non-linear power-law spectrum are often called multifractal. That is why Gaussian multiplicative chaos (and other possible extensions) are sometimes called Multifractal Random Measures (MRM for short) in the literature. It is also natural to wonder if some specific choice of the covariance kernel K may lead to replacing the symbol in Theorem 2.4 by the symbol =. It turns out that this question is related to some specific scaling relation, which we describe below. 4

15 Exact stochastic scale invariance The notion of exact stochastic scale invariance relies on the additive properties of the logarithm function. Roughly speaking, in order to produce kernels with exact scaling relations, kernels K of the type (2.5) with g = 0 (and ρ the Euclidian distance) must be considered. So we first focus on the σ-positive type of such kernels: Proposition 2.5. For d 2 and T > 0, the function is of σ-positive type. x R d λ 2 ln + T x (2.3) Proof. A straightforward computation yields: ln + T x = + 0 (t x ) + ν T (dt) where ν T is the measure (δ T is the Dirac mass at T ): Hence for any µ > 0, we have: ν T (dt) = [0,T ] (t) dt t 2 + T δ T (dt). T ln + x = µ ln T µ + x µ = + 0 (t x µ ) + ν T µ(dt). By using a Chasles relation in the integral of the right-hand side, proving that this kernel is of σ-positive type thus boils down to considering the possible values of µ > 0 such that the function ( x µ ) + is of positive type: this is the Kuttner-Golubov problem (see [77]). For d =, it is straightforward to see that ( x ) + is of positive type (compute the inverse Fourier transform). In dimension 2, Pasenchenko [2] proved that the function ( x /2 ) + is of positive type on R 2. We can thus write with K n (x) = T ln + x = K n (x) n n n (t x µ ) + ν T µ(dt) with µ = in dimension and µ = /2 in dimension 2. Theorem 2.6. Exact stochastic scale invariance. [Bacry, Muzy 2003]. Let K be the covariance kernel given by (2.3) in dimension d = or d = 2. The associated Gaussian multiplicative chaos M γ is exactly stochastic scale invariant: λ ]0, ], (M(λA)) A B(0,T/2) law = λ d e Ω λ 2 E[Ω2 λ ] (M(A)) A B(0,T/2), (2.4) where Ω λ is a Gaussian random variable, independent of the measure (M(A)) A B(0,T/2), with mean 0 and variance γ 2 ln λ. 5

16 We stress that the above equality in law is to be understood in the sense ( M(λA ),..., M(λA p ) ) law = ( λ d e Ω λ 2 E[Ω2 λ ] M(A ),..., λ d e Ω λ 2 E[Ω2 λ ] M(A p ) ) for all possible choice A,..., A p of Borelian subsets of the ball B(0, T/2). Notice that such a scaling property makes obvious several computations related to the measure M γ. The reader may, for instance, observe that the heuristic proof of Theorem 2.4 becomes rigorous for such a measure. It is also not difficult to see that this scaling relation is bound to be valid only locally (over a ball) and cannot hold on the whole space: the logarithm is not of positive type over the whole of R d. Exact scaling relations were introduced in [2] in dimension together with further generalizations in the case of log-infinitely divisible random measures. It is natural to wonder how to construct such measures in dimension higher than 3. The procedure is somewhat complicated by the following observations: for d = 3, it is an open question to know whether the kernel (2.3) is of σ-positive type and for d 4, it is even not of positive type. Another approach has therefore been suggested in [27]. The main ideas are the following: what only matters to construct exactly stochastic scale invariant measures is that the covariance kernel must be the logarithm function over a ball centered at 0. The way of truncating the logarithm does not matter in order to obtain the scaling relation (2.4) but may be sensitive to the dimension when regarding positive definiteness: truncating the logarithm with the function ln + does not resist increasing the dimension. In [27], another truncation is suggested: we can find an isotropic function g : R d R that is constant on a neighborhood of 0 and such that the kernel T K(x) = ln + + g(x) (2.5) x is of σ-positive type. Briefly, the construction is the following: let us denote by S the sphere of R d and by σ the unique uniform measure on the sphere such that σ(s) =. This (probability) measure is invariant under rotations. Let us define the function K(x) = S ln + T σ(ds), (2.6) x, s where, stands for the canonical inner product of R d. Since σ is invariant under rotations, the function K is isotropic. Fix x R d such that x T and write x = x e where e S. Then we have K(x) = S ln T x e, s σ(ds) = γ2 ln T x + γ2 ln S e, s σ(ds). By invariance under rotations of σ, the second term in the right-hand side does not depend on x and turns out to be finite: this can be seen by noticing that, under σ, the random variable e, s has the law of the first entry of a Haar vector. K thus coincides with the logarithm over a neighborhood of 0, up to an additive constant. Gaussian multiplicative chaos with associated kernel (2.6) are exactly stochastic scale invariant in the sense of (2.4). It remains an open question to know to which extent the notion of exact stochastic scale invariance uniquely determines the covariance structure of the associated kernel K. 6

17 Star scale invariance As explained above, the notion of exact stochastic scale invariance is a local notion (valid only over a ball). We now present a global notion of stochastic scale invariance, called star scale invariance in [5] in reference of earlier works by Mandelbrot in the case of multiplicative cascades on trees. It stems from the need of characterizing Gaussian multiplicative chaos with functional equations. Indeed, tractable functional equations may provide efficient tools in identifying Gaussian multiplicative chaos as, for instance, scaling limits of discrete models. Definition 2.7. Log-normal -scale invariance. A random Radon measure M on R d is said lognormal -scale invariant if for all 0 < ε, M obeys the cascading rule ( M(A) ) A B(R d ) law = ( A e ωε(x) M ε (dx) ) A B(R d ) (2.7) where ω ε is a Gaussian process, which is assumed to be stationary and continuous in probability, and M ε is a random measure independent from X ε satisfying the scaling relation ( M ε (A) ) A B(R d ) law = ( M( A ε )) A B(R d ). (2.8) Notice that the process ω ε is unknown. Roughly speaking, we look for random measures that scale with an independent lognormal factor on the whole space. This property is shared by a large class of Gaussian multiplicative chaos. And for those Gaussian multiplicative chaos that do not share this property, they are very close to satisfying it. If the reader is familiar with branching random walks (BRW), here is an explanation that may help intuition. If we consider a BRW the reproduction law of which does not change with time (i.e. is the same at each generation), the law of the branching random walk will be characterized by a discrete version of the above -scale invariance called fixed point of the smoothing transform (in the lognormal case of course, see [60, 29, 08]). If the reproduction law evolves in time, then we have to change things a bit to adapt to this time evolution. The same argument holds for the log-normal -scale invariance: it characterizes these Gaussian multiplicative chaos that do not vary along scales. It is proved in [5] that E[e ωε(r) ] = ε d as soon as the measure M possesses a moment of order + δ for some δ > 0. Furthermore, up to weak regularity conditions on the covariance kernel of the process ω ε summarized in the definition below, all the log-normal -scale invariant random measures with enough moments can be identified. Definition 2.8. We will say that a stationary random measure M satisfies the good lognormal -scale invariance if M is lognormal -scale invariant and for each ɛ <, the covariance kernel K ɛ of the process ω ɛ involved in (2.7) is continuous and satisfies: k ɛ (r) 0 as r +, (2.9) r, r R d \ {0}, k ɛ (r) k ɛ (r ) C ɛ θ ( min( r, r ) ) r r (2.20) for some positive constant C ɛ and some decreasing function θ :]0, + [ R + such that + θ(u) ln(u) du < +. (2.2) 7

18 Lognormal -scale invariant random measures are then characterized as: Theorem 2.9. [Allez, Rhodes, Vargas, 20] Let M be a good lognormal -scale invariant stationary random measure. Assume that E[M([0, ] d ) +δ ] < + for some δ > 0. Then M is the product of a nonnegative random variable Y L +δ and an independent Gaussian multiplicative chaos: A B(R d γ2 γx(x) ), M(A) = Y e 2 E[X(x)2] dx (2.22) with associated covariance kernel given by the improper integral A x R d \ {0}, K(x) = + k(xu) u du (2.23) for some continuous covariance function k such that k(0) = and γ 2 2d +δ. Conversely, given some datas k and Y as above, the relation (2.22) defines a lognormal -scale invariant random measure M with finite moments of order +β for every β [0, δ). It is plain to see that the covariance structure (2.23) can be rewritten as x R d \ {0}, K(x) = ln + + g(x) (2.24) x for some continuous bounded function g, thus making the connection with Kahane s theory presented in subsection 2.. The first star-scale invariant kernel appearing in the literature goes back to Kahane s T original paper [85] in order to approximate the kernel ln + x y in dimension 3 but he did not study the related scaling relations. Kahane chose the kernel K(x) = + e u x Another example was exhibited in [2]. It corresponds to the choice K(x) = + k(ux) u u du. du with k(x) = ( x T ) [0,T ]( x ) and is based on a one dimensional geometric construction. Also, the authors in [2] made the connection with star scale invariance. To sum up, the solutions of the star-equation are now well identified in what is called the subcritical regime, which can be characterized by the fact that associated solutions possess enough moments (larger than ). It remains to investigate the characterization of solutions with only few moments (smaller than ): we will see in section 6 that new types of solutions are then involved, giving rise to quite new structures of the chaos. 8

19 3 Extensions of the theory 3. Limitations of Kahane s theory In view of natural applications, Kahane s theory appears unsufficient. Here are a few points that the theory does not address: A kernel K of σ-positive type is nonnegative and positive definite. Is the reciprocal true? What happens if one works with convolutions of K instead of nondecreasing approximating series? Kahane s theory is a theory which ensures equality in distribution. Can one build a theory that ensures almost sure equality? The first point raised above is still an open question. The two following points have been addressed recently as we will now describe in the next two subsections. 3.2 Generalized Gaussian multiplicative chaos Kahane s theory of Gaussian multiplicative chaos relies on the notion of kernels of σ- positive type. However, on the one hand it is not always straightforward to check such a criterion and on the other hand it seems natural to think that the theory should remain valid for kernels of positive type. A way of getting rid of σ-positive typeness has been developed in [3, 32]. The idea is to make a convolution product of the covariance kernel with a sequence of continuous functions approximating the Dirac delta function in order to smooth down the singularity of the kernel. More precisely, consider a positive definite function K in R d (or a bounded domain of R d ) such that K(x) = ln + T + g(x) (3.) x and g(x) is a bounded continuous function. Let θ : R d R be some continuous function with the following properties:. θ is positive definite, 2. θ has compact support and is α-hölder for some α > 0, 3. R d θ(x)dx =. Here is the main theorem of [32]: Theorem 3.. [Robert, Vargas, 2008] Let γ 2 < 2d. For all ɛ > 0, we consider the centered Gaussian field (X ɛ (x)) x R d defined by the convolution: E[X ɛ (x)x ɛ (y)] = (θ ɛ K)(y x), where θ ɛ = θ(. ɛ d ɛ ). Then the associated random measure A B(R d ), M ɛ,γ (A) = 9 A γ2 γxɛ(x) e 2 E[Xɛ(x)2] dx

20 converges in law in the space of Radon measures (equipped with the topology of weak convergence) as ɛ goes to 0 towards a random measure M γ, independent of the choice of the regularizing function θ with the properties., 2., 3. above. This theorem is useful to define a Gaussian multiplicative chaos associated to the kernel K(x) = ln + T x in dimension 3 and hence give a rigorous meaning to the Kolmogorov- Obukhov model. Indeed, remind that this kernel is of σ-positive type in dimension and 2. In dimension d = 3, the function ln + T x is positive definite (see [32]) but it is an open question whether it is of σ-positive type. Therefore, we are bound to apply Theorem 3. instead of Kahane s theory to define the associated chaos in dimension 3. In dimension greater than 4, the kernel ln + T x is no more positive definite. Let us also mention that it is proved in [32] that in dimension, 2, 3 the random variable M γ (B(0, r)) (for r < T ) possesses a C density with respect to the Lebesgue measure. Remark 3.2. Starting with a Gaussian field (X(x)) x R d whose covariance is given by (3.), one could also state theorem 3. in an equivalent way in terms of the fields X ɛ = θ ɛ X where θ is a function of average (but not necessarily of positive type) and which satisfies a decreasing condition at infinity. Extension to open domains and non-stationary smooth fields In what follows, we explain why the generalized Gaussian multiplicative chaos theory straightforwardly extends to the situation of bounded domains and possibly non stationary covariance kernels. Let D be an open set of R d. For all δ > 0, we set: D (δ) = {x D; dist(x, D) δ} By convention, if D = R d, we set D (δ) equal to the Euclidean ball of center 0 and radius. We consider a positive definite kernel K (non necessarily stationary) satisfying: δ x, y D, K(x, y) = ln + T + g(x, y) (3.2) x y where γ 2 < 2d and g is a bounded continuous function over D (δ) for all δ > 0. Let X be a random centered Gaussian distribution with covariance given by (3.2). We introduce the following notion of smooth approximation which will play a central role in the rest of the review: Definition 3.3. Smooth Gaussian approximations. We say a that a sequence of centered Gaussian fields (X ɛ ) ɛ>0 is a smooth Gaussian approximation of K if: for all x, y D, E[X ɛ (x)x ɛ (y)] converges to K(x, y) as ɛ goes to 0. for all δ > 0, there exists some constant C > 0 and α > 0 such that for all ɛ > 0: x, y D (δ), E[(X ɛ (x) X ɛ (y)) 2 ] C x y α ɛ α. Remark 3.4. In the above definition, the second point is a technical assumption. By standard results on Gaussian processes (see [04] for example), this assumption implies 20

21 the following useful property: for all δ > 0 and A > 0, there exists C > 0 such that for all x D (δ) : P( sup X ɛ (y) X ɛ (x) t) Ce C t2 2, t 0 y; y x Aɛ It is in fact this property implied by the second point which is crucial in the proofs of the results where smooth Gaussian approximations appear. For example, the convolution constructions of the previous section are smooth Gaussian approximations. It is not very difficult to see that convolutions of X with circles, balls or smooth bounded domains are smooth Gaussian approximations of K. In fact, the techniques of [32] can be straightforwardly adapted to give: Theorem 3.5. Assume that we are given two smooth Gaussian approximations (X ɛ ) ɛ>0 and ( X ɛ ) ɛ>0 of K such that: ) for some γ 2 γ2 γxɛ(x) < 2d, the random measure M ɛ,γ (dx) := e 2 E[Xɛ(x)2] dx converges almost surely (possibly along some subsequence) to some random Radon measure M γ in the sense of weak convergence of measures. 2) for all δ > 0 and A > 0, sup x,y D (δ), x y Aɛ E[X ɛ (x)x ɛ (y)] E[ X ɛ (x) X ɛ (y)] where C A > 0 is some constant independent from ɛ. 3) for all δ > 0, C A = lim sup ɛ 0 x,y D (δ), x y Aɛ C A. E[X ɛ (x)x ɛ (y)] E[ X ɛ (x) X ɛ (y)] goes to 0 as A goes to infinity. Under the above assumptions, the random measure M ɛ,γ (dx) = e γ X ɛ(x) γ2 2 E[ X ɛ(x) 2] dx converges in law as ɛ goes to 0 in the space of Radon measures on D (equipped with the topology of weak convergence) towards the random measure M γ. Let us mention here a simple and straightforward consequence of this theorem. Let X be a random centered Gaussian distribution with covariance given by (3.2). We consider a function θ of class C such that R d θ(x)dx = (3.3) and compactly supported in the ball B(0, ). Then we define the ɛ-approximation of X by: X ɛ (x) = X(x) θ ɛ (x y) dy. 2

22 Of course, the field X ɛ only makes sense for x D (ɛ). Then the law of the limiting measure M γ = lim ɛ 0 e γxɛ(x) 2 E[Xɛ(x)2] dx does not depend on the choice of the regularizing function θ. We will see in subsection 5.2 another application of this theorem in the case of the so-called Liouville measure. 3.3 Duplantier-Sheffield s approach in dimension 2 After the two aforementioned theories (Gaussian multiplicative chaos [85] and its generalized version [3, 32]), the authors of [56] came up with another contribution to Gaussian multiplicative chaos theory in the special case when the Gaussian distribution X in (2.4) is the GFF in a bounded domain and σ the Lebesgue measure. The situation can be roughly summarized as follows. Assume for instance that you are given a two-dimensional GFF X on a bounded domain D and you want to define the approximations (X n ) n in (2.2) almost surely as measurable functions of the whole distribution X. For instance, you may define X n as the projections of X along the first n vectors of an orthonormal basis of the Sobolev space H0 (D) and obtain a first way of defining almost surely the associated multiplicative chaos, called Liouville measure, as a measurable function of the whole GFF distribution. This falls under the scope of Kahane s theory since you are adding independent Gaussian fields. To be precise, the vectors of the basis are not necessarily nonnegative so that the construction of the measure works but you do not get the uniqueness property of Theorem 2.3 (except in some special cases: for instance, the Haar basis is composed of nonnegative functions and therefore in this case you are strictly within Kahane s framework). A second way of defining cutoff approximations of X that are measurable functions of the whole GFF distribution is to use convolution techniques: you may also define X ɛ (x) as the average value of X over the circle of radius ɛ centered at x (or any other regularizing function) and plug this quantity in place of X n in (2.2). The approximating measures γ2 γxɛ(x) M ɛ,γ (dx) = e 2 E[Xɛ(x)2] dx are then measurable functions of the GFF. Though convergence in law follows from the techniques of the generalized Gaussian multiplicative chaos theory [3, 32], it is important to focus on almost sure convergence of these measures. Theorem 3.6. [Duplantier, Sheffield, 2008]. The approximating measures (M ɛ,γ ) ɛ almost surely converge as ɛ 0 for the topology of weak convergence of measures along suitable deterministic subsequences. Another important point is the following. We have now at our disposal two ways to produce the Liouville measure as the almost sure limit of suitable approximations: H 0 (D) expansions or circle average. In this formulation, it is implicitly assumed that these two approximations yield the same limiting object. Is it really the case? Let us first stress that the two above constructions have the same law, and they also have the same law as the Liouville measure based on a white noise decomposition of the GFF introduced in [28] (see Theorem 5.5 below). There is another important contribution in [56] to the theory of Gaussian multiplicative chaos 22

23 Theorem 3.7. [Duplantier, Sheffield, 2008]. The Liouville measures constructed via an expansion along an orthonormal basis of H0 (D) or circle average approximations are almost surely the same measures. Though carried out in the case of GFF, this uniqueness result makes sense in a more general context: to which extent can we prove that Gaussian multiplicative chaos constructed with different approximations almost surely defined as measurable functions of the whole Gaussian distribution coincide almost surely? 4 Multifractal analysis of the measures The purpose of this section is to give a brief insight into multifractal analysis. More precisely, Kahane proved that the carrier of a Gaussian multiplicative chaos has Hausdorff dimension greater or equal to d γ2 2. It is natural to wonder whether further pieces of information can be given about the structure of this carrier. We will relate this question to the Parisi-Frisch formalism [20]. Consider a smooth cutoff approximation (X ɛ ) ɛ of the field X with covariance K given by (3.2) on a domain D and such that X ɛ X ɛ is independent from σ{x u ; u ɛ} for all ɛ < ɛ. We suppose that there exists a constant C > 0 such that ln ɛ C E[X ɛ(x) 2 ] ln ɛ + C for all x D. In this context, Kahane introduced the set of points { x D; limɛ X ɛ (x) ln ɛ = γ} and showed that it gives full measure to the chaos M γ, thus showing that it has Hausdorff dimension greater or equal to d γ2 2. This set of points has been called thick points of the GFF in [82] when the field X is the GFF and X ɛ corresponds to circle averages (hence, we are not exactly in the framework of Kahane as circle averages of the GFF do not correspond to adding independent fields; nonetheless, the two frameworks are very similar). We stick to this terminology as we feel that it is well-sounding. The Hausdorff dimension of thick points is derived in [82] in the context of circle averages of the GFF. Let us further mention that exact equality of Hausdorff dimension is well known in the closely related context of Mandelbrot s multiplicative cascades, see [4] and references therein. In the next subsection, we will generalize the results in [82, 85] to a large class of log-correlated Gaussian fields in all dimensions. After we will discuss general multifractal analysis. 4. The Peyrière probability measure and the thick points Let us first set the framework of this subsection. Consider a bounded domain D R d and a positive definite kernel K (non necessarily stationary) satisfying: x, y D, K(x, y) = ln + T + g(x, y) (4.) x y where g is a bounded continuous function over D. We work on a fixed probability space (Ω, F, P). On this space, we consider a centered Gaussian field (X ɛ (x)) ɛ>0,x D such that: 23

24 . (ɛ, x) X ɛ (x) is almost surely continuous on ]0, [ D, 2. (X ɛ ) ɛ>0 is a smooth Gaussian approximation of K (in the sense of Definition 3.3), 3. for each fixed x D, ɛ X ɛ (x) has independent increments, 4. There exists a constant C > 0 such that for all η < ɛ: ln + x y + ɛ C E[X ɛ(x)x η (y)] ln + + C, x, y D. (4.2) x y + ɛ 5. for all γ 2 < 2d, there exists a random Radon measure M γ such that M ɛ,γ (dx) = e γxɛ(x) 2 E[Xɛ(x)2] dx converges almost surely (possibly along some subsequence) to M γ on D. For γ 2 < 2d and q ]0, 2d γ [, we set: { ln M γ (B(x, ɛ)) K γ,q = x D; lim = d + ( } { } ɛ 0 ln ɛ 2 X ɛ (x) q)γ2 x D; lim 2 ɛ 0 ln ɛ = γq. We can now state the following theorems (see appendix): Theorem 4.. For all γ 2 < 2d, almost surely, the set K γ,q gives full mass to the measure M qγ, i.e. M qγ ( c K γ,q ) = 0. Theorem 4.2. For all γ 2 < 2d, almost surely, the set K γ, has Hausdorff dimension d γ2 2. A quite simple and concise proof of the above theorem is gathered in the appendix and, to our knowledge, is new in such generality (though much more is known in dimension : see [4, 24]). In dimension 2 and in the context of circle average approximations of the GFF, another paper [82] focused on the thick points of the GFF, improving Kahane s result by giving the upper bound for the Hausdorff dimension. Let us also stress here that the place of the almost sure is important in Theorem 4. and Corollary 4.2. In the case of multiplicative cascades, a stronger statement is proved in [5], where the almost sure is valid simultaneously for all γ 2 < 2d. In case the reader wishes to skip the whole proof of Theorem 4., we sketch here Kahane s argument about the lower bound for the Hausdorff dimension of the thick points of X. We focus on the case where q = and D has Lebesgue measure. The key point is to introduce the so-called Peyrière probability measure (also called rooted measure in [56, 57, 58]): Q(F (ω, x)) := E[ F (ω, x)m γ (dx)]. (4.3) D This measure was introduced by Peyrière in the context of discrete cascades (see [87]) and was used by Kahane in his work on Gaussian multiplicative chaos. It is obvious that under 24

25 this measure the process (x, t) X e t(x) is a Brownian motion with drift γ, hence we have M γ ( c Kγ )) = 0 where: { } X ɛ (x) K γ = x D; lim ɛ 0 ln ɛ = γ. (4.4) In particular, it is straightforward from Theorem 2.6 that K γ has Hausdorff dimension at least d γ2 2. Let us further stress that similar ideas are used in [82]. 4.2 General multifractal formalism: a heuristic introduction In this section, we discuss a bit of general multifractal formalism for the measures M γ. This discussion is essentialy based on heuristics though it should not be very difficult to make it rigorous mathematically. Let us mention that the study of multifractal formalism represents a very wide domain of mathematics and physics (in particular turbulence); therefore, being exhaustive in this field is way beyond the scope of this review where we only mention the case of measures of the form e X(x) dx for some log correlated gaussian field X. Multifractal analysis is based on the calculation of the L q -spectrum of the measure M γ, defined as q R τ M (q) = lim inf r 0 + { } log sup i M γ(b(x i, r)) q, log(r) where the supremum is taken over all the centered packing of [0, ] d by closed balls of radius r. This study is achieved in [4, 24] in dimension and, with a bit of effort, it should not be difficult to generalize this result to higher dimensions: ( d + γ 2 ) 2 q if q 2d τ Mγ (q) = ξ(q) d if q [ 2d γ ( d γ 2 ) 2 q if q γ, 2d γ., 2d γ ], The so-called multifractal formalism holds for M γ : we can relate the L q -spectrum of M γ to the local regularity of M γ. More precisely, define E δ = { x [0, ] d ; lim inf r 0 + ln M γ (B(x, r)) ln(r) = δ } (δ 0). The singularity spectrum of M γ, i.e. the mapping δ 0 dim E δ (dim meaning Hausdorff dimension), is given by the celebrated Parisi-Frisch formula (see [20]): δ 0 τ M γ (δ) = inf{δq τ Mγ (q) : q R} 0. (4.5) Note that the above formula is a Legendre transform since the formula comes in fact from a large deviation argument (which can be made rigorous in terms of box counting dimensions by using the Gartner-Ellis theorem). Multifractal analysis is essentially focused on studying the valididity of (4.5) in the broad context of all (random or deterministic) 25

26 measures. From this and the explicit expression on τ Mγ, one can deduce the following dimension result: dim E δ = { d 2 ( d γ + γ 2 δ γ )2, if δ [( d γ 2 ) 2, ( d + γ 2 ) 2 ]. 0 elsewhere. (4.6) In essence, multifractal formalism and computing the dimension of thick points is the same thing (and hence we get the same associated dimensions) if we admit the following commonly used heuristic: M γ (B(x, r)) C(x)r d e γxr(x) γ 2 2 ln r (4.7) r 0 where C(x) is some random constant of order (nearly independent of r but very dependent on x). Note that the notion of star scale invariance is nothing but a rigorous formulation of the heuristic (4.7). Swithching from multifractal formalism to thick point formalism rigorously implies handling the C(x) term which is of order but fluctuates wildely as a function of x. As an example of application of heuristic (4.7), let us recover formula (4.6) from corollary 4.2. By using (4.7), we get the following equivalence: M γ (B(x, r)) r 0 r δ X r (x) r 0 ( d γ + γ 2 δ γ ) ln r, hence leading to formula (4.6) thanks to corollary Applications of Gaussian multiplicative chaos In this section, we review some applications in direct relation with Kahane s theory: some of them are well known, some of them are new. Some further applications, rather related to recent generalizations of the theory, will be given in Section Volatility of a financial asset or boundary Liouville measure The main application of the theory is to give a meaning to the limit-lognormal model introduced by Mandelbrot in []. The limit-lognormal model corresponds to the choice of a stationary kernel K on R R given by: K(x, y) = ln + T + g(x, y) (5.) x y where T is a positive parameters, g is a bounded continuous function and σ is chosen to be the Lebesgue measure on R. This model has many applications: part of them are discussed below. If (X(t)) t 0 is the logarithm of the price of a financial asset, the volatility M γ of the asset on the interval [0, t] is by definition equal to the quadratic variation of X: M γ [0, t] = lim n k= n (X(tk/n) X(t(k )/n)) 2 26

27 The volatility M can be viewed as a random measure on R. The choice for M γ of multiplicative chaos associated to the kernel K(x, y) = ln + T x y satisfies many empirical properties measured on financial markets: lognormality of the volatility, long range correlations (see [4] for a study of the SP500 index and components and [42] for a general review). With this kernel, the measure is called the lognormal multifractal random measure (MRM) and is a particular case of the log-infinitely divisble multifractal random measures (see [2, 33] for the log-poisson case and [2] for the general case). Note that K is indeed of σ-positive type so M γ is well defined. In the context of finance, γ 2 is called the intermittency parameter in analogy with turbulence and T is the correlation length. Volatility modeling and forecasting is an important field of finance since it is related to option pricing and risk forecasting: we refer to [49] for the problem of forecasting volatility with this choice of M γ. Given the volatility M γ, the most natural way to construct a model for the (log) price X is to set: X(t) = B Mγ[0,t] (5.2) where (B t ) t 0 is a Brownian motion independent of M γ. Formula (5.2) defines the Multifractal Random Walk (MRW) first introduced in [0] (see [] for a recent review of financial applications of the MRW model). (a) SP500 returns (b) Returns simulated with Black-Scholes (c) Returns simulated with MRW (d) Returns simulated with MRW Figure 2: Intermittency in financial markets 27

28 In Figure 2, observe the burst of activity (intermittency) for the SP500 index or the MRW model. This reflects the Parisi-Frisch formalism (or multifractal formalism): the (strong) variations of regularity are related to the non-linearity of the power law spectrum. MRM with infinite correlation length or boundary Liouville measure Motivated by financial applications where the correlation length T is too big to be measured on markets, the authors of [49] adressed the issue of forecasting volatility M γ in the limit T. More precisely, they start by forecasting the log volatility, the /f noise which lives in the quotient space of distributions defined up to some additive constant. In this context, taking the exponential, the associated random measure M γ is then defined up to a multiplicative constant and is the boundary Liouville measure in the upper half plane considered in [56]. 5.2 Liouville Quantum Gravity and KPZ Let us first roughly explain the original motivations coming from the physics literature. We want to define a random distribution (g, X), the partition function of which formally writes Z = DgDX e S M (X,g) λv g(σ) (5.3) where S M is some conformally invariant action for matter fields coupled to a compact simply connected two dimensional surface Σ with metric g, λ is a constant (we do not discuss its value), V g is the volume form of g and X is an embedding from Σ into a c- dimensional spacetime. Here we adopt standard path integral notations: the above integral just means that we sum over all possible embeddings and metrics. Example 5.. For the free bosonic string, we consider the Polyakov action S M = g X g X dv g 8π where X specifies the embedding of Σ into flat D-dimensional space-time. Example 5.2. For the massive Ising model, we consider S M = X ( g X + m XX ) dv g. Notice that every metric g on Σ can be decomposed as f g = e ϕ g 0, where g 0 is a fixed metric on Σ, f is a g 0 -diffeomorphism and f g is the pullback metric of g along f. So we may perform the path integral by gauge-fixing: we choose a gauge defining an equivalence class over all the metrics and perform the above integral over a slice that cuts through once each gauge equivalence class. In view of the above factorization property, a natural choice of the gauge is the conformal gauge. We choose a family (ĝ) of representatives of each equivalence class of conformally equivalent metrics and we perform a sum over (ĝ) and over the equivalence class of ĝ for all ĝ. The Jacobian of such a change of variables is the so-called Faddeev-Popov determinant F P (ĝ). We do not detail this here but the reader is referred to [25, 24, 47] for further details and references. Let us just say that once this 28

29 determinant has been computed, it remains to make sure that the quantity resulting from these computations does not depend on the choice of the family of representatives (ĝ): in physics language, we have to compute the Weyl anomaly. Performing these computations lead to considering the Liouville action (R g is the Ricci tensor of the metric g) S L (ϕ, g) = 2πγ 2 ( g ϕ g ϕ + QR g ϕ + µe ϕ )dv g and the matter action S M (X, ĝ) in such a way that DXDg e S(X,g) = DXDfDϕ F P (ĝ) e S M (X,ĝ) S L (ϕ,ĝ). For c ], ], the value of γ [0, 2] is related to the central charge of the matter by γ = 25 c c 6. When the cosmological constant µ is set to 0, the Liouville action reduces to that of a free massless boson. Mathematically speaking the corresponding field ϕ is a Gaussian Free Field. In critical Liouville quantum gravity, we are therefore led to considering random metrics of the form e γϕ ĝ and an area measure e γϕ dvĝ, where ϕ is a Free Field in the background metric ĝ. Therefore, the world sheet Σ may be equipped with two metrics: the background metric ĝ and the quantum metric e γϕ ĝ. Furthermore, conditionally on a fixed background metric ĝ (just discarding ĝ from the randomness), the metric e γϕ ĝ and the matter field are independent as may be seen from the resulting partition function. Knizhnik, Polyakov and Zamolodchikov have derived in [94] a relation between the scaling exponents of the background metric and the quantum metric e γϕ ĝ, the so-called KPZ formula. This is a very rough description of Liouville quantum gravity and the reader may consult [25, 94, 24, 47] for further details. Let us now explain why this KPZ formula may be of interest in the study of models of statistical physics at their critical point. Physicists understood a long time ago (see [6, 7, 94, 47, 48] and certainly many others) that this continuum model of quantum gravity admits a discretized counterpart via random triangulations (or other p-angulations) of surfaces. The prototype of such p-angulations is the Brownian map studied in (see [05, 06, 07, 4]) and corresponds to the pure gravity case c = 0. But we may also couple a model of statistical physics (for instance random walks, percolation, Ising model, Potts model,...) to discrete quantum gravity, i.e. by considering a model of statistical physics on the p- angulation in such a way that, as in the continuum case, the partition function involves both the p-angulation and that of the model of statistical physics. The point is that, at their critical point, these models in two dimensions should behave as a conformal field theory and may be thought of as the matter field described above. By taking the limit as the discretization step goes to 0, these models of discrete quantum gravity should converge towards Liouville quantum gravity. Interestingly, the independence of the Liouville field and the matter field (conditionally on the background metric) suggests that the same phenomena should occur when taking the limit in the discrete model. Therefore the KPZ formula may be applied to this discrete matter field, which are roughly independent of the fluctuating quantum metric: it becomes particularly useful when the scaling exponents 29

30 of a particular model can be more easily computed in its quantum gravity form than its background one, or vice-versa. For instance, Duplantier [5] have used these techniques to conjecture the exact values of the Brownian intersection exponents, which were finally rigorously derived in [0, 02, 03] via Schramm-Loewner-Evolution (SLE) techniques. The reader may consult [45, 46, 66] for pure string models with c = 0 or c = 2, and [30, 53, 89, 90] for critical systems on random p-angulations like Q-Potts model, percolation or tree like polymers. Understanding Liouville quantum gravity from a mathematical rigorous angle is a wide task, which mathematicians have tackled only recently and may take on various aspects, some of them have obviously connections with Gaussian multiplicative chaos theory. As explained above, the mathematical formulation of the problem of constructing (critical) 2d-Liouville quantum gravity could be roughly summarized as follows: construct a random metric on a two dimensional Riemannian manifold D, say a domain of R 2 (or the sphere) equipped with the Euclidean metric dz 2, which takes on the form e γx(z) dz 2 (5.4) where X is a Gaussian Free Field (or possibly other Free Fields) on the manifold D and γ [0, 2) is a coupling constant. Figure 3: Simulation of a GFF on the unit torus. The issue of constructing the distance associated to the metric remains unsolved. Yet recent progress are made in [7, 72, 30] concerning the Brownian motion, Laplace- Beltrami operator or heat kernel of 2d-Liouville quantum gravity. Nevertheless, we focus below on the volume form, which obviously falls under the scope of Gaussian multiplicative 30

31 chaos theory. This theory allows us to construct a random measure of the type: γ2 γx(x) M γ (dx) = e 2 E[X(x)2] dx, (5.5) which will be called Liouville measure. The points that we will address below are the following. In order to apply Gaussian multiplicative chaos theory and define the above measure, we have to choose a cutoff approximation of the GFF and there are several possible choices, which we will discuss. Furthermore we will explain why these cutoff approximations lead to the same limiting measure M γ. Recall that the GFF over a bounded simply connected domain D with for instance Dirichlet boundary condition is a centered Gaussian distribution with covariance kernel given by the Green function G of the Laplacian, i.e. G(x ) = 2πδ x, with Dirichlet boundary condition. Actually, other types of boundary conditions may be imposed but it suffices to detail the Dirichlet boundary conditions to draw a clear picture of the techniques involved. Decomposition of the GFF via eigenfunctions of the Laplacian Let us consider the eigenfunctions (e n ) n of the Laplacian with Dirichlet boundary conditions. They form an orthonormal basis of L 2 (D) with negative associated eigenvalues (λ n ) n. A natural choice of decomposition of the GFF is to write (formally): X(x) = n Y n (x) where Y n is a smooth Gaussian field defined by Y n (x) = β n λn e n (x). Here (β n ) n is a sequence of i.i.d standard Gaussian random variables. Observe that the sequence (β n ) n can be chosen to be measurable with respect to the whole GFF distribution: it suffices to choose β n = ( λ n ) 2 X(x)e n (x) dx. The covariance kernel of Y n matches D k n (x, y) = ( λ n ) e n (x)e n (y). It is well known that the eigenfunctions are smooth so that Y n is a smooth Gaussian field. The important point here is that the approximating sequence n n, X n (x) = Y k (x) is almost surely defined as a function of the whole GFF distribution X. Furthermore, the sequence (Y n ) n are independent Gaussian processes. Though elegant and simple, this decomposition also possesses drawbacks because the covariance kernel of each Y n is not nonnegative and, generally speaking, it is hard to get a tractable expression of k n (or rather their partial sums), excepted maybe in terms of lattice approximations (discrete GFF, see [34]). 3 k=

32 Remark 5.3. Actually, any orthonormal basis of H0 (D) produces a decomposition of the GFF function à la Kahane, i.e. a sum of independent Gaussian processes with continuous covariance kernels. Another very important decomposition relying on an H0 (D)-basis is the projection of the GFF onto the Haar basis. In that case, the corresponding kernels k n are continuous and positive. White noise decomposition of the GFF Another possible decomposition of the Green function is based on the formula: G D (x, y) = π 0 p D (t, x, y)dt. where p D is the (sub-markovian) semi-group of a Brownian motion B killed upon touching the boundary of D, namely p D (t, x, y) = P x (B t dy, T D > t) with T D = inf{t 0, B t D}. Note that the π term ensures that: Hence we can write: G D (x, y) ln x y 0 x y. G(x, y) = 2 K n (x, y) with K n (x, y) = π n p D (t, x, y)dt, n (5.6) n 0 2 n and K 0 (x, y) = π p D (t, x, y)dt. The continuity of p D implies that K n is continuous. The symmetry of p D implies that K n is positive definite. Indeed, for each smooth function ϕ with compact support in D, we have for n : D D ϕ(x)k n (x, y)ϕ(y) dx dy = 2π D D = π D = π 0. D 2 n 2 n 2 n 2 n ϕ(x)p D (t, x, y)ϕ(y) dt dx dy 2 n ϕ(x)p D (t/2, x, z)p D (t/2, z, y)ϕ(y) dt dz dx dy D 2 n D ( D ϕ(x)p D (t/2, x, z) dx) 2 dt dz Since K n is obviously positive, we can apply Kahane s theory of Gaussian multiplicative chaos to define the Liouville measure (5.5). We further stress that this argument implies a white noise decomposition of the underlying GFF: the most direct way to construct a GFF is then to consider a white noise W distributed on D R + and define X(x) = π p D ( s, x, z)w (dz, ds). D R

33 One can check that E[X(x)X(y)] = π 0 p D (s, x, y) ds = G D (x, x ). One can even work with a continous parameter ɛ and define the Liouville measure as the almost sure limit as γ2 γxɛ(x) ɛ 0 of M ɛ (dx) = e 2 E[Xɛ(x)2] dx where the corresponding cut-off approximations X ɛ are given by: X ɛ (x) = π p D ( s, x, z)w (dz, ds). D [ɛ 2, [ 2 Indeed, within this framework introduced in [28], the sequence (M ɛ (A)) ɛ>0 is a positive martingale for all compact set A. Note the following expression for the covariance of X ɛ : E[X ɛ (x)x ɛ (y)] = π ɛ 2 p D (s, x, y) ds Once we define the Liouville measure with this white noise construction, it is not hard to see that we fall in fact under the scope of theorem 3.5 (see our theorem 5.5 below). In particular, we claim: Lemma 5.4. The sequence X ɛ is a smooth Gaussian approximation of G D. We will sketch a proof of this point in the appendix. Circle average As explained in subsection 3.3, the authors in [56] have suggested a slightly different approach: instead of using the σ-positivity of the covariance kernel of the GFF to construct an approximating sequence (2.2) that is a martingale, they regularize the GFF along circles to construct their approximating sequence. More precisely, consider a GFF X and define X ɛ (z) as the mean value of X along the circle centered at z with radius ɛ, formally understood as: X ɛ (x) = 2π X(x + ɛe iθ ) dθ. 2π 0 The covariance kernel is given by G ɛ (z, z ) = G D (x, y) µ x ɛ (du)µ y ɛ (dv) where µ x ɛ (du) stands for the uniform probability measure on the circle centered at x with radius ɛ. This expression can be given a rigorous sense [56]. The main advantage of this construction is that it is well fitted to play with the spatial Markov property of the GFF. Nevertheless, the increments (X ɛ X ɛ ) ɛ<ɛ are not independent, getting trickier the proof of the almost sure convergence of the chaos. On the other hand, this circle average construction falls under the scope of the regularization procedures developed in [32] in order to get convergence and uniqueness in law. Equivalence of the constructions The first question that you must have in mind is: To which extent do the above cut-off approximations yield the same limiting multiplicative chaos?. We claim: 33

34 Theorem 5.5. The law of the limiting chaos does not depend on the cutoff approximations listed above, namely white noise decomposition, eigenvalues of the Laplacian, H 0 (D) expansions or circle average (or more generally convolution by C functions or averages on smooth domains, like ball-averages...). Before proving this theorem, let us make some further comments. In dimension 2, the Lebesgue measure is obviously in the class R 2 ɛ + for all ɛ > 0 and the Green function can be rewritten as G D (x, y) = ln + + g(x, y) (5.7) x y for some bounded continuous function g. Therefore, all the Kahane machinery applies. In particular, Theorem 2.5 ensures that the Liouville measure is non trivial if and only if γ 2 < 4, whatever the choice of the cut-off approximation. Proof of Theorem 5.5. Almost sure equivalence between a given H0 (D) expansion and circle average is already proved in [56], and therefore equivalence in law holds. First proof: Therefore, it suffices to prove equivalence in law between the white noise decomposition and the circle average construction to get theorem 5.5. In view of theorem 3.5, one could for instance establish: For all δ > 0 and A > 0, sup G D (x, y) µ x ɛ (du)µ y ɛ (dv) π x,y D (δ), x y Aɛ where C A > 0 is some constant independent from ɛ. ɛ 2 p D (s, x, y) ds C A. For all δ > 0, C A = lim ɛ 0 sup x,y D (δ), x y Aɛ G D (x, y) µ x ɛ (du)µ y ɛ (dv) π ɛ 2 p D (s, x, y) ds goes to 0 as A goes to infinity. in order to prove equivalence between circle average and white noise decomposition. These two estimates are not very difficult to obtain but we will not detail this point since a second very direct proof is possible. Second proof: Therefore, it suffices to prove equivalence in law between the white noise decomposition and one H0 (D) expansion. Just notice that the white noise decomposition and the expansion along the Haar basis correspond to two σ-finite decompositions (2.) of the Green function. Hence, by theorem 2.3, the two constructions are equivalent in law. 34

35 KPZ formula: almost sure Hausdorff version As explained above, the KPZ formula has been introduced in Liouville quantum gravity and can be thought of as a bridge between the values of the scaling exponents computed with the quantum metric e γx(z) dz 2 and the scaling exponents computed with the standard Euclidian metric. Dealing with metrics is convenient to have a direct definition of the scaling exponents but, as previously explained, a rigorous construction of the quantum metric has not been achieved yet. Nevertheless, a definition of scaling exponents via measures instead of metrics is also possible, and this is what we discuss below. We consider the Liouville measure over a bounded domain D R 2 : M(dx) = γ2 γx(x) e 2 E[X(x)2] dx (5.8) where γ 2 < 4 and X is a GFF on D, say with Dirichlet boundary conditions. If M were a random metric, we could associate a notion of (random) Hausdorff dimension to this metric. Since M is only a measure, the associated notion of Hausdorff dimension is not straightforward, excepted maybe in dimension [27, 28]. We can nevertheless associate to the measure M a notion of Hausdorff dimension: this just consists in replacing carefully quantities related to distances in the standard definition of Hausdorff dimension with similar quantities defined in terms of measures. This yields: given a Radon measure µ on R d and s [0, ], we define for a Borelian set A of R d : Hµ s,δ (A) = inf { µ(b k ) s} where the infimum runs over all the covering (B k ) k of A with open Euclidean balls with radius r k δ. Since the mapping δ > 0 Hµ s,δ (A) is decreasing, we can define the s-dimensional µ-hausdorff metric outer measure: H s µ(a) = lim δ 0 H s,δ µ (A). The limit exists but may be infinite. Since H s,δ µ is metric, all the Borelian sets are H s µ- measurable. The µ-hausdorff dimension of the set A is then defined as the value k dim µ (A) = inf{s 0; H s µ(a) = 0}. (5.9) Notice that dim µ (A) [0, ]. When µ is diffuse (without atoms), the µ-hausdorff dimension of a set A can also be expressed as: dim µ (A) = sup{s 0; H s µ(a) = + }. (5.0) Therefore, when µ is diffuse, the above relations allow us to characterize the µ-hausdorff dimension of the set A as the critical value at which the mapping s H s µ(a) jumps from + to 0. For a given compact set K of D (or a random compact set independent of M), the KPZ formula establishes a relation between the Hausdorff dimension of K computed with µ = M, call it dim M (K), and the Hausdorff dimension of K computed with µ equal to the Lebesgue measure, call it dim Leb (K). We claim (see [28] for a proof of this statement, or also [7]): 35

36 Theorem 5.6. KPZ formula [Rhodes, Vargas, 2008] Let K be a compact set of D. Almost surely, we have the relation: dim Leb (K) = ( + γ2 4 )dim M(K) γ2 4 dim M(K) 2. We develop below a heuristic to understand what is behind the KPZ formula, mainly the power law spectrum of the measure as explained in Theorem 2.4. To begin with, we recall the definition of the s-dimensional M-Hausdorff metric outer measure H s M(K) = lim δ 0 inf { M(B(x n, r n )) s ; n K n } B(x n, r n ), r n δ. Take the expectation and perform an outrageous inversion of limits: { E[HM(K)] s = lim inf E[M(B(x n, r n )) s ]; K } B(x n, r n ), r n δ. δ 0 n n Now compute the expectations via Theorem 2.4 to get: { E[HM(K)] s C s lim inf rn ξ(s)/2 ; K δ 0 n n =C s H ξ(s)/2 Leb (K). } B(x n, r n ), r n δ Because ξ(s)/2 = ( + γ2 ) γ 2 s 4 4 s2, we recover at least heuristically the KPZ formula. Nevertheless, we draw attention to the fact that we do not claim that the relation E[HM(K)] s H ξ(s)/2 Leb (K) is true. There are possibly logarithmic corrections in the choice of the gauge function involved in the definition of the s-dimensional Hausdorff measures for such a relation to be true. KPZ formula: expected box counting version In this subsection, we summarize the KPZ statements proved in [56]. The KPZ theorem of [56] relies on the notion of expected box counting dimension as a definition of scaling exponents. To state the theorem, one must introduce the following definition: Definition 5.7. isothermal quantum ball. For any fixed measure µ on D, let B δ (z) be the Euclidean ball centered at z with radius given by µ(b δ (z)) = δ. If there does not exist a unique δ with this property, take the radius to be sup{ε : µ(b ε (z)) δ}. When µ is the measure M of (5.8), the ball B δ (z) is called the isothermal quantum ball of area δ centered at z. When µ is the Lebesgue measure then B δ (z) is nothing but the Euclidean ball centered at z and radius ɛ where δ = πɛ 2, denoted by B ɛ (z). 36

37 Given a subset K D, the ɛ-neighborhood of K is defined by: B ɛ (K) = {z : B ε (z) K }. The isothermal quantum δ-neighborhood of K is defined by: B δ (K) = {z : B δ (z) K }. Finally, the authors of [56] introduce the notion of scaling exponent. Fix γ [0, 2) and let λ denote Lebesgue measure on D. A fractal subset K of D has Euclidean expectation dimension 2 2x and Euclidean scaling exponent x if the expected area of B ɛ (K) decays like ɛ 2x = (ɛ 2 ) x, i.e., log Eλ(B ε (X)) lim ε 0 log ε 2 = x. The set K has quantum scaling exponent if we have log EM(B δ (X)) lim =. δ 0 log δ Theorem 5.8. [Duplantier, Sheffield, 2008] Fix γ [0, 2) and a compact subset K of D. If K has Euclidean scaling exponent x 0 then it has quantum scaling exponent, where is the non-negative solution to x = γ ) ( γ2. (5.) 4 This theorem also extends to the case where K is a random compact set independent of M. In the physics litterature, the KPZ relation is usually stated under this form (5.) in which case x and are the weights of conformal operators. To get a formulation in terms of dimensions, one must make the correspondence 2 2x dim Leb (K) and 2 2 dim M (K). Let us finally mention that in [56] is also proved a one dimensional boundary version of KPZ. This corresponds to proving the theorem with the lognormal MRM measure of section 5.. Remark 5.9. Further comments and references on KPZ. The KPZ formula has been proved in [27] in the case of multiplicative cascades in dimension (see also [3] for a multidimensional version), in [56] in the case where X is a 2-dimensional GFF, and in [28] (see also [7]) in the case where X is a log-correlated infinitely divisible field in any dimension. Roughly speaking, infinitely divisible fields are to the family of random distributions what Lévy processes are to the family of stochastic processes. Log-correlated Gaussian fields, like two dimensional Free Fields, are a subclass of log-correlated infinitely divisible fields. Therefore, the main point here is to draw attention to the fact that the KPZ formula is a property specific to log-correlated fields: it is neither specific to the dimension, nor to the conformal invariance of the 2d-GFF, nor to the Gaussian nature: the only point that makes 2d-Liouville quantum gravity (i.e. Gaussian multiplicative chaos with respect to the 2d-GFF) satisfy a KPZ relation is the fact that the Green function of the Laplacian in dimension 2 (and in dimension 2 only) has a logarithmic singularity. Also, it may be interesting to know if a Liouville quantum gravity picture can be drawn for log-correlated infinitely divisible fields instead of Gaussian Free Fields. 37

38 Liouville quantum gravity and KPZ on Riemannian surfaces One may wonder what becomes Liouville quantum gravity and the KPZ formula on a n- dimensional Riemannian manifold (S, g) where g is the Riemannian tensor of the manifold. By Liouville quantum gravity, we mean here a Gaussian multiplicative chaos with respect to a Gaussian distribution X defined on the manifold. As long as the random Gaussian distribution X possesses a kernel of σ-positive type, Kahane s theory allows to define a Gaussian multiplicative chaos associated to this Gaussian distribution. If the covariance kernel of the Gaussian distribution is of the type (2.5) (where ρ is the distance associated to the Riemannian metric) and the measure σ in (2.4) is the volume form on S, then the non-degeneracy conditions of the chaos is γ 2 < 2n (Theorem 2.5). Since a n-dimensional Riemann surface is locally isometric to the unit ball of R n, we deduce from [28] (or [7]) that the KPZ formula holds for the Gaussian multiplicative chaos M on this Riemann surface: it reads Theorem 5.0. KPZ formula on Riemann manifolds [[28], 2008] Let K be a compact set of S. Almost surely, we have the relation: dim σ (K) = ( + γ2 2n )dim M(K) γ2 2n dim M(K) 2. In particular, we see that the curvature of the surface does not affect the KPZ relation. For instance, we can consider a 2-dimensional Riemann surface, like a sphere or an hyperbolic half-plane, and the GFF on a domain of this surface with appropriate boundary conditions in order to define the associated Liouville measure. As explained above, in dimensions different from 2, the GFF does not possess logarithmic correlations so that it does not make sense to look for KPZ relations based on the GFF. Nevertheless, in dimensions different from 2, it is plain to construct other log-correlated Gaussian distributions X: various examples of log-correlated Gaussian fields are described in the present manuscript but also in [59]). Another situation of interest is to consider massive or generalized Free Fields. On a domain D R 2, the Massive Free Field (MFF) with Dirichlet boundary conditions is defined as a standard Gaussian in the Hilbert space defined as the closure of Schwartz functions over D with respect to the inner product (f, g) h = m 2 (f, g) L 2 (D) (f, g) L 2 (D). The real m > 0 is called the mass. Its action on L 2 (D) can be seen as a Gaussian distribution with covariance kernel given by the Green function G D m of the operator m 2, i.e.: (m 2 )G D m(x, ) = 2πδ x with Dirichlet boundary conditions. When D is the whole plane, the massive Green kernel is a star-scale invariant kernel in the sense of [5]. We may also consider Generalized Free Fields as defined in [75]. Definition 5.. A Generalized Free Field with Dirichlet boundary conditions over a domain D R 2 is defined as a random centered Gaussian distribution (say on the space of Schwartz functions on D) the covariance kernel of which is given by G D ϱ (x, y) = G D m(x, y) ϱ(dm).

39 where G D m is the massive Green function on D with mass m and ϱ is a Radon measure on R +, called the Källen-Lehmann weight, satisfying and 0 k N, ln m ϱ(dm) < m k ϱ(dm) < +. The construction of such a field can be straightforwardly adapted from the previous white noise decomposition of the free field. Denote by p D (t, x, y) the transition densities of a Brownian motion killed upon touching the boundary of D. Consider a Gaussian noise, white in space and time and ϱ-colored in mass, i.e. a Gaussian random measure W (dx, ds, dm) distributed on D R + R + such that for all Borel sets A, A D, B, B, C, C R + E[W (A, B, C)W (A, B, C )] = A A B B ϱ(c C ), where A stands for the Lebesgue measure of A. The generalized free field can then be defined X(x) = m π 2 4 s p D ( s, x, y) W (dy, ds, dm), 2 and cutoff approximations: ɛ 0, X ɛ (x) = π D R + R + e D [ɛ,+ [ R + e m 2 4 s p D ( s, x, y) W (dy, ds, dm). 2 For all these fields, theory of Gaussian multiplicative chaos and Theorem 5.0 apply. Let us just stress that the result of the KPZ formula remains unchanged for Massive Free Fields whereas γ 2 must be replaced with γ 2 ϱ(r + ) for Generalized Free Fields. A last example of interest is the boundary Liouville measure introduced in [56]. The authors suggests to consider a smooth domain D of R 2 together with a GFF X on D with free boundary conditions and to define the boundary Liouville measure on D as ν(dx) = lim e ɛ 0 γ γ2 Xɛ(x) 2 4 ln ɛ dx where dx stands for the length measure on D and X ɛ (x) is the mean value of X over B(x, ɛ) D. Actually, the boundary Liouville measure is nothing but a Gaussian multiplicative chaos over a -dimensional Riemannian manifold (in fact C is enough). By noticing that the correlations are given by E[ γ 2 X ɛ(x) γ 2 X ɛ(y)] = γ2 2 ln + + g(x, y) (5.2) x y for some continuous bounded function g, Theorem 2.5 ensures that the boundary Liouville measure is non-degenerate provided that γ 2 < 4. Furthermore, the boundary KPZ holds: just take care of replacing γ 2 in Theorem 5.0 by γ 2 /2 because of the unusual normalization in 5.2 (the term in front of the log is γ2 2, averaging along semi-circles yields an extra factor 2). The expected box counting KPZ formula for the boundary Liouville measure is proved in [56] whereas the Hausdorff dimension version of the KPZ formula for boundary Liouville measure is proved in [28]. 39

40 5.3 Convergence of discrete Liouville measures on isoradial graphs Here we consider a planar (isoradial) graphs. We weight the vertices of this graph by the exponential of the Discrete Gaussian Free Field (DGFF for short) to obtain a discrete Liouville measure: the measure having as density exponential of the DGFF with respect to the discrete canonical volume measure on the graph. We will prove that this measure weakly converges in law towards the Liouvile measure as the mesh size of the graph converges to 0. Before stating a clear theorem, we need to explain the framework in further details. The reader may find in [37] all the basic tools (and much more) about isoradial graphs described below. We stick to the notations used in [37]. A planar graph Γ embedded in C is called isoradial iff each face is inscribed into a circle of common radius ɛ. If all circle centers are inside the corresponding faces, then one can naturally embed the dual graph Γ in C isoradially with the same ɛ, taking the circle centers as vertices of Γ. The name rhombic lattice is sometimes used for these graphs because all the quadrilateral faces of the corresponding bipartite graph Λ (having vertices Γ Γ ) are rhombi with sides of length ɛ. We will make the following assumption (see [37]) the rhombi half-angles are uniformly bounded away from 0 and π 2. (5.3) Roughly speaking, Λ does not possess too flat rhombi. This entails that the Euclidean distance between the vertices of Γ is comparable to the graph distance. For such graphs, the canonical volume measure µ Γ (A) of a subset A Γ is given by µ Γ (A) = z A W z where the weight W z of z Γ is the Lebesgue measure of the face of the dual graph containing z. We consider a sequence (Γ n ) n of isoradial graphs as indicated above with radius (ɛ n ) n. For all the quantities defined above, the subscript n means that it is related to the graph Γ n. For instance, µ n stands for the volume measure of the graph Γ n. We also assume that the radius ɛ n of Γ n goes to 0 as n. Let us now consider a bounded simply connected domain D of C. Let us denote by D n the graph D Γ n, i.e. we keep the vertices and edges that entirely lie in D. We consider a DGFF X n on the vertices of D Γ n. Recall that the DGFF is a collection (X n (z)) z Dn of centered Gaussian random variables with covariance kernel given by 2πG Dn, where G Dn denotes the Green function on D n with 0-boundary condition (see [37, Definition 2.6]). Then, for γ [0, 2[, we define the discrete Liouville measure on D n : γ2 γxn(z) M n,γ (dz) = e 2 E[Xn(z)2] µ n (dz). (5.4) Theorem 5.2. Let D be a bounded simply connected open domain in C. For γ [0, 2[, the discrete Liouville measures (M n,γ (dz)) n on D n weakly converge in law towards the Liouville measure on D (see subsection 5.2), that is the Gaussian multiplicative chaos: γ2 γx(x) M γ (dx) = e 2 E[X(x)2] dx where X is a GFF on D with 0-boundary condition. 40

41 The proof can be found in Appendix B. In fact, we will prove the result for the discrete GFF on the square lattice with mesh ɛ going to 0. It turns out that the proof of Theorem 5.2 also works at criticality (see subsection 6. for further details and references), yielding the following: Theorem 5.3. For γ = 2, the discrete critical Liouville measures ( ln ɛ n M n,2 (dz)) n on D n weakly converge in law towards the critical Liouville measure on D (see [57, 58] or subsection 6. below), that is the Gaussian multiplicative chaos: M (dx) = 2 π (2E[X(x)2 ] X(x))e 2X(x) 2E[X(x)2] dx where X is the GFF on D with Dirichlet boundary condition. 5.4 Kolmogorov-Obhukov model in turbulence We refer to [65] for an introduction to the statistical theory of 3 dimensional turbulence. Consider a stationary flow at high Reynolds number, that is when the velocity of the fluid is large in comparison with the viscosity forces. It is believed that at small scales the velocity field of the flow is homogeneous and isotropic in space. By small scales we mean scales much smaller than the integral scale R characteristic of the time stationary force driving the flow. In the works [96] and [8], Kolmogorov and Obukhov proposed to model the mean energy dissipation per unit mass in a ball B(x, l) of center x and radius l << R by a random variable ɛ l such that ln(ɛ l ) is normal with variance σl 2 given by: σ 2 l = λ 2 ln( R l ) + A where A is a constant and λ 2 is the intermittency parameter. As noted by Mandelbrot ([]), the only way to define such a model is to construct a random measure ɛ by a limit procedure. Then, one can define ɛ l by the formula: ɛ l = 3 < ɛ > 4πl 3 ɛ(b(x, l)) where < ɛ > is the average mean energy disspation per unit mass. Formally, one is looking for a random measure ɛ such that: A B(R d γ2 γx(x) ), ɛ(a) = e 2 E[X(x)2] dx (5.5) where (X(x)) x R d is a Gaussian field whose covariance kernel K is given by (2.2). Therefore, one can give a rigorous meaning to energy dissipation (5.5) by using Gaussian multiplicative chaos theory. Let us mention here that the objective of describing a stochastic representation of the velocity field is a much more ambitious task (see subsection 6.3 for more comments and perspectives on this topic). A 4

42 5.5 Decaying Burgers turbulence Consider the Burgers equation t v + v v = ν v + f(x, t) with initial condition v(x, 0) = v 0 (x). (5.6) ν is the viscosity parameter. When f 0, this equation is called (randomly) forced Burgers equation. When f = 0 and v 0 (x) 0, this equation is called decaying Burgers turbulence. We consider here the case f = 0 with random initial data v 0. The solution of decaying Burgers turbulence, via Hopf-Cole transform, is given by v(t, x) = ( 2ν ln Z(t, x)), where Z(t, x) = R e y x 2 2ν 2t 2ν V (y) dy. (5.7) 4πνt where v 0 = V. If we choose V as a log-correlated Gaussian random potential, i.e. with covariance kernel of the type (2.2) then Z appears as a Gaussian multiplicative chaos of the type (2.4) and integrating measure σ given by the standard heat kernel on R. This situation corresponds to a power law correlated random profile of Gaussian distributed initial velocities: E[v 0 (x)v 0 (y)] x y 2. Gaussian distributed initial velocities with power law correlations is a subject of interest in Burgers turbulence. The reader may consult [69] for an account of physical motivations, further references and a study of the present situation. It is shown in [69] that this equation exhibits a phase transition at a critical viscosity parameter ν c, which is related to the phase transition of Gaussian multiplicative chaos (see Theorem 2.5). For ν < ν c, freezing phenomena occur, highlighting a glassy phase conjecturally as that explained in subsection 6.2 (frozen phase). 6 Generalizations of the theory In this section, we review some generalizations of Kahane s theory while making connections with possible applications. 6. Gaussian multiplicative chaos at criticality Kahane s construction of Gaussian multiplicative chaos makes sense when the factor γ appearing in (2.) satisfies γ 2 < 2d. This gives rise to the issue of constructing random measures in the same spirit for γ 2 2d. Necessarily these measures will present a different structure. The case γ 2 = 2d is of special interest since it corresponds to a phase transition. We call this situation the critical case. Kahane s theory ensures that the associated martingale (M n ) n appearing in (2.2) almost surely converges towards 0. Several approaches are possible in order to make sense of a suitable random measure corresponding to criticality γ 2 = 2d. It turns out that all these approaches are conjecturally the same. The first approach is based on the convergence of a fundamental object that is called derivative martingale. Such an object has been intensively studied in the case of multiplicative cascades, branching random walks [28, 98] or branching Brownian motions [6]. The convergence is achieved in [57] in the context of Gaussian multiplicative chaos associated 42

43 to star scale invariant kernels (2.23) and may be achieved for other type kernels provided that one can use a white-noise decomposition cutoff (see [58]), for instance including the GFF in a bounded domain (see previous subsection White noise decomposition of the GFF). The limiting measure can formally be written as M (A) = A ( γ E[X 2 (x)] X(x) ) γ2 γx(x) e 2 E[X2 (x)] dx with γ = 2d. (6.) The reader may object that this construction should be possible for every value γ 2 2d and therefore argue that the interest for the only value γ 2 = 2d is not natural. It turns out that there is an abrupt change in the behaviour of this object at the critical value: the exponential term penalizes those points x where the value of the process X(x) lies above the expectation term γe[x(x) 2 ] with a strength that depends on γ: for γ 2 < 2d the measure may be well defined but the penalization term is not strong enough so that the sign of the term ( γ E[X 2 (x)] X(x) ) alternates. For γ 2 = 2d, the exponential term is strongly penalizing, forcing the measure M to be nonnegative, which is not straightforward at first sight. Furthermore, in the situation when the random distribution X possesses a star scale invariant kernel, the limit of the derivative martingale M yields a solution to the star equation for the only value γ 2 = 2d. Indeed, when trying to derive this scaling relation for M, you are left with an unusual extra term proportional to the standard chaos M. For γ 2 < 2d, this term does not vanish whereas for γ 2 = 2d Kahane s theory ensures that this extra term disappear, making the measure M star scale invariant. It is also proved in [57] that the random measure M has almost surely full support and no atom. Further improvements are made in [20]: the authors determine the exact asymptotics of the right tail of the distribution of the total mass of the measure, and an almost sure upper bound for the modulus of continuity of the cumulative distribution function of the measure. A lower bound for the increments of the measure is also investigated, showing that the measure is supported on a set of Hausdorff dimension 0. Another possible approach is to find a suitable renormalization of (2.2) in order to get a non trivial limit. This approach is carried out in [58], based on the works [99] for branching Brownian motions or [, 83] for branching random walks. The renormalization turns out to be the square root of the variance of the Gaussian field X n, the convergence holds in probability and the limit measure is the derivative measure M, up to a perfectly explicit multiplicative factor. This approach is very convenient in order to use Kahane s convexity inequalities with the measure M, which is not obvious at first sight when looking at the expression (6.). In particular, a complete description of the moments of M is achieved in [57, 58] as well as the computation of the power-law spectrum of M. Let us point out that another approach of criticality via a KPP equation has been investigated in the case of multiplicative cascades [9] or branching Brownian motion [99]. No Gaussian multiplicative chaos counterpart has ever been derived rigorously. The main reason is the structure of correlations, which are more intricate than in the discrete framework. But this is clearly a perspective of interest. In Liouville quantum gravity, the Liouville measure at criticality γ = 2 corresponds to a central charge c =. It is expected, for instance, to be the scaling limit of the O(n = 2) loop model or the Q = 4- Potts model (see [7] for a description of these models). Let us mention that, in the standard physics literature about c = Liouville field theory 43

44 Figure 4: Height landscape of the derivative martingale measure plotted with a logarithmic scale color-bar, showing that the measure is very peaked (for t = 2, a multiplicative factor of about 0 8 stands between extreme values, i.e., between warm and cold colors). [32, 73, 74, 79, 88, 92, 9, 23], the Tachyon field ϕe 2ϕ presents an unusual dependence on the Liouville field (or GFF) ϕ (compare with the usual form e γϕ for c < ). This is clearly in relation with the mathematical formulation (6.). The KPZ formula (Theorem 5.6) has been proved in [58] at criticality, i.e. for γ = Atomic Gaussian multiplicative chaos We know review some recent progress motivated by the super-critical regime γ 2 > 2d. They have been mainly inspired by the seminal paper [60] in the context of multiplicative cascades. Typically, a new class of multiplicative chaos emerges, called atomic multiplicative chaos, which can be split in two parts. Dual phase Before coming to mathematical considerations, let us first briefly outline the physics motivations for considering the dual phase. The name dual phase comes from physics literature 44

45 about Liouville quantum gravity. Remind that the coupling constant γ appearing in (5.4) or (5.5) is determined by the underlying model coupled to gravity. It is related to the so-called central charge of the model by the relation [94] c = 6 ( γ 4 ) 2. 4 γ The central charge belongs to ], ]. The special case c = (and γ = 2) corresponds to criticality. Otherwise, the central charge belongs to ], [ yielding two possible γ > 0. The first solution γ belongs to ]0, 2[ and is given γ = 25 c c 6. (6.2) It corresponds to the standard branches of Liouville gravity detailed in subsection 5.2. The second solution, call it γ, belongs to ]2, + [: γ = 25 c + c 6. (6.3) Observe that γγ = 4. In a series of papers [9, 44, 52, 9, 92, 93, 97], physicists have investigated what they called the other/non-standard/unconventional/dual branches of gravity initially via modified random matrix models generating random surfaces in order to interpret these other possible values of the coupling constant γ ]2, + [. They have noticed several intriguing relations between the standard and dual branches of the Liouville action, laying the foundations of what they called Duality of Liouville quantum gravity. More recently, duality of Liouville quantum gravity has been digged up in [55, 50] at an heuristic yet interesting level. The purpose of what follows is to present mathematical results and conjectures on the mechanisms involved in this dual phase, which turns out to be a very rich area from the mathematical or physics angles. To understand mathematically how to handle the dual phase, let us continue this pedagogical introduction in the context of Liouville quantum gravity. So we consider a bounded domain D of R 2 and a GFF X on D with Dirichlet boundary conditions. The first observation is that Theorem 2.0 tells us that the standard chaos M γ (dx) = e γx(x) γ2 2 E[X(x)2] dx reduces to 0 because γ 2 > 4. Therefore another construction has to be found. To understand which construction is involved, let us assume for a while that the chaos M γ does not reduce to 0 (though it does). Then, as explained in Theorem 2.4 (and in its proof), the chaos M γ would then satisfy the scaling relation γ2 2 M γ (B(x, r)) = r e γωr 2 E[Ω2 r ] M γ (B(x, )) (6.4) for some Gaussian random variable Ω r independent of M γ and variance ln r. A straightforward computation shows that this relation may be rewritten as: M γ (B(x, r)) = (r ) 2 γ2 /α γωr e 2 E[Ω2 r] M γ (B(x, )) (6.5) 45

46 where α = 4 γ (recall that γ γ = 4). This relation suggests that if 2 M γ were non degenerate, it should scale for small r like ( /α. M γ (B(x, r)) M γ (B(x, r))) (6.6) The convenient fact of the above relation is that the Gaussian multiplicative chaos M γ in the right-hand side is non trivial since γ 2 < 4. Therefore, we are heuristically looking for a measure that may be interpreted as the α -th root of the standard chaos M γ. From the mathematical angle, this α-th root perfectly makes sense in terms of independently scattered random measures with prescribed control measure: roughly speaking, this consists in throwing a stable point process over a landscape described by the standard Gaussian multiplicative chaos M γ. More precisely, the law of this measure can be generated as follows: sample the standard Gaussian multiplicative chaos γ2 γx(x) M γ (dx) = e 2 E[X(x)2] dx. sample a random measure M γ whose law, conditionally to M γ, is that of an independently scattered random measure characterized by q 0, E[e qm γ(a) M γ ] = e qα M γ(a). This construction is called subordination procedure. Interestingly, it is a purely atomic random measure as suggested in physics literature. Until now, we have identified, at least in law, what kind of object we are looking for to model the dual branch of Liouville quantum gravity. Let us also mention here that we recover the suggestion made in [50] concerning the dual measure. It seems that the arguments used in [50] to guess the exact form of the dual Liouville measure are based on the knowledge of the scaling exponents involved in the dual branch of gravity. From the theoretical angle, this is the beginning of interesting questions. First, we mention that the scaling relation (6.6) has an exact equivalent in the context of branching random walks [60, 29, 08], the solutions of which are called the fixed points of the smoothing transform. Interestingly, it is proved that the corresponding scaling relation has a unique solution in law, which is the exact equivalent for branching random walks of the measure M γ described above. This is a strong argument to validate the law of the above measure. The question of uniqueness in the context of Gaussian multiplicative chaos theory can be formulated in terms of star scale invariance and will be discussed in subsection 6.4. Second, the above construction of M γ is not quite satisfactory. Indeed, the measure M γ is not a measurable function of the only field X: the measure M γ also depends on the randomness originating from the sampling of the point process. It is not a measure intrinsically generated by the only field X. Third, the construction of M γ by subordination does not appear as an almost sure limit of a suitably renormalized sequence in the spirit of (2.2) for instance. We will see that this does not allow us to understand the mechanisms involved in duality and that this is also related to the latter objection. Finally, a KPZ formula for the dual measure has to be proved rigorously. 46

47 A different way of constructing the measure M γ has been carried out in [7]. It is important to understand here that the main differences with the subordination procedure take place at the level of almost sure properties of the measure M γ. Basically, the construction consists in replacing the dx measure in (2.) by an independently scattered infinitely divisible random measure, call it n α for some α ]0, [, whose law is characterized by q 0, E[e qnα(a) ] = e qα A for all Borelian set A of D. The random measure n α is assumed to be independent of the free field X. For γ 2 > 4, it then makes sense to define an approximating sequence in the spirit of (2.2) (with the same notations): M n, γ (dx) = e γxn(x) α γ2 2 E[Xn(x)2] n α (dx). (6.7) Theorem 6.. Convergence of the atomic Gaussian multiplicative chaos [[7], 202] Almost surely in n α, the sequence of random measures (M n, γ ) n converges in P X - probability towards a non trivial limiting measure M γ (dx) = e γx(x) α γ2 2 E[X(x)2] n α (dx), called atomic Gaussian multiplicative chaos, which is purely atomic. We point out that the above construction is not a Gaussian multiplicative chaos in the usual sense. Indeed the lognormal weight is not normalized to have expectation. The expectation blows up, giving rise to a situation that qualitatively deeply differs from standard Gaussian multiplicative chaos theory. Therefore, the simplicity of Kahane s construction is lost in this construction: ) we do not deal with martingales because the weights are not normalized, 2) these measures are not integrable because the measure n α is not. Surprisingly, the atomic Gaussian multiplicative chaos is independent of the measure n α used in the construction: Theorem 6.2. The random measure M γ obtained as the limit in probability of the sequence (M n, γ ) n defined by (6.7) is independent of the stable measure n α. It is therefore a measurable function of the only field X. Proof. It is standard that the stable measure n α can be decomposed as n α (A) = z m α (dx, dz) A R + where m α is a Poisson random measure on R d R + with intensity given by dx dz. Let z +α + B e γ α X(x) γ 2 us first prove that 2α E[X(x)2] B c z m α (dx, dz) = 0 for all c > 0 and all balls B. To this purpose, observe that this term is the limit in probability as n of + e γ γ2 Xn(x) α 2α E[Xn(x)2] z m α (dx, dz), c 47

48 Figure 5: Chaos and dual chaos for the value γ 2 = 0, 25 (and then α = 0.05). Left: density of the standard chaos. Middle: position and weights of the atoms of the dual measure. Right: position and weights of the atoms of the dual measure with a logarithmic ordinate scale. which is nothing but a finite sum of terms like Z i e γ α Xn(x i) γ2 2α E[Xn(x i) 2], where Z i are random variables independent of the field X. By the strong law of large numbers, for all x R, almost surely in X, we have γ lim n α X n(x) γ2 2α E[X n(x) 2 ] =. Our claim follows. The Kolmogorov 0 law then proves the result. Furthermore, the measure M γ may be seen as the volume form associated to a Riemannian tensor of the form α γ2 e γx(x) 2 E[X(x)2] (n α (dx)) 2 This can be seen by regularizing the measure n α with a sequence of mollifiers (ρ n ) n : M γ appears as the limit as n of the volume forms associated to the smooth metric tensors obtained by mollifying the measure n α with the sequence (ρ n ) n. For instance, consider the setup drawn in subsection 5.3. Instead of considering the canonical volume form µ n on D n, meaning each vertice z of D n has a weight corresponding to the Lebesgue measure of the face of the dual graph containing z, we assign to each vertic z of D n the weight W α z of the n α -measure of the dual face of D n containing z. We define the measure µ n,α (A) = z A W α z. This is some kind of Bouchaud trap model on isoradial graphs. We now consider this kind of Bouchaud trap model in the gravitational dressing, i.e. we consider the discrete GFF 48

49 X n on D n and define the discrete dual Liouville measure for γ ]2, + [ M n, γ (dz) = e γxn(z) α γ2 2 E[Xn(z)2] µ n,α (dz). It is straightforward to see that the sequence of measures (M n, γ ) n converges in law towards the measure M γ. The same conclusion holds if, instead of assigning a stable law on each face of the dual graph, we assign a Random Energy Model at low temperature on the collection of faces of the dual graph of D n. (a) γ 2 = 0.0 (b) γ 2 = (c) γ 2 = 3.6 Figure 6: Influence of γ on the spatial localization of the atoms of the dual measure. Let us further mention that the construction of atomic Gaussian multiplicative chaos carried out in [7] is more general than the situation explained here: it is valid in any dimension d for log-correlated Gaussian fields X, and also for all possible values of α ]0, [ and γ 2 > 2d, meaning that we do not impose α = 2d (= 4 γ in dimension d = 2). Also, a γ 2 2 rigorous proof of the KPZ formula for dual measures appears in [7]. Special care must be taken here to handle purely atomic measures. In Figure 6, we simulate for different value of γ a few atoms of the dual measure (the biggest ones). The colored background stands for the height profile of the associated subcritical measure M plotted with a logarithmic intensity scale: red for areas with large mass 49

50 and blue for areas with small mass. Localization of atoms is plotted in black. The larger γ is, the more localized on areas with large potential the atoms are. Frozen phase The second part of the super-critical regime consists in throwing a stable point process over a landscape described by the critical measure M of subsection 6.:. sample the critical measure M M ( (dx) = γ E[X 2 (x)] X(x) ) γ2 γx(x) e 2 E[X2 (x)] dx with γ = 2d, A 2. sample a point process N α whose law, conditionally to M, is that of an independently scattered random measure characterized by q 0, E[e qn α (A) M] = e qα M (A). This phase is called frozen due to a linearization of the free energy involved in models converging towards these measures. Nevertheless, the term frozen may also be understood by the fact that the landscape on which points are thrown is frozen and matches the critical measure M : only the height of the atoms varies with the parameter α. Furthermore, these family of measures is conjectured to be involved in the glassy phase and freezing phenomena observed in log-correlated random potentials. The reader may consult [33, 67, 68] for an account of physics motivations and results, [9, 23, 09, 38] for rigorous results in the case of discrete models and [7, 57] for precise conjectures in the context of Gaussian multiplicative chaos theory. More precisely, the glassy phase of log-correlated Gaussian potentials is concerned with the renormalization of measures beyond the critical value γ 2 > 2d. For γ 2 > 2d, consider the measure: γ2 γxn(x) M n,γ (dx) = e 2 E[Xn(x)2] dx, where (X n ) n is your favorite cutoff approximation of the Gaussian distribution X with covariance kernel of the type (.2). The limiting measure, as n, vanishes as shown by Theorem 2.0. Therefore, it is natural to wonder how to renormalize the sequence of measures (M n,γ ) n in order to get a non trivial limiting measure. As pointed out in [57], the sequence ( c 3γ 2 2d n e cn ( γ 2 d) 2 M n,γ (dx) ) t 0 is tight and every converging subsequence is non trivial, where c n = Var(X n ). This argument is based on the results in [09, 38]. Let us stress that the definition of c n is clear when X n is stationary. If X n is not stationary, there is usually a clear candidate for the definition of c n, since the behaviour of the variance usually does not depend too much on the spatial localization, like in the case of the GFF. Of course, if one chooses the cutoff approximations in a very bad way, it may happen that the renormalization constant c n does not straightforwardly make sense. We do not discuss here in further details these technical considerations. Based on heuristics on star scale invariance, it is conjectured in [7, 57] that this renormalized sequence actually admits only one possible limit: 50

51 Conjecture 6.3. Prove that: ( γ 2 c e cn ) 2 d M n,γ(dx) law d γ N α (dx), as t (6.8) 3γ 2 2d n where d γ is a positive constant depending on γ and the law of the random measure N α is that described above with α = 2d γ. In particular, physicists are interested in the behaviour of the Gibbs measure associated to M n,γ (dx) on a ball B. It is the measure renormalized by its total mass: From (6.8), we deduce G γ n(dx) = M n,γ(dx) M n,γ (B). G γ law t (dx) N α(dx), as n. (6.9) N α (B) The size reordered atoms of the latter object form a Poisson-Dirichlet process as conjectured by physicists [33]. This Poisson-Dirichlet approach has recently made important progress [8], which is presently the more accurate mathematical result to describe the glassy phase of log-correlated random potentials. A further step is to prove Conjecture 6.3 as it offers a more complete picture of the underlying phenomena than the Poisson- Dirichlet approach. The main reason is that it makes precise the spatial localization of the atoms together with their heights whereas the Poisson-Dirichlet approach only focuses on the heights of the atoms. We stress that this conjecture has been proved in the context of lognormal Mandelbrot s cascades in [9, 38] and in the context of branching random walks in [23], based on results appearing in [09]. (a) α = 0.2 (b) α = 0.5 Figure 7: Localization of atoms in the frozen phase for different values of α. We further stress that the same remarks as in the dual phase remain valid here concerning an almost sure construction of the frozen phase: though the techniques developed in [7] should apply, this remains to be done rigorously. In Figure 7, we simulate for different value of γ a few atoms of the frozen measure N α (the biggest ones). The colored background stands for the height profile of the derivative 5