The Tutte Embedding of the Poisson–Voronoi Tessellation of the Brownian Disk Converges to 8/3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{8/3}$$\end{document}-Liouville Quantum Gravity

Recent works have shown that an instance of a Brownian surface (such as the Brownian map or Brownian disk) a.s. has a canonical conformal structure under which it is equivalent to a 8/3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{8/3}$$\end{document}-Liouville quantum gravity (LQG) surface. In particular, Brownian motion on a Brownian surface is well-defined. The construction in these works is indirect, however, and leaves open a basic question: is Brownian motion on a Brownian surface the limit of simple random walk on increasingly fine discretizations of that surface, the way Brownian motion on R2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb {R}^2$$\end{document} is the ϵ→0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\epsilon \rightarrow 0$$\end{document} limit of simple random walk on ϵZ2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\epsilon \mathbb {Z}^2$$\end{document}? We answer this question affirmatively by showing that Brownian motion on a Brownian surface is (up to time change) the λ→∞\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda \rightarrow \infty $$\end{document} limit of simple random walk on the Voronoi tessellation induced by a Poisson point process whose intensity is λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda $$\end{document} times the associated area measure. Among other things, this implies that as λ→∞\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda \rightarrow \infty $$\end{document} the Tutte embedding (a.k.a. harmonic embedding) of the discretized Brownian disk converges to the canonical conformal embedding of the continuum Brownian disk, which in turn corresponds to 8/3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{8/3}$$\end{document}-LQG. Along the way, we obtain other independently interesting facts about conformal embeddings of Brownian surfaces, including information about the Euclidean shapes of embedded metric balls and Voronoi cells. For example, we derive moment estimates that imply, in a certain precise sense, that these shapes are unlikely to be very long and thin.


1.1.
Overview. This paper concerns relationships between several different topics, including Liouville quantum gravity and the Brownian map. Let us begin by briefly reviewing the objects under consideration.
A planar map is a graph together with an embedding into the plane so that no two edges cross. Two planar maps are considered to be equivalent if they differ by an orientation preserving homeomorphism of the plane. The study of planar maps goes back to work of Tutte [Tut68] and Mullin [Mul67] from the 1960s. A planar map can be viewed as a metric measure space by equipping it with the graph distance and assigning each vertex one unit of mass. In recent years, there has been considerable progress in studying the large scale metric behavior of planar maps chosen uniformly at random. Of particular relevance to the present article are the scaling limit results which give the convergence of uniformly random planar maps towards a continuous object in the Gromov-Hausdorff-Prokhorov topology. The first results of this type were due to Le Gall [Le13] and Miermont [Mie13] which are both focused on random planar maps with the sphere topology. In this case, the limiting object is a random metric measure space with the topology of the sphere [LP08,Mie08] called the Brownian map, which was first defined (in different forms) in [MM06,Le07], building on [CV81,Sch97,CS04]. The works [Le13,Mie13] have since been extended to uniformly random planar maps with several other topologies, including the disk [BM17,GM19d], the plane [CL14], and the half-plane [BMR16,GM17b]. The limiting objects that one obtains are collectively known as Brownian surfaces.
The aforementioned limit theorems are concerned with the metric measure space structure of large uniformly random planar maps, but not how they are embedded into the plane. However, it is an important problem to understand scaling limits for canonical embeddings of random planar maps. This is motivated in part by a desire to better understand the relationship between statistical mechanics models (e.g., percolation, selfavoiding walks, the Ising model) on random planar maps and their counterparts on planar lattices. It has long been believed that if the embedding is of a conformal type (e.g., Riemann uniformization, circle packing, or the Tutte embedding we consider here), then the large scale geometry of the statistical mechanics model should be the same as if it were considered on a planar lattice. This idea underlies the famous KPZ relation [KPZ88,DS11], which serves to convert critical exponents computed on a random geometry to the corresponding exponents on a deterministic geometry, and has been used numerous times to give predictions for exponents for critical lattice models which were later verified using SLE techniques (e.g., [Dup98,LSW01a,LSW01b,LSW02]).
Liouville quantum gravity (LQG) is another theory of random surfaces which was introduced by Polyakov [Pol81a,Pol81b] in the 1980s in the context of string theory. To define LQG, one starts with a (form of) the Gaussian free field (GFF) h on a domain D and then considers the random two-dimensional Riemannian manifold with metric tensor e γ h(z) (dx 2 + dy 2 ) (1.1) where γ ∈ (0, 2] is a parameter. This definition does not make rigorous mathematical sense since h is a distribution and not a function. Making rigorous sense of various aspects of LQG has been a major topic of research in recent years. Some of this work builds on [DS11], which constructs the volume form associated with (1.1), which is a random measure μ h on D (see [Kah85,RV14] for a more general construction of random measures of this type). One can construct various kinds of LQG surfaces by varying the precise definition of h (what domain it lives on, how boundary conditions are chosen, whether one "weights" the law of h in some locally absolutely continuous way). Generally, one writes γ -LQG to refer to LQG surfaces with parameter γ . The special case γ = √ 8/3 has long been known to be special: √ 8/3-LQG surfaces, like Brownian surfaces, are related to "undecorated" planar maps, and are also called pure LQG surfaces.
In fact, a recent series of works by the second two authors has shown that the Brownian map and the so-called √ 8/3-LQG sphere are in some sense equivalent (and similar statements can be made for disk, whole-plane, or half-plane variants) [MS15b,MS16a,MS16b]. Although both √ 8/3-LQG spheres and Brownian maps come with a natural measure, they also have additional structure: a Brownian surface has a metric, and a √ 8/3-LQG surface has a conformal structure (i.e., an embedding into a flat domain, defined up to conformal automorphism of that domain). The works [MS15b,MS16a,MS16b] show that each of these objects can be canonically endowed with the other one's structure, and that once this is done the objects agree in law.
We have made the present work as self-contained as possible, which in particular means that the reader is not required to have read the series [MS15b,MS16a,MS16b] relating Brownian and √ 8/3-LQG surfaces, or any papers about quantum Loewner evolution. The results we do need will be recalled and restated.
Very roughly, the argument in [MS15b,MS16a,MS16b] proceeds as follows. First, [MS15b] uses h to construct a metric D h on D using a growth process called quantum Loewner evolution (QLE) [MS16d]. We will not need to recall the precise construction of D h or the definition of QLE here. Section 2.4 contains all of the background on D h that is necessary to understand this paper. Second, [MS16a] shows that in the case of the √ 8/3-LQG sphere the corresponding metric measure space agrees in law with the Brownian map. Thus, by sampling h (which determines the √ 8/3-LQG sphere), and then generating the corresponding Brownian map, one obtains a coupling of the √ 8/3-LQG sphere and the Brownian map. Finally, [MS16b] uses certain "welding and resampling" arguments to show that in this coupling, the latter object almost surely determines the former. In other words, one can recover the √ 8/3-LQG instance (i.e., the distribution h) as a measurable function of the corresponding Brownian surface instance.
This measurable function is obtained in a remarkably non-explicit way. Although the argument tells us that an instance of a Brownian map (or disk, plane, or half-plane) has a canonical conformal structure, i.e., a canonical embedding into C (or a domain in C) defined up to Möbius transformation, it gives us no way to compute or even approximate that embedding using only information about the "metric measure space" structure of the Brownian surface. Similarly, since a Brownian surface has a canonical conformal structure, we know that Brownian motion on a Brownian surface is well-defined, at least modulo a monotone reparameterization of time. 1 But the argument tells us nothing about how to construct that Brownian motion: in particular, it does not tell us whether Brownian motion on a Brownian surface is the limit of simple random walk on some natural increasingly fine discretizations of that surface, the way Brownian motion on R 2 is the → 0 limit of simple random walk on Z 2 .
The purpose of the present work is to construct the conformal structure of a Brownian surface in an explicit manner. To this end, we will start with an instance of a Brownian surface and then approximate it with the graph of cells associated with the Poisson-Voronoi tessellation. More precisely, we will fix λ > 0 and then pick a Poisson point process with intensity measure given by λ times the area measure on the Brownian surface. The cell corresponding to a given point x in the Poisson point process consists of those points that are at least as close to x as they are to any other point of the Poisson point process (w.r.t. the metric on the Brownian surface). Cells are considered to be adjacent if they have non-empty intersection.
Our main result (Theorem 1.1) says that as λ → ∞, the random walk on the adjacency graph of Voronoi cells converges modulo parameterization to a limiting continuous path. Moreover, this path is a Brownian motion (modulo parameterization) when it is embedded into C using the identification of Brownian surfaces with √ 8/3-LQG. This gives an intrinsic way of constructing Brownian motion (modulo time parameterization) on a Brownian surface without any reference to LQG: indeed, the Brownian motion on the surface is simply the λ → ∞ limit of the random walk on the Voronoi cells.
One can define the Tutte embedding of the adjacency graph of cells in terms of hitting probabilities for the simple random walk (in the same way that Riemann uniformization can be defined using the hitting probabilities for Brownian motion, i.e., harmonic measure). Our results show that this Tutte embedding converges to √ 8/3-LQG in an appropriate sense as λ → ∞ (Theorem 1.2). This in particular gives a new, more explicit, proof of the main result of [MS16b], which says that a √ 8/3-LQG surface is a.s. determined by its metric measure space structure.
Finally, we remark on some related works. Voronoi tessellations of Brownian surfaces have also been considered in other contexts: see, e.g., [Cha16,Gui17]. These papers raise interesting questions about the law of the partitioning of volume among the Voronoi cells, but they do not consider the adjacency graph of cells as we do here.
The first two authors in [GM19b], building on [DDDF19,DFG+19,GM19c,GM19a], recently constructed a metric on a γ -LQG surface for general γ ∈ (0, 2) using a completely different construction from the one in [MS15b,MS16a,MS16b]. It was shown in [GM19b] that the two constructions give the same metric. This paper will make no use of [GM19b], but many of our estimates for the √ 8/3-LQG metric also work for general γ ∈ (0, 2); see the discussion just after Theorem 1.2 and also Remark 2.5.

Main results.
Let (X , D, μ) be an instance of the Brownian map, disk, plane, or half-plane, equipped with its metric D and its natural area measure μ. Conditional on 1 We will not consider the time parameterization in this paper, but we note that there is a canonical way to parameterize Brownian motion on a √ 8/3-LQG surface, called Liouville Brownian motion [Ber15,GRV16]. We expect, but do not prove, that Liouville Brownian motion is the scaling limit of the parameterized walk in the setting of Theorem 1.1; see Problem 5.2.
(X , D, μ), let P λ for λ > 0 be a Poisson point process on X with intensity measure λμ. In the case of the Brownian map or disk (when μ is a finite measure) with area equal to A we can equivalently define P λ as follows: sample N ∼ Poisson(λA), independently from (X , D, μ), and then conditional on N and (X , D, μ) sample N points uniformly and independently from μ.
For z ∈ P λ , let H λ z be the Voronoi cell which is the closure of the set of points in X which are D-closer to z than to any other point of P λ . There is a natural graph structure on P λ whereby z, w ∈ P λ are connected by an edge if and only if the cells H λ z and H λ w intersect along their boundaries. Equivalently, z, w ∈ P λ are joined by an edge if and only if there exists u ∈ X such that D(u, z) = D(u, w) and D(u, x) ≥ D(u, z) for each x ∈ P λ \{z, w}. In the case of the Brownian disk or half-plane, we define ∂P λ to be the set of z ∈ P λ for which the corresponding cell intersects ∂X . Our first main result (from which all of our other main results will follow) says that the simple random walk on P λ converges modulo time parameterization and that the limiting process coincides with Brownian motion under the embedding of X into C which arises from its identification with a √ 8/3-LQG surface. In fact, the convergence is uniform over all choices of starting point in any given compact set K ⊂ X (in the case of the Brownian map or disk, we can just take K = X ). This gives us an explicit, intrinsic definition of Brownian motion on a Brownian surface which coincides with the definition which comes from LQG theory.
For z ∈ X , let Y z,λ be the simple random walk on P λ started from the point of P λ whose corresponding cell contains z (if there is more than one such point, we choose one in an arbitrary manner). We extend Y z,λ from the integers to [0, ∞) by declaring that for each j ∈ N, the path Y z,λ | [ j− 1, j] traverses the D-geodesic from Y z,λ j−1 to Y z,λ j at constant speed. 2 In the case of the Brownian disk or half-plane, we stop Y z,λ at the first time it hits a point of ∂P λ . Theorem 1.1 (Brownian motion on a Brownian surface). If K ⊂ X is a compact set chosen in a manner which is measurable with respect to (X , D, μ), then as λ → ∞ the conditional law of the walk Y z,λ given (X , D, μ) converges in probability as λ → 0, uniformly over all z ∈ K , with respect to the topology on curves viewed modulo time parameterization (see Sect. 2.1.3 for a review of this topology). If we identify X with the Riemann sphere, unit disk, complex plane, or upper half-plane using the correspondence between Brownian and √ 8/3-LQG surfaces, then the limit of the conditional law of Y z,λ given (X , D, μ) is standard planar Brownian motion started from z (and stopped upon hitting the domain boundary in the case of a surface with boundary), viewed modulo time parameterization.
In light of Theorem 1.1, one can define Brownian motion on X to be the limit of the processes Y z,λ as λ → ∞. See Theorem 3.3 for a more general version of Theorem 1.1 which also applies to other √ 8/3-LQG surfaces. Theorem 1.1 allows us to give intrinsic 2 It is easy to see that for any compact set K ⊂ X , the following is true. The supremum over all j such that Y z,λ j ∈ K of the D-diameter and the Euclidean diameter of Y z,λ | [ j −1, j] converges to zero in probability as λ → ∞. This is the only property of the continuous extension of the walk which is needed for our purposes. To see why this property holds, we observe that the D-diameter of Y z,λ | [ j−1, j] equals D(Y z,λ j−1 , Y z,λ j ). Since Y z,λ j−1 and Y z,λ j are contained in adjacent Voronoi cells and the maximum D-diameter of all of the cells which intersect K tends to zero in probability as λ → ∞ (by Lemma A.3), we get the desired statement for the D-diameters. The desired property for Euclidean diameters follows since D induces the same topology as the Euclidean metric.
constructions of other objects on Brownian surfaces which can be derived from random walk. For example, by combining it with the main result of [YY11] we obtain an intrinsic definition of SLE 2 on a Brownian surface as the limit of the loop-erased random walk on Poisson-Voronoi tessellations. Our techniques also imply variants of Theorem 1.1 in which the edges of P λ are assigned random conductances in a sufficiently nice way. For example, Theorem 1.1 is still true if the conductances are i.i.d. and the conductances and their reciprocals have finite expectation. Essentially, this is because the key step in the proof of Theorem 1.1 uses the main result of [GMS18], and the result in [GMS18] allows for variable conductances. We will not state a general variable conductance theorem here, but we remark that if the reader wants to extend Theorem 1.1 to a particular variable conductance model, the key step will be checking that the hypotheses of [GMS18] remain satisfied.
Theorem 1.1 also gives us an explicit construction of the embedding of a Brownian surface into C. More precisely, one can define an explicit way of embedding the graphs P λ into C-called the Tutte embedding-under which the metric and area measure on P λ inherited from (X , D, μ) converge to the √ 8/3-LQG metric and area measure, respectively, as λ → ∞. For concreteness, let us focus on the case when (X , D, μ) is the Brownian disk (we do this since our embedding has a simpler definition when our surface has a boundary).
We will now define the Tutte embedding λ into the closed unit disk D of the Poisson point process P λ (equipped with the graph structure discussed above). Recall that ∂P λ denotes the points in P λ whose corresponding cells intersect ∂X . We first specify the points which will be sent (approximately) to 0 and 1. Let z be sampled uniformly from μ and let z 0 be the a.s. unique (see Lemma A.5) element of P λ such that z ∈ H λ z 0 . Also let x be sampled uniformly from the canonical length measure on ∂X (which can be defined, e.g., by taking a limit of the re-scaled μ-mass of D-neighborhoods of boundary arcs [LG19]) and let x 0 be the a.s. unique element of ∂P λ for which x ∈ H λ x 0 . Enumerate the other points in ∂P λ as x 1 , . . . , x m so that if j, k ∈ {0, . . . , m}, then j < k if and only if we hit the cell H λ x j before the cell H λ x k when we traverse ∂X in the counterclockwise direction started from x.
For j ∈ {0, . . . , m}, we declare that λ (x j ) = e 2πi p j ∈ ∂D, where p j is the probability that a simple random walk on P λ started from z 0 first hits ∂P λ at one of the points x 0 , . . . , x j . Note that this makes it so that the harmonic measure as seen from z 0 approximates the uniform measure on ∂D. This defines an embedding of ∂P λ . We then extend this embedding to be discrete harmonic on the rest of P λ , equivalently we require that the position of each interior vertex of P λ under our embedding is the average of the positions of its neighbors.
Note that under this embedding, the center point of the cell x 0 containing x maps to e 0 = 1. Moreover, the P λ -harmonic measure on ∂P λ as viewed from z 0 approximates the uniform measure on ∂D, the average of which is zero. So, λ (z 0 ) is close to zero (but not necessarily exactly equal to zero) when λ is large. Theorem 1.2 (Tutte embedding convergence). Let (X , D, μ) be a Brownian disk (with any fixed positive choice of area and boundary length) and let h be the GFF-type distribution on D which parameterizes the √ 8/3-LQG surface 3 corresponding to (X , D, μ) under the correspondence of [MS15b,MS16a,MS16b]. Also let μ h and D h be the √ 8/3-Liouville quantum gravity area measure and metric, respectively, induced by h, so that Although the cells differ greatly in Euclidean size, they appear to have moderate "length-to-width ratios" in the sense that a typical cell contains a Euclidean disk of diameter comparable to its own diameter. A mathematical version of this observation (see Propositions 4.4 and 4.5) plays a role in the proof that simple random walk on the adjacency graph of cells approximates Brownian motion (X , D, μ) = (D, D h , μ h ) as metric measure spaces. If we identify P λ with its image under the Tutte embedding, then we have the following convergence in probability as λ → ∞.
• The measure which assigns to each vertex of P λ a mass equal to the μ-mass of its corresponding cell converges to μ h . The same is true of the counting measure on vertices of P λ , scaled by the factor λ −1 .
• The maximum over all pairs of embedded vertices z, w ∈ P λ of the quantity |D h (z, w) − D(z, w)| converges to zero.
• The simple random walk on vertices of P λ started from z 0 converges modulo time parameterization to Brownian motion started from 0 and stopped upon hitting ∂D in the quenched sense (i.e., its conditional law given P λ and (X , D, μ) converges weakly in probability as λ → ∞).
As we will explain in Sect. 3.3, Theorem 1.2 is a consequence of Theorem 1.1. Indeed, to prove Theorem 1.2 we just need to show that when λ is large, the Tutte embedding of P λ is close to the image of P λ under the a priori embedding of (X , D, μ) which comes from its identification with (D, D h , μ h ). Due to the manner in which the Tutte embedding is defined, this, in turn, follows from the statement that under the a priori embedding, the random walk on P λ converges to Brownian motion modulo time parameterization, as asserted in Theorem 1.1.
We now briefly outline the proof of Theorem 1.1. We want to show that under the a priori embedding, the random walk on P λ converges to Brownian motion modulo time parameterization. This is a random walk in random environment (RWRE) problem: we have a random walk on P λ -viewed as a random graph drawn in D-and we want to show that it approximates Brownian motion. However, this problem falls outside of the usual RWRE or random conductance model framework (as surveyed, e.g., in [BAF16,Bis11]) because the Euclidean sizes of the Voronoi cells under the a priori embedding vary dramatically from one location to another, so the environment is highly spatially inhomogeneous (see Fig. 1) and in particular its law is not stationary with respect to spatial translations.
Nevertheless, as explained in Sect. 3, if we consider a certain special √ 8/3-LQG surface called a 0-quantum cone (which does not correspond to one of the "standard" Brownian surfaces) then the associated adjacency graph of Poisson-Voronoi cells is in a certain sense "translation invariant modulo a global rescaling." The paper [GMS18] gives conditions under which random walk converges to Brownian motion modulo time parameterization in a random environment which is only required to satisfy this weaker form of translation invariance. We re-state the particular theorem from [GMS18] which we will use as Theorem 2.2. Once certain properties of our Voronoi cells have been established, this theorem shows that random walk on the Poisson-Voronoi tessellation of the 0-quantum cone converges to Brownian motion modulo time parameterization. One can then transfer to other √ 8/3-LQG surfaces (such as the ones corresponding to the Brownian map, disk, plane, and half-plane) via local absolute continuity considerations.
The key quantitative condition needed to apply the above RWRE theorem is that the quantity diam(H 0 ) 2 deg(H 0 )/area(H 0 ) has finite expectation, where H 0 is the Voronoi cell containing the origin for the γ -quantum cone with λ = 1 and diam, deg, and area denote its Euclidean diameter, degree (in the adjacency graph of Voronoi cells), and Lebesgue measure, respectively. We will show in Sect. 3.4, using a mass-transport principle, that in fact it suffices to prove that where B H is the smallest D h -metric ball centered at the center point of H (i.e., the point of P which is in H ) which contains H . This will be important for our purposes since it is easier to lower-bound the Lebesgue measure of an LQG metric ball than a Voronoi cell. In order to verify (1.2), we need to establish a number of estimates for √ 8/3-LQG metric balls and Voronoi cells which are of independent interest (see Sect. 4).
In particular, we show in Proposition 4.4 that a √ 8/3-LQG metric ball is extremely unlikely to be "long and skinny" in the sense that it typically contains a Euclidean ball of radius comparable to its Euclidean diameter. This is done using a percolation argument for the GFF, similar to ones in [DD19,DG16,DZZ18,DG18,DF18,DD18,DDDF19]. We also show in Proposition 4.8 that the √ 8/3-LQG mass of an LQG metric ball is highly concentrated around the fourth power of its √ 8/3-LQG radius. This is done by starting with estimates for the Brownian map [Le10], then using the local independence properties of the GFF to establish a suitable concentration bound.
The high-level strategy used in this paper (especially, the application of [GMS18]) is similar to the strategy used in [GMS17] to prove the convergence to LQG of the Tutte embedding of the so-called mated-CRT map. The mated-CRT map is a random planar map built by mating a pair of continuum random trees, which has an a priori embedding into C due to the results of [DMS14]. However, the proof of the needed bound for diam(H 0 ) 2 deg(H 0 )/area(H 0 ) in this paper, including the reduction to (1.2) and the estimates for √ 8/3-LQG metric balls, is very different from the proof of the analogous bound in [GMS17]. The proof of this estimate comprises most of the technical work in this paper. All of the arguments in this paper carry over verbatim to the γ -LQG metric for general γ ∈ (0, 2), as defined in [GM19b], except for the proof of the ball volume concentration bound in Proposition 4.8 (which uses estimates for the Brownian map, so only works for γ = √ 8/3). If we had an analog of Proposition 4.8 for general γ ∈ (0, 2), we could immediately extend our results to Poisson-Voronoi tessellations of γ -LQG surfaces for all γ ∈ (0, 2). Remark 1.3 (Embeddings of random planar maps). Theorem 1.2 implies a scaling limit result for certain "coarse-grained" embeddings of random planar maps toward √ 8/3-LQG, as we now explain. Suppose {M n } n∈N is a sequence of random planar maps with boundary which converge in law to the Brownian disk in the following sense. There are scaling constants a n , b n , c n > 0 such that if we view the planar maps as curvedecorated metric measure spaces equipped with a −1 n times the graph distance, b −1 n times the counting measure on vertices, and the path which traces the boundary according to the natural ordering in such a way that each edge is traversed in c −1 n units of time, then the maps converge in law w.r.t. the Gromov-Hausdorff-Prokhorov-uniform (GHPU) topology, the analog of the Gromov-Hausdorff topology for curve-decorated metric measure spaces introduced in [GM17b].
For λ > 0, we can define a Poisson-Voronoi tessellation of M n using a Poisson point process with respect to λb −1 n times the counting measure on vertices of M n . We can then define a "coarse-grained" Tutte embedding of M n using this Poisson-Voronoi tessellation in exactly the same manner as in Theorem 1.2. It can be seen from the GHPU convergence of M n to the Brownian disk that for each fixed λ > 0, the adjacency graph of Voronoi cells on M n converges in the total variation sense to the adjacency graph P λ defined above (here we emphasize that for fixed λ the typical number of Voronoi cells is a tight random variable as n → ∞). Hence if we send λ n → ∞ sufficiently slowly as n → ∞, we get an analog of Theorem 1.2 for the λ n -coarse-grained Tutte embedding of M n .
More details regarding the above appeared in an earlier arXiv version of this paper, but were cut from the current version for brevity.

Outline.
The remainder of this article is structured as follows. In Sect. 2, we will fix some notation, state the scaling limit result from [GMS18] which is used in the proof of our main results, and recall some facts about metric spaces and √ 8/3-LQG surfaces. In Sect. 3, we prove Theorems 1.1 and 1.2 assuming that a certain moment bound for Voronoi cells is satisfied. In Sect. 4 we will prove the required moment bound, along with a number of estimates for  For z ∈ C and r > 0, we write B r (0) for the open Euclidean ball of radius r centered at z. We write diam(·) and area(·) for Euclidean diameter and Lebesgue measure on C, respectively.
For a metric space (X, d), x ∈ X , and r > 0, we write B r (x; d) for the open d-metric ball of radius r centered at x.
If We similarly define events which occur with polynomially, superpolynomially, and exponentially high probability as a parameter tends to ∞. We will often specify any requirements on the dependencies on rates of convergence in O(·) and o(·) errors, implicit constants in , etc., in the statements of lemmas/ propositions/theorems, in which case we implicitly require that errors, implicit constants, etc., appearing in the proof satisfy the same dependencies.
where the supremum is over all partitions P : Note that the d-length of a curve may be infinite. A d-geodesic between two points x, y ∈ X is a path from x to y of minimal d-length.
is the infimum of the d-lengths of paths joining z and w.

Metric on curves modulo time parameterization.
Our scaling limit results for random walk on embedded planar maps are with respect to the topology on curves modulo time parameterization, which we now recall. If β 1 : [0, T β 1 ] → C and β 2 : [0, T β 2 ] → C are continuous curves defined on possibly different time intervals, we set where the infimum is over all increasing homeomorphisms φ : [0, T β 1 ] → [0, T β 2 ] (the CMP stands for "curves modulo parameterization"). It is shown in [AB99, Lemma 2.1] that d CMP induces a complete metric on the set of curves viewed modulo time parameterization.
In the case of curves defined for infinite time, it is convenient to have a local variant of the metric d CMP . Suppose β 1 : [0, ∞) → C and β 2 : [0, ∞) → C are two such curves. For r > 0, let T 1,r (resp. T 2,r ) be the first exit time of β 1 (resp. β 2 ) from the ball B r (0) (or 0 if the curve starts outside B r (0)). We define so that d CMP loc (β n , β) → 0 if and only if for Lebesgue a.e. r > 0, β n stopped at its first exit time from B r (0) converges to β stopped at its first exit time from B r (0) with respect to the metric (2.2). Note that the definition (2.2) of d CMP β 1 | [0,T 1,r ] , β 2 | [0,T 2,r ] makes sense even if one or both of T 1,r or T 2,r is infinite, provided we allow d CMP β 1 | [0,T 1,r ] , β 2 | [0,T 2,r ] = ∞ (this does not pose a problem due to the definition of the integrand in (2.3)).
If (X, d, x 0 ) is a general metric space with a marked point, one can similarly define the metric on curves modulo time parameterization on X but with d-distances in place of Euclidean distances and d-metric balls centered at x 0 in place of Euclidean balls centered at 0.

Scaling limit for random walk on graphs of cells.
In this subsection we state a version of the main result of [GMS18] which gives general conditions under which random walk on the adjacency graph of a random collection of cells (e.g., Voronoi cells) on C converges to Brownian motion. This same result is also used in [GMS17] to prove an embedding convergence result for a different discretization of LQG. Let us first describe what we mean by an "adjacency graph of cells". Definition 2.1. A cell configuration on C consists of the following objects.
1. A locally finite collection H of compact connected subsets of C ("cells") with nonempty interiors whose union is all of C and such that the intersection of any two elements of H has zero Lebesgue measure. 2. A symmetric relation ∼ on H×H ("adjacency") such that if H ∼ H , then H ∩ H = ∅ and H = H .
We will typically slightly abuse notation by making the relation ∼ implicit, so we write H instead of (H, ∼). We view H as a weighted graph whose vertices are the cells of H and whose edge set is In [GMS18], one also allows for a conductance function on the edges of H. Here we will only consider cell configurations with unit conductances. We note that the intersections of the cells of H are required to have zero Lebesgue measure. We will check this condition for Voronoi cells in Lemma A.6. We define a metric on the space of cell configurations by where each of the infima is over all homeomorphisms f r : C → C such that f r takes each cell in H which intersects B r (0) to a cell in H which intersects B r (0) and preserves the adjacency relation between these cells, and f −1 r does the same with H and H reversed.
In [GMS18], we proved that the simple random walk on a random cell configuration H which satisfies the following hypotheses converges to Brownian motion. Here, for C > 0 and z ∈ C we write C(H − z) for the cell configuration obtained by translating all of the cells by −z then scaling all of the cells by C.
1. Translation invariance modulo scaling. There is a (possibly random and Hdependent) increasing sequence of open sets U j ⊂ C, each of which is either a square or a disk, whose union is all of C such that the following is true. Conditional on H and U j , let z j for j ∈ N be sampled uniformly from Lebesgue measure on U j . Then the shifted cell configurations H − z j converge in law to H modulo scaling as j → ∞, i.e., there are random numbers C j > 0 (possibly depending on H and z j ) such that C j (H − z j ) → H in law with respect to the metric (2.5).

Ergodicity modulo scaling. Every real-valued measurable function F = F(H)
which is invariant under translation and scaling, i.e., F(C(H − z)) = F(H) for each z ∈ C and C > 0, is a.s. equal to a deterministic constant. 3. Finite expectation. With H 0 the cell in H containing 0, where diam, area, and deg denote Euclidean diameter, Lebesgue measure, and vertex degree in H, respectively. 4. Connectedness along lines. Almost surely, for each horizontal or vertical line segment L ⊂ C, the subgraph of H induced by the set of cells which intersect L is connected.
The combination of hypotheses 1 and 2 is referred to as ergodicity modulo scaling in [GMS18]. Several equivalent formulations of hypothesis 1 are given in [GMS18, Definition 1.2] (we will use a different formulation, in terms of a "mass transport principle" in Sect. 3.4). The version of hypothesis 3 given here is slightly simpler than the version in [GMS18] since we are assuming unit conductances. Hypothesis 4 is automatically satisfied if any two cells which intersect are considered to be adjacent. This will always be the case for cell configurations considered in this paper. The following is [GMS18, Theorem 3.10].
Theorem 2.2 [GMS18]. Let H be a random cell configuration satisfying the above four hypotheses. For z ∈ C, let Y z denote the simple random walk on H started from H z (with conductances c). For j ∈ N 0 , let Y z j be an arbitrarily chosen point of the cell Y z j and extend Y z from N 0 to [0, ∞) by piecewise linear interpolation. There is a deterministic covariance matrix with det = 0 such that the following is true. For each fixed compact set A ⊂ C, it is a.s. the case that as → 0, the maximum over all z ∈ A of the Prokhorov distance between the conditional law of Y z/ given H and the law of Brownian motion started from z with covariance matrix , with respect to the topology on curves modulo time parameterization (as defined in Sect. 2.1.3), tends to 0.

Liouville quantum gravity surfaces.
Fix γ ∈ (0, 2) (in fact, we will always take γ = √ 8/3). For k ∈ N 0 , a γ -Liouville quantum gravity surface with k marked points is h is a distribution on U (which we will always take to be a realization of some variant of the GFF on U ), and z 1 , . . . , We think of two equivalent (k + 2)-tuples as above as corresponding to different parameterizations of the same surface. We refer to the distribution h corresponding to an LQG surface as the embedding of the surface. The above definitions first appeared in [DS11], and also play an important role, e.g., in [She16,DMS14].
If the law of the field h is locally absolutely continuous with respect to the law of the Gaussian free field on U , then we can define the γ -LQG area measure μ h on U , which is the a.s. limit of regularized versions of e γ h(z) dz as well as the γ -LQG boundary length measure ν h on ∂U (in the case when U has a boundary). There are several equivalent ways to construct these measures: see, e.g., [Kah85,DS11,RV14]. By [DS11, Proposition 2.1], if h and h are related by a conformal map as in (2.7), then This means that μ h and ν h can be viewed as measures on the LQG surface.
In the special case when γ = √ 8/3, an LQG surface also admits a metric D h , as shown in [MS15b,MS16a,MS16b], and this metric is compatible with coordinate changes of the form (2.7). We will review the basic properties of this metric in Sect. 2.4.

Quantum cones
The LQG surface which we will work with must frequently is the α-quantum cone for α ∈ (−∞, Q), which is defined in [DMS14, Definition 4.10]. The α-quantum cone is a doubly marked surface (C, h, 0, ∞) whose γ -LQG measure μ h has infinite total mass, but assigns finite mass to every bounded subset of C. One way to obtain an α-quantum cone is to start with a whole-plane GFF plus α log(1/| · |) then "zoom in" near the origin and re-scale (i.e., add a constant to the field) so that the LQG area of a fixed set remains of constant order. See [DMS14, Proposition 4.13(i)] for a precise statement.
Quantum cones with parameter α ∈ {0, √ 8/3} are especially natural, and these will be the main types of quantum cones which we will consider. The case α = √ 8/3 is special since a GFF a.s. has a − √ 8/3-log singularity at a point sampled from its √ 8/3-LQG measure [DS11, Section 3.3], so this surface can be thought of as describing the behavior of a general √ 8/3-LQG surface near such a point. Moreover, the √ 8/3quantum cone, equipped with its √ 8/3-LQG metric and area measure, is equivalent to the Brownian plane as defined in [CL14] (see Sect. 2.4). In a similar vein, the 0quantum cone describes the local behavior of a √ 8/3-LQG surface near a Lebesgue typical point. The 0-quantum cone can be used to construct cell configurations which satisfy the translation invariance modulo scaling condition of Theorem 2.2, which says that the origin is in some sense "Lebesgue typical".
We will need some properties of quantum cones which follow from the definition in [DMS14, Definition 4.10], so we now recall this definition. Let α < Q and let A : R → R be the process such that A t = B t + αt for t ≥ 0, where B is a standard linear Brownian motion; and for t < 0, let A t = B −t + αt, where B is a standard linear Brownian motion conditioned so that B t + (Q − α)t > 0 for all t > 0 and taken to be independent of B. We define h to be the random distribution such that if h r (0) denotes the average of h on ∂ B r (0) (see [DS11, Section 3.1] for the definition and basic properties of the circle average), then t → h e −t (0) has the same law as the process A; and h − h |·| (0) is independent from h |·| (0) and has the same law as the analogous process for a whole-plane GFF.
Since a quantum cone has only two marked points, one can get a different choice of h corresponding to the same LQG surface (i.e., a different embedding of the quantum cone) by re-scaling space and applying the LQG coordinate change formula (2.7). We will almost always consider the particular choice of distribution h defined just above. This choice of h is called the circle average embedding, and is characterized by the fact that 1 = sup{r > 0 : h r (0) + Q log r = 0}. If h is the circle-average embedding of an α-quantum cone, then h| D agrees in law with the corresponding restriction of a whole-plane GFF plus −α log | · |, normalized so that its average over ∂D is 0.
The α-quantum cone possesses a certain special scale invariance property, which we now describe. Let α < Q, let h be the circle-average embedding of an α-quantum cone, and let {h r (z) : r > 0, z ∈ C} be its circle average process. We define where here Q is as in (2.7). That is, R b gives the largest radius r > 0 so that if we scale spatially by the factor r and apply the change of coordinates formula (2.7), then the average of the resulting field on ∂D is equal to 1 √ 8/3 log b. Note that R 0 = 0 in the case of the circle average embedding. It is easy to see from the above definition of h (and is shown in [DMS14, Proposition 4.13(i)]) that for each fixed b > 0, It is immediate from the definitions of μ h and the √ 8/3-LQG metric D h that adding 1 √ 8/3 log b to the field scales √ 8/3-LQG areas by b and √ 8/3-LQG distances by b 1/4 (in the case of the metric, see [MS16a, Lemma 2.2] or Lemma 2.3). By the √ 8/3-LQG coordinate change formulas for μ h and D h , we therefore see that (2.9) implies that where here we mean equality in law as metric measure spaces. We note that the property (2.10) is not true with, say, a whole-plane GFF in place of a quantum cone. This property is a major reason for considering quantum cones.

Quantum disks, spheres, and wedges.
We will also have occasion to consider other special quantum surfaces besides just quantum cones. We will not need as many properties of these, so we just briefly mention their definitions and refer to the cited references for more details.
A quantum disk is a quantum surface (D, h) defined in [DMS14,Definition 4.21] which behaves locally like a free-boundary GFF on D, but is defined in a slightly different way. Its √ 8/3-LQG area measure and boundary length measure each have finite total mass, and one can consider quantum disks with specified boundary length (and random area) or with specified boundary length and area. One can also define a quantum disk with any number of marked boundary points and interior points sampled uniformly from its √ 8/3-LQG boundary length and area measures, respectively. A quantum sphere is a quantum surface (C, h) introduced in [DMS14, Definition 4.21] with μ h (C) < ∞. Typically one considers a unit-area quantum sphere, which means we fix μ h (C) = 1. Quantum spheres with other areas are obtained by re-scaling (equivalently, adding a constant to h). As in the case of the quantum disk, one can consider quantum spheres with one or more marked points sampled uniformly from the √ 8/3-LQG area measure. For α ≤ Q, an α-quantum wedge is a quantum surface (H, h, 0, ∞) defined in [DMS14,Definition 4.5] which has finite mass in every neighborhood of 0 but infinite total mass. It is the half-plane analog of the α-quantum cone considered above and satisfies the same scaling property (2.10) as the α-quantum cone.

The
√ 8/3-Liouville quantum gravity metric. Suppose that U ⊂ C is a connected open set and h is a random distribution on U which is locally absolutely continuous with respect to the GFF on U , in the sense that for every z ∈ U , there is a neighborhood V of z such that the law of h| V is absolutely continuous with respect to the corresponding restriction of the GFF. The papers [MS15b,MS16a,MS16b] show that one can define a √ 8/3-LQG metric D h associated with h. We will not need the precise definition of this metric here. Rather, we will only use a small number of basic properties of D h , which we now record.
1. Bi-Hölder with respect to Euclidean metric. The identity map from U , equipped with the Euclidean metric, to (U, D h ) and its inverse are each a.s. locally Hölder continuous with a (non-explicit) Hölder exponent. In particular, D h induces the same topology on U as the Euclidean metric. 2. Existence of geodesics. Almost surely, for any z, there is a D h -geodesic from z to w, i.e., a path from z to w of minimal D h -length. If the law of h is absolutely continuous with respect to that of a free-boundary GFF in a neighborhood of every point of ∂U , then D h extends to a metric on U and a.s. for each z, w ∈ U , there is a D h -geodesic in U from z to w. 3. LQG coordinate change formula. If U, U ⊂ C and φ : U → U is a conformal map, then In particular, metric balls are locally determined by h. We will also use the fact that certain special LQG surfaces (U, h) are equivalent to Brownian surfaces, in the sense that the metric measure space (U , D h , μ h ) agrees in law with a Brownian surface. 4 See [MS16a, Corollary 1.5] for the sphere, disk, and plane cases and [GM17b, Proposition 1.10] for the half-plane case.
• The quantum sphere is equivalent to the Brownian map.
• The quantum disk is equivalent to the Brownian disk.
• The √ 8/3-quantum cone is equivalent to the Brownian plane. • The √ 8/3-quantum wedge is equivalent to the Brownian half-plane.
We will now explain another elementary property of D h which allows us to define D h+ f whenever f : U → R is a random continuous function coupled with h, even if the law of h + f is not locally absolutely continuous with respect to the GFF.

Lemma 2.3. Suppose h is a random distribution on a connected open set U ⊂ C and let f : U → R be a random continuous function (not necessarily independent from h).
If the laws of h and h + f are both locally absolutely continuous with respect to the GFF on U , then a.s.
In fact, it is a.s. the case that for each z, w ∈ U , where the infimum is over all simple paths from z to w parameterized by their D h -length.
Proof. If f ≡ c is constant, then by [MS16a, Lemma 2.2], one has D h+ f = e c/ √ 6 D h . In fact, the proof of [MS16a, Lemma 2.2] shows that a.s. (2.12) holds. We will now deduce (2.13) from (2.12). To this end, fix > 0. Since f is continuous, we can find a (possibly random) Now fix z, w ∈ U and let γ : [0, T ] → U be a path from z to w whose D h+ f -length is at most D h+ f (z, w) + . Choose finitely many times 0 = t 0 < t 1 < · · · < t n = T such that max j∈ [1,n] in place of U ) and our choice of δ, we have (in the notation (2.1)) This shows that the D h -length of γ is finite and, if γ is parameterized by D h -length, that the D h+ f -length of γ and the integral len(γ ;D h ) 0 e f (γ (t))/ √ 6 dt differ by a factor of at most e / √ 6 . Sending → 0 shows that D h+ f (z, w) is at least the right side of (2.13). We similarly get the reverse inequality.
Let h be a random distribution on a connected open set U ⊂ C whose law is locally absolutely continuous with respect to the GFF on U and let f : U → R be a random continuous function (not necessarily independent from h). If D h is defined, we define D h+ f by the formula (2.13). We need to make sure that D h+ f is well-defined (i.e., we get the same metric if we make a different choice of h and f with h + f = h + f ) and that it is a measurable function of h + f (a priori we only know that D h+ f is a measurable function of (h, f )). The following lemma is an easy consequence of Lemma 2.3. We will give the proof just below.

Lemma 2.4.
Let h and f be as above and define D h+ f by the formula (2.13). Then D h+ f is a.s. determined by h + f . Moreover, if (h , f ) is another pair consisting of a random distribution on a connected open set U ⊂ C whose law is locally absolutely continuous with respect to the GFF on U and a random continuous function such that Once Lemma 2.4 is established, it follows from the above properties of D h that D h+ f is locally bi-Hölder continuous with respect to the Euclidean metric in the sense of property 1 above (so induces the same topology on U as the Euclidean metric) and satisfies the LQG coordinate change formula (2.7) and the locality property 4. We also note that Lemma 2.4 implies that for a given choice of h, a.s. D h+ f can be defined via the formula (2.12) simultaneously for every choice of continuous function f : U → R. (2.14) Indeed, this follows since the two metrics in (2.14) induce the same length measure on each path in U .
Suppose now that we are given two couplings (h, f ) and (h , f ) of a distribution whose law is locally absolutely continuous w.r.t. the GFF and a random continuous and (h , f ) are conditionally independent given h + f . By (2.14) applied to the functions f and f − f together with Lemma 2.3 applied to the distributions h and h = h + f − f , we get that a.s.
Moreover, since D h+ f and D h + f are conditionally independent given h + f (by our choice of coupling), it follows that D h+ f is a.s. determined by h + f . Remark 2.5. All of the properties of D h discussed in this section also hold for the γ -LQG metric for general γ ∈ (0, 2) from [GM19b], except that 1/ √ 6 is replaced by γ /d γ , where d γ is the Hausdorff dimension of the γ -LQG metric (which is not known explicitly). In fact, [GM19b] shows that a list of properties similar to the ones discussed in this section uniquely characterize the γ -LQG metric up to a deterministic multiplicative constant.

Proof of Main Results, Assuming Finite Expectation Hypothesis
(3.1) Extending the notation above, for w ∈ D we write H λ h,w for the (a.s. unique for w deterministic, by Lemma A.6) Voronoi cell which contains w. We define We will often omit the subscript h and/or the superscript λ when these objects are clear from the context. In this subsection, we will prove all of our main results conditional on the following proposition.
h as above with λ = 1 and let H 0 be the cell which contains the origin. Then where here deg(H 0 ) denotes the degree of H 0 as a vertex of H.
Proposition 3.1 is used in Sect. 3.1 to check the finite expectation hypotheses of Theorem 2.2 for the Voronoi cell configuration associated with 0-quantum cone. The proof of Proposition 3.1 is given in Sect. 4. This proof is the most difficult step in the proofs of our main results, and requires us to establish several estimates for √ 8/3-LQG metric balls which are of independent interest.
The rest of this section is structured as follows. In Sect. 3.1, we explain why Proposition 3.1 together with Theorem 2.2 implies a scaling limit result for random walk on the adjacency graph of Voronoi cells associated with a 0-quantum cone. In Sect. 3.2, we transfer this result to random walk on Voronoi cells on other types of √ 8/3-LQG surfaces (including the ones corresponding to the Brownian map, disk, plane, and halfplane) using local absolute continuity and thereby prove Theorem 1.1. In Sect. 3.3, we deduce Theorem 1.2 from our scaling limit result for random walk. The arguments in these three subsections are similar to the analogous arguments in [GMS17, Section 3]. In Sect. 3.4, we give a re-formulation of Proposition 3.1 which involves bounds for √ 8/3-LQG metric balls instead of Voronoi cells, and which turns out to be easier to prove than Proposition 3.1 itself.
Throughout this section, we will use several elementary properties of Voronoi cells whose proofs are collected in "Appendix A" to avoid interrupting the main argument.

3.1.
Cell configuration corresponding to a 0-quantum cone. Let (C, h, 0, ∞) be a 0quantum cone and write P = P 1 h and H = H 1 h for its associated Poisson point process and collection of Voronoi cells with λ = 1. Proposition 3.2. The conclusion of Theorem 2.2 holds for the cell configuration H above. Moreover, the covariance matrix of the limiting Brownian motion is a positive scalar multiple of the identity matrix. Proof. By Lemma A.4, the cells of H are a.s. compact with non-empty interior and H is locally finite. By Lemma A.6, a.s. the intersection of any two cells of H has zero Lebesgue measure, and by definition any two cells which are adjacent in H intersect. Therefore H satisfies the conditions of Definition 2.1. We will now check the conditions of Theorem 2.2.
Translation invariance modulo scaling. For j ∈ N, let R j be the largest r > 0 for which h r (0) + Q log r = γ −1 log j, where h r (0) denotes the circle average, as in (2.8). We will check the needed resampling property for U j = B R j (0). By (2.9), the field h j := h(R j ·) + Q log R j − γ −1 log j agrees in law with h. In particular, by the discussion just after [DMS14, Definition 4.10], h j | D agrees in law with the corresponding restriction of a whole-plane GFF, normalized so that its circle average over ∂D is 0. Consequently, if we sample z j uniformly from Lebesgue measure on B R j (0), then the proof of [DMS14, Proposition 4.13(ii)] along with the translation invariance of the law of the whole-plane GFF, modulo additive constant, shows that the there is a sequence of random constants C j → ∞ such that the law of h(C j (· − z j )) + Q log C j , restricted to any compact subset K ⊂ C, converges to the law of h| K in the total variation sense as j → ∞. By the LQG coordinate change formula, this implies that the joint law Ergodicity modulo scaling. It is easily checked that R>0 σ h| C\B R (0) is the trivial σ -algebra (see, e.g., [HS18, Lemma 2.2] for the case of the whole-plane GFF; the case of h can be treated in an identical manner due to [DMS14, Definition 4.10]). From this, it follows that also the intersection over all R > 0 of the σ -algebra generated by h| C\B R (0) and the set of points of P which are contained in To deduce condition 2 from this, consider a real-valued function In general, the conditional law of F given H(B R (0)) must be deterministic by the preceding sentence, so F is independent from H(B R (0)), whence the above claim implies that F is equal to a deterministic constant a.s.
Finite expectation. This is the content of Proposition 3.1, which will be proven in Sect. 4.

Connectedness along lines.
This follows since by definition two cells of H are connected by an edge of EH if and only if they intersect and the collection of cells is locally finite. The covariance matrix is a scalar multiple of the identity since the law of h, and therefore the law of H, is invariant under rotations around the origin.

Random walk on cells converges to Brownian motion.
The following theorem is a generalization of Theorem 1.1 (recall the correspondence between Brownian and √ 8/3-LQG surfaces as described in Sect. 2.4).
Theorem 3.3. Suppose that we are in one of the following situations.
is a doubly marked quantum sphere, e.g., with fixed area. For each deterministic compact set K ⊂ D, the supremum over all z ∈ K of the Prokhorov distance between the conditional law of Y z,λ given (h, P λ h ) and the law of a standard two-dimensional Brownian motion started from z (and stopped when it hits the boundary in the case of a quantum wedge or quantum disk), with respect to the metric on curves viewed modulo time parameterization (i.e., the metric (2.2) in the disk or half-plane case or the metric (2.3) in sphere or whole-plane case) converges to 0 in probability as λ → ∞.
We note that in Theorem 3.3, the walk is extended by piecewise linear interpolation whereas in Theorem 1.1 it follows D-geodesics between the points of P λ . This does not affect the conclusion of the theorem: indeed, by Lemma A.3 and the fact that D h induces the Euclidean topology, for any fixed compact set K ⊂ D, the maximum over all adjacent pairs of vertices z, w ∈ P λ h ∩ K of the Euclidean diameter of every D h -geodesic from z to w tends to zero in law as λ → ∞. The same is true with D h -diameters in place of Euclidean diameters and/or line segments in place of D h -geodesics.
We first prove Theorem 3.3 in the case of the 0-quantum cone, using Proposition 3.2. This is the step in the proof where we go from a.s. converges to convergence in probability.
Lemma 3.4. Theorem 3.3 is true with a 0-quantum cone in place of a γ -quantum cone.
Proof. By Brownian scaling the statement of the lemma is invariant under the operation of changing the embedding h (i.e., replacing h by h(r ·) + Q log r for some possibly random r > 0), so we can assume without loss of generality that h has the circleaverage embedding, as described in Sect. 2.3 and [DMS14, Definition 4.10] (we could also, e.g., embed so that μ h (D) = 1).
Proposition 3.2 together with Theorem 2.2 tells us that a.s. the conditional law given H 1 h of the random walk on H 1 h converges in law as → 0 to standard two-dimensional Brownian motion modulo time parameterization, and the convergence is uniform over all starting points in any fixed compact subset of C. For Since R λ → ∞ as λ → ∞, we now get the desired convergence in probability from Proposition 3.2 and Theorem 2.2.
Using local absolute continuity, we can now transfer to other quantum surfaces, starting with the case of quantum cones.
Lemma 3.5. Theorem 3.3 is true in the case of the α-quantum cone for α < Q.
Proof. As in the proof of Lemma 3.4, we work with the circle-average embedding of the α-quantum cone, which has the property that h| D agrees in law with the corresponding restriction of a whole-plane GFF plus −α log | · |, normalized so that its circle average over ∂D is 0. We also let h be the circle-average embedding of a 0-quantum cone in The statement of the lemma is essentially a consequence of Lemma 3.4 and local absolute continuity (in the form of [MS16c, Proposition 3.4]), but a little care is needed since we only have local absolute continuity between the laws of a h and h on domains at positive distance from 0 (due to the α-log singularity of h) and from ∂D (due to our choice of embedding). Throughout the proof, the Prokhorov distance is always taken with respect to the metric on curves viewed modulo time parameterization.
For ρ > 0 and z ∈ B ρ (0), let J z,λ ρ for n ∈ N be the exit time from B ρ (0) of the embedded walk Y z,λ on H λ h . Also let B z be a standard two-dimensional Brownian motion started from z and let τ z ρ be its exit time from B ρ (0). We need to show that for each ρ > 0, the supremum over all z ∈ B ρ (0) of the Prokhorov distance between the conditional laws of Y z,n | [0,J z,n ρ ] and B z | [0,τ z ρ ] given h, P λ h converges to zero in probability as λ → ∞.
We first consider a radius ρ ∈ (0, 1) and deal with the log singularity at 0. For δ ∈ (0, ρ), choose ζ = ζ(δ) ∈ (0, δ) such that the probability that a Brownian motion started from any point of D\B δ (0) hits B ζ (0) before leaving D is at most δ. By Lemma 3.4 and local absolute continuity it holds with probability tending to 1 as λ → ∞ that for each z ∈ B ρ (0)\B δ (0), the Prokhorov distance between the conditional laws of Y z,λ | [0,J z,λ ρ ] and B z | [0,τ z ρ ] given h is at most δ. Since the law of B z | [0,τ z ρ ] depends continuously on z, the Prokhorov-distance diameter of the set of laws of the curves B z | [0,τ z ρ ] for z ∈ B δ (0) tends to 0 as δ → 0.
By the last two sentences of the preceding paragraph and the strong Markov property of Y z,λ and of B z , it holds with probability tending to 1 as λ → ∞ that for each z ∈ B δ (0), the Prokhorov distance between the conditional laws of Y z,λ | [J z,λ , at a deterministic rate depending only on ρ. The distance between the curves Y z,λ | [J z,λ δ ,J z,λ ρ ] and Y z,λ | [0,J z,λ ρ ] , viewed modulo time parameterization, is at most 2δ. Sending δ → 0 now gives the theorem statement in the case ρ < 1.
The case when ρ ≥ 1 follows from the case when ρ ∈ (0, 1) and the scale invariance property of the α-quantum cone [DMS14, Proposition 4.13(i)], applied similarly as in Proposition 3.2.
Lemma 3.6. Theorem 3.3 is true in the case of the quantum sphere.
Proof. This is immediate from Lemma 3.5 and local absolute continuity.
Lemma 3.7. Theorem 3.3 is true in the case of the quantum disk.
Proof. Let (C, h, 0, ∞) be a doubly marked quantum sphere conditioned on the event that the D h -distance from 0 to ∞ is at least 1 and let U be the connected component of C\B 1 (∞; D h ) which contains 0. Then the conditional law of the quantum surface (U, h| U , 0) given ν h (∂U ) is that of a quantum disk with one marked point in its interior, with given boundary length (this follows, e.g., from the construction of D h using QLE in [MS15b]).
If we let P λ h be a Poisson point process with intensity measure λμ h , then P λ h ∩ U is a Poisson point process on U with intensity measure λμ h| U . Let H λ h| U be the configuration of Voronoi cells defined using the set of points P λ h ∩ U and the metric D h| U . Then each cell of H λ h|U which does not intersect ∂U is identical to the corresponding cell of H λ h with the same center point. It therefore follows from Lemma 3.6 that the maximum over all z ∈ U of Prokhorov distance between the following two laws, with respect to the topology on curves viewed modulo time parameterization, tends to 0 as λ → ∞: • The conditional law given (h, H λ h ) of the random walk on H λ h| U stopped upon hitting a cell which intersects ∂U , embedded into U and linearly interpolated as in Theorem 3.3.
• The law of Brownian motion a started from 0 and stopped upon hitting ∂U .
By the conformal invariance of Brownian motion and the first paragraph, this gives the statement of the lemma for the quantum disk with random boundary length ν h (∂U ). By scale invariance, this implies the statement of Theorem 3.3 for a doubly marked quantum disk with any fixed boundary length. By conditioning on the area of such a quantum disk, we also get the statement for a quantum disk with fixed area and boundary length.
Proof of Theorem 3.3. Lemmas 3.5, 3.6, and 3.7 give the theorem statement in the quantum cone, quantum sphere, and quantum disk cases, respectively. The case of the quantum wedge follows from the case of the quantum disk and the same argument as in the proof of Lemma 3.5.

Proof of Tutte embedding convergence result.
Proof of Theorem 1.2. Let (D, h, 0, 1) be a quantum disk with fixed boundary length and area, with one marked boundary point and one marked interior point. Let D h and μ h be the √ 8/3-LQG metric and area measure and let ξ h be the path which traverses ∂D counterclockwise from 1 to 1 in such a way that it traverses one unit of √ 8/3-LQG length in one unit of time. By [MS16a, Corollary 1.5], we know that the curve-decorated metric measure space (D, D h , μ h , ξ h ) is a Brownian disk with unit area and boundary length. Furthermore, by the definition of a marked quantum disk, if we condition on this curve-decorated metric measure space then the marked point 0 is a uniform sample from μ h .
For λ > 0, define the Poisson point process P λ , the Voronoi tessellation H λ , and the Tutte embedding λ : P λ → D as in the discussion just above Theorem 1.2 for the Brownian disk (X , D, μ, ξ) = (D, D h , μ h , ξ h ). Note that here the space X is identified with D, so in particular P λ ⊂ D.
(3.5) Indeed, Theorem 3.3 implies that the maximum over all vertices z ∈ P λ h of the Prokhorov distance between the Euclidean harmonic measure on ∂D as viewed from z and the P λharmonic measure on ∂P λ as viewed from z tends to zero in probability as n → ∞. From this and the definition of λ , we get (3.5).
The first two convergence statements in the theorem statement are immediate from (3.5) (for the convergence of re-scaled counting measure, we use that λ −1 times the counting measure on P λ converges in probability to μ h since the intensity measure of P λ is λμ h ). The convergence statement for the random walk on P λ follows from (3.5) and Theorem 3.3.

3.4.
A reformulation of the finite expectation hypothesis. As in Sect. 3.1, let (C, h, 0, ∞) be a 0-quantum cone and write P = P 1 h and H = H 1 h for its associated Poisson point process and collection of Voronoi cells with λ = 1. Proving Proposition 3.1 (i.e., the finite expectation hypothesis in Theorem 2.2 for H) directly turns out to be difficult since Voronoi cells depend on the field in a rather delicate way, so it is not clear how to lower-bound the Lebesgue measure of the origin-containing cell H 0 . Instead, we will use the following lemma which allows us to lower-bound the Lebesgue measure of an LQG metric ball instead. Lemma 3.8 is essentially a consequence of the "mass transport" definition of translation invariance modulo scaling in [GMS18, Definition 1.2]. However, the balls B H are not functions of the cell configuration H itself (they depend on additional randomness from the field) so we will need a trivial reformulation of the mass transport condition which allows for this.
A decorated cell configuration is a cell configuration H together with a compact set K H ⊂ C associated with each cell H ∈ H. We can define a topology on the space of decorated cell configurations by the obvious extension of (2.5): where d Haus denotes the Hausdorff distance and each of the infima is over all homeomorphisms f r : C → C such that f r takes each cell in H(B r (0)) to a cell in H (B r (0)) and preserves the adjacency relation, and f −1 r does the same with H and H reversed. increasing sequence of open sets U j ⊂ C, each of which is either a square or a disk, whose union is all of C such that the following is true. Conditional on (H, {K H } H ∈H ) and U j , let z j for j ∈ N be sampled uniformly from Lebesgue measure on U j . Then there are random numbers C j > 0 (possibly depending on (H, {K H } H ∈H ) and z j ) such that in law with respect to the metric (3.7).
Exactly as in [GMS18, Definition 1.2], one can formulate various equivalent definitions of translation invariance modulo scaling for cell configurations and prove that these definitions are equivalent via exactly the same arguments as in [GMS18, Appendix A]. For our purposes, we will need the "mass transport" formulation of translation invariance modulo scaling for decorated cell configurations.

Lemma 3.10. (Mass transport condition) A random decorated cell configuration is (H, {K H } H ∈H ) is translation invariant modulo scaling in the sense of Definition 3.9 if and only if it satisfies the following condition. Suppose that F(H, {K H } H ∈H , x, y) is a non-negative measurable function on the space of decorated cell configurations with two marked points in C such that F is covariant with respect to dilations and translations of the plane in the sense that for each C > 0 and z ∈ C,
(3.8) Then In light of Lemma 3.8, we only need to prove that the expectation on the right side of (3.6) is finite. Actually, we will prove the following much stronger statement.
Proposition 3.11. For each p > 0, we have We emphasize that Proposition 3.11 does not imply that diam(H 0 ) 2 area(H 0 ) deg(H 0 ) has finite moments of all positive orders (rather, we only get that it has a finite first moment) since Lemma 3.8 does not allow us to compare moments of order greater than 1. The rest of the paper is devoted to the proof of Proposition 3.11.

Estimates for LQG Metric Balls
The goal of this section is to establish Proposition 3.11, which together with Lemma 3.8 will conclude the proof of our main results. Along the way, we will establish a number of estimates for LQG metric balls which are of independent interest (see in particular Propositions 4.3, 4.4, 4.5 and 4.8).
We start out in Sect. 4.1 by introducing two random distributions on C which are defined using the white noise decomposition of the GFF and which have certain nice properties that the GFF itself does not. The first of these distributions, which we call h, possesses exact scale and translation invariance properties (not just scale and translation invariance modulo additive constant, like the whole-plane GFF). The second, which we call h tr , has the property that its restrictions to two sets at distance at least 1/5 from each other are independent. We then state a general lemma (Lemma 4.1) which allows us to compare h, h tr , and the whole-plane GFF.
In Sect. 4.2, we prove an upper bound for the LQG distance across a Euclidean annulus (Proposition 4.3) as well as estimates to the effect that an LQG metric ball B typically contains a Euclidean ball of radius comparable to the Euclidean diameter of B (Propositions 4.4, 4.5). These are proven using percolation arguments which rely crucially on the local independence property of h tr . In Sect. 4.3, we prove that the LQG area of an LQG metric ball of radius r is tightly concentrated around r 4+o r (1) (Proposition 4.8). This is proven by starting with known estimates for metric balls in the Brownian map, then transferring to the GFF using the equivalence of Brownian surfaces and LQG surfaces [MS16a], and finally using the local independent property of h tr to go from events of high probability to events of superpolynomially high probability.
In Sect. 4.4, we transfer the estimates of the preceding subsections from the wholeplane GFF to the 0-quantum cone. In Sect. 4.5, we conclude the proof of Proposition 3.11.

White-noise approximation of the Gaussian free field.
In this subsection we will introduce various white-noise approximations of the GFF which are often more convenient to work with than the GFF itself. Similar approximations to the ones used here were also studied in [DG16,DZZ18,DG18]. Let W be a space-time white noise on and Borel measurable sets A ⊂ C and I ⊂ [0, ∞), we slightly abuse notation by writing For an open set U ⊂ C, we write p U (s; z, w) for the transition density of Brownian motion killed upon exiting U , so that for s ≥ 0, z ∈ C, and A ⊂ U , the integral A p U (s; z, w) dw gives the probability that a standard planar Brownian motion B started from z satisfies B([0, s]) ⊂ U and B s ∈ A. We also write We define the centered Gaussian process We set h := h 0 . By [DG16, Lemma 3.1] and Kolmogorov's criterion, each h t for t ∈ (0, 1] admits a continuous modification. Henceforth whenever we work with h t we will assume that it has been replaced by such a modification. The process h does not admit a continuous modification, but its integral against any smooth compactly supported test function has finite variance, so it makes sense as a distribution. We record for reference the formula which is immediate from (4.1). The distribution h is often more convenient to work with than the GFF thanks to the following symmetries, which are immediate from the definition.
• Rotation/translation/reflection invariance. The law of h is invariant with respect to rotation, translation, and reflection of the plane.
One property which h does not possess is spatial independence. To get around this, we will sometimes work with a truncated variant of h where we only integrate over a ball of finite radius. We define and we interpret h tr as a random distribution. The key property enjoyed by h tr is spatial independence: if A, B ⊂ C with dist(A, B) ≥ 1/5, then h tr | A and h tr | B are independent. Indeed, this is because h tr | A and h tr | B are determined by the restrictions of the white noise W to the disjoint sets B 1/10 (A) × R + and B 1/10 (B) × R + , respectively. Unlike h, the distribution h tr does not possess any sort of scale invariance but its law is still invariant with respect to rotations, translations, and reflections of C.
The following lemma is proven using elementary calculations for the transition density p U (t; z, w) together with the Kolmogorov continuity theorem (see, e.g., [DG18, Lemma 3.1]). It will allow us to use h tr or h in place of the GFF in many of our arguments.  is a coupling (h, h, h tr ) of a wholeplane GFF normalized so that h 1 (0) = 0 and the distributions from (4.1) and (4.3) such that the following is true. For h 1 , h 2 ∈ {h, h, h tr }, the distribution (h 1 − h 2 )| K a.s. admits a continuous modification and there are constants c 0 , c 1 > 0 depending only on K such that for A > 1, (4.4) In fact, in this coupling one can arrange so that h and h tr are defined using the same white noise W .
The existence of continuous modifications in Lemma 4.1 allows us to define the √ 8/3-LQG metrics D h and D h tr . Indeed, this is because we know how to define D h+ f when h is a GFF and f is a continuous function (see the discussion just after Lemma 2.3). Moreover, we get that D h and D h tr each a.s. induces the Euclidean topology on C. Due to Lemma 2.3, the estimate (4.4) will allow us to compare D h and D h tr to the √ 8/3-LQG metrics induced by a GFF.
The following lemma will be used when we apply the scaling property of h.

Comparing LQG metric balls and Euclidean balls.
Throughout this subsection, we let h be a whole-plane GFF normalized so that h 1 (0) = 0. We will prove the following three propositions, which relate D h -metric balls and Euclidean balls. Our first estimate implies in particular that a D h -metric ball is extremely unlikely to have an unusually large Euclidean diameter. This estimate is related to the fact that μ h (D) has negative moments of all orders (see [DS11,Lemma 4.5] or [RV14, Theorem 2.12]) but is proven in a very different way. We next state two closely related estimates to the effect that a D h -metric ball typically contains a Euclidean ball of radius comparable to its Euclidean diameter. These are analogs for D h -metric balls of estimates for space-filling SLE cells from [GHM15, Section 3] and [GMS17, Section 4], and will be used to control the ratio diam(H 0 ) 2 /area(B H 0 ) appearing in Proposition 3.11. For ρ > 0, define the square annulus A ρ and its inner and outer boundaries by and ∂ out A ρ := ∂([−2ρ, 2ρ] 2 ). (4.7) The main step in the proof of the above propositions is Lemma 4.6. For n ∈ N, we consider the restriction to A n of the truncated white noise field h tr of (4.3). The reason for considering A n instead of A 1 , say, is that two sets need to be at Euclidean distance at least 1/5 from each other for the restrictions of h tr to be independent and we want to define lots of independent events. Basic properties of the √ 8/3-LQG metric show that for each 1 × 1 square S ⊂ A n , it holds with probability tending to 1 as C → ∞ that the D h tr -distance from the boundary of the 1/2-neighborhood ∂ B 1/2 (S) to S is at least 1/C and each Euclidean ball of radius e −Cn 2/3 which intersects S has D h tr -diameter at most e −n 2/3 (this last condition would also hold with e −Cn 2/3 and e −n 2/3 replaced by, e.g., n −C and n −1 or 1 and 1/C, but we use e −n 2/3 since we will get an error of order e n 1/2 in Lemma 4.7). The restrictions of h tr to squares which lie at distance at least 1/5 from one another are independent, so the adjacency graph of "good" squares which satisfy the above properties looks like a very supercritical percolation on Z 2 when C is large. Hence with exponentially high probability in n, there is path of such good squares which separates the inner and outer boundaries of A n . This implies analogs of Propositions 4.3 and 4.4 for h tr | A n . In Lemma 4.7, we set n (log −1 ) 3/2 and transfer from h tr | A n to h| A ρ (and thereby to h| A ρ ) using Lemma 4.1 and the scale invariance properties of the field h.
The square annulus A n of (4.7) is shown in light green. To prove Lemma 4.6, we use a percolation argument (based on the local independence property of h tr ) to show that with extremely high probability, we can find a collection of 1 × 1 squares S (light blue) whose union disconnects the inner and outer boundaries of A n and such that each path which crosses one of these squares has to have D h tr -length at least 1/C and has to hit a Euclidean ball of radius e −Cn 2/3 and D h tr -diameter at most e −n 2/3 . After re-scaling by 1/n, this statement is used to prove Propositions 4.3, 4.4 and 4.5 Lemma 4.6. Define A n for n ∈ N as in (4.7). There are universal constants a 0 , a 1 > 0 and C > 1 such that for each n ∈ N, it holds with probability at least 1 − a 0 e −a 1 n that the following is true.
1. The D h tr -distance from ∂ in A n to ∂ out A n is at least 1/C. 2. Each path from ∂ in A n to ∂ out A n intersects a Euclidean ball with Euclidean radius e −Cn 2/3 and D h tr -diameter at most e −n 2/3 . See Fig. 2 for an illustration of the statement and proof of Lemma 4.6. We will eventually apply the lemma with n (log −1 ) 3/2 , so that e −n is smaller than any power of and e −n 2/3 ≈ .
Proof of Lemma 4.6. Let p ∈ (0, 1) be a small universal constant to be chosen later, in a universal manner. For n ∈ N, let S(A n ) be the set of unit side length squares with corners in Z 2 whose Euclidean 1-neighborhood satisfies B 1 (S) ⊂ A n . For S ∈ S(A n ) and C > 1, let E S (C) be the event that the following is true. Henceforth fix such a C. View S(A n ) as a graph with two squares considered to be adjacent if they share an edge. We claim that if p is chosen sufficiently small, then for appropriate constants a 0 , a 1 > 0 as in the statement of the lemma, it holds for each n ∈ N that with probability at least 1 − a 0 e −a 1 n , we can find a path P of squares in S(A n ) which disconnects ∂ in A n from ∂ out A n such that E S (C) occurs for each S ∈ P.
Assume the claim for the moment. If a path P as in the claim exists, then each Euclidean path from ∂ in A n to ∂ out A n must pass through one of the squares S ∈ P. Since E S (C) occurs for each such square, each such path must hit a Euclidean ball of radius e −Cn 2/3 centered at a point of S which has D h tr ≤ D h tr | B 1/2 (S) -diameter at most e −n 2/3 , i.e., condition 2 in the lemma statement holds. Furthermore, since B 1/2 (S) ⊂ A n for each S ∈ S(A n ), any path from ∂ in A n to ∂ out A n must cross one of the annuli B 1/2 (S)\S for some S ∈ P. Since E S (C) occurs for each such S, condition 1 in the lemma statement holds.
It remains only to prove the claim. Let S * (A n ) be the graph whose squares are the same as the squares of S(A n ), but with two squares considered to be adjacent if they share a corner or an edge, instead of only considering squares to be adjacent if they share an edge. We define the inner and outer boundaries of S * (A n ) to be the set of squares which lie at Euclidean distance 1 from ∂ in A n and ∂ out A n , respectively (recall that squares in S(A n ) satisfy B 1 (S) ⊂ A n ). By planar duality, it suffices to show that if p, a 0 , a 1 are chosen appropriately, then it holds with probability at least 1 − a 0 e −a 1 n that there does not exist a simple path in S * (R n ) from the inner boundary of A n to the outer boundary of A n consisting of squares for which E S (C) does not occur. This will be proven by a standard argument for subcritical percolation.
By the definition (4.3) of h tr , the event E S (C) is a.s. determined by the restriction of the white noise W to B 1 (S) × R + . In particular, E S (C) and E S (C) are independent whenever B 1 (S) ∩ B 1 ( S) = ∅. For each fixed deterministic simple path P in S * (R n ), we can find a set of at least |P|/100 squares hit by P for which the neighborhoods B 1 (S) are disjoint. Since the events E S (C) for these |P|/100 squares are independent and each has probability at least 1 − p, the probability that E S (C) fails to occur for every square in P is at most p |P|/100 . We now take a union bound over all simple paths P in S * (A n ) connecting the inner and outer boundaries. For k ∈ [n, 16n 2 ] Z , the number of such paths with |P| = k is at most 4n8 k since there are 4n possible initial squares along the inner boundary of A n and 8 choices for each step of the path. Combining this with the estimate in the preceding paragraph, we find that the probability of an inner-outer crossing in S * (A n ) consisting of squares for which E S (C) does not occur is at most 4n 16n 2 k=n p k/100 8 k+1 , which is bounded above by an exponential function of n provided we take p < 8 −100 .
We now transfer from h tr | A n to h| A ρ .
Lemma 4.7. There is a universal constant C > 1 such that for each ρ ∈ (0, 1), it holds with superpolynomially high probability as → 0 that the following is true.

intersects a Euclidean ball with Euclidean radius at least C and D h -diameter at most .
Proof. We will apply Lemma 4.6 with n (log −1 ) 3/2 together with a scaling argument. We first establish an estimate for D h -distances in A n . By Lemma 4.1 and a union bound over O n (n 2 ) Euclidean balls of unit radius which cover A n , we can find constants c 0 , c 1 > 0 and a coupling of h and h tr such that If max z∈A n |( h− h tr )(z)| ≤ A, then by the scaling property of LQG distances (Lemma 2.3), Setting A = √ 6n 1/2 and applying Lemma 4.6, we see that there is a universal constant C > 1 such that with exponentially high probability as n → ∞, the following is true.
Choose a finite set Z n of O n (n 4 ) points z ∈ A n such that A n ⊂ z∈Z n B ρ/n (z). By Lemma 4.2 (applied with δ = ρ/n and C = (1/2)n 1/2 ), the Gaussian tail bound applied to the O n (n 4 ) centered Gaussian random variables h ρ/n (z) for z ∈ Z n , each of which has variance log(ρ/n), and a union bound, we can find constants c 0 , c 1 > 0, depending only on ρ, such that We know that ( h − h ρ/n )((ρ/n)·) d = h, so by combining (4.8) and the conclusion of the first paragraph with ( h − h ρ/n )((ρ/n)·) in place of h, we get that (after possibly shrinking c 0 and c 1 ) it holds with probability at least 1 − c 0 n 4 e −c 1 n that the following is true. 1. The D h -distance from ∂ in A ρ to ∂ out A ρ is at least e −O n (n 1/2 ) . 2. Each path from ∂ in A ρ to ∂ out A ρ intersects a Euclidean ball with Euclidean radius (ρ/n)e −Cn 2/3 and D h -diameter at most e −n 2/3 +O n (n 1/2 ) .
We now choose n = (log −1 ) 3/2 . This makes it so that n 4 e −c 1 n/ log n decays faster than any positive power of , e −n 1/2 decays slower than any positive power of , (ρ/n)e −Cn 2/3 = C+o (1) , and e −n 2/3 +O n (n 1/2 ) = 1+o (1) . Making this choice of n and possibly slightly adjusting C and concludes the proof.
Proof of Proposition 4.3. By Lemma 4.1, we can couple h and h in such a way that (h − h)| D is a continuous function and with superpolynomially high probability as → 0, one has max z∈D |(h − h)(z)| ≤ (log −1 ) 2/3 . Combining this with Lemma 4.7 and the scaling property of LQG distances shows that for each fixed square annulus A ⊂ D such that the ratio of its inner and outer side lengths is 4, it holds with superpolynomially high probability as → 0 (at a rate depending on A) that the D h -distance between the inner and outer boundaries of A is at least . We can find finitely many such square annuli contained in D\B ρ (0) such that the union of their inner boundaries disconnects the inner and outer boundaries of D\B ρ (0). Each path between the inner and outer boundaries of D\B ρ (0) must cross between the inner and outer boundaries of one of these square annuli, so applying the preceding estimate once to each such annulus and taking a union bound concludes the proof.
Proof of Proposition 4.4. Via the same argument as in the proof of Proposition 4.3, Lemma 4.7 implies that there is a universal constant C > 1 such that with superpolynomially high probability as → 0, each path from B ρ (0) to ∂ B (1+ρ)/2 (0) intersects a Euclidean ball with Euclidean radius at least C and D h -diameter at most . In particular, each D h -ball B which intersects both ∂ B ρ (0) and ∂D intersects a Euclidean ball of radius at least C and D h -diameter at most which is contained in B (1+ρ)/2 (0). On the other hand, Proposition 4.3 shows that with superpolynomially high probability as → 0, the D h -distance from B (1+ρ)/2 (0) to ∂D is at least 2 , in which case the aforementioned Euclidean ball is contained in the D h -metric ball B. Replacing C by concludes the proof.
Proof of Proposition 4.5. Observe that the conclusion of Proposition 4.4 does not depend on the choice of additive constant for h. By the scale and translation invariance of the law of h, modulo additive constant, we see that Proposition 4.4 implies that for each ρ ∈ (0, 1), r > 0, and z ∈ C, it holds with superpolynomially high probability as → 0, at a rate which is uniform in r and z, that each D h -metric ball which intersects both ∂ B ρr (z) and ∂ B r (z) contains a Euclidean ball of radius at least r .

Volume estimates for LQG metric balls.
The goal of this subsection is to establish the following estimate for the LQG mass of LQG metric balls.
Proposition 4.8. Let h be a whole-plane GFF normalized so that h 1 (0) = 0. For each ζ ∈ (0, 1), it holds with superpolynomially high probability as → 0 that (4.9) We will extract Proposition 4.8 from known ball volume estimates for the Brownian map, which say that a.s. the volume of every ball of radius s in the Brownian map simultaneously is bounded above and below by constants times s 4−ζ and s 4+ζ (stated as Lemma 4.9). These estimates together with the equivalence of the Brownian map and the quantum sphere do not immediately imply (4.9) since we are working with a whole-plane GFF instead of a quantum sphere. One could attempt to transfer the estimates using some sort of quantitative local absolute continuity, but we instead take a different approach which we find to be easier. We note that Proposition 4.8 has not previously appeared in the Brownian map literature, although closely related results have been established (see the proof of Lemma 4.9).
Local absolute continuity (without any quantitative Radon-Nikodym derivative bound) shows that a.s. the μ h -mass of every D h -ball of radius s ∈ (0, 1) which is contained in D is bounded above and below by constants times s 4−ζ and s 4+ζ . To turn this into a bound which holds with superpolynomially high probability instead of just a.s., we first use Lemma 4.1 and scale invariance considerations to transfer from h| D to the restriction of the truncated white-noise field h tr of (4.3) to B R (0) for a large value of R (Lemma 4.10). The restrictions of h tr to radius-1 Euclidean balls contained in B R (0) which lie at distance at least 1/5 from one another are independent. Hence, the fact that an event (in our setting, bounds for the μ h -mass of D h -balls contained in the Euclidean ball) holds simultaneously for all such Euclidean balls with high probability shows that in fact the event for a single Euclidean ball has to hold with extremely high probability (Lemma 4.11). We then transfer back to h to conclude the proof.
Let us first record what we get from Brownian map estimates.
The same is true with the field h of (4.1) in place of h.
Proof. We will use estimates for the Brownian map, so we need to work with a quantum sphere due to the equivalence of the Brownian map and quantum sphere [MS16a, Corollary 1.4]. Let h Sph be an embedding into C of the quantum sphere (say, conditioned to have LQG area at least 1), normalized so that two marked points sampled uniformly from μ h Sph are sent to 0 and ∞ and so that 1 = sup{r > 0 : h Sph r (0) + Q log r = 0}, provided h Sph r (0) + Q log r = 0 for some r > 0. This choice of normalization makes it so that the laws of h Sph | D\B 1/2 (0) and h| D\B 1/2 (0) are mutually absolutely continuous on the event {h Furthermore, since the continuum label process (the "head of the Brownian snake") used to define the Brownian map is a.s. Hölder continuous with any exponent less than 1/4, a.s. inf By the local absolute continuity between h and h Sph , we obtain (4.10) with D\B 1/2 (0) in place of D. By the translation invariance of the law of h, modulo additive constant, we get (4.10). By Lemma 4.1, the same is true with h in place of h.
We will now transfer to an estimate for h tr restricted to a large ball.
Lemma 4.10. Fix p > 1 and ζ ∈ (0, 1). There is a universal constant c > 0 and a random set Z r p ⊂ B r p (0) ∩ (3Z 2 ) independent from h tr such that with probability tending to 1 as r → ∞, one has #Z r p ≥ cr 2 p and for each z ∈ Z r p , (4.11) The reason why we need to restrict to the set Z r p in Lemma 4.10, instead of looking at all points in B r p (0) ∩ (3Z 2 ), is as follows. To transfer from an estimate on D to an estimate on B r p (0), we will use the scale invariance property of the white noise field h, which says that h(r − p ·)− h r − p (r − p ·) has the same law as h and is independent from h r − p . We will restrict attention to the set of points where h r − p is not unusually large, which Proof of Lemma 4.10.
Step 1: re-scaling from D to B r p (0). We will re-scale with the eventual aim of transferring Lemma 4.9 to an estimate with B r p (0) in place of D and h in place of h. This will lead to the definition of Z r p . If we set h r p : then h r p d = h and h r p is independent from h r − p . Let (4.12) We emphasize that Z r p is determined by h r − p , so is independent from h r p .
We will now argue that there is a universal constant c > 0 such that P #Z r p ≥ cr 2 p → 1 as r → ∞.
(4. 13) To see this, we observe that each h r − p (r − p z) for z ∈ 3Z 2 is Gaussian with variance log r p . By the Gaussian tail bound, P[| h r − p (r − p z)| ≤ ζ 2 2 log r ] tends to 1 as r → ∞, uniformly over all z ∈ 3Z 2 . By Markov's inequality, it holds with probability tending to 1 as r → ∞ that the number of z ∈ B r p (0) ∩ (3Z 2 ) for which | h r − p (r − p z)| ≤ ζ 2 2 log r is at least 1 2 # B r p (0) ∩ (3Z 2 ) , say. This last quantity is at least cr 2 p for some universal constant c > 0. By Lemma 4.2, it holds with probability tending to 1 as r → ∞ that By combining these estimates, we get (4.13).
Step 2: estimate for LQG balls centered at points of Z r p . By the LQG coordinate change formula (2.11) and Lemma 2.3, for each z ∈ Z r p and each x, y ∈ B 1 (z), Moreover, the analogous properties for the LQG measure show that Combining this with Lemma 4.9 (with h in place of h) shows that with probability tending to 1 as r → ∞, it holds for each z ∈ Z r p that (4.14) where here a > 0 is a universal constant. Note that we used that 4/ √ 6 = √ 8/3 to cancel two large powers of r and we used that r ζ 2 ≤ s −ζ for s ≤ r −ζ to absorb a small power of r into a power of s.
Step 3: transferring from h to h tr . By Lemma 4.1 and a union bound over O r p (r 2 p ) Euclidean balls of radius 1 which cover B r p (0), we can couple h r p and h tr in such a way that with probability tending to 1 as r → ∞, we have max z∈B r p +1 (0) |( h r p − h tr )(z)| ≤ (log r ) 2/3 . Since h r p is independent from Z r p , we can take h tr to be independent from Z r p in this coupling. Our choice of coupling together with (4.13) and (4.14) shows that with probability tending to 1 as r → ∞, it holds for each z ∈ Z r p , each w ∈ B 1 (z), and each s ≤ r −ζ e After adjusting ζ appropriately, this gives (4.11).
We can now go from events with probability tending to 1 to events with superpolynomially high probability.
Lemma 4.11. For each ζ ∈ (0, 1), it holds with superpolynomially high probability as → 0 that (4.15) Proof. Fix p > 1, which we will eventually send to ∞. For z ∈ 3Z 2 and r > 1, let Note that the event in (4.15) is the same as E −1/ζ (0). Furthermore, if we let Z r p ⊂ B r p (0) ∩ (3Z 2 ) be the random set independent from h tr from Lemma 4.10, then that lemma tells us that with probability tending to 1 as r → ∞, we have #Z r p ≥ cr 2 p and E r (z) occurs for every z ∈ Z r p .
The fields h tr | B 1 (z) for different choices of z ∈ 3Z 2 are independent and the law of h tr is invariant with respect to spatial translations. Consequently, the events E r (z) for different choices of z ∈ 3Z 2 are independent. Since h tr is independent from Z r p , we get that If r is large enough that this 1 − o r (1) is at least e −c , then re-arranging gives where here we have used the elementary inequality 1 − e −x ≤ x. Since p > 1 can be made arbitrarily large, we get that E r (0) occurs with superpolynomially high probability as r → ∞. Setting r = −1/ζ now concludes the proof.
Proof of Proposition 4.8. By Lemmas 4.1 and 4.11 together with the scale invariance of the law of h, modulo additive constant, and the fact that the law of h 2 (0) is Gaussian with constant-order variance, it holds with superpolynomially high probability as → 0 that By Proposition 4.3, it holds with superpolynomially high probability as → 0 that D h (∂D, ∂ B 2 (0)) ≥ 1/2 . Combining these estimates gives (4.9).

Estimates for the 0-quantum cone.
We now want to shift attention from the wholeplane GFF to the 0-quantum cone, with a view toward proving Proposition 3.11. To this end, we will transfer the main results of the preceding subsections to the case of a 0-quantum cone. We start with estimates for the LQG areas of LQG metric balls which follows from Proposition 4.8.
1. With superpolynomially high probability as C → ∞, one has C −1 ≤ μ h (B 1 (0; D h )) ≤ C. 2. For ζ ∈ (0, 1), it holds with superpolynomially high probability as C → ∞ that By the scaling property (2.10) of the 0-quantum cone, the law of (C, D h , μ h ) as a metric measure space is invariant under scaling distances by b 1/4 and areas by b, for any b > 0. We will often use this fact in conjunction with Proposition 4.12 without comment.
Proof of Proposition 4.12. The proposition statement does not depend on the choice of embedding for h, so we can assume without loss of generality that h is given the circle-average embedding. Recall that h| D agrees in law with the corresponding restriction of a whole-plane GFF normalized so that its circle average over ∂D is zero. By Proposition 4.3, it holds with superpolynomially high probability as C → ∞ that B C −1 (0; D h ) ⊂ D. Hence Proposition 4.8 (applied with ζ = 1, say) shows that with superpolynomially high probability as C → ∞, (4.17) By the scale invariance property of the 0-quantum cone (2.10), we have (C, D h , μ h ) d = (C, C D h , C 4 μ h ) as metric measure spaces. Therefore, (4.17) implies that with superpolynomially high probability as C → ∞, one has C −1 ≤ μ h (B 1 (0; D h )) ≤ C, which is assertion 1.
We now prove assertion 2 via a similar argument. By Proposition 4.8 and our above description of the law of h| D , it holds with superpolynomially high probability as C → ∞ that Furthermore, by Proposition 4.3, it holds with superpolynomially high probability as Hence with superpolynomially high probability as C → ∞, We now scale distances by C 2 and areas by C 8 and apply (2.10) as above to deduce assertion 2 from (4.18).
We next record an estimate to the effect that D h -metric balls have to contain Euclidean metric balls of radius comparable to their Euclidean diameters.
Proposition 4. 13. Let (C, h, 0, ∞) be a 0-quantum cone and let ζ ∈ (0, 1). With superpolynomially high probability as C → ∞, each D h -ball B which is contained in B C (0) and which has D h -radius at least C −1 contains a Euclidean ball of radius at least C −ζ diam(B).
We will deduce Proposition 4.13 from Propositions 4.3 and 4.5. Before we can do so, however, we need some basic polynomial tail estimates for the minimal and maximal radii of Euclidean balls. This is because Proposition 4.13 only holds for LQG balls with sufficiently small Euclidean diameter and because diam(B) 1+ζ can be much smaller than Proof. By standard estimates for the √ 8/3-LQG measure (see, e.g., the proof of [DG18, Lemma 3.7] with γ = √ 8/3), for p > 2 √ 8/3, it holds with probability at least 1 − δ 3 p 2 /16−2 that each Euclidean ball centered at a point of D with radius at least δ has μ h -mass at most δ 10/3− p . We now fix ζ ∈ (0, 1), which we will eventually send to 0. Applying the above estimate with p = 10/3 − (4 + ζ )/q and δ = q = (4+ζ )/(10/3− p) and shows that with probability at least 1− α(q)+o ζ (1)+o (1) (with the o ζ (1) deterministic and independent of ), each Euclidean ball centered at a point of D with radius q has μ h -mass at most 4+ζ .
By Proposition 4.8 (applied with B 2 (0) in place of D and 1/2 in place of ), we see that with superpolynomially high probability as → 0, μ h (B (z; D h )) < 4−ζ for each z ∈ D. By Proposition 4.5, it holds with superpolynomially high probability as → 0 that each D h -ball centered at a point of D which has Euclidean diameter at least q contains a Euclidean ball of radius at least q . By the preceding estimates, with probability at least 1 − β( q)+o ζ (1)+o (1) , each such D h -ball has μ h -mass at least 4−ζ and hence D h -radius strictly larger than . Sending ζ → 0 and then q → q concludes the proof.
Proof of Proposition 4.13. Let q ∈ 0, We will eventually send q to 0 and q to ∞. By Lemmas 4.14 and 4.15, it holds with probability at least 1 − C −(3α(q))∧β(q)+o C (1) that each D h -ball contained in D with D h -diameter between C −3 and C −1 has Euclidean diameter between C −3q and C −q . By Proposition 4.5, it holds with superpolynomially high probability as C → ∞ that each D h -ball contained in D with Euclidean diameter at most C −q contains a Euclidean ball of radius at least diam(B) 1+ζ /(2q) . By Proposition 4.3, it holds with superpolynomially high probability as C → ∞ that B C −1 (0; D h ) ⊂ D. Hence with probability at least to h (i.e., scale distances by C 2 ) to get that the event in the statement of the lemma holds with probability at least 1 − C −(3α(q))∧β(q)+o C (1) . Since α(q), β(q) → ∞ as q → ∞ and q → 0, this concludes the proof. Proof. Let Z be a maximal collection of points in B C (0; D h ) such that the balls B 1/2 (z; D h ) for z ∈ Z are disjoint. By Proposition 4.12, it holds with superpolynomially high probability as C → ∞ that min z∈Z μ h (B 1/2 (z; D h )) ≥ C −ζ /2 and μ h (B C (z; D h )) ≤ C 4+ζ /2 , which implies that #Z ≤ C 4+ζ . By the maximality of Z, The following lemma shows that Voronoi cells are extremely unlikely to have a larger quantum diameter than one would expect. To prove Lemma 4.17, we will use the following lemma to build a "wall" of Voronoi cells which separate the boundaries of two concentric D h -balls.
Lemma 4.18. For ζ ∈ (0, 1), it holds with superpolynomially high probability as C → ∞ that the following is true. For each z ∈ B C (0; D h ), we can find a finite collection of at most C ζ D h -metric balls of radius 1/2 which are contained in B 3 (z; D h )\B 1 (z; D h ) and whose union disconnects B 1 (z; D h ) from C\B 3 (z; D h ).
Proof. By Proposition 4.12, it holds with superpolynomially high probability as C → ∞ that Henceforth assume that (4.22) holds and fix z ∈ B C (0; D h ). We will construct a collection of D h -balls as in the statement of the lemma.
Let C be a maximal collection of points in ∂ B 2 (z; D h ) such that the balls B 1/4 (w; D h ) for w ∈ C are disjoint. The union of the balls B 1/2 (w; D h ) for w ∈ C covers ∂ B 2 (z; D h ) (otherwise, we could find a point in ∂ B 2 (z; D h ) which lies at distance at least 1/2 from each w ∈ C, which contradicts the maximality of C).
Finally, since the balls B 1/4 (w; D h ) for w ∈ C are disjoint and contained in B 3 (z; D h ), we see from (4.22) that #C ≤ C ζ .
Proof of Lemma 4.17. By Lemma 4.16, with superpolynomially high probability as C → ∞, we can find a collection Z C of at most C 4+ζ points of B C (0; D h ) such that the union of the balls B 1 (z; D h ) for z ∈ Z C covers B C (0; D h ). By Lemma 4.18 and the scaling property (2.10) of the 0-quantum cone, with superpolynomially high probability as C → ∞ we can find for each z ∈ Z C a finite collection C(z) of at most C ζ D h -balls of radius 1 2 C ζ which are contained in B 3C ζ (z; D h )\B C ζ (z; D h ) and whose union disconnects B C ζ (z; D h ) from C\B 3C ζ (z; D h ). By Proposition 4.12, it holds with superpolynomially high probability as C → ∞ that the μ h -mass of each ball in each of the collections C(z) is at least C 4ζ (1−ζ ) By the formula for the Poisson distribution, if this is the case then the conditional probability given h that each of the balls in z∈Z C C(z) contains a point of P is at least 1 − exp −C 4ζ(1−ζ ) . By a union bound over the at most C 4+2ζ balls in z∈Z C C(z), we find that with superpolynomially high probability as C → ∞, each of the balls this collection contains a point of P. Similarly, it holds with superpolynomially high probability as C → ∞ that B C ζ /4 (z; D h ) contains a point w z ∈ P for each z ∈ Z C . Since each point of B 1 (z; D h ) lies within D h -distance C ζ /4 + 1 of w z , this means that the center point of each Voronoi cell which intersects If the events described in the preceding paragraph are satisfied, then for z ∈ Z C , each point of C\B 3C ζ (z; D h ) is D h -closer to a point of P which is contained in one of the balls in C(z) than it is to any point of B C ζ /2 (z; D h ). This means that no such point can be contained in a cell whose center point is in B C ζ /2 (z; D h ), hence no such point can be contained in a cell which intersects B 1 (z; D h ). Hence each Voronoi cell which intersects B 1 (z; D h ) is contained in B 3C ζ (z; D h ) with superpolynomially high probability as C → ∞. Since the union of the balls B 1 (z; D h ) for z ∈ Z C covers B C (0; D h ), this gives the statement of the lemma with 3C ζ in place of C ζ , which is sufficient.
We can now prove our main moment estimates. Recall that B H denotes the smallest D h -metric ball containing the cell H which is centered at the center point of H . Proof. Fix ζ ∈ (0, 1), which we will eventually send to 0. Lemma 4.17 shows that with superpolynomially high probability as C → ∞, each H ∈ H with 0 ∈ B H has D h -diameter at most C ζ , so is contained in B 2C ζ (0; D h ). We will now argue that with probability at least 1− O C (C −4/ζ +ζ ), the D h -radius of B H for each such cell H is at least C −1/ζ . Indeed, Proposition 4.12 shows that with superpolynomially high probability as then the D h -radius of B H must be at least C −1/ζ , as required. By Proposition 4.13, it holds with superpolynomially high probability as C → ∞ that each D h -ball B centered at a point of B C (0; D h ) with D h -radius at least C −1/ζ contains a Euclidean ball of radius at least C −1/2 diam(B). Combining this with the preceding paragraph shows that (4.23) holds with probability at least 1 − O C (C −4/ζ +ζ ). Sending ζ → 0 concludes the proof.
To bound the number of cells with 0 ∈ B H and their degrees, we will need the following lemma.
we see that the probability that B 5C (0; D h ) contains more that C 4+ζ points of P decays superpolynomially in C. This gives (4.24). Consequently, this sum has finite moments of all positive orders.

Open Problems
Perhaps the most natural question to ask about Brownian motion on the Brownian map is the following.
Problem 5.1. Show that random walk on uniform random planar maps (e.g., uniform quadrangulations or triangulations) converges to Brownian motion on the Brownian map with respect to the Gromov-Hausdorff-Prokhorov-uniform topology, the natural topology for curve-decorated metric measure spaces [GM17b].
It is known that self-avoiding walk and percolation interfaces on uniform random planar maps converge to SLE 8/3 and SLE 6 , respectively [GM16a,GM17a]. In contrast to the case of Z 2 , however, random walk on a random planar map seems harder to analyze than SAW or percolation interfaces since the random walk can re-trace its past, so one cannot explore the curve and the planar map simultaneously using peeling.
Our results only concern Brownian motion on Brownian surfaces viewed modulo time parameterization. The natural way to parameterize Brownian motion on a √ 8/3-LQG surface, equivalently a Brownian surface, is called Liouville Brownian motion and is constructed in [Ber15,GRV16].
Problem 5.2. Show that in the setting of Theorem 1.1, the random walk on P λ , parameterized so that it traverses one edge in one unit of time, converges to Liouville Brownian motion with respect to the uniform topology as λ → ∞.
There are other natural types of random walks that one can consider on Brownian surfaces which one would expect to converge to Brownian motion in the scaling limit. For example, one can generate a random walk which at each step moves to a point sampled uniformly at random from the metric ball of radius centered at its current position. As a second example, the Brownian snake construction of the Brownian map involves describing the Brownian map as a gluing of the tree of geodesics back to the root together with a dual tree (and instance of the CRT) rooted at the dual root. The peanocurve which "snakes between these two trees" is a space-filling curve, which one may use to give a graph approximation analogous to the mated-CRT map considered in [GMS17,GMS19]. In particular, if (X , D, μ) denotes the unit-area Brownian map, the Brownian snake construction gives a quotient map p : [0, 1] → X . We then fix > 0 and consider the random walk on the adjacency graph of μ-mass cells One of the appeals of this construction is that one can sample from it in linear time (one just needs to generate an instance of the Brownian snake) and then compute its Tutte embedding efficiently using a sparse matrix package (c.f. [GMS17, Remark 1.2]). The analog of the property of Problem 5.4 for Brownian motion on certain random planar maps (like the UIPT or the UIPQ) follows from the so-called spatial Markov property, a.k.a. peeling; see, e.g., [BC13]. It is known that the complementary connected components of SLE 6 on a Brownian surface are Brownian disks (this follows from the results of [DMS14,MS15c] and the equivalence of Brownian and √ 8/3-LQG surfaces), so it may be possible to solve Problem 5.4 using the relationship between SLE 6 and Brownian motion. An alternative approach to Problem 5.4 is via Theorem 1.1. Indeed, we know that the complementary connected components of a metric ball on the Brownian map are Brownian disks conditional on their boundary lengths [MS15a,LG19]. Moreover, a Poisson-Voronoi cell is determined by the metric ball centered at its center point whose radius equals twice the distance from the center point to the boundary of the cell, together with the points of the Poisson point process which intersect this ball. It is possible that one could apply this property at the cells hit by the walk to solve Problem 5.4.
Problem 5.4 might have some relevance to Problem 5.1. Indeed, if one can show that Brownian motion on the Brownian map is uniquely characterized by the Markov property of Problem 5.4 together with the law of the boundary lengths of the complementary connected components, then potentially this could be used to identify a subsequential scaling limit of random walk on random planar maps (one would also have to establish tightness). A similar strategy is used to prove the convergence of percolation on random planar maps to SLE 6 in [GM17a]. Theorem 1.1 together with the result of Yadin and Yehudayoff [YY11] allow us to give an intrinsic definition of SLE 2 on a Brownian surface as the limit of the loop-erased random walk on Poisson-Voronoi tessellations.
Problem 5.5. Does the perspective of this paper lead to any insights about SLE 2 on a Brownian surface (concerning, e.g., the law of the surface parameterized by its complement or its relationship to random planar maps)?
Problem 5.5 would be very interesting to solve since currently very little is known about the behavior of SLE κ curves on a γ -LQG surface for κ / ∈ {γ 2 , 16/γ 2 }. "Appendix A" includes several basic properties of Voronoi cells which are needed in the proofs of our main results. However, there are many questions about such cells which have not been answered, for example the following.
Problem 5.6. Is the boundary of a Voronoi cell a.s. given by the union of finitely many disjoint simple curves? What is the Hausdorff dimension of this boundary (with respect to the Euclidean or √ 8/3-LQG metric)? Is the collection of Voronoi cell boundaries a.s. conformally removable?
The simulation in Fig. 1 seems to suggest that the answer to the first part of Problem 5.6 is affirmative. Although we will not explain this in detail here, we expect that the Hausdorff dimension w.r.t. the √ 8/3-LQG metric should be 2. (Roughly speaking, this is because one expects that on a Brownian surface, the set of points equidistant to generic points z 1 and z 2 should be a curve that has the same local structure as a branch of the dual of the tree of geodesics drawn toward a fixed root.) We do not have a conjecture for the Euclidean Hausdorff dimension of the cell boundaries. and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

A. Basic Properties of Voronoi Cells
In this section we will prove a number of a.s. properties of Voronoi cells which are intuitively obvious from, e.g., the simulations (see Fig. 1). In particular, we will show that such cells are connected (Lemma A.2), they are compact and locally finite (Lemma A.4), and their boundaries have zero LQG measure and zero Lebesgue measure (Lemmas A.5, A.6). Throughout, we will consider the setup described at the beginning of Sect. 3, so that for a GFF-type distribution h and λ > 0, P λ h and H λ h denote the associated Poisson point process and collection of Voronoi cells, respectively. The arguments in this section do not use any of the results of Sects. 3 and 4, so can be read independently of those sections.
Remark A.1. We will often consider the Voronoi tessellations H λ h for fields h for which there is a choice in the way that the additive constant is fixed. Different choices scale μ h , and hence the intensity measure for P λ h , by a constant factor. However, for any fixed choice of additive constant for h, the conditional law of P λ h given h is well-defined. Furthermore, if μ h is a.s. finite and we condition on h and the event {#P λ h = n} for n ∈ N, then the conditional law of P λ h is that of a collection of n i.i.d. samples from μ h , so this conditional law does not depend on the choice of additive constant for h (or on λ). If μ h is only locally finite, one can apply the preceding sentence with h replaced by its restriction to a compact set. As a particular consequence of this, the a.s. statements for Voronoi cells which we consider in this section do not depend on the choice of additive constant, so we will not specify it.
We start with the following elementary deterministic fact.
Lemma A.2. Suppose h is a whole-plane GFF. For each z ∈ P λ h and each u in the corresponding cell H z (resp. each u in the interior of H z ), each D h -geodesic from z to u is contained in H z (resp. the interior of H z ). In particular, H z and its interior are both connected. The same is true with h replaced by a free-boundary GFF on a Jordan domain or an embedding of a quantum cone, sphere, disk, or wedge.
Proof. For each > 0, a.s. K can be covered by finitely many D h -balls of radius at most and each of these balls has positive μ h -mass. The conditional probability given h that each such ball contains a point of P λ h tends to 1 as λ → ∞. On this event, each cell which intersects K has D h -diameter at most 2 .
Lemma A.4. Suppose h is a whole-plane GFF and λ > 0. Almost surely, each cell in H λ h has non-empty interior and is compact. Furthermore, a.s. each compact subset of C intersects only finitely many cells of H λ h . The same is true with h replaced by a free-boundary GFF on a Jordan domain or an embedding of a quantum cone, sphere, disk, or wedge.
The proof may seem harder than the reader expects. This is because we need to rule out cells which have extremely long (perhaps even infinitely long) "tentacles" which could make the cell unbounded or cause the cell to intersect a compact set very far from its center point (which may cause difficulties with local finiteness). This requires some basic control on how "spread out" the points of H λ h can be. We will prove much more quantitative estimates for cells in Sect. 4.

Proof of Lemma A.4.
Since μ h is a.s. a locally finite measure, a.s. we can find a D h -ball centered at any given point of P λ h which does not contain any other points of P λ h . Since D h induces the Euclidean topology on the domain of h a.s. each D h -ball contains a Euclidean neighborhood of its center point. Hence a.s. every cell in H λ h has non-empty interior. By the continuity of (z, w) → D h (z, w) with respect to the Euclidean topology (which follows from the fact that D h induces the Euclidean topology on its domain), we see that the cells of H λ h are a.s. closed. We will now argue that a.s. all of the cells of H λ h are bounded and that a.s. each compact subset of C intersects only finitely many cells of H λ h in the case when h is a whole-plane GFF normalized so that its circle average over ∂D is zero. We first claim that a.s. for each large enough k ∈ N, there is a point of P λ h in B 2 k (0)\B 2 k−1 (0). To see this, we first observe that by standard estimates for the Since t → h e t (0) evolves as a standard linear Brownian motion [DS11, Section 3.1], by the Borel-Cantelli lemma it follows that a.s. μ h B 2 k (0)\B 2 k−1 (0) grows exponentially in k as k → ∞. Since P λ h is a Poisson point process with intensity measure λμ h , our claim now follows.
For k ∈ N, let E k be the event that max z,w∈B 2 k (0)\B 2 k−1 (0) i.e., the D h -diameters of the two dyadic annuli on either side of B 2 k−1 (0)\B 2 k−2 (0) are each strictly smaller than the D h -distance across this annulus and the D h -distance across ∂ B 2 k−3 (0)\B 2 k−4 (0). We claim that a.s. E k occurs for arbitrarily large values of k ∈ N. It is easily seen from the local absolute continuity properties of the GFF (see, e.g., [GM16b, Lemma 4.2]) that P[E 1 ] > 0. By the scale invariance of the law of h, viewed modulo additive constant, we see that there is a universal constant p > 0 such that P[E k ] ≥ p for each k ∈ N. The event E k is determined by h| C\B 2 k−2 (0) , so since the tail σ -algebra r >0 σ h| C\B 2 k−2 (0) is trivial (see, e.g., [HS18, Lemma 2.2]), we obtain our claim.
Combining the preceding two paragraphs shows that a.s. there exists arbitrarily large values of k ∈ N for which E k occurs and B 2 k (0)\B 2 k−1 (0) and B 2 k−2 (0)\B 2 k−3 (0) each contain a point of P λ h . If k is one of these values, then each point of C\B 2 k−1 (0) is D hcloser to one of the points of P λ h in B 2 k (0)\B 2 k−1 (0) than it is to any point of B 2 k−2 (0), so cannot be contained in a cell centered at a point of B 2 k−2 (0). Hence every cell centered at a point of B 2 k−2 (0) is contained in B 2 k−1 (0), so in particular is bounded. Furthermore, each point of B 2 k−2 (0) is D h -closer to a point of P λ h which is in B 2 k−2 (0)\B 2 k−3 (0) than it is to any point of C\B 2 k−1 (0). Hence each cell which intersects B 2 k−2 (0) must be centered at a point in B 2 k−2 (0). Since μ h (B 2 k−2 (0)) < ∞, a.s. B 2 k−2 (0) contains at most finitely many points of P λ h , so B 2 k−3 (0) intersects only finitely many cells of H λ h . Since this happens for arbitrarily large values of k, we conclude the proof in the case of the whole-plane GFF.
The case when h is a quantum cone or a quantum wedge can be deduced from the case of a whole-plane GFF using local absolute continuity, or alternatively be treated similarly to the case of a whole-plane GFF using the radii R b of (2.8) with b = 2 k and the relation (2.9).
From basic properties of the Brownian disk, we get the following.

Lemma A.5. Suppose h is a whole-plane GFF (with some choice of additive constant). Almost surely, the boundary of each of the cells in H λ
h has zero √ 8/3-LQG area measure. In fact, it is a.s. the case that for any z, w ∈ P λ h , the set of u ∈ C for which D h (u, z) = D h (u, w) has zero √ 8/3-LQG measure. The same is true with h replaced by a freeboundary GFF on a Jordan domain or an embedding of a quantum cone, sphere, disk, or wedge.
Proof. We will prove the lemma in the case of the quantum disk. The general case follows from this and the local absolute continuity of the fields mentioned in the lemma with respect to an embedding of the quantum disk (this local absolute continuity holds away from the boundary of the domain in the case of the whole-plane GFF or the quantum cone or sphere and up the domain boundary in the case of a free-boundary GFF or quantum wedge). Since the quantum disk is equivalent to the Brownian disk, general Brownian disk theory shows that if v is sampled uniformly from μ h , then a.s. for each r > 0 one has μ h (∂ B r (v; D h )) = 0: indeed, this follows, e.g., from the fact that the set of times for which the head of the Brownian snake takes any particular value r > 0 has zero Lebesgue measure. Since P λ h is a Poisson point process with intensity measure λμ h , this shows that the probability that two points of P λ h lie on the boundary of the same D h -ball centered at v is zero. That is, the probability that D h (u, v) = D h (u, w) is zero. Since v is sampled uniformly from μ h , the statement of the lemma follows.
We next check that cell boundaries have zero Lebesgue measure. The basic idea of the proof is as follows. If φ is a smooth bump function and a ∈ R, then the laws of h +aφ and h are mutually absolutely continuous. Furthermore, a certain "good" generic event for h + aφ occurs for Lebesgue-a.e. choice of a. By taking a to be random according to a distribution with a density with respect to Lebesgue measure (e.g., sampled from the standard Gaussian distribution) this implies that the desired generic behavior holds with probability 1, which will then allow us to conclude that the cells have zero Lebesgue measure. Similar arguments can be used to obtain that cells of H λ a.s. have "generic" behavior in certain senses. Some results to this affect appeared in an earlier arXiv version of this paper. Lemma A.6. Suppose h is a whole-plane GFF (with some choice of additive constant). Almost surely, the boundary of each of the cells in H λ h has zero Lebesgue measure. In fact, it is a.s. the case that for any z, w ∈ P λ h , the set of u ∈ C for which D h (u, z) = D h (u, w) has zero Lebesgue measure. The same is true with h replaced by a free-boundary GFF on a Jordan domain or an embedding of a quantum cone, sphere, disk, or wedge.
Proof. We will prove the lemma in the case when h is a free-boundary GFF on D. This implies the lemma in general by local absolute continuity. Note that the lemma statement does not depend on the choice of additive constant for h by Remark A.1 and since adding a constant to h scales D h by a constant factor. If we condition on h and the total number of points in P λ h , then the elements of P λ h are i.i.d. samples from μ h , normalized to be a probability measure. It therefore suffices to show that if z, w ∈ D are independent samples from μ h , normalized to be a probability measure, then a.s. the set of u ∈ D with D h (u, z) = D h (u, w) has zero Lebesgue measure. For this purpose, it suffices to show that for any fixed u ∈ D, we have P[D h (u, z) = D h (u, w)] = 0.
Step 1: reducing to an event with deterministic sets. Henceforth fix u ∈ D and for open sets V ⊂ V ⊂ D with V ⊂ V and V = D, let E h = E h (u, V, V ) be the event that the following is true.
1. z ∈ V and D h (u, z) = D h (u, w). 2. There is a D h -geodesic from w to u which does not enter V .
Since D h -geodesics have zero μ h -mass, the probability that every D h -geodesic from z to u passes through w is zero, and the same is true with z and w interchanged. Consequently, on the event {D h (u, z) = D h (u, w)} there a.s. exists deterministic open sets V ⊂ V ⊂ D with V ⊂ V such that E h (u, V, V ) occurs, and we can take V and V to be finite unions of Euclidean balls with rational centers and radii.
It therefore suffices to show that for any fixed deterministic choice of V and V as above, one has P[E h ] = 0. Henceforth fix such a deterministic choice of V and V . Since adding a constant to h does not effect the occurrence of the event E h , we can assume without loss of generality that the additive constant for h is fixed so that μ h (U ) = 0 for some deterministic open set U ⊂ C which is disjoint from V .
Step 2: randomly perturbing the field. Consider a smooth bump function φ : D → [0, 1] with φ| V ≡ 1 and φ| D\V ≡ 0. For a ∈ R, the laws of the fields h + aφ and h are mutually absolutely continuous (this follows from, e.g., [MS16c,Lemma 3.4] and the fact that φ| U ≡ 0, so adding aφ does not affect the choice of normalization for the field). It follows from this that the joint laws of (h, z, w) and (h + aφ, z, w) are mutually absolutely continuous. We emphasize that this still holds even though we are always sampling z and w from μ h , rather than from μ h+aφ , since adding a smooth function to h results in an LQG measure which is absolutely continuous with respect to μ h .
Let E h+aφ be defined in the same manner as above but with (D h+aφ , z, w) in place of (D h , z, w). Also let A be a standard Gaussian random variable, independent from everything else. The laws of (h + Aφ, z, w) and (h, z, w) are mutually absolutely continuous, so it suffices to show that P[E h+Aφ ] = 0. For this purpose, it is enough to show that if E h+aφ occurs for some a ∈ R, then E h+bφ does not occur for any b ∈ R\{a} (since this implies that P[E h+Aφ |(h, z, w)] = 0).
Let us therefore suppose that E h+aφ occurs and b = a. We seek to show that E h+bφ does not occur, so we can assume that all of the conditions in the definition of E h+bφ occur except possibly for the condition that D h+bφ (u, z) = D h+bφ (u, w). We seek to show that D h+bφ (u, z) = D h+bφ (u, w). Since there is a D h+aφ -geodesic and a D h+bφgeodesic from w to u which do not enter V and φ is supported on V , we see that D h+aφ (u, w) = D h+bφ (u, w). On the other hand, since z ∈ V Lemma 2.3 implies that D h+aφ (u, z) < D h+bφ (u, z) if b > a, and one has the reverse inequality if b < a. This shows that D h+bφ (u, z) = D h+bφ (u, w), as required.