## Dimers on Riemann surfaces I:

## Temperleyan forests

Nathana¨ el Berestycki ∗

Benoit Laslier †

Gourab Ray ‡

July 24, 2024

## Abstract

This is the first article in a series of two papers in which we study the Temperleyan dimer model on an arbitrary bounded Riemann surface of finite topological type. The end goal of both papers is to prove the convergence of height fluctuations to a universal and conformally invariant scaling limit. In this part we show that the dimer model on the Temperleyan superposition of a graph embedded on the surface and its dual is well posed, provided that we remove an appropriate number of punctures. We further show that the resulting dimer configuration is in bijection with an object which we call Temperleyan forest, whose law is characterised in terms of a certain topological condition. Finally we discuss the relation between height differences and Temperleyan forest and prove that the convergence of the latter (which is the subject of the second paper in this series) implies the existence of a conformally invariant scaling limit of height fluctuations.

## Contents

| 1   | Introduction                              | Introduction                                                        |   2 |
|-----|-------------------------------------------|---------------------------------------------------------------------|-----|
|     | 1.1                                       | Background and main result . . . . . . . . . . . . . . . . . . .    |   2 |
|     | 1.2                                       | Relation with previous work . . . . . . . . . . . . . . . . . . .   |   3 |
|     | 1.3                                       | Overview of the proof and organisation of the paper . . . . .       |   5 |
| 2   | Background and setup                      | Background and setup                                                |   8 |
|     | 2.1                                       | Riemann surfaces and embedding . . . . . . . . . . . . . . . .      |   8 |
|     | 2.2                                       | Universal cover . . . . . . . . . . . . . . . . . . . . . . . . . . |   9 |
|     | 2.3                                       | Intrinsic and topological winding . . . . . . . . . . . . . . . .   |   9 |
|     | 2.4                                       | Discretization setup . . . . . . . . . . . . . . . . . . . . . . .  |  11 |
|     | 2.5                                       | Height function and forms . . . . . . . . . . . . . . . . . . . .   |  15 |
| 3   | Temperley's bijection on Riemann surfaces |                                                                     |  17 |
|     | 3.1                                       | Notion of Temperleyan cycle rooted spanning forest; bijection       |  17 |
|     | 3.2                                       | Criterion for a wired CRSF to be Temperleyan . . . . . . . .        |  21 |
|     | 3.3                                       | Wilson's algorithm to generate wired oriented CRSF . . . . .        |  22 |

∗ University of Vienna. Email: nathanael.berestycki@univie.ac.at

† Universit´ Paris Cit´ and Sorbonne Universit´, CNRS, Laboratoire de Probabilit´s, Statistique et Mod´lisation. e e e e e Email: laslier@lpsm.paris

‡ University of Victoria. Email:gourabray@uvic.ca

| 4   | Winding and height function              | Winding and height function                                       |   23 |
|-----|------------------------------------------|-------------------------------------------------------------------|------|
|     | 4.1                                      | Winding field of embedded trees and choice of reference flow      |   24 |
|     | 4.2                                      | Winding of CRSF and height function . . . . . . . . . . . .       |   29 |
|     | 4.3                                      | Mapping back to the manifold . . . . . . . . . . . . . . . . .    |   32 |
| 5   | Local coupling                           | Local coupling                                                    |   34 |
|     | 5.1                                      | Cutset exploration . . . . . . . . . . . . . . . . . . . . . . .  |   35 |
|     | 5.2                                      | Full coupling . . . . . . . . . . . . . . . . . . . . . . . . . . |   39 |
| 6   | Convergence of height function and forms | Convergence of height function and forms                          |   40 |
|     | 6.1                                      | Precise statement of the result . . . . . . . . . . . . . . . .   |   41 |
|     | 6.2                                      | Some a priori tail estimates on winding . . . . . . . . . . .     |   42 |
|     | 6.3                                      | From partial path to the full spine. . . . . . . . . . . . . . .  |   47 |
|     | 6.4                                      | Convergence of height function . . . . . . . . . . . . . . . .    |   50 |
| A   | Geometry of spines                       | Geometry of spines                                                |   56 |

## 1 Introduction

## 1.1 Background and main result

This paper is the first part of a series of two papers in which we show universality and conformal invariance of the Temperleyan dimer model on Riemann surfaces.

Given a finite graph G (with an even number of vertices) and nonnegative weights w = ( w e ) e ∈ E on the edges, a dimer configuration is a perfect matching of the vertices along the edges of G . The dimer model is the random perfect matching obtained by sampling a dimer configuration m on G with probability proportional to the product of the weights of the edges in the matching m :

$$\mathbb { P } ( \mathbf m ) = \frac { \prod _ { e \in \mathbf m } w _ { e } } { Z },$$

where Z is the partition function. If all the weights are equal to 1, P is simply a uniformly picked perfect matching. Following Thurston, a particularly convenient way to encode the dimer model on a planar bipartite graph is through a notion of height function. Informally, it is a real valued function defined on the faces of the graph, which is uniquely determined by the dimer configuration. This turns the dimer model into a model of discrete random surfaces. A central question in the study of the dimer model is to describe the fluctuations around the mean of the height function. It turns out that the height function of the dimer model is extremely sensitive to boundary conditions (see for example the quote of Kasteleyn at the beginning of [14]). In the so-called liquid phase , the Kenyon-Okounkov conjecture predicts that the fluctuations of the height function around its mean are asymptotically given by a Gaussian free field in a certain nontrivial conformal gauge (see [27] section 2.3 for details). Despite remarkable progress on this conjecture over the last decade, this remains an outstanding problem in all but a handful of cases: see [4, 19, 26, 34, 35] for examples of models where this is already shown, and see [7, 37] for an introduction to the Gaussian free field. We emphasise that, given the extreme sensitivity of the model to its boundary conditions, the occurrence of a universal limit such as the Gaussian free field is to some extent very surprising.

The goal of this paper is to describe the fluctuations of the height function for certain graphs embedded on Riemann surfaces of finite topological type, i.e., with finitely many handles and holes. Our assumptions on these graphs ensure that the dimer model is liquid on the entire Riemann surface (as opposed to having frozen or gaseous regions, say as in the Aztec diamond)

and the conformal gauge is simply the identity; we call the resulting dimer model Temperleyan , and we say that the sequence of graphs G # δ is a Temperleyan discretisation of the surface. Even though this assumption simplifies the problem, conformal invariance in the scaling limit is by no means obvious; for instance, in the simply connected setting and for the square lattice, this is in fact precisely the content of Kenyon's landmark paper [25].

Before we give a simplified statement of the main theorem, we point out that in the context of a Riemann surface M , the height function becomes a height one-form h # δ : that is, the gradient ∇ h # δ is well defined but this gradient does not derive itself from a single-valued function on the surface. Instead it can be viewed as the gradient of a (single-valued) function defined on the universal cover of the surface. For concreteness we will denote by h # δ the associated function on the universal cover of the surface. Readers who prefer it can however equivalently think of h # δ as the random one-form on M associated to the dimer configuration. (See Sections 2.5 and 6.1 for precise definitions.) The function h # δ is the central object of this sequence of papers, and we still refer to it as the height function.

Our main theorem, stated informally below, confirms that the dimer model admits a universal scaling limit on Riemann surfaces. This scaling limit is furthermore conformally invariant.

Theorem 1.1. Let M denote a bounded Riemann surface, possibly with a boundary, which is not the sphere or a simply connected domain. Let G # δ be a sequence of Temperleyan discretisations of M satisfying the Invariance principle and a crossing estimate for random walk described in Section 2.4. Let h # δ denote the height function of the dimer model on G # δ . For all smooth compactly supported test functions f on the universal cover of M ,

$$\int ( h ^ { \# \delta } - \mathbb { E } ( h ^ { \# \delta } ) ) f$$

converges as δ → 0 in law and the limit is conformally invariant and independent of the sequence G # δ .

The precise assumptions on the surface M will be given in Section 2.1, the meaning of Temperleyan discretisations and our precise assumptions on the sequence G # δ can be found in Section 2.4. Finally, note that the smooth test function f is taken on the universal cover ˜ M of M , so that the integration takes place on ˜ M (equipped e.g. with Lebesgue measure). See Sections 2.5 and 6.1 for details. We refer to Theorem 6.1 for the complete version of the theorem.

Although we cannot find what the limiting object is, as an artefact of our proof, we can show that the limit is non trivial. Indeed, we prove later in Proposition 5.10 that the limiting field restricted to a small enough neighbourhood of M away from the punctures can be coupled to the restriction of a Gaussian free field in a slightly bigger neighbourhood, with positive probability; in particular this field is nontrivial. We conjecture that the limiting field is a compactified Gaussian free field with appropriate parameters depending only on the surface M and the positions of the punctures ( x , . . . , x 1 k ) on M .

## 1.2 Relation with previous work

One of the first rigourous works on the dimer model on Riemann surfaces is the inspiring paper by Boutillier and de Tili`re [9], who were motivated by nonrigourous physics predictions based e on the Coulomb gas formalism (see e.g. [17] and the references in [9]). They described the topological part of the fluctuations in the dimer model on the honeycomb lattice on the torus. This was complemented some years later by the fundamental paper of Dub´dat [19], in which the e full law of the limiting height fluctuations was described for isoradial graphs on the flat torus, and identified as a so-called compactified Gaussian free field. Dub´ edat also conjectured a greater form of universality beyond the isoradial setting; such a universality statement was subsequently

proved by Dub´ edat and Gheissari [20] but only for the topological part of fluctuations on the torus. Note that in particular our Theorem 1.1 settles this conjecture of Dub´dat. e

We point out that all techniques in the paper mentioned above, based on Kasteleyn theory (in contrast with our approach), are very specific to the case of the torus and cannot easily be adapted to general Riemann surfaces. Cimasoni [10] shows that for general surfaces of genus g (and no boundary), the partition function of the dimer model becomes a weighted (in fact, signed) sum of 2 2 g determinants of Kasteleyn matrices in which the orientation has been reversed along some of the 2 g cycles forming a basis of the first homology group of the surface. The signs in front of the determinants themselves form a nontrivial algebraic topological invariant of the surface which Cimasoni, remarkably, was able to identify in terms of the so-called Arf invariant. (In the case of the torus and the honeycomb lattice, this goes back to the work of Boutillier and de Tili`re [9].) e This illustrates the difficulty in generalising the Kasteleyn approach to general surfaces. See also [1, 11, 12, 15] for related results.

The above difficulty illustrates the well known fact that the dimer models becomes 'less integrable' as the genus of the surface increases, or at any rate integrability becomes harder to exploit. In this paper, we will therefore follow the general strategy of our previous work [4], which relies instead on Temperley's bijection between the dimer model and uniform spanning trees (in fact, see the overview of the proof in Section 1.3 for a discussion of this strategy and the relation between the two papers). The increased robustness of this method accounts in large part for the fact that we are able to tackle generic Riemann surfaces (in contrast with methods based on Kasteleyn theory or interacting particle systems). This is not to say that the difficulties caused by the increased complexity of the surface vanish entirely, but they manifest themselves in a different manner, which boils down to the following. As we will describe in more details in Section 2, a Temperleyan discretisation of the surface is the graph resulting from the superposition of a graph Γ embedded on the surface and its dual. The resulting graph is bipartite, with its black vertices being given by the vertices and faces of Γ, and the white vertices being given by the edges of Γ. However, a straightforward calculation using Euler's formula shows that this does not admit a dimer covering unless we remove 2 g + b -2 = | χ | many white vertices from the graph (where g is the genus and b the number of boundary components, and χ is the Euler characteristic); these can be thought of as punctures in the surface, or monomers. Note that the number of punctures increases with the complexity of the surface, and when the surface is a torus or an annulus no such puncture needs to be removed. Roughly speaking, we will see that the locally tree-like random subgraphs appearing when we generalise Temperley's bijection, which we will call Temperleyan self-dual pair , have one topological constraint for each puncture. In the case of a torus or an annulus where there is no puncture, we will see that the Temperleyan self-dual pair boils down to the simpler notion of cycle rooted spanning forest studied by Kassel and Kenyon in [24].

The necessity to remove a certain number of punctures from the surface in order to formulate the Temperleyan dimer model has a geometric flavour which is reminiscent of Liouville conformal field theory , in which correlation functions are only well defined with the appropriate number of singularities (known in that context as insertions). The required number of insertions is specified by the so-called Seiberg bounds, see e.g. [16] on the sphere and [21] for Riemann surfaces (see also [7] for an introduction). As is the case here, the number of punctures is entirely specified by the topology of the surface, but in a way that is in some sense opposite to ours: indeed the required number of insertions is an increasing function of χ in Liouville CFT, whereas it is decreasing here.

This difference can be understood informally as follows. To first order, Liouville CFT may be viewed as yielding a random metric which is approximately hyperbolic (this can for instance be made rigourous in the semiclassical limit where the parameter γ → 0, see [29]). However, on the Riemann sphere, for which the theory is of particular interest, and on other surfaces of

non-negative Euler characteristic, it is impossible to equip the surface with a metric of constant negative scalar curvature without first adding conical singularities at a certain number of points, as follows from the Gauss-Bonnet theorem. By contrast, the Temperleyan dimer model is trying to equip the surface with a flat metric , a fact which will be particularly apparent when we relate the height function to the winding of the Temperleyan forest in Section 4 (note that the winding is computed on the universal cover, where we ignore the effects of curvature, see Section 2.2). Yet, on a hyperbolic surface (i.e., when χ &lt; 0) no flat metric exists, again by the Gauss-Bonnet theorem. In summary, Liouville CFT may be viewed as trying to build a hyperbolic metric on a surface which may be parabolic or elliptic, whereas the Temperleyan dimer model tries to impose a flat metric on a surface which may be hyperbolic. In both cases the difficulty is resolved by puncturing the surface at an appropriate number of singularities. (The above parallel between Liouville CFT and dimer model may perhaps also be compared with the two couplings between SLE and Gaussian free fields given respectively by the so-called forward and reverse couplings; the forward coupling describes Liouville conformal field theory whereas the reverse coupling describes imaginary geometry [32, 33]).

## 1.3 Overview of the proof and organisation of the paper

As mentioned above, our general approach to study height fluctuations on Riemann surfaces is motivated and inspired by our earlier work [4] connecting the dimer model to imaginary geometry in the simply connected setting. Let us recall the main idea of that work briefly in the original simply connected setting, before explaining how this relates with and differs from the case of Riemann surfaces.

- · The starting point of [4] is Temperley's bijection, which relates the dimer model on a Temperleyan discretisation of a simply connected domain to a pair of dual uniform spanning trees with respectively wired and free boundary conditions. Furthermore, the height function of the dimer model can be determined relatively simply from the associated pair of spanning trees: roughly speaking, the height difference between two points is simply the total winding (sum of turning angles) of the unique path connecting these points in either tree.
- · The second observation is that, as follows from the landmark paper of Lawler, Schramm and Werner [30], paths between points in the wired uniform spanning tree (UST) have a scaling limit which may be described in terms of Schramm-Loewner Evolutions with parameter κ = 2 or SLE 2 for short.
- · Finally, in the continuum, there is a coupling known as imaginary geometry between SLE κ curves and an appropriate multiple of the Gaussian free field. In this coupling, the values of the GFF may informally be identified as the winding of the associated SLE curves (note that this requires careful interpretation, since the Gaussian free field is not defined pointwise and SLE κ paths are not smooth). Nevertheless, this coupling can be thought of as a continuous form of Temperley's bijection in the case κ = 2.

The strategy of [4] is to exploit Temperley's bijection on the one hand, and imaginary geometry on the other hand, to show that the winding of paths in a UST is indeed given by the appropriate multiple of the Gaussian free field asymptotically (essentially, one has to show that we can exchange the order of taking limits, or that a diagram commutes). This strategy is very robust; in particular it does not appeal to the more solvable aspects of the dimer model. Indeed, the convergence of loop-erased random walks (which describe the branches of a wired UST at the discrete level) to SLE 2 is known in great generality.

Roughly speaking, we wish to implement a similar approach in the case of Riemann surfaces. However, it should be clear that each part will be substantially different in this new setup.

- 1. To begin with, we observe that Temperley's bijection is local and so may be applied on a Riemann surface to output a pair of random subgraphs that are dual to one another and locally tree-like, and which we call Temperleyan self-dual pair . While in the simply connected case, these are just dual spanning trees, the combinatorial structure (i.e. state space) of the Temperleyan self-dual pair on a Riemann surface is a priori not entirely clear, and its law even less so. Our first result will be a characterisation of the Temperleyan selfdual pair. We also describe its law explicitly, which has a simple and concise expression. It turns out that one of the two components, which we call the Temperleyan cycle rooted spanning forest or Temperleyan CRSF for short, determines the other (and hence the entire dimer configuration) almost uniquely (up to some choice of orienting cycles), and also has a tractable law. Essentially, a Temperleyan CRSF is a random spanning subgraph in which every cycle is non-contractible, and satisfying a certain further topological condition. The introduction and study of Temperleyan CRSFs take up a substantial part in both the articles in this series.
- 2. A second problem is that the connection between height function and winding is not as straightforward as in the simply connected case; in particular, in addition to the already mentioned fact that only the gradient of the height function is well defined (leading to a multivalued height function), there may not be a path connecting two given points in the Temperleyan CRSF, even when they are neighbours. We thus need to develop a systematic method for computing the height function given the Temperleyan CRSF.
- 3. The next difficulty is to obtain a scaling limit result for our Temperleyan CRSF. Let us immediately note here that there are specific difficulties with even defining variants of SLE on surfaces since the simply connected uniformisation theorem lies at the heart of the standard definition of SLE.
- 4. Finally if we were to follow the blueprint of [4] one would also need to establish a version of imaginary geometry for Riemann surfaces, an endeavour which seems out of reach with current techniques. Fortunately, we are in fact able to entirely bypass this step because it turns out that the (discrete) 'error' estimates needed anyway to exchange limits are strong enough to imply convergence regardless of any information in the continuum. Clearly, the downside of this approach that we cannot hope to describe the law of the limit in this manner.

The topological and geometric structures of the surface M raise substantial specific difficulties at each of the steps of this roadmap. In this first paper we develop the relevant discrete/combinatorial arguments corresponding to the first two items above and present the conclusion according to the last item. The second paper in this sequence, [5], is devoted to the third bullet point, i.e the proof of convergence of the Temperleyan CRSF in the scaling limit. This is the most involved step which requires several new ideas, but can be read essentially independently.

At this point, it is also useful to discuss a bit more precisely the role that the different structures associated to a Riemann surface will play in our paper. In order to even discuss convergence on the surface M one needs in principle to specify a Riemannian metric on M , and it is natural to require that this metric be consistent with the conformal structure of M . However, ultimately, we want to prove convergence to a conformally invariant object so in principle only the conformal class of this metric matters (here we mean that two metric tensors are equivalent if we obtain one from the other by multiplying by a smooth function e φ .)

It turns out however that, after lifting to the universal cover (which in most of the interesting cases for this paper will be the unit disc, see below for more details), it will be more convenient for us to equip this cover with the standard Euclidean metric instead of the a priori more natural hyperbolic metric.

This impacts the way we formulate some of our assumptions in Section 2.4. Unfortunately, some of our technical estimates involve indirectly the Euclidean metric on the unit disc: for instance, through our assumptions that the graph G has a bounded density and each edge has a bounded winding (items i and ii in Section 2.4). Nevertheless, changing the metric to another one in the same class only changes the constants involved in these assumptions, and so does not impact the main results. We explicitly prove this fact in Lemma 2.5.

## Organisation of the paper

- · In Section 2, we recall some basic facts about dimers on surfaces and in Section 2.3, we recall the notions of windings of curves which were developed in our previous article [4]. Sections 2.1 and 2.2 contain small reminders about Riemann surfaces. We introduce precisely our discrete setup in Section 2.4. As one of the main differences with the classical case is that the height function now becomes multivalued we recall in Section 2.5 the language of height one-forms on graphs which is a classical way to handle such multivalued functions. In Section 2.5, we introduce another way to deal with multivalued functions by lifting them to single-valued functions on the universal cover. We also illustrate how an application of the Hodge decomposition theorem allows us to decompose the (multivalued) height function on the surface into a single-valued function (scalar component) and a canonical representative of the 'multivalued (or topological) part', the so-called instanton component.
- · In Section 3, we introduce the generalisation of Temperley's bijection to Riemann surfaces, essentially covering point 1 above. Section 3.1 contains the bijection itself. We also provide in Section 3.2 a simple criterion for a CRSF to be Temerleyan which is the basis of the second paper in the series [5].
- · In Section 4, we carefully develop the relation between the height differences and the winding of Temperleyan forests, i.e point 2 of the overview. This is in the spirit of the work of Kenyon-Propp-Wilson [28], who treated the (planar) simply connected case. In that case of course the bijection is with a uniform spanning tree. Compared to that work there are two additional points that we need to handle. The first one is that the edges of the graph cannot be assumed to be straight lines as in [28]. The second and more significant one is that there are additional terms coming from the fact that the forest is not connected: given points x and y on the universal cover, the height difference h x ( ) -h y ( ) (which is unambiguously defined) is essentially given by the intrinsic winding of any path connecting x to y , plus additional discontinuities every time the path jumps between components. The resulting key formula is stated in Theorem 4.10 (see also Lemma 4.6).
- · In Section 5 we extend our local and global coupling results from [4] to the framework of Riemann surfaces. This coupling is a key ingredient of our approach, which allows us to show that given a finite number of points ( z i ) 1 ≤ ≤ i k on the surface (or, rather, their lifts to the universal cover), the respective geometries of the CRSF in neighbourhoods of these k points are almost independent. This is crucial because certain types of 'error terms' in our approach vanish in expectation but are large in absolute value.
- · In Section 6, we conclude by proving our main convergence result: Theorem 6.1. Sections 6.2 and 6.3 contain some technical a priori estimates on winding of the CRSF branches on the universal cover. Finally in Section 6.4, we finish the proof of Theorem 6.1.
- · Finally, appendix A contains some geometric facts about spines in the CRSF (that is, the lifts to the universal cover of a cycle of the CRSF). It is given as a appendix because the proofs relies on slightly more involved elements of the classical theory of Riemann surfaces than the rest of the paper and we wanted to avoid overloading the already lengthy setup section.

## Acknowledgements

The authors are grateful to Antoine Dahlqvist and Julien Dub´dat for a number of useful discuse sions. NB's is supported by FWF grant P33083, 'Scaling limits in random conformal geometry'. GR's research is supported by NSERC 50311-57400. BL's research is supported by ANR-18CE40-0033 'Dimers'. The paper was finished while NB was in residence at the Mathematical Sciences Research Institute in Berkeley, California, during the Spring 2022 semester on Analysis and Geometry of Random Spaces , which was supported by the National Science Foundation under Grant No. DMS-1928930. Finally, we thank the anonymous referee for valuable feedback on an earlier version of the manuscript.

## 2 Background and setup

## 2.1 Riemann surfaces and embedding

In this article, we work with a Riemann surface M (a connected, one dimensional complex manifold) satisfying the following properties.

- · M is of finite topological type, meaning that the fundamental group π 1 ( M ) is finitely generated. In other words, we assume that the surface has finitely many 'handles' and 'holes'.
- · M can be compactified by specifying a boundary ∂M . We denote by M the compactified Riemann surface with the boundary. More precisely, every point is either in the interior and hence has a local chart homeomorphic to C or is on the boundary and has a local chart homeomorphic to the closed upper half plane ¯ . H Also there are finitely many such charts which cover the boundary. Note that this condition implies that M has no punctures. (However for future reference, we note here that we will later introduce punctures on M which will align with the removed vertices of Temperleyan graphs. The resulting punctured manifold will be denoted by M ′ .)

We say M is nice if M satisfies the above properties.

Given such a Riemann surface M , it is possible to equip it with a Riemannian metric g (with associated distance function denoted by d M ) which extends continuously to the boundary, and turns M into a smooth Riemannian surface. This Riemannian metric can be constructed by considering the hyperbolic metric on the universal cover of the surface if it has no boundary, see [23, Section 2.4] in combination with the uniformisation theorem ([23, Theorem 4.4.1], also recalled and discussed below). If the surface has boundary components, we may apply the same result to the surface obtained by gluing it to itself along the boundary - also called the double.

Note that the choice of the metric above exclude surfaces such as the hyperbolic plane. This will simplify certain topological issues later when we deal with the Schramm topology.

Classification of surfaces. Riemann surfaces can be classified into the following classes depending on their conformal type (see e.g. [18] for an account of the classical theory):

- · Elliptic: this class consists only of the Riemann sphere, i.e., M ≡ ˆ C
- · Parabolic: this class includes the torus, i.e., M ≡ T = C / ( Z + τ Z ) where ℑ ( τ &gt; 0), the cylinder M ≡ C \{ } 0 , and the complex plane itself (or the Riemann sphere minus a point), M ≡ C .
- · Hyperbolic: this class contains everything else. This includes examples such as the twotorus, the annulus, as well as proper simply connected domains in the complex plane, etc.

The proofs in this paper (and its companion [5]) are concerned with the hyperbolic case (subject to the above conditions) as well as the case of the torus. So from now on we always assume that M is such a manifold. We note that the case of simply connected proper domain in C is covered in our previous work [4].

## 2.2 Universal cover

The universal cover ˜ M of the Riemann surface M will play an important role in our analysis. Recall that any Riemann surface admits a (regular) covering map by a covering space which is simply connected. This covering space may furthermore be endowed with a conformal structure, with respect to which the covering map is then analytic. Being simply connected, this covering space is (by the uniformisation theorem, see Theorem 4.4.1 and Theorem 2.4.3 in [23]) conformally equivalent to a Riemann surface ˜ , which is either the unit disc M D ⊂ C (hyperbolic case), the whole plane (parabolic case), or the Riemann sphere ˆ (elliptic case). C Altogether we obtain a map 1 p : ˜ M → M which is both analytic and a regular covering map of M . By definition of a covering map, p is a local homeomorphism: for every z ∈ M , there exists a neighbourhood N containing z so that p is injective in every component of p -1 ( N ).

Recall further that in each of the three cases (elliptic, parabolic, hyperbolic), ˜ M can be endowed with a metric (compatible with the complex structure) of constant Gaussian curvature equal to 1 0 and , -1 respectively. Furthermore, there is a subgroup F of conformal automorphisms on ˜ M acting properly discontinuously and freely - i.e., without fixed points - (see, e.g., [23, Section 2.4] for definitions), such that p descends to a conformal equivalence between ˜ M/F and M . Put it more simply, F is a discrete subgroup of the corresponding set of M¨bius transo formations describing the conformal automorphisms of ˜ , M and M is conformally equivalent to ˜ M/F (i.e., to either ˆ C /F , C /F or D /F ).

In case of the torus, this discrete subgroup is isomorphic to Z 2 and the generators specify translations in the two directions of the torus. In the hyperbolic case, this class of subgroups is much bigger and are known as Fuchsian groups (see e.g. [18] for a general account). This particular representation of a hyperbolic Riemann surface, sometimes called a Fuschian model, will be particularly convenient because it allows us to describe the scaling limits of Temperleryan or cycle rooted spanning forests in the universal cover , rather than on the surface itself. This will allow us to import directly a number of the ideas and result from [4] on the simply connected case.

Concretely, one of the main advantages of the Fuchsian model of a hyperbolic Riemann surface - and more generally of a lift to the universal cover which inherits the Euclidean metric - is that, equipped with this Euclidean metric , we have a flat embedding in which the sum of angles around a point is 2 π , the total winding of a simple loop is 2 π (see also Lemma 2.1 below). The Fuchsian structure itself plays a more technical role, see for instance appendix A. We rely on this for a canonical description of the height function in terms of the extended Temperley bijection.

To conclude, moving to the universal cover allows us in a sense to trade a natural law on a graph possessing a nontrivial topology, for a law with complicated periodicity-like constraints but on a planar graph.

## 2.3 Intrinsic and topological winding

The goal of this section is to recall several notions of windings of curves drawn in the plane, which we use in this paper. We refer to [4] for a more detailed exposition. A self-avoiding (or simple) curve in C is an injective continuous map γ : [0 , T ] → C for some T ∈ [0 , ∞ ].

1 For concreteness one can fix a point x 0 in M and consider p such that p (0) = x 0 . In the hyperbolic case fix an arbitrary choice of rotation as well.

The topological winding of such a curve γ around a point p / γ ∈ [0 , T ] is defined as follows. We first write

$$\gamma ( t ) - p = r ( t ) e ^ { i \theta ( t ) },$$

where the function θ : [0 , ∞ → ) R is taken to be continuous, which means that it is unique modulo a global additive constant multiple of 2 π . We define the winding for an interval of time [ s, t ], denoted W γ s,t , p ( [ ] ), to be

$$W ( \gamma [ s, t ], p ) = \theta ( t ) - \theta ( s )$$

(note that this is uniquely defined). Notice that if the curve has a derivative at an endpoint of γ , we can take p to be this endpoint by defining

$$W ( \gamma [ 0, t ], \gamma ( 0 ) ) \coloneqq \theta ( t ) - \lim _ { s \to 0 } \theta ( s )$$

and similarly

$$W ( \gamma [ s, T ], \gamma ( T ) ) \coloneqq \lim _ { t \to T } \theta ( t ) - \theta ( s ).$$

With this definition, winding is additive: for any 0 ≤ s ≤ t ≤ T

$$W ( \gamma [ 0, t ], p ) & = W ( \gamma [ 0, s ], p ) + W ( \gamma [ s, t ], p ). & ( 2. 2 )$$

The notion of intrinsic winding we describe now, also discussed in [4], is perhaps a more natural definition of windings of curves. This notion is the continuous analogue of the discrete winding of non backtracking paths in Z 2 which can be defined just by the number of anticlockwise turns minus the number of clockwise turns. Notice that we do not need to specify a reference point with respect to which we calculate the winding, hence our name 'intrinsic' for this notion.

̸

We call a curve smooth if the map γ is smooth (continuously differentiable). Suppose γ is smooth and for all t, γ ′ ( ) = 0. t We write γ ′ ( ) = t r int ( ) t e iθ int ( ) t where again θ int : [0 , ∞ → ) R is taken to be continuous. Then define the intrinsic winding in the interval [ s, t ] to be

$$W _ { \text{int} } ( \gamma [ s, t ] ) \coloneqq \theta _ { \text{int} } ( t ) - \theta _ { \text{int} } ( s ).$$

The total intrinsic winding is again defined to be lim t → T W int ( γ, [0 , t ]) provided this limit exists. Note that this definition does not depend on the parametrisation of γ (except for the assumption of non-zero derivative). The following topological lemma from [4] connects the intrinsic and topological windings for smooth curves.

Lemma 2.1 (Lemma 2.1 in [4]) . Let γ [0 , 1] be a smooth self-avoiding curve in C then,

$$W _ { \text{int} } ( \gamma [ 0, 1 ] ) = W ( \gamma [ 0, 1 ], \gamma ( 1 ) ) + W ( \gamma [ 0, 1 ], \gamma ( 0 ) ).$$

Note that the right hand side in the above expression makes sense as soon as the endpoints of γ are smooth. From now on we therefore extend the definition of intrinsic winding to all such curves using the equality of Lemma 2.1. Since the topological winding is clearly a continuous function on the curve away from the point p , this extension is consistent with any regularisation procedure. We also recall the following deformation lemma from [4] (see Remark 2.5).

Lemma 2.2. Let D be a domain and γ ⊂ ¯ D a simple smooth curve (or piecewise smooth with smooth endpoints). Let ψ be a conformal map on D and let arg ψ ′ ( D ) be any realisation of argument on ψ ′ ( D ) . Then

$$W _ { \text{int} } ( \psi ( \gamma ) ) - W _ { \text{int} } ( \gamma ) = \arg _ { \psi ^ { \prime } ( D ) } ( \psi ^ { \prime } ( \gamma ( 1 ) ) ) - \arg _ { \psi ^ { \prime } ( D ) } ( \psi ^ { \prime } ( \gamma ( 0 ) ) ).$$

## 2.4 Discretization setup

We will require the graphs we are working with to be embedded on M in a way that preserves its topology and geometry, we make this precise now.

Let M be a Riemann surface with b holes and g handles satisfying the assumptions of Section 2.1. We say that a graph Γ is faithfully embedded in M if it satisfies the following (see Figure 1):

- · the embedding is proper, i.e. edges do not cross;
- · Γ is connected;
- · Every face of Γ not intersecting ∂M has the topology of a disc.
- · Γ has a marked vertex (called a boundary vertex, the set of which is denoted by ∂ Γ) for each component of ∂M ;

Note that with this convention, we can think of each boundary vertex as being 'delocalised' along the whole boundary component and indeed later we will consider cycle rooted spanning forests (CRSFs) with a wired condition on each boundary component. This 'delocalised' boundary vertex is best illustrated to be sitting inside the hole, see Figure 1.

Let Γ † denote the dual graph of Γ, defined as follows. Every face of Γ contains a single vertex of Γ † and two vertices of Γ † are connected by an edge if their corresponding faces share an edge of Γ. We also assume the following for Γ : †

- · The embedding is proper, i.e., the edges do not cross.
- · Γ † is connected.
- · Every face of Γ † corresponding to a non-boundary vertex of Γ is a topological disc.
- · The vertices of Γ † corresponding to the faces of Γ which contain a given boundary vertex form a simple cycle homotopic to the corresponding boundary component of ∂M . We call this the boundary cycle , and we denote this by ∂ Γ . †

Figure 1: The graphs Γ Γ , † , G ˆ , and G . The graph Γ is pictured in red and Γ † in black.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00832v3-figure-e1c80e48-5f75-49f6-b94e-9bbf2c2d296c.jpg)

We assign an oriented weight w e ( ) = w Γ ( e ) ≥ 0 to every oriented edge of Γ. This defines a continuous time Markov chain on Γ where for every edge ( x, y ) in Γ, the chain jumps from x to y at rate w x, y ( ), and with a slight abuse of terminology we will from now on refer to this chain as the random walk on Γ. Note that we defined it in continuous time for simplicity but we will in fact always be interested only by the geometry of its paths so its time parametrisation is mostly irrelevant. With a slight abuse of notations, we will often identify the trajectory of the random walk (or other processes such as loop-erased paths) with the continuous path obtained from the set of edges used by the walk, either seen as a path or as a subset of M . Similarly, with

a slight abuse of notation, we will use direct image notations to denote sub-paths of the random walk, i.e if s &lt; t are two times, we will denote by X s, t [ ] the path followed by the random walk between times s and t . Finally we will always stop the random walk whenever it reaches a boundary vertex so we assume that all oriented edges starting at such a vertex have weight 0. We emphasize that there are no weights on the dual graph Γ † and indeed we will never consider the random walk on the dual graph.

Remark 2.3. To explain the asymmetry between Γ and Γ , we point out that Γ can be thought † of as a graph that is wired at the boundary of ∂M whereas Γ † has free boundary conditions.

In order to apply locally the bijection between dimers and (a variant of) spanning trees, we need to define the superposition of Γ and Γ . † We introduce a new set of vertices W at each point where some edge e and dual edge e † intersect. Define the graph ˆ G whose vertices are V = V (Γ) ∪ V (Γ ) † ∪ W and whose edges are obtained by joining every vertex in W which is on the edge e with the endpoints of e and e † (in other words, each pair of edges e and e † corresponds to four edges in ˆ ). G We call ˆ the G superposition graph . Then we let G denote the graph where we remove from ˆ the boundary vertices G ∂ Γ and its corresponding half edges. See Figure 1 for an illustration. Clearly, G is a bipartite graph where we take W to be the set of white vertices and the rest to be black. Also every non-outer face of G is bounded by 4 edges (i.e. every non-outer face is a quadrangle). We can naturally define un oriented weights on G from the ones on Γ. If w is the white vertex in the middle of the edge ( x, y ) of Γ, we set w G ( x, w ) = w Γ ( x, y ) and w G ( y, w ) = w Γ ( y, x ) while w G ( e ) = 1 if e is part of an edge of Γ . † It will be useful later for the definition of the height (see Section 4.1) to also choose in each face f of ˆ a G diagonal d f ( ) which is just a curve connecting the primal and dual vertices adjacent to f without exiting f , together with a midpoint m f ( ) which is just a point in the interior of the diagonal. If the primal vertex adjacent to f is a boundary point choose the diagonal as a path connecting the boundary and the dual vertex. It will be clear from Section 4.1 that the choice of these diagonals and midpoints only affects some of the arbitrary conventions needed to define the height function and that any statement on the centered height h -E ( h ) is independent of this choice.

We will prove in Section 3 that if the manifold M is not a torus or an annulus, one must introduce punctures in the graph G for it to even admit dimer coverings. More precisely we will remove finitely many white vertices from G and call a removed vertex and the edges associated with it a puncture (equivalently, the puncture may be thought of as forcing the vertex to be a monomer instead of in a dimer edge). Note that each puncture creates an octagon G consisting of four white vertices, two black vertices from Γ and two black vertices from Γ † (see Figure 2). We will call the resulting graphs respectively G , ′ Γ ′ , (Γ ) † ′ and we call the manifold obtained by removing the white vertices M ′ . We will see later in Section 3 that the number of punctures we remove must be k = | χ | where χ is the Euler characteristic of M . When this is the case we will call G ′ a Temperleyan graph on M or Temperleyan discretisation of M .

Figure 2: Removing a white vertex (puncture) to create the dimer graph G ′ .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00832v3-figure-066596fb-b263-47f3-bd0e-59879b5f23c8.jpg)

When considering scaling limits, we will of course consider a sequence of graphs Γ # δ , (Γ ) † # δ , G # δ from the setup above. We now introduce our assumptions on such a sequence. Let p : ˜ M → M be a map which is both analytic and a regular covering of M by its universal

cover ˜ , M which is either the unit disc D ⊂ C or the whole plane C (recall the discussion in Section 2.2 for the existence of this map.) In the end the choice of this covering map does not affect the following assumptions, see Lemma 2.5 for a precise discussion about this covering map and also about the covering map for the punctured manifold.

Recall from Section 2.1 that we equipped M with a distance function d M derived from a Riemannian metric, which extends to the closure ¯ ; M we still denote this extended distance function by d M .

- (i) (Bounded density) There exists a constant C independent of δ such that for any x ∈ M , the number of vertices of Γ # δ in the ball { z ∈ M d : M ( x, z ) &lt; δ } is smaller than C .
- (ii) (Good embedding) The edges and diagonals of the graph are embedded as smooth curves and for every compact set K ⊂ ˜ , M the intrinsic winding of every edge or diagonal in the lift ˜ Γ # δ intersecting K is bounded by a constant C = C K depending only on K . (Note that this allows edges to wind quite a bit near holes.)
- (iii) (Invariance principle) As δ → 0, the continuous time random walk { ˜ X t } t ≥ 0 on Γ ˜ # δ started from a nearest vertex to 0 satisfies:

$$( \tilde { X } _ { t / \delta ^ { 2 } } ) _ { t \geq 0 } \xrightarrow { ( d ) } ( B _ { \phi ( t ) } ) _ { t \geq 0 },$$

where ( B , t t ≥ 0) is a two dimensional standard Brownian motion in ˜ M (killed when it leaves ˜ , M if ˜ M = D ) started from 0, and ϕ is a nondecreasing, continuous, possibly random function satisfying ϕ (0) = 0 and ϕ ( ∞ ) = ∞ . The above convergence holds in law in Skorokhod topology.

We remark that the above condition is equivalent to asserting that simple random walk from some fixed vertex converges to Brownian motion on the Riemann surface itself up to time parametrisation (see e.g. [22]).

Figure 3: An illustration of the crossing condition.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00832v3-figure-d37b52a6-2976-4502-9eff-d8d0951b7d3d.jpg)

- (iv) (Uniform crossing estimate). Let R be the horizontal rectangle [0 , . 3] × [0 , . 1] and R ′ be the vertical rectangle [0 , . 1] × [0 , . 3]. Let B 1 := B (( 05 . , . 05) 0 025) be the , . starting ball and B 2 := B (( 25 . , . 05) , . 025) be the target ball (see Figure 3).

The uniform crossing condition is the following. There exist universal constants δ 0 , α 0 &gt; 0 such that for every compact set K ⊂ M , there exists a δ K such that for all δ ∈ (0 , δ K ) the following is true. Let ˜ K = p -1 ( K ) be the lifts of K . Let R ′′ be a set of the form cR + , z where c ≥ δ/δ 0 and z ∈ R 2 (i.e. a scaling and translate of R ), and is such that R ′′ ⊂ ˜ . K Let B ′′ 1 = cB 1 + z and B ′′ 2 = cB 2 + . For all z v ∈ ˜ Γ # δ ∩ B ′′ 1 ,

$$\mathbb { P } _ { v } ( \tilde { X } \text{ hits } B _ { 2 } ^ { \prime \prime } \text{ before exiting } R ^ { \prime \prime } ) > \alpha _ { 0 }.$$

We emphasise that this crossing condition is defined in the Euclidean metric in the disc, and not the (perhaps more standard) hyperbolic metric when dealing with the universal cover of a hyperbolic surface.

In what follows, sometimes for a compact set S ⊂ ˜ , M we will write δ S to mean δ p S ( ) as defined above.

- (v) (punctures) We remove k = | χ | many points from (the interior of) M and call the resulting surface M ′ , where χ is the Euler characteristic of M . We assume that each of the ( k many) monomers (white vertices removed from ˆ G to obtain G ′ ) converge to a unique puncture as δ → 0 which is a given point in the interior of the Riemann surface. We also assume without loss of generality that for all δ , the punctures are sufficiently far (say at graph distance at least some large fixed constant) from the boundary.

Finally, we also suppose that the following holds. Let u , v i i be the endpoints of the edge of Γ corresponding to the white monomer removed for 1 ≤ i ≤ k . Then there exist paths γ u i , γ v i (viewed as a ordered collection of adjacent edges (( x , x 1 2 ) , ( x , x 2 3 ) , . . . , ( x r -1 , x r )) or as a continuous path depending on context) satisfying the following. The paths start with x 1 = u i or x 1 = v i respectively for γ u i , γ v i , and we assume

$$M \, \bigvee _ { 1 \leq i \leq k } \left ( \gamma _ { u _ { i } } \cup \gamma _ { v _ { i } } \cup \left ( u _ { i }, v _ { i } \right ) \right )$$

is topologically a disjoint union of annuli.

Figure 4: An illustration of the assumption in Item v. The paths γ u 1 , γ v 1 consists of p and q , which decompose the surface into a number of disconnected annuli: three in the first case, two in the second. The dotted edges are the ones removed from G to get G ′ .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00832v3-figure-a0a787a3-abb7-4d0e-95ac-8cc784fa0fc2.jpg)

Remark 2.4. Typically the existence of the paths γ u i and γ v i for δ small enough is guaranteed by the crossing condition and the pants decomposition of M . In words, the paths decompose the surface into disjoint annuli, see Figure 4. However, it is possible to construct pathological examples (for large δ ) where there is no paths such that (2.5) holds, so this needs to be made as a separate assumption. In fact, we will see that (v) is equivalent to the graph G admitting a dimer cover (see Figure 1), so we could equivalently assume the latter instead - this would be more elegant but less concrete.

Lemma 2.5 (Conformal invariance of assumptions) . The assumptions (ii), (iii), (iv) on the sequence (Γ # δ ) δ&gt; 0 are invariant with respect to the choice of the covering map p , in the sense that if p ′ is any other analytic and regular covering of M then these assumptions are satisfied for some (possibly different) choices of δ 0 , α 0 , ( δ K K ) ⊂ M , ( C K K ) ⊂ ˜ M .

Furthermore, if (Γ # δ ) δ satisfies the assumptions (i)-(v) for a choice of covering map p and constants C, δ 0 , α 0 , ( δ K K ) ⊂ M , ( C K K ) ⊂ ˜ M then these assumptions are invariant with respect to conformal transformations in the sense that if ψ : ( M,x ,x ,... , x 1 2 k ) → ( N,x ,x ,. . . , x ′ 1 ′ 2 ′ k ) is a conformal bijection mapping M to N and x i to x ′ i for 1 ≤ i ≤ k , then ( ψ (Γ # δ )) δ&gt; 0 also satisfies assumptions (i)-(v) for the choice of the covering map p ψ ◦ and constants C, δ 0 , α 0 , ( δ ψ -1 ( K ) ) K ⊂ M ′ , and ( C K K ) ⊂ ˜ M .

Proof. We start with the second of these two assertions, for which items (ii)-(v) are trivial. For (i), if M has no boundary then d M ( x, y ) = d N ( ψ x , ψ y ( ) ( )) (essentially because M¨bius maps o

are isometries with respect to the hyperbolic metric on the unit disc). Thus (i) holds. If M has a boundary, then note that the conformal bijection ψ between M and N extends to a conformal bijection ¯ between the doubles of ψ M,N ; and we are back to the previous case.

Let us now prove the first assertion. Note that by the uniformisation theorem, there exists a conformal bijection (M¨bius map) o ϕ : ˜ M → ˜ M such that p ◦ ϕ = p ′ . Note ϕ is a M¨bius map o from the unit disc to itself in the hyperbolic case and from the complex plane to itself in the torus case. Assumption (ii) is trivially preserved in the torus case as M¨bius maps are translations o and rotations. For the hyperbolic case, we can employ the change in winding under conformal map formula (Lemma 2.2), and note that derivatives of M¨bius maps are bounded on compact o subsets of the open unit disc. In assumption (iii), we only require convergence of the random walk up to time change which is preserved under M¨bius maps. o Assumption (iv) is easily seen to be preserved in the torus case as a M¨bius map maps a rectangle to another rotated and o translated rectangle, which can be crossed by concatenating bounded number of vertical and horizontal rectangles of smaller scales. In the hyperbolic case, a rectangle is mapped to a domain bounded by four circular arcs, and the starting and target discs are mapped to discs inside this domain. Furthermore, compact sets are mapped to compact sets. Thus concatenating domains of this type, it is easy to see that uniform crossing estimate holds for M , perhaps for different positive constants δ 0 , α 0 , ( δ K K ) ⊂ M .

Remark 2.6. The Invariance principle assumption (item iii) actually implies something stronger in combination with the other assumptions: for any point x in ˜ , the random walk started from M a vertex x # δ nearest to x converges to a Brownian motion started from x up to a time change as above. This is a consequence of the fact that random walk from 0 comes close to x with uniformly positive probability using the crossing estimate and the strong Markov property of Brownian motion.

̸

In case ∂M = ∅ , recall that the set of boundary cycles ( ∂ Γ ) † # δ corresponds to the connected components of ∂M . One consequence of the Invariance principle assumption (item iii) is that each boundary cycle converges in the Hausdorff metric (induced by d M ) to the associated component of ∂M .

Sometimes we drop the superscript δ from Γ # δ , (Γ ) † # δ for clarity, when there is no possibility of a confusion.

Remark 2.7. In the hyperbolic case, note also that while we have stated the assumptions (ii), (iii) and (iv) on the universal cover of M , these assumptions are also valid for the sequence (Γ # δ ) δ&gt; 0 on the universal cover of the punctured surface M ′ as well if we endow it with the metric it inherits from M via inclusion. This can be checked using the fact that there is a map from the universal cover of M ′ to ˜ M p \ -1 ( { x , . . . , x 1 k } ) (where x , . . . , x 1 k is the set of punctures) which is locally a conformal bijection. Also note that (i) is also satisfied in M ′ trivially by the choice of the metric.

## 2.5 Height function and forms

A flow ω is a real valued antisymmetric function on the oriented edges ⃗ E of G , i.e., for every oriented edge ( u, v ), ω u, v ( ) = -ω v, u ( ). The total flow out of a vertex v is defined to be ∑ w ∼ v ω v, w ( ). Similarly, the total flow into a vertex v is defined to be ∑ w ∼ v ω w,v ( ). A flow f is a closed 1-form if the sum over any oriented contractible cycle is 0: i.e., for any oriented cycle ( v , v 0 1 , . . . , v n = v 0 ) in G so that the embedding of ∪ n -1 i =0 ( v , v i i +1 ) in M forms a contractible loop,

$$\sum _ { i = 0 } ^ { n - 1 } \omega ( v _ { i }, v _ { i + 1 } ) = 0.$$

It is clear that if M is simply connected, then there exists a function f on the vertices of G (uniquely defined up to a global constant) such that f ( v ) -f ( u ) = ω u, v ( ) for all u ∼ v .

We now associate to any dimer configuration m on G a closed 1-form on ⃗ E . Let m be a dimer configuration on G , and let ⃗ e = ( w,b ) be an oriented edge, where w is a white vertex and b a black vertex. We define the flow ω m by setting ω m ( ⃗ e ) = 1 { ∈ e m } . Also, ω m is defined in an antisymmetric way: ω m (( b, w )) = -ω m (( w,b )). Note that the total flow out of a white vertex is 1 and that out of a black vertex is -1.

To any flow ω on oriented edges, one can associated a dual flow ω † defined on the oriented edges of the dual graph G † , where if e † crosses the edge e = ( w,b ) with w on its right and b on its left, then we set ω † ( e † ) = ω e ( ). Note also that if ω is divergence free (i.e., the flow out of every vertex is 0), then ω † is a closed 1-form on ⃗ E † .

Consider any reference flow ω 0 which has total flow out of white vertex equal to 1 and total flow out of a black vertex equal to -1. Then ω = ω m -ω 0 defines a divergence free flow on ⃗ E . We call ω † the height 1-form corresponding to m with reference flow ω 0 .

↦

When G is embedded on a simply connected domain (so that, in particular, no cycle in G † is non-contractible), every closed 1-form ω on ⃗ E † becomes exact: i.e., there exists a function on the faces F G ( ) of G h , : F G ( ) → R , so that for any two adjacent faces f, f ′ ,

$$h ( f ^ { \prime } ) - h ( f ) = \omega ( f, f ^ { \prime } ).$$

Observe further that this function is defined only up to a global constant. The function h is then called the height function of the dimer m , admitting an abuse of terminology.

We recall the following simple but useful observation. A path in G (or G † ) is a sequence of vertices ( v , . . . , v 0 n ) (or faces ( f , . . . , f 0 n )) of G so that v i is adjacent to v i +1 (or f i is adjacent to f i +1 in G † ) for all 0 ≤ i ≤ n -1.

Lemma 2.8 (Unique path lifting property) . Let γ = ( f , f 0 1 , . . . , f n ) be a path (not necessarily simple) in G † . Let ˜ f 0 be a lift (i.e. one pre-image) of f 0 to ˜ M . Then there exists a unique path ˜ = ( ˜ γ f , f 0 ˜ 1 , . . . , f ˜ ) n in ˜ M so that ˜ f i is the lift of f i . Further, ∑ n -1 i =0 ω f , f ( ˜ i ˜ i +1 ) = h f ( ˜ n ) -h f ( ˜ ) 0 .

We now turn to the definition of height function in the more complicated case when M is no longer assumed to be simply connected. In that case, when we sum the values of the height 1-form (defined above) along any non-contractible cycle, we may get a nonzero value. One can use the Hodge decomposition theorem , to isolate out the part of the height 1-form which is encoded by the topology of the underlying surface. The Hodge decomposition theorem works in great generality, but in the present context, it takes the following simple form. For any function f on the vertices of G we define df to be the closed 1-form defined on ⃗ E as

$$d f ( u, v ) = f ( v ) - f ( u ).$$

A harmonic 1-form h is a closed 1-form which is divergence free, so that ∑ v ∼ u h ( u, v ) = 0.

Theorem 2.9 (Hodge decomposition [2, 8, 31]) . For any closed 1-form ω on G (or G † ), there exist a function f on the vertices of G and a harmonic 1-form h defined on ⃗ E such that

$$\omega = d f + \mathfrak { h },$$

and f is unique up to an additive global constant, and h is unique. Furthermore, h is completely determined by summing ω over a finite set of oriented non-contractible cycles which forms the basis of the first homology group of M .

In this paper, we will analyze ( f, h ) corresponding to the divergence free flow ω m -ω 0 , where m is dimer configuration chosen from the law (1.1) subject to certain natural conditions and ω 0 is a carefully chosen reference flow (in fact, we will consider the height function see Section 6

for a precise statement). We will call h the instanton component . We remark that changing the reference flow changes ( f, h ) by a deterministic additive factor, and in particular does not affect the fluctuations of ( f, h ) around their mean. To be more precise, our main theorem (Theorem 6.1) will be stated in terms of the single-valued function associated to ω m -ω 0 on the universal cover of M ′ . See Remark 6.16 for a statement concerning both the scalar and instanton parts of the height function.

Throughout the paper, rather than working with the scalar and instanton components of the height 1-form, it will be more convenient to lift the height 1-form ω to the universal cover of M . Since the latter is always simply connected, this allows us to work with actual functions without having to worry about the Hodge decomposition Theorem 2.9. We will then check that the convergence of height function on the universal cover implies convergence of each of the components in the Hodge decomposition.

Our assumptions on the graph G where the dimer model lives are such that ˜ = G p -1 ( G ) is a planar graph embedded on ˜ . M Moreover, the height 1-form ω on the dual edges of G lifts to a height 1-form ˜ ω on the dual edges of ˜ . G Since ˜ M is simply connected, and since ˜ ω is a closed one-form on the dual edges of ˜ G (this is a local property, so remains true when we lift to ˜ ), G we can define a height function h = h ( m , G ) (up to a global constant) on the dual graph ˜ G † . The instanton component h can be related to the height function h on the universal cover by summing up the value of ˜ along any path in the dual graph of ω ˜ corresponding to a G non-contractible loop in the dual edges of G . This is easier to explain on an example.

Example 2.10. If M is the flat torus T := C / ( Z + τ Z ) for some complex number τ with ℑ ( τ ) &gt; 0, then the universal cover is the complex plane C . The universal cover can be thought of as many copies of the fundamental domain (a parallelogram determined by 1 and τ ).

Fix v 0 in the fundamental domain. Then by periodicity of dh , the height function h on ˜ G evaluated at a point v = v 0 + m + τn (where m,n ∈ Z ) is given by

$$h ( v ) = h ( v _ { 0 } + m + i n ) = h ( v _ { 0 } ) + a m + \mathfrak { b } n,$$

for some a b , ∈ R which do not depend on either v , m, n 0 . Let us describe what a b , are. Consider the two loops on the torus described by L 1 := ( +1 t / 2 τ : t ∈ [0 , 1]) and L 2 := (1 / 2+ tτ : t ∈ [0 , 1]) in the fundamental domain (i.e., L 1 and L 2 are the two non-contractible loops in the torus which form the basis of the homology group). Then a is the sum of the values of ω along any loop in ⃗ E † which is homotopic to L 1 , whereas b is the sum of the values of ω along any loop which is homotopic to L 2 . Clearly, the choice of these curves in G † do not matter since the height 1-form is closed. Furthermore, in the Hodge decomposition of Theorem 2.9, the harmonic 1-form h is uniquely determined by the numbers a and b .

## 3 Temperley's bijection on Riemann surfaces

## 3.1 Notion of Temperleyan cycle rooted spanning forest; bijection

The goal of this section is to extend the classical Temperley bijection to graphs embedded in surfaces of higher genera. The bijection is between dimer covers of certain natural class of graphs introduced in Section 2.4 and certain objects which we call Temperleyan forests . The main result is formulated in Theorem 3.3.

Before defining the objects in question, let us briefly recall the classical Temperley bijection on simply connected surfaces [14]. The bijection map is local: given a dimer cover, one 'extends' (in the direction from black to white) the dimer to an oriented edge as depicted in Figure 6. This yields two oriented trees which are dual to each other. In the opposite direction, given a pair of trees which are dual to each other, one can naturally orient them towards a pre-assigned 'root vertex'; the dimer cover then consists of the first half of each oriented edge in the two dual

trees. The primary issue that needs to be dealt with in higher genus is that the relevant objects are no longer pairs of trees that are dual to each other, but pairs ( t, t † ) of arborescences (i.e., collection of oriented edges such that there is a single outgoing edge out of every vertex, except for boundary vertices) on Γ and Γ † respectively, dual to each other, and where each component contains a non-contractible oriented cycle. Thinking only about the arborescence t on the primal graph, which we aim to ultimately understand via a variant of Wilson's algorithm, the requirement that its planar dual t † can be endowed with the orientation of an arborescence imposes nontrivial topological constraints on t which will be elucidated below.

Let Γ be a graph faithfully embedded on a surface with a certain specified set of boundary vertices. Assume every edge e of Γ comes with a specified weight w e ( ) &gt; 0. Before introducing the notion of Temperleyan forest, we start with the simpler notion of cycle rooted Spanning Forest.

Definition 3.1. A wired oriented cycle rooted spanning forest (which we abbreviate: wired oriented CRSF) of Γ with the specified boundary is an oriented subgraph t of Γ where

- · Every non-boundary vertex of Γ has exactly one outgoing edge in t . Every boundary vertex has no outgoing edge. (As a result, any cycle of t must be consistently orientated).
- · Every cycle of t is non-contractible.

To each CRSF t , we may associate the weight ∏ e ∈ t w e ( ) . We will consider the probability measure on CRSFs such that the probability of t is by definition proportional to its weight.

This is equivalent to the notion of essential CRSF on a graph with wired boundary introduced by Kassel and Kenyon [24]. Ignoring the orientation of t gives an unoriented graph, its connected components will simply be called the connected components of t without any additional precision. Note that if t is a wired oriented CRSF, every connected component of t contains at most one cycle: more precisely, every boundary component must have zero cycles, while every non-boundary component contains exactly one cycle.

We will refer to the set of all non-contractible cycles of a wired oriented CRSF to mean the set of unique cycles corresponding to each (non-boundary) component of the wired oriented CRSF.

Let us come back to the setup of Section 2.4 and recall that we had the graph Γ, its dual Γ † and the superposition graph G all embedded nicely in a Riemann surfaces with g handles and b holes. Compared to the planar setting, there is an immediate topological difficulty, which is that this graph G in general does not admit a dimer cover. Indeed by Euler's formula, if G has v vertices, e edges and f many contractible faces, then v -e + f = χ := 2 -2 g -b (where χ is the Euler's characteristic). On the other hand if G admits a dimer cover, one must have e = v + f since e is the number of white vertices which must match v + f (the number of vertices and contractible faces of Γ ). † Thus we need to remove k = | χ | many edges from the superposition graph for it to have any chance of having a dimer cover (see Figure 2). These removed edges can be thought of as creating punctures in the surface. Call the graph with these punctures G ′ . (Recall that, as per item v, we assume throughout that the punctures are macroscopically far apart in M , and are away from the boundary as well, and converge as δ → 0 to a set of pairwise distinct points in M ). Note that if M is a torus or an annulus then χ = 0 so no punctures need to be removed. See also Ciucu-Krattenthaler [13] and Dub´dat [19] for other situations where e punctured dimers arise.

To every wired oriented CRSF t of Γ ′ with boundary ∂ Γ, one we can associate a natural dual free cycle rooted spanning forest t † (abbreviated free CRSF) of (Γ ) † ′ as follows. The vertices of t † are given by the vertices of (Γ ) † ′ (i.e., it spans (Γ ) ) and an edge † ′ e † is present in t † if and only if its dual e is absent in t . Note that a priori t † does not come with an orientation and that its cycles can overlap (or equivalently a component might contain several cycles), see Figure 5 for an example. In fact, it is not hard to see that in cases similar to Figure 5, it is

Figure 5: A non-Temperleyan CRSF (in blue). The surface M is the 'pair of pants': a domain of the plane with two holes (in grey in the picture). In this example, t does not contain any cycle, and hence any connected component flows to the boundary of M . Its dual t † must contain a component with two cycles which overlap. The cycles go around each of the two holes, and must be connected, as otherwise there would have to be a path in t separating them; however this is impossible as such a path would have to connect two distinct boundary points. So t is not Temperleyan. Note that χ = 2 -2 g -b = -1 here. See Lemma 3.5 for a more general argument.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00832v3-figure-c1de5850-d507-4710-a72e-a65e0b6f10e4.jpg)

not possible to consistently orient the free CRSF, so there is no chance to apply any version of Temperley's bijection. This motivates the introduction of the following definition.

Definition 3.2. We say that the wired oriented CRSF t is Temperleyan if each connected component of t † contains exactly one cycle.

An example of a wired oriented CRSF t that is not Temperleyan is provided in Figure 5. Temperleyan forests also come with a natural law, which is simply the law of the CRSF conditioned on the event that each component of the dual contains exactly one cycle. We will soon formulate a more explicit criterion for a CRSF t to be Temperleyan (see Section 3.2 and in particular Theorem 3.7). For now, observe that if t is a Temperleyan wired oriented CRSF and t † is its dual, we can assign an orientation to each cycle in each component of t † arbitrarily from one of the two possible choices. Then we orient all other edges of t † towards the cycle of that component. We call ( t, t † ) a Temperleyan pair or Temperleyan self-dual pair if t is a Temperleyan CRSF and t † is its dual with a choice of orientation as above. We warn the reader that the notation slightly obscures the choice of orientation needed to go from t to t † in order to highlight the duality. Hopefully this should not introduce any confusion since we will always want to think of our CRSF as being oriented.

We now state the Temperley bijection for general surfaces. Recall that we assign oriented weights w e to each edge e in Γ , no weight (or unit weight) to edges of Γ ′ † and that this turns G ′ into a weighted un oriented graph. Indeed, if e = ( x, y ) is an oriented edge of Γ , let ′ w denote the white vertex in the middle of e . Then we assign to the unoriented edge { x, w } of G ′ the weight w ( x,y ) and to the edge { w,y } of G ′ the weight w ( y,x ) .

We define the measure P Temp as the measure on Temperleyan pairs with weights inherited from the weights of Γ and we will call ( T , T † ) a generic associated random variable, i.e

$$\mathbb { P } _ { \text{Temp} } ( ( \mathcal { T }, \mathcal { T } ^ { \dagger } ) = ( t, t ^ { \dagger } ) ) = \frac { 1 } { Z _ { \text{Temp} } } 1 _ { \{ ( t, t ^ { \dagger } ) \text{ Temperleyan pair} \} } \prod _ { e \in t } w ( e ), \quad \quad ( 3. 1 )$$

where Z Temp is the partition function.

Figure 6: Illustration of the local transformation in Temperley's bijection. Left: a dimer configuration on a superposition graph. Right: a collection of oriented edges forming dual oriented CRSF.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00832v3-figure-a431784a-3245-4243-99a1-8d3ad772204e.jpg)

Theorem 3.3 (Temperley bijection on general surfaces) . Let M, Γ ′ , (Γ ) † ′ , G ′ be as in Section 2.4. Then there exists a bijection ψ between the set of Temperleyan pairs and the set of dimer configurations on G ′ . Furthermore if ( T , T † ) has the law (3.1) then m = ψ (( T , T † )) has law (1.1) with unoriented weights on G ′ described above.

We emphasise that it is unclear at this point whether the measure Z Temp P Temp (and consequently the measure on dimers) is not trivially zero as we have only deduced that the construction of G ′ is necessary, but have not deduced that it is sufficient. We prove this later in Proposition 3.4.

Proof of Theorem 3.3. We assume that the set of Temperleyan CRSF and the set of dimer covers are both nonempty (this assumption is validated in Proposition 3.4). Given a Temperleyan pair ( t, t † ), we obtain a configuration of edges m = ψ (( t, t † )) as follows: for every oriented edge ⃗ e ∈ t (resp. ⃗ e ∈ t † ), we can write ⃗ e = e 1 ∪ e 2 where e , e 1 2 are the first and second halves of ⃗ e and set e 1 ∈ m (see Figure 6). The resulting m is a matching on G ′ because every (non boundary) vertex has a unique outgoing edge in either t or t † . Furthermore, since t ∪ t † spans the black vertices of G ′ , the matching is a perfect matching.

Also ψ is injective: if ( t 1 , t † 1 ) and ( t 2 , t † 2 ) are distinct, then there must a black vertex v on G ′ (i.e., a vertex of Γ ′ or (Γ ) ) such that the unique outgoing edge from † ′ v in t 1 or t † 1 is different from the unique outgoing edge from v in t 2 or t † 2 . Hence v will be matched to two distinct white vertices in ψ (( t 1 , t † 1 )) and ψ (( t 2 , t † 2 )).

We now check ψ is onto. Given a matching m of G ′ , we can obtain a pair ( t, t † ) by extending the matched edges: for every non-boundary black vertex b of G ′ , let w be the white vertex matched with b in m . By construction ( bw ) is a part of a unique edge e = ( bb ′ ) (either of Γ of Γ ) and we include the edge ( † bb ′ ) (oriented from b to b ′ ) in t or t † as appropriate. The fact that neither t nor t † contain contractible cycles follows from the same argument as in the standard, planar case: if say t has a contractible cycle C , then since Γ \ C has a dimer cover, an elementary counting argument then shows that v -e + f = 0 where v, e, f are the number of vertices, edges and faces of the contractible component Γ \ C . On the other hand, Euler's formula implies v -e + f = 1 (again excluding the outer face). This shows that all cycles are non-contractible. Since every vertex v of t † has a unique outgoing edge, t † must have exactly

one cycle per component and so t is Temperleyan. It also easily follows from the definition of the weights of the Temperleyan tree and the dimer covers that the bijection is weight preserving. This concludes the proof.

## 3.2 Criterion for a wired CRSF to be Temperleyan

In this Section, we prove the following proposition:

Proposition 3.4. The graph G ′ obtained in Section 2.4 has a dimer cover. In particular, the measure P Temp defined in eq. (3.1) is a probability measure.

As a by product of the proof of Proposition 3.4, we derive a simple criterion for a CRSF to be Temperleyan (Theorem 3.7). We start with a lemma about the χ = 0 case.

Lemma 3.5. Let M be nice with g handles and b boundary components and Γ Γ , † , G be embedded as above (i.e. without punctures). A wired Temperleyan oriented CRSF of Γ exists if and only if M has the topology of either a torus or an annulus. Furthermore in these cases, all oriented CRSF are Temperleyan.

Proof. Note that if M has the topology of an annulus (with all the nice properties of Section 2.1), then finding a Temperleyan oriented CRSF is straightforward. Indeed, any wired spanning forest in the annulus (where both boundaries are wired) is Temperleyan: the dual is a graph containing a single cycle separating the two components touching each boundary. Also notice that for a torus, every oriented CRSF is Temperleyan [20]: essentially an oriented CRSF must contain a cycle (since there is no boundary on a torus) and cutting along this cycle gives us a (bounded) cylinder or equivalently an annulus.

For the converse, observe that the extended Tempereley's bijection (Theorem 3.3) allows us to construct a dimer configuration from the Temperleyan CRSF and its dual (by endowing each cycle with arbitrary orientation and orienting every other edge towards the unique cycle of its component). Notice that although Theorem 3.3 is written for the punctured graph G ′ , since in the case of torus and annulus G = G ′ , the same proof goes through. Now recall that a dimer configuration exists without removing any punctures only if χ = 0 or equivalently 2 g + b = 2 (as discussed in Section 3.1). This equation has only two feasible solutions: g = 1 , b = 0 (i.e. a torus) and g = 0 , b = 2 (i.e. an annulus). This completes the proof of the 'only if' part.

For the last assertion, note that every oriented CRSF in the annulus divides the annulus into several smaller disjoint annuli. The rest follows easily from the arguments in the first paragraph above.

Proof of Proposition 3.4. In light of Lemma 3.5, we assume M is neither a torus, nor an annulus. Using assumption Item v, recall that there exist a collection of paths γ u i , γ v i on Γ ′ for 1 ≤ i ≤ k starting respectively from u i and v i , such that

$$M \, \sum _ { 1 \leq i \leq k } \left ( \gamma _ { u _ { i } } \cup \gamma _ { v _ { i } } \cup \left ( u _ { i }, v _ { i } \right ) \right )$$

is topologically a disjoint union of annuli when we view the paths as continuous curves on M and ( u , v i i ) denotes the edge from u i to v i . We start building a dimer configuration m on G ′ as follows. First, for each path γ in this collection, which by definition starts from either u i or v i for some 1 ≤ i ≤ k , we orient it away from the puncture from which it emanates, i.e. away from u i or v i . We then consider the collection of dimer edges naturally associated with γ formed by the first halves of each of the oriented edges of γ , as in Figure 6 and the proof of Theorem 3.3. Then include all these edges in m .

The remaining disjoint annuli are nice in the sense of Section 2.1, so we can apply Lemma 3.5 to obtain a Temperleyan forest on each of these annuli, and hence by Theorem 3.3 a dimer cover

on the part of G ′ corresponding to these annuli. Patching these together gives us a dimer configuration m on all of G ′ as desired.

Remark 3.6. Note that this implies that the assumption in Item v is, as already mentioned earlier, necessary in order for G ′ to admit a dimer cover. Indeed if there is a dimer cover then by applying the local mapping in Figure 6, starting from the vertices u i and v i adjacent to the punctures, we can obtain paths as described in Item v. Each component of the complement of these paths also has a dimer cover, but no puncture. According to Lemma 3.5, this is only possible if (2.5) is satisfied.

We now deduce from the above an extremely convenient criterion for a wired oriented CRSF to be Temperleyan. Let us suppose M is not the torus or an annulus, whence k = 2 g + b -2 &gt; 0 (we already know that every CRSF is Temperleyan otherwise). Define the branch starting from a primal vertex v of Γ ′ to be the path obtained by going along the unique outgoing edge from each vertex (which necessarily ends when a loop is formed or a boundary is hit). Recall that, at a puncture (i.e., a white monomer), there are exactly two vertices u , v i i of Γ ′ on either side of the puncture. Let B i 1 , B i 2 be the branches in Γ ′ of u , v i i for 1 ≤ i ≤ k = 2 g + b -2. We call the union of these branches, together with the removed edges e i , the skeleton s of t . That is,

$$\mathfrak { s } = \bigcup _ { i = 1 } ^ { k } ( \mathfrak { B } _ { i 1 } \cup e _ { i } \cup \mathfrak { B } _ { i 2 } ).$$

Simply put, the skeleton of t is the union of the branches emanating out of the punctures.

Theorem 3.7. Suppose M is not a torus. A wired oriented CRSF is Temperleyan if and only if every component of M \ s has the topology of an annulus.

(Of course, in the case of an annulus, as already noted there are no punctures and so s = , ∅ so the condition is automatically satisfied, as we already know.) The criterion for a CRSF to be Temperleyan is thus just to say that the skeleton cuts the surface into disjoint annuli. Figure 4 gives two examples on a surface M (the 'pair of pants') where M \ s consists of topological annuli, where s is the paths p, q in that figure. (We get three annuli in the first example, and two in the second example).

Proof. Let t be a Temperleyan CRSF and let t † be its dual with a choice of orientation. Note that the vertices (in G ′ ) of s are all matched with each other in the dimer configuration associated to ( t, t † ). Therefore in each component of M \ s , all vertices are also matched with each other. By Lemma 3.5, this implies that these components are annuli.

Conversely, note that by definition t † cannot cross s , and so t can be restricted to each component of M \ s to form a wired CRSF. By the other implication in Lemma 3.5, it follows that t is Temperleyan in each such component and so is globally Temperleyan.

## 3.3 Wilson's algorithm to generate wired oriented CRSF

Recall the measure P Temp from eq. (3.1), i.e. the law of ( T , T † ) on Temperleyan pairs which is naturally associated with the dimer model. We will not study directly P Temp but rather a slightly altered version which can be sampled through Wilson's algorithm and is defined as follows: we first sample a wired oriented CRSF of Γ with law

$$\mathbb { P } _ { \text{Wils} } ( \mathcal { T } = t ) = \frac { 1 } { Z _ { \text{Wils} } } 1 _ { \{ t \text{ Temperleyan} \} } \prod _ { e \in t } w ( e ) ;$$

and given T = , we pick t T † an oriented dual uniformly among all possibilities of orientation of the dual.

Clearly the only difference between the two is due to the fact that any cycle of the dual t † of a Temperleyan oriented CRSF t can be oriented in two possible ways to determine a dual pair ( t, t † ). We deduce the following relationship:

Lemma 3.8. Let ( t, t † ) be a Temperleyan pair such that t † contains exactly k † non-contractible cycles. Then the Radon-Nikodym derivative satisfies

$$\frac { d \mathbb { P } _ { \text{Temp} } } { d \mathbb { P } _ { \text{Wils} } } ( t, t ^ { \dagger } ) = \frac { Z _ { \text{Wils} } } { Z _ { \text{Temp} } } 2 ^ { k ^ { \dagger } }$$

In particular, conditioned on having k † non-contractible cycles for the dual forest, the law P Temp and P Wils coincide.

Let Γ Γ , † be faithfully embedded on a nice Riemann surface M . We now describe Wilson's algorithm to generate a wired (but not necessarily Temperleyan) oriented CRSF on Γ. We prescribe an ordering of the vertices ( v , v 0 1 , . . . ) of Γ.

- · We start from v 0 and perform a loop-erased random walk until a non-contractible cycle is created or a boundary vertex (i.e., a vertex in ∂ Γ) is hit.
- · We start from the next vertex in the ordering which is not included in what we sampled so far and start a loop-erased random walk from it. We stop if we create a non-contractible cycle or hit the part of vertices we have sampled before.

There is a natural orientation of the subgraph created since from every non-boundary vertex there is exactly one outgoing edge through which the loop erased walk exits a vertex after visiting it. Let ˜ P Wils be the law of the resulting wired oriented CRSF.

## Proposition 3.9. We have

$$\tilde { \mathbb { P } } _ { \text{Wils} } ( t ) = \frac { 1 } { \tilde { Z } _ { \text{Wils} } } \prod _ { e \in t } w ( e ).$$

In particular, ˜ P Wils generates a wired oriented CRSF of Γ as described by Definition 3.1 with law given by (3.2) . Furthermore, conditionally on being Temperleyan, ˜ P Wils coincides with the first marginal of P Wils .

$$P r o o f. \text{ This follows from Theorem 1 and Remark 2 of Kassel-Kenyon [24].}$$

As a consequence of the result of this section, we obtain a relatively usable description of the law P Temp of interest: start with ˜ P Wils , which can be sampled via Wilson's algorithm, and condition on the topological event described in Theorem 3.7 (in the second paper [5] we denote this conditioning event by A M ). Finally bias the law by the Radon-Nikodym derivative coming from Lemma 3.8. One of the main difficulties in [5] stems from the fact that A M is asymptotically degenerate (i.e., the probability of A M for a wired CRSF tends to zero as the 'mesh size' δ tends to zero). The Radon-Nikodym derivative in Lemma 3.8 on the other hand is relatively benign and does not introduce much further complication.

## 4 Winding and height function

In this section, we explain the connecting between winding of the Temperleyan forest and height function in our setup. In the process we develop a systematic way to relate these two notions in a very general setting which we believe is of interest even in simply connected domains.

A first difficulty is that measuring angles and computing the total angle along say a simple loop is well known to be more complicated on a surface than on the plane. In the following, we will sidestep completely this issue by only establishing the connection for the lifts of the Temperleyan forest and the height function where we will be able to use the usual planar theory

of angles and winding. This will also have the advantage that the dimer height become an actual function instead of a one form.

Using this approach, we are left with a situation analogous to the original generalised Temperley's bijection in [28]: there is a fixed graph with a given planar embedding and we have to relate the dimer height to the winding of the associated 'tree' at a deterministic level. (We also point out that a version on a torus already appeared in [38].) In particular, the fact that in our current setting the graphs (and trees) arise as lifts and therefore have many symmetries, should not play a significant role (recall that [28] is valid in great generality).

There are however still significant technicalities to work out. First, we cannot assume that our embedding uses only straight lines (since in the hyperbolic case, the lift of an edge cannot be a straight line simultaneously in different copies of the fundamental domain). Second, our 'tree' objects are actually forests so the winding between vertices belonging to different components must be properly defined. Third, we need to check that the objects defined on the lift can be reasonably mapped back to the surface and state any simplification of the general theory that would arise specifically in our setting.

Finally, the way the connection is phrased in [28] is relatively unwieldy so it would be nice to make it more 'user-friendly', for example to make it symmetric with respect to the primal and dual forests. Recall that the winding field is naturally defined on either primal or dual vertices while the dimer height function is naturally defined on faces of the superposition graph, so we need some natural local rule to move between faces and vertices. Our solution for this technicality will be to extend both the primal and dual trees along the diagonal of each face, up to the midpoint of that face, a construction which results in what we call augmented trees. We will see that with this construction, the primal and dual trees play a completely symmetric role, as desired.

Our strategy in this section will be to proceed by successive generalisations from the definition of a winding field on a single deterministic embedded tree to dual pairs of spanning tree and the relation to dimer height function in Section 4.1. Section 4.2 covers the situation specific to our Riemann surface setting, with the main result of this section stated as Theorem 4.10. Finally Section 4.3 treats the relation between the instanton component and the set of non-contractible loops.

## 4.1 Winding field of embedded trees and choice of reference flow

Embedded tree. We temporarily forget about dimers and Temperleyan forests, and focus on how to compute the winding field of a deterministic tree.

Consider a graph (finite or infinite) which is a tree oriented towards either a leaf or one of its ends. Embed it on C with smooth edges (although having points of non-differentiability at the vertices, which we call corners , is allowed, so that the branches are only piecewise smooth). Call the embedding T , we want to think of it as a union of smooth curves in C . Recall that, using Lemma 2.1, the intrinsic winding of any finite self avoiding path in this tree is well defined. From every point x on T , we let γ x be the unique path from x following the orientation (so it either stops at a leaf or at the end towards which the tree is oriented) and call it the branch of x . Suppose for now that

$$| W _ { \text{int} } ( \gamma _ { x } ) | < \infty \text{ for all } x \in \mathcal { T }.$$

We can then define a winding field { h T ( x ) : x ∈ T } of the tree simply as h T ( x ) := W int ( γ x ).

The first generalisation step is to extend the definition of h T even if (4.1) is not satisfied which follows from the following elementary lemma. If x / γ ∈ y and y / γ ∈ x , notice that γ x and γ y eventually merge and γ y merges either to the right or to the left of γ x (since the tree is embedded in C , this makes sense). Let γ xy be the unique path connecting x and y in T .

Lemma 4.1. In the above setup, if γ y merges with γ x to its right then

$$h _ { \mathcal { T } } ( x ) - h _ { \mathcal { T } } ( y ) = W _ { \text{int} } ( \gamma _ { x y } ) + \pi.$$

If γ y merges with γ x to its left, then

$$h _ { \mathcal { T } } ( x ) - h _ { \mathcal { T } } ( y ) = W _ { \text{int} } ( \gamma _ { x y } ) - \pi.$$

If y ∈ γ x and y is not a corner (or x ∈ γ y ( ) and x is not a corner), then

$$h _ { \mathcal { T } } ( x ) - h _ { \mathcal { T } } ( y ) = W _ { \text{int} } ( \gamma _ { x y } ).$$

Proof. Notice that the last assertion follows simply from additivity of intrinsic winding. Indeed, for example if y ∈ γ x , there is no discontinuity of intrinsic winding at y since y is not a corner. For the rest, take m ∈ γ x ∩ γ y . Notice that by additivity of winding,

$$h _ { \mathcal { T } } ( x ) - h _ { \mathcal { T } } ( y ) & = W _ { \text{int} } ( \gamma _ { x m } ) - W _ { \text{int} } ( \gamma _ { y m } ) \\ & = W _ { \text{int} } ( \gamma _ { x m } ) + W _ { \text{int} } ( \gamma _ { m y } ) \\ & = W _ { \text{int} } ( \gamma _ { x y } ) + \varepsilon \pi$$

where ε = +1 (or ε = -1) if γ y merge with γ x to its right (or left). This is clear since we need to do a half turn to move from γ xm to γ my at m and the turn is clockwise if γ y merges to the right of γ x and anticlockwise otherwise.

We remark that the last assertion of the above lemma could be made to work even at corners by adding the appropriate angle, however in what follows, we avoid using winding field at corners so we don't need to introduce these additional difficulties. Finally the formulas for h x ( ) -h y ( ) described in Lemma 4.1 can be taken to be the definition of the winding field (or rather its gradient) for any tree embedded with piecewise smooth edges, even if (4.1) is not satisfied. This will be the typical situation for our setup.

Dual pairs and dimers. Take an infinite planar graph Γ, properly embedded in C with smooth edges and suppose that the embedding is locally finite in the sense that every vertex and face has finite degree, every face is a topological disc, and any bounded set intersects only finitely many edges. Let T be a one ended spanning tree of Γ (we emphasise that we are still considering a spanning tree for the moment and not yet a forest). Since Γ is locally finite, T admits a dual T † which is also a one-ended spanning tree of the dual graph Γ . † Let G be the superposition graph (in C ), as introduced in Section 2.4 (with the obvious modifications to take into account that our graph is infinite and without creating any puncture). It will be useful to augment the trees T and T † as follows. Recall that in Section 2.4, for each face f of G , we fixed a diagonal d f ( ) which is just a smooth simple curve connecting the primal and dual vertex staying inside the face and a midpoint m f ( ) which was just a point in d f ( ). For each face, add to T (resp. T † ) the portion of the diagonal d f ( ) connecting the point m f ( ) to the unique primal (dual) vertex touching d f ( ). This way the primal and dual trees meet in each face f at the (smooth) point m f ( ) on the diagonal d f ( ) of that face. With a small abuse of notation, we will still denote by T and T † these augmented trees.

Note that we can orient T , T † towards their respective ends. This allows us to define two winding fields as above (one for T and one for T † ). Having oriented T and T † we can also apply the local operation of Figure 6 as described in Section 3.1 to obtain a dimer covering of G . We will first need to define a suitable reference flow on G , which will then allow us to speak of the height function associated to the dimer configuration and then show the relation between this dimer height function and the two above winding fields.

Figure 7: The definition of the height function in terms of winding.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00832v3-figure-8d4126d7-474f-4c1d-94f9-1f727d56ad47.jpg)

Definition 4.2 (reference flow) . Let w,b be two adjacent white and black vertices and let f , f l r be the faces to the left and right of the oriented edge ( bw ) . Define

$$\omega _ { r e f } ( w b ) = \frac { 1 } { 2 \pi } \left ( W _ { \text{int} } \left ( ( m ( f _ { r } ), b ) \cup ( b, m ( f _ { l } ) ) \right ) + \pi \right ),$$

where ( m f ( l ) , b ) is the portion of d f ( l ) joining m f ( l ) and b , and similarly ( m f ( r ) , b ) . Define ω ref ( bw ) = -ω ref ( wb ) (see Figure 7).

Figure 8: Proof of Lemma 4.3.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00832v3-figure-04a1d061-4c95-457e-b798-8628d801c1c1.jpg)

Lemma 4.3. ω ref defined above is a valid reference flow. That is, the total mass sent out of any white vertex w , ∑ b ∼ w ω ref ( wb ) , is equal to 1; and the total mass received to any black vertex b , ∑ w ∼ b ω ref ( wb ) , is also 1.

Proof. Let w be a white vertex. Notice that in ∑ b ∼ w ω ref ( wb ), the oriented diagonals form a clockwise loop whose total winding is -2 π . Adding + π for each of the surrounding black vertices, and dividing by 2 π , we see that the total flow out of w is indeed 1 as desired.

For a black vertex b , the argument is better explained by considering a picture (see Figure 8). Fatten the 'star' formed by the half-diagonals incident to b into a star shaped domain. Then notice that the total flow out of b is simply the limit of the total winding, divided by 2 π , of the boundary of this domain (again in the clockwise orientation this time), as the domain thins into the 'star'. Indeed, the -π term in the definition of ω ref in (4.2) counts the half-turn as we move from the left side to the right side of a half-diagonal. Since the total winding of such a curve is -2 π , it follows that the total flow out of b is -1, as desired.

We are now ready to relate the three notions of height function defined by the pair of dual spanning trees T , T † . To do so, note that we can extend the definition of the winding field of both T and T † from Lemma 4.1 to the augmented trees.

Proposition 4.4. In the above setup, let h T and h T † be the winding field of T and T † respectively. Let h dim be the height function corresponding to the dimer configuration obtained from ( T , T † ) with reference flow ω ref . For any two faces f, f ′ ,

$$h _ { \mathcal { T } } ( m ( f ^ { \prime } ) ) - h _ { \mathcal { T } } ( m ( f ) ) = h _ { \mathcal { T } ^ { 1 } } ( m ( f ^ { \prime } ) ) - h _ { \mathcal { T } ^ { 1 } } ( m ( f ) ) = 2 \pi ( h _ { \dim } ( f ^ { \prime } ) - h _ { \dim } ( f ) ).$$

Remark 4.5. The winding fields in T and T † now play a completely symmetric role.

Proof. For a face f , we define γ f (respectively γ † f ) to be the branch of the augmented tree T (resp. T † ) starting from m f ( ). Fix faces f and f ′ and assume without loss of generality that γ f ′ is to the right of γ f when they merge. Note that this means that γ † f ′ is to the left of γ † f . Let γ ff ′ and γ † f f ′ denote the paths in T and T † connecting respectively m f ( ) to m f ( ′ ) and m f ( ′ ) to m f ( ). With the above convention on the relative position of γ f and γ f ′ , the concatenation of γ ff ′ with γ † f f ′ has a to a simple clockwise loop and there is no jump of the winding at m f ( ) and m f ( ′ ) because by assumption the midpoint m f ( ) is a smooth point of d f ( ) so

$$W _ { \text{int} } ( \gamma _ { f f ^ { \prime } } ) + W _ { \text{int} } ( \gamma _ { f ^ { \prime } f } ^ { \dagger } ) = - 2 \pi$$

The first equality easily follows by applying the definition of winding field from Lemma 4.1.

For the last equality let f , f l r be adjacent faces. Let the common (oriented) edge be ( bw ) with b being the black and w the white vertex and let f r lie to its right. We assume without loss of generality that b is a primal vertex. From the definition of ω ref and recalling the sign convention of the flow defining the height function (Section 2.5)

$$2 \pi ( h _ { \dim } ( f _ { r } ) - h _ { \dim } ( f _ { l } ) ) & = \omega _ { r e f } ( w b ) - 2 \pi 1 _ { ( b w ) \, \text{occupied by dimer} } \\ & = W _ { \text{int} } ( ( f _ { r }, b ) \cup ( b, f _ { l } ) ) + \pi - 2 \pi 1 _ { ( b w ) \, \text{occupied by dimer} }.$$

Note that that γ f r and γ f l merge at b . Also note from Temperley bijection that if ( bw ) is occupied by a dimer, then γ b starts by using the (half) edge bw which implies that γ f l lies to the left of γ f r . Otherwise, γ f l lies to the left of γ f r . We conclude using the definition of winding field from Lemma 4.1 and the above equation.

From tree to forest. One can think of Proposition 4.4 as giving a way to compute the (gradient of the) winding field by following the path between two vertices in the the tree. Before dealing with forests, we generalize the above by showing how to compute a height difference along an arbitrary path in Γ, which might use edges not in the tree.

We suggest the reader to look at Figure 9 while reading the following definitions. Fix f , f ′ two faces of G and let P be a self avoiding path starting at m f ( ), going along d f ( ) up to the vertex of Γ adjacent to f , then moving along edges of Γ up to f ′ and finally following d f ( ′ ) up to m f ( ′ ). Partition this path into connected components belonging to T (called segments from now on) separated by edges not belonging to T , allowing for the segments to be reduced to a single vertex (so any two segments are joined by a single edge). We call these segments and edges respectively ( P , . . . , P 1 k ) and ( e , . . . , e 1 k -1 ) and we let ( x , y i i ) be the starting and ending points of P i and z i be the white intersection vertex in e i .

̸

If x i = y i , let ε i = +1 (resp. -1) if γ y i lies to the right (resp. left) of γ x i for 1 ≤ i ≤ k . If x i = y i (i.e. when P i is a single vertex), then let ε i = +1 (resp. -1) if γ x i starts to the left (resp. right) of P . Note that by definition the edges e i are not in T so z i are all in T † . We let δ i = +1 (resp. -1) if γ † z i starts to the left (resp. right) of P .

Lemma 4.6. Under the above setup,

$$h _ { \mathcal { T } } ( m ( f ) ) - h _ { \mathcal { T } } ( m ( f ^ { \prime } ) ) = W _ { \text{int} } ( P ) + \pi ( \sum _ { i = 1 } ^ { k } \varepsilon _ { i } + \sum _ { i = 1 } ^ { k - 1 } \delta _ { i } ) = 2 \pi ( h _ { \dim } ( f ) - h _ { \dim } ( f ^ { \prime } ) ).$$

Figure 9: An example for the computation of the height along an arbitrary trajectory. The path P in Γ is the union of solid blue and dashed blue edges and here it is composed of k = 3 partitions. The solid blue denotes the segments in the tree T with arrows giving the orientation of each edge. The union of solid red and solid blue edges is the modified path. The orange dotted paths denote the direction in which the primal tree goes off to ∞ and the green dotted path denotes the same for the dual tree. This determines the ε s ′ i and δ ′ i s. Here ε 1 = -1 , ε 2 = 1 , ε 3 = -1 and δ 1 = 1 , δ 2 = -1 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00832v3-figure-7f7b3926-ac46-49b1-93d2-1cba393e39ad.jpg)

Proof. Modify P as follows. Observe that the oriented edges e i = ( y , x i i +1 ) has two faces of G to their left. Call the one incident to y i , f -i and the one incident to x i +1 , f + i and m ,m -i + i the respective midpoints of the diagonals of these faces. By convention, also call x 1 = m + 0 and y k = m -k . We also join y i to x i +1 using the diagonal segments of f -i and f + i . Finally we delete the edges e i . This completes the modification of the path P and note that this modification does not change the intrinsic winding, so we are still going to denote the modified path P (see Figure 9).

Note that after this modification, connected components connecting diagonal midpoints staying in T alternate with components connecting diagonal midpoints staying in T † so we can apply Lemma 4.1 and Proposition 4.4 for all such connected components. The terms ε i and δ i match the ± π terms in Lemma 4.1 so we obtain the result.

We now define the winding field given by a dual pair of spanning forests. We are in the same setup as above for the graphs Γ Γ , except , † T , T † are not necessarily spanning trees but spanning forests, and they are dual to each other.

Definition 4.7. Suppose we are in the above setup and ( T , T † ) is a dual pair of spanning forests where each component is oriented towards one of its ends. The winding field of ( T , T † ) is a function defined on all diagonal midpoints which satisfies the first equality from Lemma 4.6. It is well defined and unique up to a global shift in R .

Proof that Definition 4.7 is consistent. Once we show consistency, uniqueness up to a global shift by R is clear. To that end, note that one can apply Temperley's bijection to ( T , T † ) to obtain a dimer configuration whose height function is always well defined up to a global shift in R and which can be computed along any path. The second equality of Lemma 4.6 shows that the definition is consistent.

Remark 4.8. Note that there is no assumption about the number of ends of ( T , T † ) in the above definition: it makes sense as long as an orientation is fixed towards a unique end in each component. Further note that in the above proof, we only need the pair ( T , T † ) to be

Figure 10: Splitting of the boundary vertex of Γ (red graph) into a cycle. The black graph is Γ † .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00832v3-figure-7208d0ae-70ff-487b-9b52-13e4a340c9e3.jpg)

associated to a dimer configuration by Temperley's bijection so Definition 4.7 also extends to the case where Γ and ( T , T † ) are obtained by lifting a graph and a Temperleyan pair ( t, t † ) (satisfying the assumptions in Section 2.4) from the manifold M ′ to its cover. Note however that we need to be a bit careful with the boundary vertices of Γ as they are somewhat special. In order not to deal with boundary components separately, we replace the boundary vertex by a cycle whose length is the same as the degree of the vertex, and with every vertex in the cycle having one neighbour in Γ (see Figure 10). It is not too hard to see that the neighbours can be chosen so that the embedding of Γ still satisfies the first three items of the definition of faithful embedding (see Section 2.4). Let us call these the boundary cycles of Γ (the boundary cycles defined in Section 2.4 belonged to Γ , so this should not cause confusion). † We now replace the boundary vertex in t by this cycle in a natural way. In this setup, every component of t, t † has a unique non contractible cycle, with the boundary component of t corresponding to the boundary cycles we just defined. It is clear from the unique path lifting property that an oriented loop in t, t † corresponds to a collection of bi-infinite simple paths in the lift, each one obtained by going along the loop infinitely many times in clockwise and anti-clockwise direction. Hence, after the above surgery of the boundary vertex, t , t † lift to a forest T , T † with every component having exactly two ends. Fixing an arbitrary orientation of the boundary cycles, we can also orient the lifts towards one of its ends. We can then readily apply Definition 4.7.

## 4.2 Winding of CRSF and height function

In this section, we give further details on the connections between winding and height functions which are specific to our setting, i.e graphs and spanning forests obtained as lifts. More precisely our main goal will be to provide for every pair of faces f, f ′ a canonical path γ ff ′ that makes the formula from Lemma 4.6 as simple as possible, see Theorem 4.10 for the final statement.

At this point, we need to spare a few words related to the removal of white vertices from the graph G to obtain a graph with a dimer cover. This operation can be interpreted as inserting certain discrete version of magnetic operators on the free field (e.g. in the sense of [19]). If we want to interpret the height function as winding, the height function would be additively multivalued where it picks up an additional ± 2 π winding when it goes around a removed white vertex. Recall that we introduced a puncture corresponding to each face obtained by removal of a white vertex (cf. Figure 2) and that the new manifold is called M ′ . According to the above discussion, it is natural to lift G ′ to the universal cover ˜ M ′ of M ′ and not the cover of M (for the rest of this section, when we speak of lift or cover, we will always refer to ˜ M ′ ) and call it ˜ . G ′ Note that punctures are mapped to the boundary of the disc and every non-outer face of G ′ is mapped to a quadrangle.

Fix a Temperleyan CRSF t of Γ ′ and its oriented dual t † , both augmented by joining the midpoints of diagonals as explained in Section 4.1. Each component of t or t † contains either a single cycle (which is oriented) or contains a boundary component. Using the surgery in Remark 4.8, each boundary component of t also has a unique non-contractible cycle. We call

↦

↦

Figure 11: The conformal maps between an annulus and its universal cover. The dotted lines separate different copies of the fundamental domain. The blue and red curves show two loops in the surface and the associated spines in the cover.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00832v3-figure-7f00f13b-9645-45f1-97fd-de6b5aabaf0c.jpg)

each component of the lift of a cycle a spine . Note that by definition, each spine is a bi-infinite path oriented towards one of its ends.

We now state a useful lemma about the geometry of spines.

Lemma 4.9. If M is hyperbolic, then any spine S is either a simple path in D connecting two points in ∂ D or a simple loop containing a unique point in ∂ D .

If M is the torus, any spine S is a bi-infinite periodic path and all spines and dual spines have the same asymptotic direction.

The proof of the first part of Lemma 4.9 uses tools from the theory of Riemann surfaces, so we postpone it to appendix A in order to not disrupt the rest of the argument too much.

By Lemma 4.9, we see that ( T , T † ) form a dual pair of spanning forest of Γ ˜ ′ and its dual so Section 4.1 applies to them. We now turn to the definition of the terms involved in our main statement, Theorem 4.10.

Given two spines S and S ′ , we have a well defined region Ω S,S ′ between them which is bounded by S , S ′ and (in the hyperbolic case) two portions of ∂ D (if S and S ′ are loops, then these portions of ∂ D are understood as only prime ends associated with the same point). We say that a spine or dual spine S ′′ separates S from S ′ if it connects these two portions of ∂ D . Since the graph ˜ G ′ only has accumulation points on ∂ D , it is easy to see that two spines can only be separated by finitely many others. Furthermore, two adjacent dual spines must be separated by a primal spine.

Now pick two faces f, f ′ of ˜ G ′ and let S, S ′ be the spines of the components of T containing m f ,m f ( ) ( ′ ). Let Ω S,S ′ be the component between them as above. Draw a simple curve connecting S and S ′ in Ω S,S ′ , and oriented from S to S ′ . For each primal spine σ ⊂ Ω S,S ′ (hence not including S and S ′ ), let ε σ be the algebraic number of times σ is crossed by this curve (with the convention of counting +1 if σ is crossed from its left to its right by the curve and -1 otherwise). If σ is not crossed, define ε σ to be 0. Define δ σ similarly for primal spines. Notice that

$$\sum _ { \sigma \subset \Omega _ { S, S ^ { \prime } } } \varepsilon _ { \sigma } + \delta _ { \sigma }$$

Figure 12: An illustration of spines in the hyperbolic case. The spines S and S ′ are separated by a dual spine drawn in dashed blue. This spine also separates S ′′ and S ′ (left). The second and third figures from the left illustrate the choice of ε S in Theorem 4.10: in the second, ε S = 1 = -ε S ′ . In the third figure (note that the orientation of S ′ is reversed), ε S = ε S ′ = 1 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00832v3-figure-17e0de7f-1076-460d-aa6f-22f8c698ed9b.jpg)

is a topological term which does not depend on the choice of the curve (where the sum is over all primal and dual spines contained in Ω S,S ′ ). Also, its only dependence on f and f ′ is through the spines S and S ′ .

For any face f , recall that γ f is the infinite oriented path in T started from m f ( ) and going off to infinity along the unique outgoing oriented edges. Let ζ be the limit point of γ f (which exists due to Lemma 4.9). Let ζ ′ be the limiting endpoint of S ′ which lies in the same connected component of ∂M ˜ ∩ Ω S,S ′ as ζ . Note that we are not defining ζ ′ to be the limit point of γ f ′ as it could be the case that ζ and ζ ′ might lie in different boundary components of ∂M ˜ ∩ Ω S,S ′ (depending on the orientation of S ′ ). In case S, S ′ are loops through the same boundary point, we want ζ and ζ ′ to be in the same prime end of Ω S,S ′ .

- · If the limit point of γ f ′ does not lies in the same component of ∂M ˜ ∩ Ω S,S ′ as ζ and f ′ is in Ω S,S ′ , then we set ε S ′ = -1,
- · in all other case, define ε S ′ = +1.

Finally define ε S = +1.

Let ( ζ, ζ ′ ) be the arc joining ζ and ζ ′ (which could be a single point if ζ = ζ ′ ) in ∂M ˜ in the hyperbolic case. Let h dim be the dimer configuration corresponding to the pair ( T , T † ) with reference flow given by Definition 4.2. We can now state the following final result relating the dimer height function to the winding of trees:

Theorem 4.10. In the hyperbolic case, let γ := γ ff ′ be the curve formed by concatenating γ f , ( ζ, ζ ′ ) and the path in T joining ζ ′ to m f ( ′ ) . Orient γ from m f ( ) to m f ( ′ ) . We have the following deterministic relation:

$$\ast$$

$$h _ { \mathcal { T } } ( m ( f ^ { \prime } ) ) - h _ { \mathcal { T } } ( m ( f ) ) & = 2 \pi ( h _ { \dim } ( f ^ { \prime } ) - h _ { \dim } ( f ) ) = W ( \gamma, m ( f ) ) + W ( \gamma, m ( f ^ { \prime } ) ) + \pi \sum _ { \sigma \in \hat { \Omega } _ { S, S ^ { \prime } } } ( \varepsilon _ { \sigma } + \delta _ { \sigma } ). \\ \intertext { r _ { \dots } } \ a k _ { \cdots } \dots \varepsilon \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 ( 4. 4 )$$

Here the sum ∑ σ ∈ ¯ Ω S,S ′ ε σ + δ σ is as in (4.3) but also includes S and S ′ In the case of the torus, we have

.

$$h _ { \mathcal { T } } ( m ( f ^ { \prime } ) ) - h _ { \mathcal { T } } ( m ( f ) ) = - 2 \pi ( h _ { \dim } ( f ^ { \prime } ) - h _ { \dim } ( f ) ) = \\ W ( \gamma _ { f }, m ( f ) ) - W ( \gamma _ { f ^ { \prime } }, m ( f ) ) + W ( \gamma _ { f }, m ( f ^ { \prime } ) ) - W ( \gamma _ { f ^ { \prime } }, m ( f ^ { \prime } ) ) + \pi \sum _ { \sigma \in \tilde { \Omega } _ { s, s ^ { \prime } } } ( \varepsilon _ { \sigma } + \delta _ { \sigma } ).$$

Proof. Let us consider the hyperbolic case first. Take two vertices x, x ′ in S, S ′ respectively. Observe that we can find a path joining x, x ′ which goes along the primal components, and

moves from one component to the other by 'jumping' over dual spines (i.e. going along edges whose dual belong to a dual spine). Also one can make sure that the path is minimal, in the sense that every component between those of x, x ′ (i.e. components which separate S, S ′ ) is visited at most once. Call this path ( x, x ′ ). Now letting x → ζ and x ′ → ζ ′ , we can also ensure that the portion of Ω S,S ′ bounded by ( x, x ′ ) , ( ζ, ζ ′ ) and the spines S, S ′ contains none of f or f ′ . Note that for this path, the choice of the ε S , ε S ′ exactly matches with the choice defined for Lemma 4.6. We refer the readers to the second figure of Figure 12 for the case when the limit point of γ f ′ does not lie in the same component of ∂M ˜ ∩ Ω S,S ′ and f ′ lies in Ω S,S ′ , in which case ε S ′ = -1.

Let ˜ be the path obtained by concatenating ( γ m f ,x , ( ) ) ( x, x ′ ) , ( x , m f ′ ( ′ )). From Lemma 4.6, it is clear that

$$\tilde { 2 } \pi ( h _ { \dim } ( f ^ { \prime } ) - h _ { \dim } ( f ) ) = W _ { \text{int} } ( \tilde { \gamma } ) + \pi \sum _ { S } ( \varepsilon _ { S } + \delta _ { S } ).$$

Indeed, notice that since the path is minimal, the ε S , δ S terms are defined so as to match with the definition of ε , δ i i in Lemma 4.6. Furthermore, W int (˜) = γ W int ( γ ) because of the above choice of x, x ′ . We finish the proof using Lemma 2.1 and the fact that the topological winding of ˜ and γ γ are equal..

For the torus case, the same proof applies by noticing that W x,x ,m f (( ′ ) ( ))+ W x,x ,m f (( ′ ) ( ′ )) converges to 0 as x and x ′ go to infinity along S, S ′ in the same asymptotic direction (the signs in eq. (4.5) come from the fact that one goes along γ f ′ in the opposite direction).

Finally, the fact that the term π ∑ S ( ε S + δ S ) converges follows simply from the fact that the number of non-contractible components is a.s. finite in the limiting CRSF (and hence so is the number of spines which separate S, S ′ into different components).

## 4.3 Mapping back to the manifold

In this subsection, we briefly explain how the notion of height function constructed above on the universal cover corresponds to a height one-form on the surface coming from the dimer cover. We also describe how the instanton component is a function of only the non-contractible loops of the Temperleyan CRSF.

Suppose γ is a path in M ′ connecting two faces f, f ′ and let h f, f ( ′ ) denote the height difference in M ′ calculated along γ . Now imagine two copies of the lift of γ to the universal cover ˜ M ′ and observe that using our formula (Equation (4.4) or (4.5)) for computing the height difference along the lifts in the universal cover, it is not clear how these quantities are related for different lifts and how these formulas map back to the surface M ′ . In our next lemma, we show that this difference among different copies can be explicitly calculated and is deterministic, so it does not contribute to the fluctuations of the height function, which is the quantity of interest in this article.

Recall from Sections 2.2 and 2.4 that the covering map p maps ˜ M /F ′ to ˜ M ′ where F is a subgroup of Z 2 in the torus case, and is a Fuchsian group (a discrete subgroup of the M¨bius o group) in the hyperbolic case. Thus one can map one lift of a path in G ′ to another copy simply by using an element from F , which is either a translation in the torus case, or a M¨bius map o in the hyperbolic case. In the following lemma we show how this mapping affects the winding. We encourage the reader to recall the formulas in Section 2.3 before reading the lemma.

Lemma 4.11. Suppose we are in the hyperbolic case, so ˜ M ′ = D . For any M¨ obius map ϕ from D to itself, mapping ˜ G ′ to itself, for any two diagonal midpoints m,m ′ ,

$$h _ { \mathcal { T } } ( \phi ( m ^ { \prime } ) ) - h _ { \mathcal { T } } ( \phi ( m ) ) = h _ { \mathcal { T } } ( m ^ { \prime } ) - h _ { \mathcal { T } } ( m ) + \arg _ { \phi ^ { \prime } ( \mathbb { D } ) } ( \phi ^ { \prime } ( m ( f ) ) - \arg _ { \phi ^ { \prime } ( \mathbb { D } ) } ( \phi ^ { \prime } ( m ( f ^ { \prime } ) ).$$

In particular, h T -arg( ϕ ′ ) can be mapped back as a one-form on the dual of G ′ , away from vertices corresponding to punctures. In the parabolic case the same identity holds with no argument terms.

Proof. We focus on the hyperbolic case since the case of the torus is trivial. If ϕ is a conformal map mapping ˜ G ′ to itself, then it must also map T to itself and in particular γ ϕ f ϕ f ( ) ( ′ ) = ϕ γ ( ff ′ ) in the notations of Theorem 4.10. Clearly the topological terms ϵ σ , δ σ from Theorem 4.10 are invariant under ϕ :

$$\sum _ { \sigma \in \Omega _ { S, S ^ { \prime } } } ( \epsilon _ { \sigma } + \delta _ { \sigma } ) = \sum _ { \sigma \in \Omega _ { \phi ( S ), \phi ( S ^ { \prime } ) } } ( \epsilon _ { \sigma } + \delta _ { \sigma } ).$$

For the term W γ,m f ( ( )) + W γ,m f ( ( ′ )), it transforms as an intrinsic winding under conformal maps so by Lemma 2.2

$$W ( \phi ( \gamma ), \phi ( m ( f ) ) ) + W ( \phi ( \gamma ), \phi ( m ( f ^ { \prime } ) ) ) & = W ( \gamma, m ( f ) ) + W ( \gamma, m ( f ^ { \prime } ) ) \\ & \quad + \arg _ { \phi ^ { \prime } ( \mathbb { D } ) } ( \phi ^ { \prime } ( m ( f ) ) - \arg _ { \phi ^ { \prime } ( \mathbb { D } ) } ( \phi ^ { \prime } ( m ( f ^ { \prime } ) ), \ \ ( 4. 6 )$$

which concludes the proof.

Remark 4.12. A consequence of Lemma 4.11 is the following. Fix any two faces f, g of G ′ and any path γ in the dual of G ′ starting at f and ending at g . Let h f, g ( ) be the sum of the height one form calculated along γ . Suppose ˜ and ˜ γ γ ′ be two copies of the lifts of γ with endpoints ˜ ˜ and f, g ˜ f , g ′ ˜ ′ respectively. Recalling the definition of height function h dim from Sections 4.1 and 4.2, observe that

$$h _ { \dim } ( \tilde { f } ) - h _ { \dim } ( \tilde { g } ) - \mathbb { E } ( h _ { \dim } ( \tilde { f } ) - h _ { \dim } ( \tilde { g } ) ) = h _ { \dim } ( \tilde { f } ^ { \prime } ) - h _ { \dim } ( \tilde { g } ^ { \prime } ) - \mathbb { E } ( h _ { \dim } ( \tilde { f } ^ { \prime } ) - h _ { \dim } ( \tilde { g } ^ { \prime } ) ).$$

since the arg terms cancel. In other words, the dimer height fluctuation does not depend on the choice of the copy of γ used in the cover, and hence can be mapped back into M ′ without any ambiguity.

Let us now consider the instanton component. Intuitively, it should be completely determined by the homology class of the non-contractible loops of the Temperleyan forests (and in particular it should be determined by the spines). To that end, take a rooted loop λ , and consider the height difference along a particular choice of the lift of the loop, starting and ending at copies of the root, call this h λ ( ). Lemma 4.11 tells us that we can compute h λ ( ) unambiguously on any choice of the lift but we should in fact be able to say more. In particular, computing h λ ( ) using eqs. (4.4) and (4.5) a priori require more precise geometric information involving windings. The following lemma shows that in fact the computation of h λ ( ) does not involve any winding term and that the instanton component is at least a 'reasonable' function of the non-contractible loops.

Take an ordered finite set of continuous simple loops which forms the basis of the first homology group of M ′ , all endowed with a fixed orientation. Let H be the finite set of numbers which denote the height change along these loops. It is well-known (see Theorem 2.9) that h is completely determined by H (at the discrete level).

Lemma 4.13. Let ( C C , † ) be the set of oriented non-contractible loops of the primal and dual Temperlayen CRSF. Let h be the instanton component of the height 1-form of the dimer configuration corresponding to this CRSF pair given by the extended Temperley bijection. Then h is a function of ( C C , † ) only.

Proof. We focus first on the hyperbolic case. Let λ be a loop from the basis of the homology group considered above, without loss of generality we can assume that λ starts in a diagonal midpoint. Choose an arbitrary lift of λ and call its starting and ending points m and m ′ and the corresponding faces f and f ′ .

Let ϕ be the Mobius map sending m to m ′ and mapping ˜ G ′ to itself. By Theorem 4.10 and Lemma 2.1, h T ( m ′ ) -h T ( m ) = W int ( γ ) + π ∑ σ ∈ ¯ Ω S,S ′ ( ϵ σ + δ σ ) with γ obtained in this case by concatenating γ f , the arc ( ζ, ζ ′ ) and the time reverse of γ f ′ (we assume for a brief moment that

γ is smooth close to ζ and ζ ′ and therefore piece-wise smooth everywhere). By additivity of the intrinsic winding we have W int ( γ ) = W int ( γ f ) + W int (( ζ, ζ ′ )) -W int ( γ f ′ ) but since γ f ′ = ϕ γ ( f ), by Lemma 2.2, W int ( γ f ) -W int ( γ f ′ ) = -arg( ϕ ′ ( ζ )) + arg( ϕ ′ ( m )). Clearly, the above argument extends by taking an arbitrary regularisation of γ if it is not smooth near ζ and ζ ′ (since the right hand side is independent of the regularisation procedure). Thus we have proved that

$$h _ { \mathcal { T } } ( m ^ { \prime } ) - h _ { \mathcal { T } } ( m ) = \arg ( \phi ^ { \prime } ( m ) ) - \arg ( \phi ^ { \prime } ( \zeta ) ) + W _ { \text{int} } ( ( \zeta, \zeta ^ { \prime } ) ) + \pi \sum _ { \sigma \in \tilde { \Omega } } ( \epsilon _ { \sigma } + \delta _ { \sigma } ). \quad ( 4. 7 )$$

The function ϕ only depends on λ and the choice of lift m for its starting point. Also clearly ζ, ζ ′ , ¯ Ω S,S ′ and the ϵ σ , δ σ only depend on the set of spines which concludes the proof of the hyperbolic case.

The proof for the parabolic case is identical and in fact even simpler since γ f and γ f ′ are translate of each other and there is no term for the winding of ( ζ, ζ ′ ) in Theorem 4.10.

We record now the following corollary regarding the weak limit of H as the mesh size goes to 0. For this corollary, we need to rely on results from the second article in the series [5] (and is not a deterministic result like the rest of this section) which will be reproduced later in the paper in Theorem 6.3. We write a superscript δ to account for the dependence in δ as in Section 2.4.

Corollary 4.14. H # δ described above converges in law as δ → 0 to some H which is measurable with respect to the limit of ( C # δ , ( C † ) # δ ) .

Proof. The term arg( ϕ ′ ( m )) -arg( ϕ ′ ( ζ ))+ W int (( ζ, ζ ′ )) is clearly a continuous function of ( C C , † ). Since ( C # δ , ( C † ) # δ ) converges to an almost surely disjoint set of loops using Theorem 6.3, ∑ σ ∈ ¯ Ω ( ϵ σ + δ σ ) in the right hand side of (4.7) also converges in law, thereby completing the proof.

## 5 Local coupling

In this section we prove a local discrete coupling result which extends ideas of [4] to the setup of Riemann surfaces. Roughly speaking, the goal of such a result is to show that the local geometry in a small neighbourhood of a Temperleyan CRSF is given by that of a uniform spanning tree in the surface or, alternatively (and more usefully), in some reference planar domain. Moreover, locally around a finite number of given points, the configurations can be coupled to independent such USTs.

Recall that the actual Temperleyan CRSF's cannot be completely sampled using the standard Wilson's algorithm. However, due to Theorem 3.7, given the skeleton s emanating from either side of the punctures, the rest of the Temperleyan CRSF can be sampled from Wilson's algorithm. The main result in this section is Proposition 5.10.

The argument follows the same line of arguments as in [4] so let us first recall the strategy there. This consists of two main steps. Consider k points v , . . . , v 1 k in (Γ ) ′ # δ . In the first step, we choose cutsets at a small but macroscopic distance around each of the k points, such that the cutsets separate the points from each other and from the rest of the graph. We reveal all the branches emanating from these cutsets. This leaves k unexplored subgraphs Γ # δ i , one around each point. In this step the key point is to make sure that the Γ # δ i are macroscopic (e.g., contain a ball of a radius roughly of the same order of magnitude as the distance to the cutsets). Clearly, the conditional law of the tree in each Γ # δ i is that of a wired UST. Moreover, these wired USTs are also independent conditionally given the cutset exploration. (Of course, unconditionally there is still some dependence). The second step is then to say that in each Γ # δ i , one can couple a wired UST with a full plane one, which shows, among other things, that the unconditional distribution is close to being independent.

To adapt this strategy to Riemann surfaces, if the Γ # δ i are sufficiently small so that it is simply connected, then the conditional law of the CRSF in each Γ # δ i will also be that of a wired UST, so that the second step can be used directly. For the first step however (cutset exploration), we will need to redo parts of the proof to take into account the possible loops in the CRSF.

We first state few useful lemmas from the second article ([5, Lemmas 2.5 and 2.6])

Lemma 5.1 (Beurling type estimate) . For all r, ε &gt; 0 there exists η &gt; 0 such that for any δ &lt; δ η ( ) and for any vertex v ∈ Γ # δ such that η/ 2 &lt; d M ( v, ∂M ) &lt; η , the probability that a simple random walk exits B M ( v, r ) := { z ∈ M d : M ( z, v ) &lt; r } before hitting ∂ Γ # δ is at most ε .

Lemma 5.2. Let K 0 ⊂ K ′ 0 ⊂ M be open connected sets and let K ⊂ K ,K ′ 0 ′ be compact sets which are the closures of K ,K 0 ′ 0 . Also assume K contains a loop which is non-contractible in M . Then there exists a δ K,K ′ &gt; 0 and α K,K ′ &gt; 0 depending only on K,K ′ such that for all δ &lt; δ K,K ′ and all v ∈ Γ # δ such that v ∈ K , simple random walk started from v has probability at least α K,K ′ of forming a non-contractible loop before exiting K ′ .

Next, we recall a lemma from [5] about the regularity of the skeleton branches away from the punctures. As mentioned above, the skeleton s cannot be directly sampled using Wilson's algorithm, however the following (technical) result from [5] guarantees that it will not be qualitatively different from a regular loop-erased random walk away from the punctures. Let ¯ λ M denote the law of s . Let η r := { η ( ) i r = ( η ( ) i r, 1 , η ( ) i r, 2 ) } 1 ≤ ≤ i k denote the portion of the branches of the skeleton until they exit a ball of radius r around the puncture, where r is some fixed small quantity. For each such path η ( ) i r,j (with 1 ≤ i ≤ k , j = 1 2) run independent simple random , walks from their tips, where each walk is conditioned to exit the ball of radius r ′ &gt; r before hitting η ( ) i r,j , for another small choice of r ′ . After hitting r ′ the walks are allowed to run unconditionally until they either hit η r or ( ∂ Γ ) ′ # δ or their loop erasure creates a non-contractible loop. Let ˜ µ r,r ′ denote the joint law of the loop erasures of the resulting pair of walks.

Proposition 5.3 (Proposition 6.5 of [5]) . For any event E which is measurable with respect to s \ η r ,

$$\underline { \lambda } _ { \mathcal { M } } ( E ) \leq C \left ( \sqrt { \sum _ { j = 1 } ^ { 2 k } \sup _ { \eta _ { r / 2 ^ { j } } } \tilde { \mu } _ { \mathcal { M } _ { r / 2 ^ { j }, r / 2 ^ { j - 1 } } } ( E ) } \right )$$

where the supremum is over all possible η r/ 2 j .

Thanks to that proposition, we will be able to use Wilson's algorithm to prove estimates about all the branches in Temperleyan CRSF, including the skeleton.

## 5.1 Cutset exploration

We now describe the construction more precisely. Fix k points v , . . . , v 1 k in (Γ ) ′ # δ and let N v i be small enough neighbourhoods around each v i such that { N v i } 1 ≤ ≤ i k do not intersect each other or the pre-specified boundary ∂ Γ # δ of Γ # δ (if the boundary is non-empty) and for each 1 ≤ i ≤ k , p -1 ( N v i ) is a disjoint union of sets in ˜ M such that p restricted to each component is a homeomorphism to its image. Also for each i , fix one pre-image ˜ v i of v i and let ˜ N ˜ v i denote the component of ˜ in v i p -1 ( N v i ). Let R v , N (˜ i ˜ ˜ v i ) be the inradius seen from ˜ v i in ˜ M with respect to the Euclidean metric, and call

$$r = \min _ { 1 \leq i \leq k } R ( \tilde { v } _ { i } ; \tilde { N } _ { \tilde { v } _ { i } } ).$$

(We point out that in our setup it is natural to define these quantities with respect to Euclidean geometry on ˜ M and project to M , because our assumption and in particular the uniform crossing

Figure 13: A schematic representation (in solid line) of the lift of the loop erasure of the random walk on the torus until a non-contractible loop is formed. Call this path ˜ γ and let γ = p γ (˜) . The dashed square denotes the fundamental domain and the dashed paths denote some other lifts of γ . In this case, ˜ γ stops inside the fundamental domain.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00832v3-figure-7269950c-05bb-4e9d-80c9-9093f203008d.jpg)

condition are stated with respect to this geometry; whereas the intrinsic geometry of M does not really appear in our setup).

Wedenote by B euc (˜ v, r ) the Euclidean ball of radius r around ˜ and for v r ′ &gt; r , let A euc (˜ v, r, r ′ ) = B euc (˜ v, r ′ ) \ B euc (˜ v, r ) be the Euclidean annulus of inradius r and outer radius r ′ . Let H i be a set of vertices in p A ( euc (˜ v , r/ i 2 , r )) ∩ (Γ ) ′ # δ which disconnects v i from ∂N v i . The cutset exploration is done in two steps

- a. First sample the skeleton according to the correct law (i.e. according to CRSF branches with the global conditioning that the branches divide the manifold into annuli).
- b. Then, reveal all the branches from H , i 1 ≤ i ≤ k using (unconditioned) Wilson's algorithm, resulting in a subgraph T # δ H i , 1 ≤ i ≤ k .

We say a vertex v j has cutset isolation radius 6 -k r at scale r if

$$p ( B _ { \text{euc} } ( \tilde { v } _ { j }, 6 ^ { - k } r ) ) \, \text{does not contain any vertex from $\mathcal{T}_{H_{i}}^{\#\delta}$}.$$

Let us define J v j as the the minimum value such that v j has isolation radius 6 -J v j r and let J = max j J v j .

We want to show that J has exponential tail (which results in polynomial tail of the isolation radius). To this end, we rely on the following bound on the distance between a loop-erased walk and a point, which is a version of Proposition 4.11 in [4].

̸

Lemma 5.4. Let γ be a loop-erased random walk starting from a vertex v until it hits ( ∂ Γ ) ′ # δ or a non-contractible loop is created. Let u = v be another vertex and let ˜ u be one of the images under p -1 of u . Let r be small enough such that U := ¯ B euc (˜ u, r ) , is contained in the pre-image of N u which contains ˜ u . Then there exist constants α, c &gt; 0 (depending only on the initial assumptions of the graph) such that for all δ &lt; δ U and all n ∈ (0 log , 2 ( Crδ 0 δ ) -1) where δ U , δ 0 are as in the crossing condition (assumption (iv)),

̸

$$\mathbb { P } ( \gamma \cap p ( B _ { e u c } ( \tilde { u }, 2 ^ { - n } r ) ) \neq \emptyset ) < c \alpha ^ { n }.$$

Proof. This is almost identical to Lemma 4.11 in [4] except we now need to take into account the topology as well. To emphasise the differences with Lemma 4.11 in [4], we recall the strategy there. The idea is that if the loop erased walk on the manifold comes inside p B ( euc (˜ 2 u, -n r )) then some lift of the random walk must necessarily come within B euc (˜ 2 u, -n r ) - call this region (exponential) scale n . Furthermore, after the last such visit, this (lift of the) random walk must cross n many annuli without making a loop around ˜ (which we called a 'full turn'). This has u a probability bounded by e -cn .

In the current situation however, the random walk might form a non-contractible loop before exiting U , and therefore its relevant lift described above does not necessarily have to cross n many annuli (see e.g. Figure 13) after coming within Euclidean distance 2 -n r of ˜ before we u stop it. Thus we cannot simply apply the argument of the previous paper. We overcome this using the following idea which we first explain informally. The random walk on the manifold may visit p B ( euc (˜ u, r )) multiple times. Consider one such time when it enters p B ( euc (˜ u, r )) (not necessarily the first time) and suppose its relevant lift ˜ which also enters X B euc (˜ u, r ) at that time. Let X ′ be the portion of the loop-erasure of X on M at that time such that, if the random walk hits X ′ later on, then a non-contractible loop in M is formed. Assume that X ′ comes inside p B ( euc (˜ 2 u, -m r )). In order for the whole loop-erasure γ to come inside p B ( euc (˜ 2 u, -n r )) the following events must take place on one such visit. First, the lift ˜ X has to cross the first m scales without performing a full turn (otherwise it would close a non-contractible loop on the manifold and stop); this has a probability bounded by e -cm . Then ˜ X needs to come within Euclidean distance 2 -n r of ˜, but we bound crudely the probability of this event by 1. u Finally ˜ X needs to come back out after the last visit to B euc (˜ 2 u, -n r ) in such a way that no full turn occurs between distances 2 -n r and 2 -m r (for the same reason as in the simply connected case). The probability of this event can be bounded by e -c n ( -m ) . The intersection of these three events gives the right overall upper bound: e -cm e -c n ( -m ) = e -cn .

Let us fill in the details. Let C i be the Euclidean circle of radius r i := 2 -i r around ˜ for u i ≥ 0. We start a random walk X from v (to emphasise again: X is on the manifold M ). Let { τ k } be a sequence of stopping times defined inductively as follows. Set τ 0 = 0. Having defined τ k to be a time when the random walk X crosses or hits some circle p C ( i ( k ) ), define τ k +1 to be the smallest time after τ k when the random walk crosses or hits either p C ( i ( k ) -1 ) or p C ( i ( k )+1 ); this defines τ k by induction for every k ≥ 1. If i ( k ) = 0, define τ k +1 to be the smallest time after τ k when the random walk crosses or hits p C ( i ( k )+1 ). Let k exit be the smallest integer such that in the interval [ τ k exit , τ k exit +1 ] either a non-contractible loop is created by the loop-erasure of the random walk or the random walk hits the boundary ∂ Γ # δ . Let κ , κ 0 1 , . . . , κ N be the set of indices k &lt; k exit +1 when X τ k has crossed or hit p C ( 0 ). Also note that the portions X τ [ κ 0 +1 , τ κ 1 ], X τ [ κ 1 +1 , τ κ 2 ] , . . . are contained inside p C ( 0 ); while the portions X τ [ κ 0 , τ κ 0 +1 ] , X τ [ κ 1 , τ κ 1 +1 ] , . . . are contained outside p C ( 1 ). Our first claim is that N has exponential tail, i.e., ∃ α 1 ∈ (0 1), , such that for all n ≥ 1,

$$\mathbb { P } ( N > n ) < \alpha _ { 1 } ^ { n }.$$

Indeed, using Lemma 5.2, every time the walk hits p C ( 0 ), there is a positive probability independent of the past that the walk creates a non-contractible loop before returning to p C ( 1 ). Therefore, N has geometric tail. This proves (5.3).

̸

Let S be the set of points { X τ ( k ) } k ≥ 1 . Note that given S , the pieces of random walk { X τ , τ [ k k +1 ] } are independent of each other. We call such pieces elementary pieces of the walk. If i ( k ) = 0, X τ , τ [ k k +1 ] is (given S ) a random walk starting from X τ k conditioned to exit the annulus p A ( euc (˜ u, r i ( k ) -1 , r i ( k )+1 )) at X τ k +1 . If i ( k ) = 0, X τ , τ [ k k +1 ] is a random walk which is conditioned to next hit p B ( euc (˜ u, r/ 2)) at X τ k +1 . (Note that this property is lost if we solely condition on { X τ ( k ) } k ≤ κ N since the conditioning on N is complicated.) By Lemma 4.7 in [4], conditionally on S , each random walk piece X τ , τ [ k k +1 ] has a uniformly positive probability to do a full turn in the annulus p A ( euc (˜ u, r i ( k ) -1 , r i ( k ) )) for δ ≤ δ U and given range of n (where δ U comes from the uniform crossing assumption). Here we define a full turn to be the event that the random walk crosses every curve joining the inner and outer boundary in the specified annulus. Indeed, although Lemma 4.7 in [4] gives the uniform positive probability estimate for crossing in the Euclidean plane, the estimate is valid here by considering the relevant lift of X τ , τ [ k k +1 ] which is inside ¯ B euc (˜ u, r ) = U and applying the uniform crossing estimate in the universal cover. Let us also point out that δ is chosen small enough so that the uniform crossing estimate is valid inside U .

We now define certain low probability events E 1 , . . . , E n such that one of them must take

place if the loop-erasure of X is to enter B n := p B ( euc (˜ 2 u, -n r )). Let us condition on S . We say the event E j occurs if:

- (i) The portion X τ [ κ j -1 +1 , τ κ j ] intersects B n . We then let λ j be the smallest k in [ τ κ j -1 +1 , τ κ j ] such that X τ k crosses (or hits) p C ( n ).
- (ii) Let X ′ j be the portion of the loop-erasure of X [0 , τ κ j -1 +1 ] such that if the walk X τ [ κ j -1 +1 , τ κ j ] intersects X ′ j , a non-contractible loop would be created. Let m be the maximum index such that X ′ intersects the circle p C ( ).
- m Then for the event E j to hold we also require, in addition to the previous point, that: for any κ j -1 +1 ≤ k ≤ λ j if i ( k ) &lt; m , then the walk X τ , τ [ k k +1 ] does not perform a full turn in p A ( euc (˜ u, r i ( k ) -1 , r i ( k ) )). If no such m exists (e.g., if X ′ j is empty), we do not require anything further.
- (iii) Let ℓ j be the last k before κ j such that X τ k crosses (or hits) p C ( n ). Let ℓ ′ j be the first time after ℓ j that the walk intersects p C ( m +1 ) . Then in addition to the previous two points, for E j to hold we require that the walk X τ , τ [ k k +1 ] does not perform a full turn for any ℓ j ≤ k ≤ ℓ ′ j .

From the discussion above, it is clear that we have the following lemma.

## Lemma 5.5. We have

̸

$$\{ \gamma \cap B _ { n } \neq \emptyset \} \subseteq \{ N \leq n, \cup _ { j = 1 } ^ { n } { \mathbf E } _ { j } \} \cup \{ N > n \}.$$

Thus all we need to show is that the event on the right hand side of Lemma 5.5 has exponential tail bound. Now we claim that

$$\mathbb { P } ( E _ { j } | \mathcal { S } ) \leq e ^ { - c ^ { \prime } n } \text{ and so } \mathbb { P } ( E _ { j } ) \leq e ^ { - c ^ { \prime } n }.$$

This is justified using the uniform positive probability of the walk X performing a full turn, even conditionally given S . Indeed, conditioned on S , we have that the event (ii) in the definition of E j has probability at most e -c m ′ , and conditioned on event in (ii), the event in (iii) has probability at most e -c ′′ ( n -m ) since conditionally given S , the random walk portions X τ , τ [ k k +1 ] are independent. Thus the overall probability given S is at most e -cn with c = c ′ ∧ c ′′ . All in all, using eq. (5.3), Lemma 5.5 and a union bound, we obtain

̸

$$\mathbb { P } ( \{ \gamma \cap B _ { n } \neq \emptyset \} ) \leq n e ^ { - c ^ { \prime } n } + e ^ { - c ^ { \prime } n } \leq c e ^ { - c ^ { \prime } n }.$$

thereby concluding the proof.

Lemma 5.6. Let s denote the skeleton. Let u be a vertex which is not one of the punctures and let ˜ u be one of the images under p -1 of u . Let s be small enough such that U := ¯ B euc (˜ u, s ) , is contained in the pre-image of N u which contains ˜ u . Then there exist constants α, c &gt; 0 (depending only on the initial assumptions of the graph) such that for all δ &lt; δ U and all n ∈ (0 log , 2 ( Csδ 0 δ ) -1) where δ U , δ 0 are as in the crossing condition (assumption (iv)),

̸

$$\mathbb { P } ( \mathfrak { s } \cap p ( B _ { e u c } ( \tilde { u }, 2 ^ { - n } s ) ) \neq \emptyset ) < c \alpha ^ { n }.$$

Proof. This is a corollary of Proposition 5.3 and Lemma 5.4. Consider the case k = 1 and recall the notations from from Proposition 5.3. It is easy to see using the Markov property of the independent random walks used to sample ˜ µ r/ ,r 2 that after hitting r , the result of Lemma 5.4 kicks in, yielding the required exponential bound for ˜ µ r/ ,r 2 ( E ). The same argument obviously applies to ˜ µ r/ ,r/ 4 2 ( E ), yielding the desired statement.

The case k &gt; 1 is analogous. We skip the details here.

We now state the result showing an exponential tail of J , which is a combination of Lemmas 5.4 and 5.6 and Schramm's lemma, and is identical to the proof of Lemma 4.20 in [4] (see also the proof of Theorem 4.21 in [4]); as it is identical we skip the proof here.

Lemma 5.7. There exist constants c, c ′ &gt; 0 such that the following holds. Let ˜ D be a compact set containing B v , r (˜ i ) for 1 ≤ i ≤ k where r is as in (5.2) . Then for all δ ∈ (0 , δ ˜ D ) and for all m ∈ (0 log , 6 ( δ 0 r/δ ) -1) ,

$$\mathbb { P } ( J > m ) \leq c e ^ { - c ^ { \prime } m }.$$

Finally, we state a lemma which says that with exponentially high probability, a branch of the CRSF, after entering an exponential scale t does not backtrack to a smaller scale. Such a lemma for SLE curves can be found in [36] (see also Lemma 3.4 in [4] in the simply connected case).

Lemma 5.8. Fix u, U, δ U as in Lemma 5.4. Suppose γ is the loop-erasure of a simple random walk in Γ # δ started from vertex u until it hits the boundary or creates a non-contractible loop. Suppose ˜ γ is the lift of γ started from ˜ u (parametrised from ˜ u to its endpoint). There exist constants c, c ′ &gt; 0 (depending only on the initial assumptions of the graph) such that for all δ &lt; δ U and all n ∈ (0 log , 2 ( Crδ 0 δ ) -1)

$$\mathbb { P } ( \tilde { \gamma } \ e n t e r s \ B _ { e u c } ( \tilde { u }, r 2 ^ { - n } ) \ a f t e r \ e x i t i n g \ B _ { e u c } ( \tilde { u }, r 2 ^ { - n / 2 } ) ) \leq c e ^ { - c ^ { \prime } n }.$$

Proof. Let E be the event that ˜ enters γ B euc (˜ u, r 2 -n ) after exiting B euc (˜ u, r 2 -n/ 2 ). Let r be as in Lemma 5.4 and assume n is even without loss of generality. The argument for this is very similar to Lemma 5.4 and in fact simpler, so we content ourselves with a sketch. Let r, C , B i i be as in the proof of Lemma 5.4. We look at the lift ˜ of the simple random walk X X started at ˜ and we stop when either u X hits the boundary or a non-contractible loop is created. Let τ k be the set of stopping times defined as in the proof of Lemma 5.4 but for the lift ˜ X instead of X and let S be the set { ˜ ( X τ k ) } k ≥ 1 . Observe that lift of the loop erasure of X is the loop erasure of ˜ X since erasing contractible loops commutes with lifting to the universal cover. If the loop erasure of ˜ X has to backtrack to scale n , the random walk has to necessarily enter scale n after leaving scale n/ 2. Let J be the largest j such that ˜ X τ j crosses or hits C n after leaving C n/ 2 and let I be the largest index i smaller than J when ˜ X τ i crosses or hits C n/ 2 . Conditioned on S , if E occurs, then ˜ X enters B n at least once after leaving B n/ 2 , and none of the elementary pieces of the walk between τ I and τ J can perform a full turn. But again, conditioned on S there is a uniformly positive probability to do a full turn for each elementary piece. Since there are at least n/ 4 such elementary pieces contained in [ τ I , τ J ], we conclude the proof of the lemma applying the upper bound on the full-turn estimate on each elementary piece.

## 5.2 Full coupling

The results of the previous section covered the first step in the proof of the coupling. As we mentioned above, the second step is identical to the simply connected case so we only recall the main statements.

We first recall the result which we will need from [4]. We use the notations and assumptions already in force. Let ˜ D ⊂ ˜ D ′ ⊂ ˜ M ′ be two simply connected compact domains and fix ˜ v ∈ ˜ D # δ . Let T ˜ D ′ , T ˜ D denote the wired UST respectively in ( ˜ D ′ ) # δ and ˜ D # δ . Let r ˜ v denote R v, N (˜ ˜ ˜ v ) as in eq. (5.2) (the largest Euclidean radius so that p is injective).

Lemma 5.9 (Theorem 4.21, [4]) . There exists c, c ′ &gt; 0 such that the following holds. Fix ˜ v, D, D ˜ ˜ ′ as above. There exists a coupling between T ˜ D and T ˜ D ′ and a random variable R &gt; ′ 0 such that

$$\mathcal { T } ^ { \tilde { D } ^ { \prime } } \cap B _ { e u c } ( \tilde { v }, R ^ { \prime } ) ^ { \# \delta } = \mathcal { T } ^ { \tilde { D } } \cap B _ { e u c } ( \tilde { v }, R ^ { \prime } ) ^ { \# \delta }.$$

Furthermore, for all δ &lt; δ ˜ D ′ and for all m ∈ (0 log , 6 ( δ 0 r ˜ v /δ ) -1) , if we write R ′ = 6 -K v r ˜ v,D ˜ where r ˜ v,D ˜ is the minimum of r ˜ v and the Euclidean distance between ˜ v and ∂D ˜ , then

$$\mathbb { P } ( K _ { v } \geq m ) \leq c e ^ { - c ^ { \prime } m }.$$

̸

We now put together the cutset exploration with the above one-point coupling from the simply connected case as follows. Recall the notations from Section 5.1. Let T # δ H i be the branches revealed in the cutset exploration around each v i . Let T # δ H = ∪ k i =1 T # δ H i . Let Γ # δ i be the connected component of (Γ ) ′ # δ \ T # δ H containing v i . Observe that given T # δ H , the law of T # δ ∩ Γ # δ i are independent wired UST in Γ # δ i (with the natural boundary), by the generalised Wilson algorithm. Applying the coupling of Lemma 5.9 in the lift of Γ # δ i to ˜ M ′ and some fixed compact ˜ D in ˜ M ′ containing some lifts ˜ v , . . . , v 1 ˜ k of points v , . . . , v 1 k , we obtain a coupling of the oriented CRSF in Γ # δ and independent wired USTs ( T ˜ D i ) 1 ≤ ≤ i k in ˜ . D Furthermore, the r in (5.2) and r ˜ v in Lemma 5.9 differs by a constant multiplicative factor which depends only on the choice of lifts of the vertices. Furthermore, we may (and will) choose the ( T ˜ D i ) 1 ≤ ≤ i k to be independent of T # δ H . In fact, note that for each 1 ≤ i ≤ k , T ˜ D i may be chosen independent of the restriction of T # δ to ∪ j = i p B v , r/ ( (˜ j 2)).

Proposition 5.10. The above coupling has the following properties: there exists random variables R , . . . , R 1 k such that

$$\mathcal { T } ^ { \# \delta } \cap p ( B ( \tilde { v } _ { i }, R _ { i } ) ) = p ( \mathcal { T } _ { i } ^ { \tilde { D } } \cap B ( \tilde { v } _ { i }, R _ { i } ) ).$$

Furthermore if we write R i = 6 -I v i r i where r i is the minimum of r as in (5.2) and the distance between ˜ v i and ∂D ˜ , then for all δ ≤ δ ˜ D , and for all 1 ≤ i ≤ k , for all n ∈ (0 , log 6 ( δ 0 r /δ i ) -1) ,

$$\mathbb { P } ( I _ { v _ { i } } \geq n ) \leq c e ^ { - c ^ { \prime } n }$$

for some constants c, c ′ &gt; 0 (depending only on the initial assumptions on the graph). In particular, P ( I v i ≥ n ) ≤ ce -c n ′ ∨ δ c ′ .

The set of non-contractible loops of T # δ is measurable with respect to T # δ H . In particular, ( T ˜ D i ) 1 ≤ ≤ i k are also independent of the set of oriented non-contractible loops in the oriented CRSF.

Observe that when I v i is very big or r i is very small, it is possible that the ball B v , R ( i i ) is reduced to a point so the statement (5.4) is trivial (that is, (5.4) holds for a single vertex). This happens with a probability which is at most δ c ′ for some c ′ .

Proof. This follows immediately from Lemma 5.7 and Lemma 5.9 once we observe that I v i is within O (1) of the sum of J v i in Lemma 5.7 and K v i in Lemma 5.9, both of which have exponential tails.

The proof of the final assertion is a topological fact. Indeed, conditioned on T # δ H , it is a deterministic fact that none of the oriented non-contractible loops pass through the portion of the CRSF in Γ # δ i (which is a wired UST). Indeed, otherwise we would have a path joining two points of the wired boundary of a wired UST which is impossible. Thus, the wired UST in each Γ # δ i is conditionally independent of the non-contractible loops and hence so are T ˜ D i s.

## 6 Convergence of height function and forms

In this section, we precisely state our main result (this is Theorem 6.1) and then prove it. Recalling the sketch from Section 1.3, we see that what remains to be done is to go from the

convergence of the Temperleyan CRSF to the convergence of its winding field using the coupling of Section 5. The global idea is the same as in [4] but some of the estimates on winding of LERW have to be redone. The primary difficulty is that spines are made of multiple copies of LERW, and hence Wilson's algorithm cannot be used to estimate the winding of all the copies. These are dealt with in Sections 6.2 and 6.3. The proof of the main result is finally concluded in Section 6.4.

## 6.1 Precise statement of the result

From now on, we work with the manifold M ′ . Recall that it is obtained by removing 2 g + b -2 points from the interior of M and that the white vertices removed from G # δ to obtain a Temperleyan graph in Theorem 3.7 will converge to these points. Denote by ˜ M ′ the universal cover of M ′ ; in the sequel 'lift' refers to lift on ˜ M ′ . Let h # δ be the height function defined as in Section 4 sampled using eq. (1.1). Recall that h # δ is a function from the the dual of ( ˜ G ′ ) # δ to R which is defined up to a global additive constant and ¯ h # δ := h # δ -E ( h # δ ) is independent of the choice of the reference flow (see Remark 4.12). Recall the probability measures P Wils and P Temp from Section 3.3.

First, we extend h # δ to a function h # δ ext : ˜ M ′ → R by defining h # δ ext ( x ) to be the value of h # δ on the face containing x . Recall the graphs (Γ # δ ) ′ , (Γ † , # δ ) ′ from Section 2.4. Let T # δ be the Temperleyan forest sampled using P Temp . Recall that we endowed M with a Riemannian metric d M as in Section 2, and we equip M ′ with the same metric via the inclusion map. This induces a volume form d µ in ˜ M ′

Theorem 6.1. Let ˜ : f ˜ M ′ → R be a smooth compactly supported function with ∫ ˜ M ′ ˜ d f µ = 0 . Suppose the assumptions in Section 2.4 hold. Then

$$\left ( \int \tilde { f } ( x ) \bar { h } _ { e x t } ^ { \# \delta } ( x ) \, \mathrm d \mu ( x ), \mathcal { T } ^ { \# \delta } \right )$$

converges jointly in law as δ → 0 . The first coordinate also converges in the sense of all moments. Furthermore, the limit of the first coordinate is measurable with respect to the limit T of T # δ , is universal (in the sense that it does not depend on the graph sequence ( G ′ ) # δ ), and is conformally invariant.

To clarify, convergence in the sense of all moments means that for all i , denoting X # δ := ∫ ˜ ( f x h ) ¯ # δ ext ( x ) d µ x ( ), E ( | X # δ i | ) converges as δ → 0. Notice also that since ∫ ˜ d f µ = 0, the fact that ¯ h # δ ext ( x ) is defined only up to a global additive constant is irrelevant.

Let us explain briefly what is meant by conformal invariance in our setting. Suppose | χ | = k and ( M,x ,x ,... , x 1 2 k ) and ( N,y , . . . , y 1 k ) are conformally equivalent in the sense that there is a map ϕ : M → N so that ϕ x ( i ) = y i for all 1 ≤ i ≤ k and ϕ is a conformal bijection between the two Riemann surfaces. Note that if p M and p N are their respective covering maps (satisfying the properties in Section 2.4), then there exists a M¨bius map o ψ such that p N = p M ◦ ψ . Let h M (resp. h N ) denote the limit field obtained in the setup of Theorem 6.1 with the punctures being x i and y i for 1 ≤ i ≤ k . Then h N = h M ◦ ψ in the sense that for any compactly supported test function f on the universal cover as in Theorem 6.1, ∫ h M f d µ has the same law as ∫ h N ( f ◦ ψ -1 ) ( | ψ -1 ) ′ | 2 d µ .

Remark 6.2. Before we start giving the details, we remark that the upcoming proof of the theorem goes through if we can establish the convergence of the Temperleyan CRSF (in the Schramm sense) to a conformally invariant limit in which all curves are a.s. simple, together with a few additional ingredients. In addition to the assumptions on the discretisation of the surface described in Section 2.4 and in particular the crossing assumption, these are the following: for each fixed vertex v ∈ (Γ ) ′ # δ , let γ v the branch of the Temperleyan forest starting from v . Then we require that γ v satisfies Lemma 5.4 (the branch does not come close to a given point)

as well as the moments on winding of Lemma 6.9 and Lemma 6.10. Note that the convergence to Brownian motion in item iii is not really needed here; instead, convergence of random walk to a process which is absolutely continuous with respect to Brownian motion would be sufficient. Stated this way, we note that Theorem 6.1 is novel even in the simply connected case. For instance, this result is applied in [3] and in the forthcoming [6].

For ease of reference, let us assemble in a single statement the convergence results from the second paper in the series that are needed in the following (see Theorem 3.1, Corollary 3.14 and Proposition 6.1 in [5]).

Theorem 6.3. Let M be a Riemann surface satisfying the assumptions of Section 2.1, which we equip with the distance d M described in Section 2.1. Let (Γ # δ ) δ&gt; 0 be a sequence of graphs with boundary ∂ Γ # δ faithfully embedded in M satisfying the assumptions of Section 2.4. Then the limit in law as δ → 0 of the Temperleyan CRSF sampled using P Wils or P Temp exists in the Schramm topology relative to d M . Both limits are independent of the sequence Γ # δ subject to the assumptions in Section 2.4. These limits are also conformally invariant.

Furthermore, let K # δ (resp. ( K † ) # δ ) denote the number of noncontractible loops in the CRSF (resp. its dual). Then for any q &gt; 1 there exists a constant C q &gt; 0 independent of δ such that

$$\mathbb { E } _ { \text{Wils} } ( q ^ { K ^ { \# \delta } } ) \leq C _ { q },$$

where E Wils denotes the expectation under P Wils . The same bound holds also with K # δ replaced by ( K † ) # δ .

Furthermore, let ( C # δ 1 , . . . , C # δ K # δ ) be the set of noncontractible cycles of the Temperleyan CRSF,

$$( C _ { 1 } ^ { \# \delta }, \dots, C _ { K ^ { \# \delta } } ^ { \# \delta } ) \xrightarrow { ( d ) } \left ( C _ { 1 }, C _ { 2 }, \dots, C _ { K } \right )$$

where C , . . . , C 1 K # δ are almost surely disjoint noncontractible simple loops in M .

## 6.2 Some a priori tail estimates on winding

The goal of this section is to obtain tail estimates on the winding of a branch of a Temperleyan forest. This will be achieved in Lemma 6.7 which is the main result of this section. Conceptually the arguments are similar to Section 4 of [4]. However, as in Section 5, there are additional difficulties linked with the fact that Wilson's algorithm can stop because of the formation of a non-contractible loop. Furthermore we ultimately want to consider the winding of entire spines, yet only one copy of the lift of a loop is directly connected to a loop erased random walk path. We first treat the part directly obtained by loop-erasure in this section and defer estimates on the copies for the next section.

Throughout this section, we deal with a CRSF sampled from P Wils (cf. Section 3.3). First we sample the skeleton (cf. Theorem 3.7) which we denote by s . Let ˜ be the lift of s s to ˜ M ′ . Recall that s decomposes M ′ into finite number of disjoint annuli and that, conditionally on s , in each of them we can sample the rest of the Temperleyan CRSF using Wilson's algorithm.

We now work conditionally on s . Let γ v be the branch of the CRSF started from v (which can thus be sampled using Wilson's algorithm). For any vertex ˜ v ∈ p -1 ( v ), let ˜ γ ˜ v be the lift of γ v starting from ˜ and up until the time when v γ v closes a non-contractible loop or hits the boundary ∂ Γ # δ ∪ s . Let τ nc be the stopping time when the simple random walk generating ˜ γ ˜ v creates a non-contractible loop or hits the boundary. Recalling the definition of spines from Section 4, note that ˜ γ ˜ v will include a part of a spine whenever the random walk does not hit ∂ Γ # δ ∪ s when it stops (see Figure 13).

Let us orient ˜ γ ˜ v starting from ˜ and going away from it in some continuous manner in [0 1]. v , For t &gt; 0, if ˜ γ ˜ v exits B v, e (˜ - -t 1 ), let t 1 be the first time ˜ γ ˜ v exits B v, e (˜ - -t 1 ) and let t 2 be the last time it exits B v, e (˜ -t ). In this case, if ˜ γ ˜ v ends in B v, e (˜ -t ) we set t 2 = 1.

We first state a simple deterministic lemma connecting the winding of curves avoiding each other.

Lemma 6.4. Consider the annulus A = { z ∈ C : r &lt; | z -x &lt; R | } where r &gt; 0 and R ∈ ( r, ∞ ] and let γ 0 be a simple curve in A connecting the outer and inner boundaries of A , assumed to be parametrised in [0 , 1] . Then for any simple curve γ in A \ γ 0 , we have

$$| W ( \gamma, x ) | \leq \sup _ { 0 \leq s _ { 1 } < s _ { 2 } \leq 1 } | W ( \gamma _ { 0 } [ s _ { 1 }, s _ { 2 } ], x ) | + 4 \pi.$$

For R = ∞ , the statement is interpreted as γ 0 connects the boundary of the B x,r ( ) to ∞ .

Note that in the above lemma we do not need any assumption on the regularity of the curves beyond continuity.

Proof. Observe that any two points in an annulus can be joined by a smooth curve lying completely in a half annulus containing those two points. Join each endpoint of γ to its nearest points (with its intrinsic metric inherited from the Euclidean plane) in γ 0 using smooth curves η , η 1 2 satisfying the above mentioned property. Assume that η 1 and η 2 intersect γ 0 at γ 0 [ s 1 ] and γ 0 [ s 2 ] respectively. Consider the loop L obtained by γ , η , η 1 2 and γ 0 [ s , s 1 2 ]. Note that L might be non-simple, but any loop which is created by the union of γ and the two geodesics lie in the annulus A slitted by γ 0 which is simply connected and does not contain x . Hence all these loops contribute winding 0 around x . Furthermore, it is easy to see that erasing all these loops (say in a chronological manner started from one of the tips of γ ) creates a simple loop. Since the winding of L after erasing these loops is 0 or ± 2 π , the total winding of L around x is also either 0 or ± 2 π . Since by definition η 1 and η 2 lie in the smaller half annulus containing the two points, its winding at most π . Using this observation, the proof of the lemma is complete.

Assume without loss of generality and to simplify notations that t is an integer. Let r ˜ v be such that B v, (˜ 10 r ˜ v ) does not intersect ˜ and s p is injective in B v, r (˜ ˜ v ) (this is a slight modification of the previous definition of r ˜ v ). Let e -t 0 = r ˜ v . We consider concentric circles C j of radius { e -t 0 -j } 0 ≤ ≤ j t +2 (with B j the ball inside it) around ˜. v Take a random walk X in (Γ ) ′ # δ starting from ˜ stopped when it hits ( v ∂ Γ ) ′ # δ ∪ s . In case ( ∂ Γ ) ′ # δ ∪ s = ∅ (i.e., in the case of the torus), the random walk continues forever. Let ˜ X be the lift of this walk starting from ˜. v Now let { τ k } k ≥ 0 be the set of stopping times as described in the proof of Lemma 5.4 for the random walk ˜ . X That is, if we hit or cross the circle C i ( k ) at time τ k , we wait until we hit or cross C i ( k ) ± 1 at time τ k +1 . If i ( k ) = 0 (or t +2), we wait until we hit or cross C 1 (or C t +1 ). Let B i denote the disc inside C i .

Let τ i 1 , τ i 2 , . . . be the successive times in the sequence ( τ k ) k ≥ 1 when the walk hits C 0 . Note that the interval [ τ i j , τ i j +1 ) is spent completely outside B 1 . We now claim that the random walk cannot wind too much outside B 1 .

Lemma 6.5. There exist c, c ′ &gt; 0 such that for all δ &lt; δ p B ( ¯ 0 ) , n ≥ 1 , j ≥ 1 and u ∈ B 1 such that P ( X τ i j +1 = u ) &gt; 0 ,

$$\mathbb { P } \Big ( \sup _ { \mathcal { Y } \subset \tilde { X } [ \tau _ { i _ { j } }, \tau _ { i _ { j } + 1 } ] } | W ( \mathcal { Y }, \tilde { v } ) | > n \Big | X _ { \tau _ { i _ { j } + 1 } } = u \Big ) \leq c e ^ { - c ^ { \prime } n }, \\ \mathbb { P } \Big ( \sup _ { \mathcal { Y } \subset \tilde { X } [ \tau _ { i _ { j } }, \tau _ { i _ { j } + 1 } ] } | W ( \mathcal { Y }, \tilde { v } ) | > n \Big | X _ { \tau _ { i _ { j } + 1 } } \in \partial \Gamma ^ { \# \delta } \cup \mathfrak { s } \Big ) \leq c e ^ { - c ^ { \prime } n }. \\.$$

$$\begin{smallmatrix}. \end{smallmatrix}$$

Here the supremum is over all continuous paths obtained by erasing portions of ˜ [ X τ i j , τ i j +1 ] .

Proof. The proof of Lemma 6.5 is very similar to Lemma 4.8 in [4]. But it needs an input from Riemannian geometry to control the winding of the spines near the boundary of the universal cover. We postpone this proof to appendix A.

Lemma 6.5 takes care of the winding of the excursions outside B 1 . For excursions inside, most of the technical work was done in [4]. Let ˜ Y j be the loop-erasure of ˜ [0 X , j ]. For any j , we parametrise ˜ Y j in some continuous way away from ˜. Fix v m ≥ 1 and let t 1 be the first time ˜ Y τ i m exits B v, e (˜ - -t 1 ) and t 2 be the last time ˜ Y τ i m exits B v, e (˜ -t ).

Lemma 6.6. There exist C, c &gt; 0 such that for all δ &lt; δ p B ( ¯ 0 ) , t ∈ ( t 0 , log( Cδ /δ 0 )) , and n ≥ 1 , m ≥ 1 ,

$$- \nu, \\ \mathbb { P } \left ( \sup _ { \mathcal { Y } \subseteq \tilde { Y } ^ { \tau _ { i _ { m } } } [ t _ { 1 }, t _ { 2 } ] \cap B ( \tilde { v }, e ^ { - t - 1 } ) ^ { c } } | W ( \mathcal { Y }, \tilde { v } ) | > n | \tau _ { i _ { m } } < \infty \right ) < C e ^ { - c n }.$$

We emphasise that in the above lemma, m is non-random as it is going to be important in what follows.

Proof. The proof of this lemma is identical to that of [4], Lemma 4.13 and the only difference is in the setup, which we point out. In [4], we were working with a simply connected domain and we waited until the random walk exited it. The argument proceeds by conditioning on the positions S := { ˜ X τ k } k ≥ 1 and arguing that, in each interval of walk between successive points in S , the winding of any continuous subpath has exponential tail. Then we proved that we only need to look at a random number of intervals which itself has exponential tail. In the present setup, we can condition on S so that τ i m &lt; ∞ is satisfied. The exponential tail of winding inside any inner annulus follows from Lemma 4.7 in [4] (this is exactly the same as the simply connected case) and for the outer annulus, we use Lemma 6.5. The number of relevant intervals to consider also has exponential tail following verbatim the proof of Lemma 4.15 in [4].

With these lemmas we can now state and prove the main result of this section, which controls the (topological) winding in an annulus.

Lemma 6.7. There exist constants C, C , c &gt; ′ 0 so that for all n ≥ 1 , for all ˜ v , for all δ &lt; δ ¯ (˜ B v,r ˜ v ) and for all 0 &lt; t &lt; log( C r δ /δ ′ ˜ v 0 ) ,

$$\mathbb { P } \Big ( \max _ { \mathcal { Y } \subseteq \tilde { \gamma _ { e } } [ t _ { 1 }, t _ { 2 } ] \cap B ( \tilde { v }, r _ { \tilde { v } } e ^ { - t - 1 } ) ^ { c } } | W ( \mathcal { Y }, \tilde { v } ) | > n \Big ) < C e ^ { - c n ^ { 1 / 3 } },$$

where the maximum is taken over all continuous segments from ˜ [ γ ˜ v t 1 , t 2 ] . Also

$$\mathbb { P } \Big ( \max _ { \mathcal { Y } \subseteq \tilde { \gamma } _ { \delta } \cap B ( \tilde { v }, r _ { \delta } e ^ { - 1 } ) ^ { e } } | W ( \mathcal { Y }, \tilde { v } ) | > n \Big ) < C e ^ { - c n }.$$

A few remarks are in order. Firstly, the stretched exponential tail is an artefact of the proof and we believe that an exponential tail bound could be proved with more care if necessary. Secondly, note that the intersection in the argument of the max above is used so that we don't need to look into what happens very close to ˜. v This is a technicality which simplifies the proof, but later it is not going to matter. In the end, we can decompose the whole path ˜ γ ˜ v into a disjoint union over t of [ t ′ 1 , t ′ 2 ] where t ′ 1 is last exit of B v, e (˜ - -t 1 ) and t ′ 2 is the last exit of B v, e (˜ -t ) which will accomplish the desired moment bound of truncated winding by exploiting the moment bound for each of these segments. We also emphasise that, although the non-contractible loop in the component containing v could be a branch from the skeleton, ˜ γ ˜ v itself does not contain any portion of the skeleton. This is simply because we have sampled the skeleton first before defining γ v .

Proof of Lemma 6.7. The main idea is to use a union bound for the winding of each excursion between annuli. One difficulty with working with ˜ γ ˜ v is that if we condition on where or when a non-contractible loop is formed then we break the independence between the pieces of random walk. Indeed, a walk conditioned on not forming a non-contractible loop will avoid certain portions of its previous trajectory, which a priori might bias it to have a lot of winding. Recall

that at the point when a non-contractible loop is formed, the walk could be very close to the starting point. Thus the idea is to run the random walk until the boundary ( ∂ Γ ) ′ # δ is hit (which could be potentially longer than dictated by Wilson's algorithm) and control the winding of its loop erasure uniformly over all scales using Lemmas 6.5 and 6.6. After this, we need to separately compare the winding made by the last excursion into the ball to the previous estimate.

Recall that t 0 satisfies e -t 0 = ˜ . Note that (6.2) follows easily from Lemma 6.5. Indeed, let r ˜ v us call τ nc the first time where a non-contractible loop is created and τ ∂ be the first hitting of the boundary. Fix N = max { j : τ i j &lt; τ nc ∧ τ ∂ } . Recall that N has exponential tail since every time we hit C 0 , conditioned on what happened before, we have a positive probability to create a non-contractible loop or hit the boundary before hitting C 1 , by Lemma 5.2. Thus we can work on the event N ≤ n at a cost which is exponentially small in n . Now Lemma 6.5 entails that on each of the n pieces, the winding has uniform exponential tail. This completes the proof of (6.2). We emphasise here that to bound portions of the loop erasure outside radius e -1 , we can run the random walk until it hits the boundary ( ∂ Γ ) ′ # δ and not worry about the case when the walk stops deep inside the ball, as the loop erasure outside radius e -1 is already taken care of.

Now we turn to (6.1). As before, let S := { ˜ X τ k } k ≥ 1 . The idea is to compare the winding of ˜ γ ˜ v to the winding of ˜ Y τ i m because there the conditioning on S preserves the independence and Lemma 6.6 controls the winding. The main work will be to do a case-by-case analysis of what can be erased and created between τ i m and the creation of a non-contractible loop.

Let N be as above and note that using the same idea of exponential tail of N and exponential tail of winding up to a fixed number of hits of the outermost circle (Lemma 6.6), we get

$$\mathbb { P } \Big ( \sup _ { y \subseteq \tilde { Y } ^ { \tau _ { i } } _ { N } [ t _ { 1 }, t _ { 2 } ] \cap B ( \tilde { v }, e ^ { - t - 1 } ) ^ { c } } | W ( \mathcal { Y }, \tilde { v } ) | > n \Big ) < C e ^ { - c ^ { \prime } n }.$$

Let λ 1 and λ 2 be respectively the first time ˜ Y τ i N exits B v, e (˜ - -t 1 ) and the last time ˜ Y τ i N exits B v, e (˜ -t ). To ease notations, from now on, all maximums in this proof come with the additional condition that the paths stay outside B v, e (˜ - -t 1 ) without writing it explicitly. We will also assume without loss of generality that the parametrisations of ˜ γ ˜ v and ˜ Y τ i N are identical up to the first point where their traces differ.

If a non-contractible cycle is created or the boundary is hit in the interval [ τ i N , τ ( i N +1) ] then, since ˜ [ X τ i N , τ ( i N +1) ] does not intersect B v, e (˜ -t 0 -1 ), only pieces of ˜ Y τ i N [0 , λ 2 ] can be erased and nothing can be added in the time interval [ τ i N , τ ( i N +1) ] to get ˜ γ ˜ v [0 , t 2 ]. In particular we see that in this case

$$\max _ { t _ { 1 } < s _ { 1 } < s _ { 2 } < t _ { 2 } } | W ( \tilde { \gamma } _ { \tilde { v } } ( s _ { 1 }, s _ { 2 } ), \tilde { v } ) | \leq \max _ { \lambda _ { 1 } < s _ { 1 } < s _ { 2 } < \lambda _ { 2 } } | W ( \tilde { Y } ^ { \tau _ { i } } ( s _ { 1 }, s _ { 2 } ), \tilde { v } ) |,$$

and we are done using (6.3).

The only other possibility is that a non-contractible cycle is created between τ ( i N +1) and τ i ( N +1) (otherwise that would contradict the maximality of N ). Since p is injective in B v, e (˜ -t 0 ), this occurs if and only if the walk hits a copy inside B v, e (˜ -t 0 ) of a portion of the walk that is further away from ˜, as illustrated schematically in Figure 14. v From now on we assume we are on this event.

First we make a topological observation. Let β be the first exit time of B v, e (˜ -t 0 ) by ˜ Y τ i N . We claim that ˜ [ X τ i N +1 , τ nc ] cannot hit ˜ Y τ i N [0 , β ]. Indeed, suppose by contradiction that it does so at some time T ∈ [ τ i N +1 , τ nc ]. Then, after erasure, we are left with a path completely contained in B v, e (˜ -t 0 ) where p is injective. No non-contractible loop can then be created before τ i N +1 , which contradicts the maximality of N as explained above.

Let k nc be the index during which the non-contractible loop happens, i.e, the unique index k such that τ nc ∈ [ τ k , τ k +1 ]. By the topological claim in the previous paragraph, no full turn can occur in any of the intervals [ τ k , τ k +1 ] for i N ≤ k &lt; k nc . However, by Corollary 4.5 in [4], a

full turn can occur in any interval [ τ k , τ k +1 ] independently with uniformly positive probability given S . Hence, there exist constants c, C depending only on the crossing estimate such that

$$\mathbb { P } ( k _ { n c } - i _ { N } > n ) \leq \mathbb { P } ( k _ { n c } - i _ { N } > n, N \leq n ) + \mathbb { P } ( N > n ) \leq C e ^ { - c n }.$$

Combining the above with Lemma 4.7 in [4] (which bounds the winding of the random walk during an interval of the type [ τ , τ i i +1 ]), we obtain the following stretched exponential tail bound:

$$\mathbb { P } ( \max _ { \mathcal { Y } \subset \tilde { X } [ \tau _ { i _ { N } }, \tau _ { n c } ] } | W ( \mathcal { Y }, \tilde { v } ) | > n ) \leq C e ^ { - c \sqrt { n } }$$

for some C, c &gt; 0 depending only on the crossing estimate and where, as before, Y is any continuous portion obtained from X τ [ i N , τ nc ] which preserves the order of the random walk path. In particular, this gives a good control of the winding of the piece of ˜ γ ˜ v added to Y τ i N .

Figure 14: The concatenation of the red and blue curves is ˜ Y τ i N and the green part is ˜ [ X τ i N , τ nc ] . The dotted path denotes a copy of ˜ Y τ i N so overall the event represented is one where a non-contractible loop is formed inside B v, e (˜ -t 0 ) .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00832v3-figure-9f3b14dd-938a-4116-80ec-b4a3d19253a4.jpg)

However, we are still not done as illustrated in Figure 14: the times t 1 , t 2 (which were first entry, last exit times for ˜ γ ˜ v ) may be different from λ , λ 1 2 (which were the first entry, last exit of ˜ Y τ i N ), so we need additional arguments. Let

$$\lambda _ { E } \coloneqq \inf \{ \lambda \, \colon \exists t \in [ \tau _ { i _ { N } }, \tau _ { n c } ], \tilde { Y } ^ { \tau _ { i _ { N } } } ( \lambda ) = \tilde { X } ( t ) \}$$

be the last time in ˜ Y τ i N which is not erased. Note that with this notation we showed above that λ E ≥ β . This implies immediately that t 1 = λ 1 . Now we need to consider several cases depending on where λ E is in relation to λ 2 .

- · If λ E ∈ [ λ , λ 1 2 ], then ˜ [ γ ˜ v t 1 , t 2 ] can be decomposed as the union of a piece of ˜ Y τ i N [ λ , λ 1 2 ] and some erasure of X τ [ i N , τ nc ], so on that event

$$\max _ { t _ { 1 } < s _ { 1 } < s _ { 2 } < t _ { 2 } } | W ( \tilde { \gamma } _ { \tilde { v } } ( s _ { 1 }, s _ { 2 } ), \tilde { v } ) | \leq \max _ { \lambda _ { 1 } < s _ { 1 } < s _ { 2 } < \lambda _ { 2 } } | W ( \tilde { Y } ^ { \tau _ { i _ { N } } } ( s _ { 1 }, s _ { 2 } ), \tilde { v } ) | + \Delta,$$

where ∆ is a variable with stretched exponential tail using (6.4).

- · If λ E ≥ λ 2 and ˜ ( γ ˜ v λ E , t 2 ] (the loop erasure after λ E ) does not enter B v, e (˜ -t ) then ˜ [ γ ˜ v t 1 , t 2 ] = ˜ Y τ i N [ λ , λ 1 2 ] and there is nothing more to prove.
- · If λ E ≥ λ 2 and ˜ ( γ ˜ v λ E , t 2 ] enters B v, e (˜ -t ) (this is the case pictured in Figure 14), then we decompose ˜ γ ˜ v as the union of ˜ Y τ i N [ λ , λ 1 2 ], ˜ Y τ i N [ λ , λ 2 E ] and ˜ [ γ ˜ v λ E , t 2 ]. The winding of any continuous portion of the first part has exponential tail by (6.3). The winding of any continuous portion of the last part has stretched exponential tail since it is bounded by ∆ as in the first case. So we only need to take care of the middle part.

To that end, we decompose ˜ Y τ i N [ λ , λ 2 E ] into excursions inside and outside B 0 as follows:

$$\tilde { Y } ^ { \tau _ { i _ { N } } } \left [ \lambda _ { 2 }, \lambda _ { E } \right ] = \tilde { Y } ^ { \tau _ { i _ { N } } } \left [ \lambda _ { 2 }, \beta _ { 0 } \right ] \cup \tilde { Y } ^ { \tau _ { i _ { N } } } \left [ \beta _ { 0 }, \beta _ { 1 } \right ] \cup \tilde { Y } ^ { \tau _ { i _ { N } } } \left [ \beta _ { 1 }, \beta _ { 2 } \right ] \cup \dots \cup \tilde { Y } ^ { \tau _ { i _ { N } } } \left [ \beta _ { k _ { 0 } }, \lambda _ { E } \right ], \quad ( 6. 5 )$$

where β 0 = β and for k ≥ 1, β 2 k is the first exit of B v, e (˜ -t 0 ) after β 2 k -1 and β 2 k -1 is the first entry into B v, e (˜ -t 0 -1 ) after β 2 k -2 . Note that any portion of ˜ Y τ i N [ β 2 k , β 2 k +1 ] is outside B v, e (˜ -t 0 -1 ) so its winding has exponential tail using Lemma 6.5. Note also that ˜ Y τ i N [ β 2 k +1 , β 2 k +2 ] lies inside the annulus bounded between C t and C 0 and never intersects the curve ˜ [ X τ i N , τ nc ] by maximality of λ E . Let Y ′ be the loop erasure of the portion of ˜ X from τ i N until its first hit of C t . Hence using Lemma 6.4

$$\max _ { \beta _ { 2 k } < s _ { 1 } < s _ { 2 } < \beta _ { 2 k + 1 } } | W ( \tilde { Y } ^ { \tau _ { i } } ( s _ { 1 }, s _ { 2 } ), \tilde { v } ) | \leq \max _ { s _ { 1 } < s _ { 2 } } | W ( Y ^ { \prime } [ s _ { 1 }, s _ { 2 } ], \tilde { v } ) | + 4 \pi.$$

Note that the right hand side has stretched exponential tail by (6.4). Thus each term in the decomposition (6.5) has stretched exponential tail and clearly the number of terms is bounded by N which itself has exponential tail. Combining these two, we can conclude that the absolute value of winding of any continuous portion of ˜ Y τ i N [ λ , λ 2 E ] has stretched exponential tail (with exponent 1 / 3), thereby completing the proof of this case.

This concludes the proof of Lemma 6.7.

## 6.3 From partial path to the full spine.

As shown in Theorem 4.10, the path ˜ γ ˜ v defined above is only a part of what is needed to compute the height function. Recall the notion of spine from Section 4.2.and that when we follow the outgoing edges in the Temperleyan forest from any given ˜, we always end up in a unique spine v (since boundary loops have been added around every hole).

Recall from Theorem 4.10 that we are interested in the winding of the path starting from ˜, v and then moving along the spine to infinity or the boundary of the disc. Observe that the initial portion of this path is ˜ γ ˜ v (Figure 13) and then it moves along copies of the non-contractible loop in the component of ˜ in the Temperleyan CRSF. In this section, we will call this path v p ˜ v (i.e., ˜ γ ˜ v followed by a semi-infinite piece of the spine). Observe that this notation is in contrast to the notation used in Theorem 4.10, where the path p ˜ v was called γ f , whereas here we emphasise that ˜ γ ˜ v denotes what is sampled from Wilson's algorithm given the skeleton. When we apply Wilson's algorithm to sample ˜ γ ˜ v , given the skeleton, we may form a new non-contractible loop, in which case we have discovered a portion of a spine (the unique spine attached to ˜). v The winding of this portion is controlled by Lemma 6.7. However, we also need to control the winding of the rest of p ˜ v . The purpose of this section is precisely to achieve this (done in Lemma 6.8). We also need to control the winding of the other semi-infinite piece of the spine attached to ˜, v which is done in Lemma 6.9. Finally, if we want to control the height gap between ˜ and ˜, we u v need to control the winding of p ˜ u around ˜ as well, which is done in Lemma 6.10. v

Let us parametrise p ˜ v in [0 , ∞ ) in any continuous manner away from ˜. v Let t 0 be as in Section 6.2, i.e., e -t 0 = ˜ . For r ˜ v t ≥ t 0 , let t 1 be the first exit time of p ˜ v from B v, e (˜ - -t 1 ) and let t 2 be the last exit time of the same from B v, e (˜ -t ). Note again that here, it is possible that the unique spine attached to ˜ could be the same as the spine corresponding to a skeleton branch. v

Lemma 6.8. For all k ≥ 1 , there exists a constant m&gt; 0 so that for all δ &lt; δ ¯ (˜ B v,e -t 0 ) and for all t 0 ≤ t &lt; log( C δ /δ ′ 0 ) ,

$$\mathbb { E } \left ( \max _ { \mathcal { Y } \subseteq p _ { \tilde { v } } [ t _ { 1 }, t _ { 2 } ] \cap B ( \tilde { v }, e ^ { - t - 1 } ) ^ { c } } | W ( \mathcal { Y }, \tilde { v } ) | ^ { k } ) \right ) \leq m,$$

where the supremum is taken over all continuous segments. Also

$$\mathbb { E } \left ( \max _ { \mathcal { Y } \subseteq p _ { \vartheta } \cap B ( \tilde { v }, e ^ { - t _ { 0 } } ) ^ { c } } | W ( \mathcal { Y }, \tilde { v } ) | ^ { k } ) \right ) \leq m.$$

Proof. In this proof, we write p ˜ v = p to lighten notation. Let us first consider the case when p does not contain a portion of the spine corresponding to the skeleton.

We parametrise ˜ γ ˜ v in [0 , 1] as before and assume that the parametrisation of ˜ γ ˜ v and p is the same until they start to differ. We drop ˜ and write ˜ v γ ˜ v = γ throughout the rest of this proof for notational clarity. To differentiate the first and last exit times for these two parametrisations we write t 1 ( p ) (resp. t 2 ( p )) for the first exit of B v, e (˜ - -t 1 ) (resp. last exit of B v, e (˜ -t )) of p . We also write t 1 = t 1 ( γ ) and t 2 = t 2 ( γ ) for clarity.

Note that it is always the case that t 1 ( γ ) = t 1 ( p ) and that t 2 ( γ ) ≤ t 2 ( p ). On the event t 2 ( p ) = t 2 ( γ ), the maximum in this lemma is over the same set as the one in Lemma 6.7 so there is nothing to prove. We focus now on the case where t 2 ( p ) &gt; t 2 ( γ ), meaning that a part of the spine comes back close to x on the cover ˜ M ′ after the time that a non-contractible loop was created.

We start with the case of the torus. Let us denote by S the spine attached to ˜ (meaning the v full bi-infinite path). Since S is a periodic path, it has a well defined direction d (say represented by a unit vector in R 2 ) and separates the plane in two sets right and left of S . Let us assume without loss of generality that the direction d is horizontal and that ˜ is below v S . Let τ be a time at which the vertical coordinate of γ reaches its maximum (for topological reasons, this has to occur on the spine S ) and let us define γ ′ by appending to γ [0 , τ ] a vertical segment going up to infinity.

Now it is easy to see that the winding of a straight segment around a point is bounded by π , so

$$\max _ { \mathcal { Y } \subseteq \gamma ^ { \prime } \cap B ( \tilde { v }, e ^ { - t - 1 } ) ^ { c } } | W ( \mathcal { Y }, \tilde { v } ) | \leq \max _ { \mathcal { Y } \subseteq \gamma \cap B ( \tilde { v }, e ^ { - t - 1 } ) ^ { c } } | W ( \mathcal { Y }, \tilde { v } ) | + \pi. \quad \quad ( 6. 6 )$$

Using Lemma 6.4, this implies that we can always bound

$$\max _ { \Psi \subseteq \mathfrak { P } _ { 0 } \cap B ( \tilde { v }, \varepsilon ^ { t - 1 } ) ^ { c } } | W ( \mathcal { Y }, \tilde { v } ) | \leq \max _ { \mathcal { Y } \subseteq \gamma ^ { \prime } \cap B ( \tilde { v }, \varepsilon ^ { t - 1 } ) ^ { c } } | W ( \mathcal { Y }, \tilde { v } ) | + 4 \pi \leq \max _ { \mathcal { Y } \subseteq \gamma \cap B ( \tilde { v }, \varepsilon ^ { t - 1 } ) ^ { c } } | W ( \mathcal { Y }, \tilde { v } ) | + 5 \pi. \\ \varepsilon ^ { \prime } \quad \varepsilon ^ { \prime } \quad \varepsilon ^ { \prime } \quad \varepsilon ^ { \prime } \quad \varepsilon ^ { \prime } \quad \varepsilon ^ { \prime } \quad \varepsilon ^ { \prime } \quad \varepsilon ^ { \prime } \quad \varepsilon ^ { \prime } \quad \varepsilon ^ { \prime } \quad \varepsilon ^ { \prime } \quad \varepsilon ^ { \prime } \quad \varepsilon ^ { \prime } \quad \varepsilon ^ { \prime } \quad \varepsilon ^ { \prime } \quad \varepsilon ^ { \prime } \quad \varepsilon ^ { \prime } \quad \varepsilon ^ { ( 6. 7 ) }$$

By Lemma 6.7 (the ˜ γ ˜ v there is γ here) applied to all scales up to scale t , the right hand side is a sum of O t ( ) variables with uniform stretched exponential tails.

Let us now bound P ( t 2 ( p ) &gt; t 2 ( γ ); p ∩ s = ). Note that on the event ∅ t 2 ( p ) &gt; t 2 ( γ ), ˜ has to v be at distance less that e -t of its spine. Using the fact that we can sample points in any order in Wilson's algorithm, we can bound that probability by first doing a cutset exploration around ˜ at scale v r ˜ v . After this cutset exploration, the probability that any of the branches sampled intersects B v, e (˜ -t ) is at most Ce -c t ( -t 0 ) for some constants C, c (see Lemma 5.7). However, after the cutset exploration around ˜, the spine attached to ˜ is necessarily fully sampled (see v v e.g., Proposition 5.10, final assertion). Hence

$$\mathbb { P } ( t _ { 2 } ( p ) > t _ { 2 } ( \gamma ) ; p \cap \mathfrak { s } = \emptyset ) \leq C e ^ { - c ( t - t _ { 0 } ) }.$$

Thus overall, either max Y⊆ p [ t 1 ,t 2 ] ∩ B v,e (˜ - -t 1 ) c | W ( Y , v ˜) | is the same as the variable in Lemma 6.7, or with probability at most Ce -c t ( -t 0 ) , it is at most a sum of O t ( ) variables with uniform stretched exponential tail. The moment bound now easily follows.

For the case where the universal cover is the unit disc, a similar argument can be done provided we can construct a path γ ′ which is nice and avoids the spine. Let us also introduce τ S as the first time γ hits the spine and τ nc as the ending time of γ . As recalled in more details in appendix A, it is a well known fact from Riemann surfaces that there exists a M¨bius transform o ϕ = ϕ γ,v ˜ such that

$$\mathfrak { p } ( v ) = \gamma [ 0, \tau _ { S } ] \cup \bigcup _ { n \geq 0 } \phi ^ { ( n ) } ( ( \tau _ { S }, \tau _ { n c } ] )$$

where ϕ ( n ) = ϕ ◦ . . . ◦ ϕ and the union is disjoint. In other words, we keep applying the same M¨ obius transform to obtain all the copies. Furthermore ϕ can be written as ϕ = Φ -1 ◦ µ ◦ Φ where Φ is a M¨bius map from the unit disc to the upper half plane and o µ is either the translation by ± 1 or a multiplication by λ &gt; 1.

If it is a translation we use the same argument as before. Otherwise, it is a scaling by λ &gt; 1 in which case ϕ S ( ) necessarily converges to infinity. Furthermore, let τ be the time at which ℑ (Φ( γ )) reaches its minimum, we define γ ′ by appending a straight vertical segment ℓ from Φ( γ τ ( )) to R (note that this segment may not intersect Φ( γ ) nor the subsequent scalings of the image of the portion of the spine Φ( γ ([ τ S , τ nc ])) and then mapping Φ( γ [0 , τ ] ∪ ℓ ) back to the disc by Φ -1 . By construction, it is trivial to check that Φ( γ ′ ) does not intersects Φ( p ( τ, ∞ )) so γ ′ does not intersect p ( τ, ∞ ). On the other hand, Φ -1 is the image of a line segment by a M¨ obius transform so it is a circular arc and therefore its winding around any point is bounded by 2 π . We can then conclude using exactly the same reasoning as in the torus case, with only the constant π replaced by 2 π on the right hand side of (6.6) and hence obtain the analogue of (6.7).

Finally, let us consider the case when p contains a portion of the spine corresponding to the skeleton. Note that the whole argument was geometric and the only place where we used the probability was in (6.8) which is provided by Lemma 5.6.

Lemma 6.9. For all k ≥ 1 , there exist a constant C &gt; 0 such that the following holds. Fix a compact set K ⊂ ˜ M ′ . Suppose we are in the setup of Lemma 6.8. Let p = p ˜ v be as above and ¯ p be the curve which starts at ˜ v , hits the spine, and then goes to infinity in the direction opposite to the orientation of the spine. Orient and parametrise both p p , ¯ from ˜ v to ∞ so that t 2 is the last time they exit B v, e (˜ -t ) . Then for all δ &lt; δ K , t ∈ [ t 0 , ∞ ] , v ˜ ∈ K ,

$$\mathbb { E } \left ( \sup _ { \mathcal { Y } \subset \mathfrak { p } [ t _ { 2 }, \infty ) } | W ( \mathcal { Y }, \tilde { v } ) | ^ { k } \right ) & \leq C ( ( 1 + t ) \wedge \log ( 1 / \delta ) ) ^ { k }, \\ \mathbb { E } \left ( \sup _ { \mathcal { Y } \subset \mathfrak { p } [ t _ { 2 }, \infty ) } | W ( \mathcal { Y }, \tilde { v } ) | ^ { k } \right ) & \leq C ( ( 1 + t ) \wedge \log ( 1 / \delta ) ) ^ { k }.$$

We emphasize here that the spine might be identical to a spine of a skeleton branch.

Proof. First we tackle the case of p . If t &lt; log( C δ /δ ′ 0 ), we need to add O t ( ) many variables given by Lemma 6.8 and hence we are done by Minkowski's inequality. If t &gt; log( C δ /δ ′ 0 ), first we apply the previous bound up to t = log( C δ /δ ′ 0 ). For the remainder, we can simply bound the winding by the volume of the ball of radius C δ/δ ′ 0 around ˜ which is v O (1) by our assumptions on the graph (see assumption (i) in Section 2.4). Indeed, this quantity is bounded by the number of times p crosses a straight line joining ˜ to the boundary of the ball, which is v simply bounded by the volume of the ball. For ¯ p , we can use Lemma 6.4 and the bound for p and the proof is complete.

Lemma 6.10. For all k there exist C, c &gt; 0 such that for all points ˜ ˜ u, v , for all compact sets ˜ K containing ˜ ˜ u, v ∈ ˜ M ′ , for all δ &lt; δ ˜ K ,

$$\mathbb { E } \left ( | W ( \mathfrak { p } _ { \tilde { u } }, \tilde { v } ) | ^ { k } \right ) \leq C ( | 1 + \log | \tilde { u } - \tilde { v } | | \wedge \log ( 1 / \delta ) ) ^ { c }.$$

Proof. Fix T = log( C δ ′ 0 δ ) with C ′ as in Lemma 6.8. We parametrise p ˜ v , p ˜ u as in Lemma 6.9. If | ˜ u - | ˜ v &lt; e -T , we use Lemmas 6.4 and 6.9 to conclude.

If | ˜ u - | ≥ ˜ v e -T , we take t = -log | ˜ u - | ˜ . v Then using Lemma 6.4, we can write

$$| W ( \mathfrak { p } _ { \tilde { u } }, \tilde { v } ) | \leq 2 \pi + \sum _ { k = t } ^ { \infty } \sup _ { \mathcal { Y } \subset \mathfrak { p } _ { \tilde { v } } [ k _ { 2 }, \infty ) } | W ( \mathcal { Y }, \tilde { v } ) | 1 _ { | \mathfrak { p } _ { \tilde { u } } - \tilde { v } | \in [ e ^ { - k }, e ^ { - k - 1 } ] }$$

where k 2 is the last exit from the B v, e (˜ -k ) by p ˜ u . Using Lemma 5.4 or Lemma 5.6 whichever is appropriate, we get an exponential bound on the expectation of the indicator event above (notice for scales k &gt; log( C δ /δ ′ 0 ), the probability actually becomes 0 so this is a finite sum). Therefore, we can again use Lemma 6.9 and Cauchy-Schwarz to conclude.

## 6.4 Convergence of height function

In this section we prove Theorem 6.1.

Let us describe informally the general structure of the proof. Most of the work will be in the universal cover to obtain the convergence of expressions of the form E ∏ i ( h z (˜ ) i -h w ( ˜ )) for i distinct points ˜ z i and ˜ w i (Lemma 6.12). By integrating this expression, we will then obtain the convergence of the height function, seen as a function on the universal cover.

For the first part, the idea is to introduce a regularised height function h # δ t which is a continuous function of the CRSF so that the convergence of h # δ t as δ → 0 is immediate. This leaves us with two issues: the comparison between h # δ and h # δ t for fixed δ and the convergence of h t as t →∞ in the limit. Both questions are actually solved simultaneously by the estimate of Lemma 6.12 which compares h # δ and h # δ t with an error term that becomes small both in δ and t independently.

Setup and notations. Recall that our goal is to prove that

$$\int _ { \tilde { M } ^ { \prime } } \bar { h } _ { \text{ext} } ^ { \# \delta } ( z ) \tilde { f } ( z ) d \mu ( z )$$

converges in law and also in the sense of moments as δ → 0. (Recall that ¯ denotes the centered X random variable ¯ X = X -E ( X ) whenever this expectation is well defined). Implicitly, h # δ is sampled according to the dimer law (1.1). In other words, by Theorem 3.3, we need to first sample a Temperleyan pair ( T , T † ) and apply the bijection ψ in that theorem.

However, it turns out to be more convenient in the proof to work with a pair ( T , T † ) sampled from P Wils . To this pair ( T , T † ) we can apply the bijection ψ from Theorem 3.3 and obtain a dimer configuration m whose centered height function can be studied. We will first prove (6.9) for P Wils and then later explain how this implies the same for P Temp .

Recall that since ∫ ˜ M ′ ˜ ( f z dµ z ) ( ) = 0, the expression in (6.9) is in fact well defined. We wish to compute the moments of this integral but only in terms of height differences since the field is a priori defined only up to constant. To do this, we use the following trick. Note that since ∫ ˜ fdµ = 0, ∫ ˜ f + dµ = ∫ ˜ f -dµ =: Z f ( ˜ ) where f + = max { f, 0 } and f -= -min { f, 0 } . Now we can write

$$\int _ { \tilde { K } } \bar { h } _ { \text{ext} } ^ { \# \delta } ( z ) \tilde { f } ( z ) d \mu ( z ) = \int _ { \tilde { K } \times \tilde { K } } ( \bar { h } _ { \text{ext} } ^ { \# \delta } ( z ) - \bar { h } _ { \text{ext} } ^ { \# \delta } ( w ) ) \frac { \tilde { f } ^ { + } ( z ) \tilde { f } ^ { - } ( w ) } { Z ( \tilde { f } ) } d \mu ( z ) d \mu ( w ),$$

which implies

$$( \int _ { \tilde { K } } \overline { h } _ { \text{ext} } ^ { \# \delta } ( z ) \tilde { f } ( z ) d \mu ( z ) ) ^ { k } = \int _ { \tilde { K } ^ { 2 k } } \prod _ { i = 1 } ^ { k } ( \tilde { h } _ { \text{ext} } ^ { \# \delta } ( z _ { i } ) - \tilde { h } _ { \text{ext} } ^ { \# \delta } ( w _ { i } ) ) \frac { \tilde { f } ^ { + } ( z _ { i } ) \tilde { f } ^ { - } ( w _ { i } ) } { Z ( \tilde { f } ) } d \mu ( z _ { i } ) d \mu ( w _ { i } ). \quad ( 6. 1 1 )$$

Therefore we are interested in the k -point function, E [ ∏ k i =1 ( h # δ ext ( z i ) -h # δ ext ( w i ))]. Pick k distinct pairs of points ( z , w 1 1 ) , . . . , ( z k , w k ) ∈ ˜ K and let ( f ( z 1 ) , f ( w 1 )) , . . . , ( f ( z k ) , f ( w k )) be the faces containing them. Let z # δ i , w # δ i be the midpoint of the diagonals of f ( z i ) , f ( w i ).

Now recall Theorem 4.10 which relates the dimer height difference between two faces f and f ′ to the winding of a specific path γ f,f ′ connecting m f ( ) and m f ( ′ ) and additional terms of the form ± π associated with jumping over components of the CRSF. Let γ # δ i be the path γ f,f ′ when f = f ( z i ) and f ′ = f ( w i ), and orient it from z # δ i to w # δ i . Then with these notations, Theorem 4.10 says that

$$h ( z _ { i } ^ { \# \delta } ) - h ( w _ { i } ^ { \# \delta } ) = W ( \gamma _ { i } ^ { \# \delta }, z _ { i } ^ { \# \delta } ) + W ( \gamma _ { i } ^ { \# \delta }, w _ { i } ^ { \# \delta } ) + \Psi _ { i } ^ { \# \delta }$$

where Ψ # δ i is the π ∑ S ( ε S + δ S ) term in Theorem 4.10. We drop the superscript δ from now on for clarity. Thus from now on we focus on proving convergence of

$$\mathbb { E } _ { W i l s } \left ( \prod _ { i = 1 } ^ { k } ( \bar { W } ( \gamma _ { i }, z _ { i } ) + \bar { W } ( \gamma _ { i }, w _ { i } ) + \bar { \Psi } _ { i } ) \right ).$$

Let

̸

X = { z , w i i : 1 ≤ i ≤ k } , and assume p u ( ) = p v ( ) for all u = v ∈ X .

̸

Clearly, (6.13) can be expanded as a sum of 2 k many terms of the form:

$$\mathbb { E } _ { W i l s } ( \prod _ { x \in S } ( \bar { W } ( \gamma _ { x }, x ) + \bar { \Psi } _ { x } / 2 ) ) = \colon \mathbb { E } _ { W i l s } ( \prod _ { x \in S } \bar { F } _ { x } ( \gamma _ { x } ) )$$

where S is a subset of vertices of X with distinct indices and γ x is γ i for some i such that x = z i or w i (of course, the products in the expansion of (6.14) has further restrictions, but we ignore that for clarity).

We are interested in estimating and proving convergence of (6.14). To that end we employ an idea similar to that in [4]: we truncate the CRSF branch at a macroscopic scale and deal with the truncated macroscopic winding and the remaining microscopic winding part separately.

We now fill in the details. Parametrise γ = γ x . Define ( γ x ( )) t t ≥ 1 to be γ x from the first entry time into the ball B x,e ( -t ) (note that at the moment, the balls B x,e ( -t ) might overlap for different x ∈ X .) We emphasise here that we parametrise γ in the opposite direction compared to Sections 6.2 and 6.3 to be consistent with [4]. With an abuse of notation, we will denote by γ x ( -∞ , t ] the whole path from the opposite end of γ up until γ x ( ). t Define the regularised term and the error term as

$$\bar { F } _ { x } ( \gamma _ { x }, t ) \coloneqq \bar { W } ( \gamma _ { x } ( - \infty, t ] ) ; x ) + \bar { \Psi } _ { x } / 2$$

$$\bar { e } _ { x } ( t ) \coloneqq \bar { W } ( \gamma _ { x }, x ) - \bar { W } ( \gamma _ { x } ( - \infty, t ] ), x ) = \bar { F } _ { x } ( \gamma _ { x } ) - \bar { F } _ { x } ( \gamma _ { x }, t ). \quad \quad ( 6. 1 6 )$$

When we want to emphasise the role of z, w , we write γ zw , Ψ zw , F z ( γ zw ) in place of the above. We start with a general moment bound for the truncated winding.

Lemma 6.11. For all m , there exists constants c = c K, m , α ( ˜ ) = α m ( ) such that for all z, w, δ &lt; δ ˜ K , t ≥ 0 ,

$$\mathbb { E } _ { \text{Wils} } ( | \bar { F } _ { z } ( \gamma _ { z w }, t ) | ^ { m } ) \leq c | ( 1 + t + | \log | z - w | | ) \wedge \log ( 1 / \delta ) | ^ { \alpha }.$$

Proof. The proof is immediate from Lemmas 6.9 and 6.10 and the fact that Ψ zw is bounded by a constant c K ( ) times the number of non-contractible loops in the CRSF. Indeed, this is clear from the fact that for a fixed compact set ˜ K ⊂ ˜ M there exists a constant c K ( ˜ ) such that any curve P ⊂ ˜ K will cross at most c K ( ˜ ) copies of a given spine (or lift of a loop). Since the number of non-contractible loops has superexponential tail by Theorem 6.3, the moment bound of Ψ zw is immediate.

̸

$$\Delta = \Delta ( \mathcal { X } ) = \frac { 1 } { 1 0 } ( r _ { \mathcal { X } } \wedge r _ { \tilde { K } } ). \quad \quad \quad$$

By compactness, choose r ˜ K small enough so that p is injective in B z, r ( ˜ K ) for all z ∈ ˜ . K Now let r X be defined as in (5.2) but for the set of vertices p ( X ) (which are all distinct by assumption on X ). We observe that r X ≥ c K ( ˜ ) min x = y ∈X | x - | y for some constant c K ( ˜ ). We set

Lemma 6.12. There exist constants c, c ′ &gt; 0 such that for all m,m ′ ≥ 1 there exists α &gt; 0 such that the following holds. Let S ⊂ X , S ′ ⊂ X be disjoint. Also assume that | S | = m, S | ′ | = m ′ ≥ 0 . Let ∆ be as in (6.17) . Then for all δ &lt; δ ˜ K and t ≥ log( r ˜ K / ∆)

$$| \mathbb { E } _ { W i s } ( \prod _ { x \in S } \bar { e } _ { x } ( t ) \prod _ { x \in S ^ { \prime } } \bar { F } _$$

Proof. For simplicity, we first assume S ′ is empty. Recall that we can sample a CRSF from P Wils by first sampling the branches s and then the rest by Wilson's algorithm. Note that Lemma 5.6 tells us that the isolation radius corresponding of each vertex in S with respect to the skeleton has polynomial tail. Since after sampling the skeleton, the rest of the branches are sampled simply by Wilson's algorithm, we can conclude that the application of Proposition 5.10 is valid with ∆ in place of r i .

We perform the coupling in Section 5.2 (in particular Proposition 5.10) with points p ( X ) and their lift X , and a compact domain D ⊂ ˜ M ′ containing all points in X so that the minimal Euclidean distance between any point in X and ∂D is at least r X . Note that we can choose r X for the r in (5.2) there. Call ( T D x ) x ∈X the resulting independent UST in D . For x ∈ X let γ D x be the branch in the UST T D x joining x to ∂D . Let γ D x ( ) be parametrised so that t γ D x enters B x,e ( -t r ˜ K ) for the first time (going from the outside to x ) at time t .

Let R x be the isolation radius of x in the application of Proposition 5.10 and write R x = e -I x r ˜ K . Let I = max x ∈X I x . (Note that for notational clarity, I here is shifted by log( r ˜ K / ∆) from the one in Section 5.2.) We now decompose

$$e _ { x } ( t ) = \underbrace { \overline { W } ( \gamma _ { x } ^ { D } ( t, \infty ), x ) } _ { \alpha _ { x } } + \underbrace { ($$

Therefore, we need to deal with expectation of products of α , ξ x x for different indices x . Let C be the sigma algebra generated by the cutset exploration. We will first compute the conditional expectation and then take the overall expectation. Note that T D x are independent for different x and also independent of C and hence any term involving α x is 0. Thus we only have to deal with terms involving ξ x , x ∈ S .

Let Λ be the last time γ x enters B x,e ( -I r ˜ K ). In the coupling, γ , γ x D x agree inside B x,e ( -I r ˜ K ) so in particular γ x (Λ , ∞ ) = γ D x (Λ , ∞ ), so

$$\xi _ { x } = \overline { W } ( \gamma _ { x } ( t, \Lambda ), x ) - \overline { W } ( \gamma _ { x } ^ { D } ( t, \Lamb$$

Let G be the event that ξ x without the bar (meaning without the expectation terms) is 0. Observe that if G does not occur then one of the following events happen. Either

$$I - \log ( r _ { \tilde { K } } / \Delta ) > ( t - \log ( r _ { \tilde { K } } / \Delta ) ) / 2$$

which has probability at most e -c t ( -log( r ˜ K / ∆)) ∨ δ c ≤ ( e -c t ( -log( r ˜ K / ∆)) + δ c ) (Proposition 5.10). Otherwise, γ x has to exit B x,e ( -I r ˜ K ) after hitting e -t r ˜ K . This also has probability at most e -c t ( -log( r ˜ K / ∆)) (Lemma 5.8). Now we bound the moments of ξ x . Notice that we can write

$$\sim \sim \sim \kappa ^ { \prime } \sim \kappa ^ { \prime } = \sim \nu ^ { \prime } ( \text{Lemma 0.5} ). \text{ now we found the moments or $\xi_{x}$}. \text{ in voice that we can write} \\ \mathbb { E } _ { \text{Wils} } ( | \xi _ { x } | ^ { k } ) & \leq \sum _ { j = t / 2 + \log ( r _ { \tilde { K } } / \Delta ) / 2 } ^ { \infty } \mathbb { E } ( | \xi _ { x } | ^ { k } 1 _ { j \leq I \leq j + 1 } ) + c e ^ { - c ( t - \log ( r _ { \tilde { K } } / \Delta ) ) } \\ & \leq \sum _ { j = t / 2 + \log ( r _ { \tilde { K } } / \Delta ) / 2 } ^ { \infty } | 1 + j | ^ { k } ( e ^ { - c ( j - \log ( r _ { \tilde { K } } / \Delta ) ) } + \delta ^ { c } )$$

where we first used Cauchy-Schwarz and then used Lemma 6.9 to bound the moment and Proposition 5.10. Thus overall

$$\mathbb { E } _ { \text{Wils} } ( | \xi _ { x } | ^ { k } ) \leq | 1 + t | ^ { k } ( e ^ { - c ( t - \log ( r _ { \tilde { K } } / \Delta ) ) } + \delta ^ { c } ).$$

We have a product of at most m terms and so using H¨lder's inequality, taking o α = m works.

Finally, if

S

′

is non-empty, the proof is exactly the same as the vertices in

S

′

are distinct from

S

and hence local independence from the coupling still holds. We then use Lemma 6.11 to conclude using Cauchy-Schwarz. Details are left to the reader.

Corollary 6.13. Let S ⊂ X containing vertices with distinct indices such that | S | = m . There exists a constant c = ( ˜ c K,m ,α ) = α m ( ) such that for all δ &lt; δ ˜ K ,

$$| \mathbb { E } _ { \text{Wils} } ( \prod _ { x \in S } \bar { F } _ { x } ( \gamma _ { x } ) ) | \leq c | 1 + \log ^ { \alpha } ( \Delta ) |.$$

Proof. Decompose

$$\mathbb { E } _ { W i l s } ( \prod _ { x \in S } \bar { F } _ { x } ( \gamma _ { x } ) ) = \mathbb { E } _ { W i l s } ( \prod _ { x \in S } ( \bar { F } _ { x } ( \gamma _ { x }, t ) + \bar { e } _ { x } ( t ) ) )$$

with t = 2log( r ˜ K / ∆). Then this is a straightforward application of Lemmas 6.11 and 6.12.

Now we prove a convergence of the height function integrated against f in the sense of moments (still for the Wilson law P Wils ). Recall that X δ converges in the sense of moments as δ → 0 if for all k , E ( X k δ ) converges as δ → 0.

Lemma 6.14. ∫ ˜ K h # δ ext ( z f ) ˜ ( z dµ z ) ( ) converges in the sense of moments under the law P Wils . Furthermore, the limit does not depend on the sequence ( G ′ ) # δ .

Proof. Using eqs. (6.11), (6.13) and (6.14), we need to prove the convergence as δ → 0 of

$$\int _ { \tilde { K } ^ { 2 k } } \mathbb { E } _ { \text{Wils} } ( \prod _ { x \in S } \bar { F } _ { x } ( \gamma _ { x } ^ { \# \delta } ) ) \frac { \tilde { f } ^ { + } ( z _ { i } ) \tilde { f } ^ { - } ( w _ { i } ) } { Z ( \tilde { f } ) } d \mu ( z _ { i } ) d \mu ( w _ { i } ).$$

̸

̸

↦

We use Fubini to bring the expectation inside the integral in the above display. We first observe that integrating (6.20) over all sets of vertices S so that p x ( ) = p y ( ) for all x = y ∈ S is enough (recall p : ˜ M ′ → M ′ is the covering map). Indeed, using Lemma 6.11 and the fact that ∥ f ∥ ∞ &lt; ∞ , we see that the integrand can be bounded by O (log α (1 /δ )). Since the volume of { S ⊂ ˜ K 2 k : p x ( ) = p y ( ) for some x, y ∈ S } is O δ ( ) (indeed the number of preimages of any vertex in ˜ K is bounded by a constant depending only on ˜ ), we see that the integral over this K set is O δ ( log α (1 /δ )).

Thus we now concentrate on the integral (6.20) over

̸

$$\mathcal { A } ( \tilde { K } ) = \{ \text{sets of vertices $S$ so that $p(x)\neq p(y)$ for all $x\neq y\in S$} \}. \quad \quad ( 6. 2 1 )$$

̸

We now use Corollary 6.13 and dominated convergence theorem. Since ∥ f ∥ ∞ &lt; ∞ and log m (∆) is integrable for any m&gt; 0, we need to prove convergence of the expectation inside the integral above. Now we claim that the regularised part

$$\mathbb { E } _ { W i l s } ( \prod _ { x \in S } \bar { F } _ { x } ( \gamma _ { x } ^ { \# \delta }, t ) )$$

converges as δ → 0. This follows from our assumption of a.s. convergence of Temperleyan CRSF and because the term inside the expectation is a.s. continuous function of the Temperleyan

CRSF. The same can be said about the skeleton using the absolute continuity statement given by Proposition 5.3. Indeed, this follows from a.s. continuity properties of SLE 2 and the fact that the CRSF branches are made a.s. from finitely many chunks of SLE , in particular for a 2 fixed t , the SLE curve is a.s. not a tangent to the boundary of circle of radius e -t (we refer to [5, Section 3]) for details. Lemma 6.11 tells us that the random variable in (6.22) is uniformly integrable which completes the proof of convergence of the regularised term (6.22).

Now given a fixed S , choose t ≥ log( r ˜ K / ∆). Now we claim that that the error term satisfies

$$| \mathbb { E } _ { W i l s } ( \prod _ { x \in S } \bar { F } _ { x } ( \gamma _ { x } ) ) - \mathbb { E } _ { W i l s } ( \prod _ { x \in S } \bar { F } _ { x } ( \gamma _ { x }, t ) ) | < c | 1 + t | ^ { \alpha } ( e ^ { - c ( t - \log ( r _ { \tilde { K } } / \Delta ) ) } + \delta ^ { c ^ { \prime } } ). \quad \ ( 6. 2 3 )$$

Indeed, (6.23) follows by writing ¯ ( F x γ x ) = ¯ ( F x γ , t x ) + ¯ ( ) and then expanding and using the e x t bounds in Lemmas 6.11 and 6.12. To finish the proof of the lemma, fix ε &gt; 0. First choose t large enough and then δ small enough so that the right hand side of (6.23) is less than ε . Next choose a smaller δ if needed so that for all δ ′ &lt; δ ,

$$| \mathbb { E } _ { \text{Wils} } ( \prod _ { x } \bar { F } ( \gamma _ { x } ^ { \# \delta ^ { \prime } }, t ) ) - \mathbb { E } _ { \text{Wils} } ( \prod _ { x } \bar { F } _ { x } ( \gamma _ { x } ^ { \# \delta }, t ) ) | < \varepsilon$$

via the convergence of the regularised term. This completes the proof as ε is arbitrary.

Observe that (6.23) completes the proof that the limit does not depend on the sequence (( G ′ ) # δ ) δ&gt; 0 since the main term (6.22) is measurable with respect to the limiting continuum Temperleyan CRSF which is universal by Theorem 6.3.

It remains to prove convergence in law (still for P Wils for now). For this, we need to alter the definition of truncation slightly. Given t &gt; 0, take a cover { B euc ( x, r ˜ K e -10 t ) : x ∈ ˜ K } of ˜ K and take a finite subcover. Let c ( z , ) c ( w ) denote the center of one of the balls (chosen arbitrarily) in the finite subcover in which z, w belong. Define F x ( ) t , e x ( ) to be the same as t F x ( ) and t e x ( ) t but we cut off the first time γ x enters B ( c ( x , e ) -t ) (so compared to the above, the centre of the cutoff ball is shifted by an amount which is at most e -10 t ). However, in this definition, we still measure winding around x not c ( x ).

Lemma 6.15. The statements of Lemmas 6.11 and 6.12 still hold if we replace F, e by F e , everywhere.

Proof. For Lemma 6.11, the proof identically follows from Lemmas 6.9 and 6.10 by noticing that the supremum over all continuous subpaths contain the portion of the branch until it first hits the shifted ball. Notice also that we shift only by an additive term in the exponential scale which is O t ( ).

For Lemma 6.12, the proof is also identical, in particular we still consider the coupling around the points in S . Because we still shift only by O t ( ) in the exponential scale, the proof readily follows.

We can now complete the proof of Theorem 6.1.

Proof of Theorem 6.1. Again, we first prove this under P Wils before explaining how to extend this and concluding the proof of Theorem 6.1 to P Temp . We write r = r ˜ K and pick a t &gt; 0 to be taken large later. Recall from eqs. (6.10) and (6.12),

$$\int _ { \tilde { K } } \bar { h } _ { \text{ext} } ^ { \# \delta } ( z ) \tilde { f } ( z ) d \mu ( z ) = \int _ { \tilde { K } \times \tilde { K } } ( \bar { W } ( \tilde { \gamma } _ { z w } ^ { \# \delta }, z ) + \bar { W } ( \tilde { \gamma } _ { z w } ^ { \# \delta }, w ) + \bar { \Psi } _ { z w } ) \frac { \tilde { f } ^ { + } ( z ) \tilde { f } ^ { - } ( w ) } { Z ( \tilde { f } ) } d \mu ( z ) d \mu ( w ).$$

Here Z f ( ˜ ) is a deterministic constant so we can assume it to be 1 from now on without loss of generality. Firstly we recall that we only need to integrate over the set A ( ˜ ) K as in (6.21) as

the integral over the remaining part is O δ ( log (1 c /δ )). Let us denote by Y = Y ( δ ) the above integral over A ( ˜ ). K Introduce

$$X ( \delta, t ) \coloneqq \int _ { \mathcal { A } ( \tilde { K } ) } ( \bar { F } _ { \epsilon ( z ) \epsilon ( w ) } ( t ) + \bar { F } _ { \epsilon ( w ) \epsilon ( z ) } ( t ) ) 1 _ { | \epsilon ( z ) - \epsilon ( w ) | > r e ^ { - t / 1 0 } } \tilde { f } ^ { + } ( z ) \tilde { f } ^ { - } ( w ) d \mu ( z ) d \mu ( w ).$$

Note that X δ, t ( ) is a sum of regularised winding of finitely many branches and hence converges in law (as δ → 0 and t is fixed) by our assumption.

Now we show that for all t &gt; 0 and δ &lt; C δ e ′ 0 -t ,

$$\mathbb { E } ( Y ( \delta ) - X ( \delta, t ) ) ^ { 2 } \leq c ( 1 + t ) ^ { \alpha } e ^ { - c ^ { \prime } t }.$$

We expand the above square to get an integral over ( z, w, z , w ′ ′ ) ∈ A ( ˜ ) K 2 . We can again without loss of generality restrict to the the set A 2 ( ˜ ) K such that the projection p of all four points ( z, w, z , w ′ ′ ) maps to pairwise distinct points on M . In that case let ∆ = ∆( z, w, z , w ′ ′ ) be as in (6.17). We now argue that this integral over the set ∆ ≤ e -t/ 11 is exponentially small in t . Indeed, we are integrating over a set which has exponentially small measure with respect to µ 4 and furthermore the moments are integrable using Corollary 6.13 and Lemma 6.15.

For the rest of the integral, we write

$$( \bar { W } ( \tilde { \gamma } _ { z w } ^ { \# \delta }, z ) + \bar { W } ( \tilde { \gamma } _ { z w } ^ { \# \delta }, w ) + \bar { \Psi } _ { z w } ) = \bar { F } _ { z w } ( t ) + \bar { F } _ { w z } ( t ) + \bar { e } _ { z w } ( t ) + \bar { e } _ { w z } ( t )$$

and we expand again the product inside the integral to get products of terms of the form

$$\bar { \mathrm e } _ { z w } ( t ) + \bar { \mathrm e } _ { w z } ( t ) ; \ \bar { F } _ { z w } ( t ) + \bar { F } _ { w z } ( t ) - \bar { F } _ { \mathrm c ( z ) \mathrm c ( w ) } ( t ) - \bar { F } _ { \mathrm c ( w ) \mathrm c ( z ) } ( t ).$$

Note here that the indicator over | c ( z ) -c ( w ) | &gt; e -t/ 10 is included in ∆ &gt; e -t/ 11 so we can get rid of the indicator. Products containing at least one e are small because of Lemma 6.15. Also note that on ∆ &gt; e -t/ 11 , | z -w &gt; c K e | ( ˜ ) -t/ 11 . Thus we need to show

$$\mathbb { E } ( | \bar { F } _ { z w } ( t ) - \bar { F } _ { \epsilon ( w ) \epsilon ( z ) } ( t ) | ^ { 2 } 1 _ { | z - w | > c ( \tilde { K } ) e ^ { - t / 1 1 } } ) \leq c ( 1 + t ) ^ { \alpha } e ^ { - c ^ { \prime } t }.$$

Observe that

$$( \bar { F } _ { z w } ( t ) - \bar { F } _ { \epsilon ( z ) \epsilon ( w ) } ( t ) ) 1 _ { | z - w | > c ( \bar { K } ) e ^ { - t / 1 1 } } = ( \bar { F } _ { z w } ( t ) - \bar { F } _ { \epsilon ( z ) w } ( t ) ) 1 _ { | z - w | > c ( \bar { K } ) e ^ { - t / 1 1 } } \\ & + ( \bar { F } _ { \epsilon ( z ) w } ( t ) - \bar { F } _ { \epsilon ( z ) \epsilon ( w ) } ( t ) ) 1 _ { | z - w | > c ( \bar { K } ) e ^ { - t / 1 1 } } ).$$

If p ( z ) and p ( c ( z )) merge before exiting B ( c ( z , e ) -5 t ), then the paths we need to consider are identical outside B ( c ( z , e ) -t ), only the centers around which we measure the topological winding are different. This difference in winding is therefore deterministically O e ( -9 t ). On the other hand, the paths we need to consider are not identical outside B ( c ( z , e ) -5 t ) if and only if a simple random walk from z does not hit p ( c ( z )) before exiting B ( c ( z , e ) -5 t ) which has probability exponentially small in t by a Beurling-type lemma (Lemma 5.1). Hence we can conclude using Cauchy-Schwarz and Lemma 6.15. (Note that here the shifted cutoff is particularly useful since the cutoff point is the same for both z and c ( z )).

The proof of the rest of the statements is a simple consequence of the fact that the main term X δ, t ( ) is an a.s. continuous function of T # δ for a fixed t . So trivially ( X δ, t , ( ) T # δ ) converges jointly in law, the limit does not depend on the sequence ( G ′ ) # δ and is conformally invariant. It is a simple exercise to show that the L 2 bound on the error term Y ( δ ) -X δ, t ( ) as shown above is enough to conclude the proof of the theorem (i.e. the proof of the same claim as in the previous sentence for the pair ( ¯ h # δ ext ( z , ) T # δ ) ).

We have proved convergence in law of ( ¯ h # δ ext ( z , ) T # δ ) for the P Wils law. We now explain how to convert the result so that it holds under the Temperleyan law P Temp . Recall that if

T is sampled by performing Wilson's algorithm (after sampling the skeleton) and T † is sampled from the uniform distribution among all oriented duals of T , then the Radon-Nikodym derivative of ( T , T † ) with respect to P Temp is Z = 2 K † / E Wils (2 K † ) where K † is the number of nontrivial cycles of T † by Lemma 3.8. Now, observe that Z is measurable with respect to T # δ so ( Z, ( ¯ h ext , f ˜ )) jointly converges in law under P Wils by the above. Furthermore all moments of Z and of ( ¯ h ext , f ˜ ) are bounded under P Wils and therefore we conclude that

$$\left ( \int _ { \tilde { M } } ( h _ { \text{ext} } ( z ) - \mathbb { E } _ { \text{Wils} } ( h _ { \text{ext} } ( z ) ) ) \tilde { f } ( z ) \mu ( d z ) ; \mathcal { T } ^ { \# \delta } \right )$$

converges in law and in the sense of moments under P Temp . In particular, taking expectation (under P Temp ) of the first quantity, we also deduce that ∫ ˜ M ( E Temp ( h ext ( z )) -E Wils ( h ext ( z ))) ˜ ( f z µ dz ) ( ) converges. Taking the difference, it follows that

$$\left ( \int _ { \tilde { M } } ( h _ { \text{ext} } ( z ) - \mathbb { E } _ { \text{Temp} } ( h _ { \text{ext} } ( z ) ) ) \tilde { f } ( z ) \mu ( d z ) ; \mathcal { T } ^ { \# \delta } \right )$$

converges in law and in the sense of moment under P Temp , as desired. Furthermore, since the limit of T # δ is independent of the chosen graph sequence (subject to the assumptions in Section 2.4), so is the limit of ∫ ˜ ( f x h ) ¯ # δ ext ( x ) d µ x ( ).

Remark 6.16. In fact, with some more work, the convergence in Theorem 6.1 can be upgraded to specifically include information about the instanton component. More precisely, using Lemma 4.13 we obtain the following.

Take an ordered finite set of continuous simple loops which forms the basis of the first homology group of M ′ , all endowed with a fixed orientation. Let H # δ ∈ R 4 g +2 b -4 denote the vector of height differences along these loops (i.e., for each such loop, record the height accumulated by going along the loop once in the prescribed orientation). Consider the one-form

$$d \bar { h } ^ { \# \delta } ( e ) \coloneqq \bar { h } ^ { \# \delta } ( \tilde { e } ^ { + } ) - \bar { h } ^ { \# \delta } ( \tilde { e } ^ { - } ).$$

It is well-known (see Theorem 2.9) that the instanton component of the above one-form is completely determined by H # δ . Then we have:

$$\left ( \int \tilde { f } ( x ) \bar { h } _ { \text{ext} } ^ { \# \delta } ( x ) \, \mathrm d \mu ( x ), \mathbf H ^ { \# \delta }, \mathcal { T } ^ { \# \delta } \right )$$

converges jointly in law as δ → 0. The first two coordinates also jointly converge in the sense of all moments. Furthermore,

$$\lim _ { \delta \rightarrow 0 } \left ( \int \tilde { f } ( x ) \bar { h } _ { \text{ext} } ^ { \# \delta } ( x ) \, \mathrm d \mu ( x ), \mathrm H ^ { \# \delta } \right )$$

is measurable with respect to the limit T of T # δ . We do not go into the details of these claims for the sake of brevity.

## A Geometry of spines

In this section we prove Lemma 4.9 and Lemma 6.5. But before getting into the proofs, we remind the readers certain basic facts from the classical theory of Riemann surfaces.

By the uniformisation theorem of Riemann surfaces (Section 2.2), recall that there exists a conformal map from the Riemann surface D / g to M ′ where g is a Fuchsian group which is a discrete subgroup of the group of M¨bius transforms. o Furthermore, such a conformal map

↦

is unique up to conformal automorphisms (i.e. M¨ obius transforms) of the unit disc. In other words, if g g , ′ are two Fuchsian groups such that M ′ is conformally equivalent to both D / g and D / g ′ then there exists a M¨bius map o ϕ : D → D such that g ′ = ϕ -1 ◦ g ◦ ϕ . Since we have fixed a canonical lift, we have defined Γ uniquely.

It is further known that g is isomorphic to the fundamental group π 1 ( M ,x ′ ) (topologically, this is also known as the group of deck transformations). This connection is described as follows: choose a base point in the manifold x 0 and a particular lift of ˜ x 0 , then any simple loop ℓ based in x 0 can be lifted to a simple curve ˜ in ℓ D starting from ˜ x 0 , with the endpoint ˜ x 1 depending only on the homotopy class of the loop. Then to ℓ we associate the map ϕ ℓ,x ˜ 0 ∈ g that sends ˜ x 0 to ˜ x 1 . Note then that ϕ ℓ,x ˜ 0 ( ˜ ) ℓ is a curve which projects via p to the same curve ℓ in M ′ since M ′ is conformally equivalent to D / g (and in particular is homeomorphic) and it does not intersect ˜ since ℓ ℓ is a simple loop. Furthermore, since the endpoint of ϕ ℓ,x ˜ 0 ( ˜ ) ℓ is ϕ ℓ,x ˜ 0 ◦ ϕ ℓ,x ˜ 0 (˜ ), by the x 0 unique path lifting property, the curve ˜ ℓ ∪ ϕ ℓ,x ˜ 0 ( ˜ ) ℓ is the unique lift obtained by going around ℓ twice in the same direction. Iterating, we obtain that the infinite path obtained by going around ℓ in the same direction is given by ∪ ∞ n =0 ϕ ( n ) ℓ,x ˜ 0 ( ˜ ) ℓ where ϕ ( n ) ℓ,x ˜ 0 is the n -fold composition of ϕ ℓ,x ˜ 0 and that the union is a disjoint union.

We can actually say more using the classification of M¨bius maps according to their trace. o Recall that M¨bius maps preserving the unit disc have the form o

$$\phi ( z ) = e ^ { i \theta } \frac { z - a } { \bar { a } z - 1 } ; \ \ | a | < 1 ; \ \ \theta \in [ 0, 2 \pi )$$

and can be classified (up to conjugation with M¨bius transforms) depending upon the behaviour o of the trace:

- · If | e iθ - | 1 2 = 4( | a | 2 -1) 2 , then ϕ z ( ) is conjugate to either z +1 or z -1 seen as maps from H to H .
- · If | e iθ - | 1 2 &gt; 4( | a | 2 -1) 2 , then ϕ z ( ) is conjugate to λz where λ &gt; 1 seen as maps from H to H .
- · If | e iθ - | 1 2 &lt; 4( | a | 2 -1) 2 , then ϕ is conjugate to a rotation of the unit disc.

Thus there is a M¨bius map Φ from o D to D or H such that ϕ = Φ -1 ◦ µ ϕ ◦ Φ where µ ϕ is either a translation by ± 1 or a scaling by λ (in case Φ maps D to H ) or a rotation (in case Φ maps D to D ).

We argue for any loop ℓ , µ ϕ ℓ,x ˜ 0 cannot be a rotation. Indeed, if a map ϕ ℓ,x ˜ 0 was conjugate to a rotation, then its iterates would be either periodic, or a dense set (on the image of a circle). Being periodic is forbidden because of path lifting property and being dense is forbidden because g is discrete. From the two remaining cases, we can complete the proof of Lemma 4.9.

Proof of Lemma 4.9. Let S be a spine, we can choose a point x 0 on p S ( ) and a lift ˜ x 0 on S . Then we see that the previous general theory applies so we can find an open path ˜ (a certain ℓ lift of the path going once around a non-contractible loop) and a map ϕ p S ,x ( ) ˜ 0 such that

$$S = \bigcup _ { n \in \mathbb { Z } } \phi _ { p ( S ), \tilde { x } _ { 0 } } ^ { ( n ) } ( \tilde { \ell } ),$$

and the map ϕ p S ,x ( ) ˜ 0 is conjugate to either a scaling or a translation. The case of a scaling clearly gives a simple path between two different points (the image of 0 and ∞ ) while the case of a translation gives a simple loop where the image of ∞ is the unique point point of ∂ D on the loop.

In the case of the torus, the map ϕ p S ( ) is simply a translation in C so the result is trivially true.

We now prove Lemma 6.5 which provides bounds on the macroscopic winding of the spines. We repeat the statement of the Lemma here for convenience.

Lemma (Restating Lemma 6.5) . Fix a compact set ˜ K ⊂ ˜ M ′ . There exist constants c, c ′ &gt; 0 such that for all ˜ v ∈ ˜ , for all K δ &lt; δ p B ( ¯ 0 ) , n ≥ 1, j ≥ 1 and u ∈ B 1 such that P ( X τ i j +1 = u ) &gt; 0,

$$\mathbb { P } \Big ( \sup _ { \mathcal { Y } \subset \tilde { X } [ \tau _ { i _ { j } }, \tau _ { i _ { j } + 1 } ] } | W ( \mathcal { Y }, \tilde { v } ) | > n \Big | X _ { \tau _ { i _ { j } + 1 } } = u \Big ) \leq c e ^ { - c ^ { \prime } n }, \\ \mathbb { P } \Big ( \sup _ { \mathcal { Y } \subset \tilde { X } [ \tau _ { i _ { j } }, \tau _ { i _ { j } + 1 } ] } | W ( \mathcal { Y }, \tilde { v } ) | > n \Big | X _ { \tau _ { i _ { j } + 1 } } \in \partial ( G ^ { \prime } ) ^ { \# \delta } \cup \mathfrak { s } \Big ) \leq c e ^ { - c ^ { \prime } n }.$$

Here the supremum is over all continuous paths obtained by erasing portions of ˜ [ X τ i j , τ i j +1 ].

Proof. We are going to borrow the notations from Section 6.2. Take a non-contractible simple loop ℓ in M ′ through v and find a continuous path ℓ # δ in G # δ which approximates this loop in the sense that it stays within distance cδ/δ 0 from ℓ (this is guaranteed to exist by the uniform crossing estimate for small enough δ depending on ℓ ). Notice that the lift starting from ˜ of a v curve which winds infinitely many times in the clockwise (resp. anti-clockwise) direction of ℓ # δ defines an infinite path ˜ ℓ + (resp. ˜ ℓ -) which converge to the boundary in the hyperbolic case (Lemma 4.9) or goes to infinity in an asymptotic direction in case of the torus. Let ˜ ℓ ′ ± denote the portion of ˜ ℓ ± from the last exit of B 1 to infinity (given an arbitrary parametrisation starting from ˜). Notice that ( ˜ v ℓ ′ -∪ ˜ ℓ ′ + ) divides the annulus ˜ M ′ \ B 1 into two simply connected domains. By compactness we can find a constant C such that the winding of ˜ ℓ ′ + around ˜ is bounded by v C , uniformly over all points ˜ v ∈ ˜ K for a suitable choice of loop ℓ . Let τ + 1 , be the first hitting time of ˜ ℓ ′ + in the interval [ τ i j , τ i j +1 ] and by induction define τ -,i to be the first hitting time of ˜ ℓ ′ -after τ + ,i and define τ + ,i the hitting time of ˜ ℓ ′ + after τ - -,i 1 . Let I + the number of τ + ,i before τ i j +1 , i.e I + = |{ | i τ + ,i ≤ τ i j +1 }| . We now use the deterministic bound

$$\sup _ { \mathcal { Y } \subset \tilde { X } [ \tau _ { i _ { j } }, \tau _ { i _ { j } + 1 } ] } | W ( \mathcal { Y }, \tilde { v } ) | \leq ( C + 2 \pi ) ( I _ { + } + 1 ).$$

and observe that conditionally on either { X τ i j +1 = u } or { X τ i j +1 ∈ ∂ G ( ′ ) # δ ∪ s } , I + has an exponential tail. The proof of this fact is essentially the same as the second item of Lemma 4.8 in [4]. Indeed we need to show that once the random walk intersects ˜ ℓ ′ + outside B 1 , there is a positive probability (uniform in δ ) for the walk to create a non-contractible loop without intersecting ˜ ℓ ′-(as in the proof of Lemma 5.2). This is intuitively clear, but a complete proof of this needs an input from Riemannian geometry. The issue at hand is that too close to the boundary, the uniform crossing ceases to hold uniformly in δ . We control the winding of this portion of this walk as follows.

Consider a compact set A ⊂ M ′ containing ℓ which forms an approximation of ℓ , in the sense that topologically K is an annulus with ℓ being non-contractible in K . Recall from the proof of Lemma 5.2 that the simple random walk has a uniform positive probability to create a noncontractible loop inside K by winding around exactly once and we stop the simple random walk if this happens. But in this process, we can assume without loss of generality on ℓ that the lift of the walk stays inside at most four consecutive copies of K (where these four copies are mapped to each other by M¨bius transforms). So if the walk is on o ˜ ℓ + at a copy of ˜ which is more than A four copies away from ˜ and v ˜ ℓ -, the lift of that walk cannot intersect ℓ -in this process. On the other hand, applying four copies of the corresponding M¨bius map to an arbitrary ˜ o v ∈ ˜ K yields by compactness of ˜ K a slightly bigger compact ˜ K ′ . Applying the uniform crossing property for the given choice of δ to ˜ K ′ , we now simply apply the argument of Lemma 4.8 in [4].

## References

- [1] A. Ahn, M. Russkikh, and R. Van Peski. Lozenge tilings and the Gaussian free field on a cylinder. arXiv preprint arXiv:2105.00551 , 2021.
- [2] J. Aru. The geometry of the Gaussian free field combined with SLE processes and the KPZ relation . PhD thesis, Ecole normale sup´rieure de Lyon, 2015. e
- [3] N. Berestycki and L. Haunschmid-Sibitz. Near-critical dimers and massive SLE. arXiv preprint arXiv:2203.15717 , 2022.
- [4] N. Berestycki, B. Laslier, and G. Ray. Dimers and imaginary geometry. The Annals of Probability , 48(1):1-52, 2020.
- [5] N. Berestycki, B. Laslier, and G. Ray. Dimers on Riemann surfaces II: conformal invariance and scaling limits. Preprint, 2022.
- [6] N. Berestycki and M. Liu. Piecewise Temperleyan dimers and a multiple SLE . 8 https://arxiv.org/abs/2301.08513 , 2022.
- [7] N. Berestycki and E. Powell. Gaussian free field, Gaussian multiplicative chaos and Liouville quantum gravity. Book in preparation, available on authors' webpage , 2015.
- [8] R. Bott and L. W. Tu. Differential forms in algebraic topology , volume 82. Springer Science &amp; Business Media, 2013.
- [9] C. Boutillier and B. De Tili`re. e Loop statistics in the toroidal honeycomb dimer model. The Annals of Probability , 37(5):1747-1777, 2009.
- [10] D. Cimasoni. Discrete Dirac operators on Riemann surfaces and Kasteleyn matrices. J. Eur. Math. Soc. , 14:1209-1244, 2009.
- [11] D. Cimasoni and N. Reshetikhin. Dimers on surface graphs and spin structures. i. Communications in Mathematical Physics , 275(1):187-208, 2007.
- [12] D. Cimasoni and N. Reshetikhin. Dimers on surface graphs and spin structures. ii. Communications in mathematical physics , 281(2):445-468, 2008.
- [13] M. Ciucu and C. Krattenthaler. A factorization theorem for lozenge tilings of a hexagon with triangular holes. Trans. Amer. Math. Soc. , 369(5):3655-3672, 2017.
- [14] H. Cohn, R. Kenyon, and J. Propp. A variational principle for domino tilings. J. Am. Math. Soc. , 14(2):297-346, 2001.
- [15] R. Costa-Santos and B. M. McCoy. Dimers and the critical ising model on lattices of genus &gt; 1. Nuclear Physics B , 623(3):439-473, 2002.
- [16] F. David, A. Kupiainen, R. Rhodes, and V. Vargas. Liouville quantum gravity on the Riemann sphere. Communications in Mathematical Physics , 342(3):869-907, 2016.
- [17] P. Di Francesco, H. Saleur, and J.-B. Zuber. Relations between the Coulomb gas picture and conformal invariance of two-dimensional critical models. J. Stat. Phys. , 49(1-2):57-79, 1987.
- [18] S. Donaldson. Riemann surfaces . Oxford University Press, 2011.
- [19] J. Dub´ edat. Dimers and families of Cauchy Riemann operators I. J. Am. Math. Soc , 28(4):1063-1167, 2015.
- [20] J. Dub´ edat and R. Gheissari. Asymptotics of height change on toroidal Temperleyan dimer models. Journal of Statistical Physics , 159(1):75-100, 2015.
- [21] C. Guillarmou, A. Kupiainen, R. Rhodes, and V. Vargas. Segal's axioms and bootstrap for Liouville theory. arXiv preprint arXiv:2112.14859 , 2021.
- [22] E. P. Hsu. Stochastic analysis on manifolds , volume 38. American Mathematical Soc., 2002.

- [23] J. Jost. Compact Riemann surfaces . Springer, 2006.
- [24] A. Kassel and R. Kenyon. Random curves on surfaces induced from the Laplacian determinant. Ann. Probab. , 2012.
- [25] R. Kenyon. Conformal invariance of domino tiling. Ann. Probab. , 28(2):759-795, 2000.
- [26] R. Kenyon. Dominos and the Gaussian free field. Ann. Probab. , 29(3):1128-1137, 2001.
- [27] R. Kenyon and A. Okounkov. Limit shapes and the complex burgers equation. Acta Math. , 199(2):263-302, 2007.
- [28] R. W. Kenyon, J. G. Propp, and D. B. Wilson. Trees and matchings. Electron. J. Combin. , 7, 2000.
- [29] H. Lacoin, R. Rhodes, and V. Vargas. Semiclassical limit of Liouville field theory. J. Funct. Anal. , 273(3):875-916, 2017.
- [30] G. F. Lawler, O. Schramm, and W. Werner. Conformal invariance of planar loop-erased random walks and uniform spanning trees. Ann. Probab. , 32(1):939-995, 2004.
- [31] C. Mercat. Discrete Riemann surfaces. In Handbook of Teichm¨ller theory. Vol. I u , volume 11 of IRMA Lect. Math. Theor. Phys. , pages 541-575. Eur. Math. Soc., Z¨rich, 2007. u
- [32] J. Miller and S. Sheffield. Imaginary geometry I: interacting SLEs. Probab. Theor. Rel Fields , (164):553-705, 2016.
- [33] J. Miller and S. Sheffield. Imaginary geometry IV: interior rays, whole-plane reversibility, and space-filling trees. Probab. Theor. Rel. Fields , (169):729-869, 2017.
- [34] L. Petrov. Asymptotics of uniformly random lozenge tilings of polygons. Gaussian free field. The Annals of Probability , 43(1):1-43, 2015.
- [35] M. Russkikh. Dimers in piecewise Temperleyan domains. Comm. Math. Phys. , 359(1):189222, 2018.
- [36] O. Schramm. Scaling limits of loop-erased random walks and uniform spanning trees. Israel J. Math. , 118:221-288, 2000.
- [37] S. Sheffield. Gaussian free fields for mathematicians. Probab. Th. Rel. Fields , 139(3-4):521541, 2007.
- [38] W. Sun. Toroidal dimer model and Temperley's bijection. arXiv:1603.00690 , 2016.
Nathanaël Berestycki Benoit Laslier Gourab Ray
Dimers on Riemann surfaces I: Temperleyan forests
This is the first article in a series of two papers in which we study the
Temperleyan dimer model on an arbitrary bounded Riemann surface of finite
topolgical type. The end goal of both papers is to prove the convergence of
height fluctuations to a universal and conformally invariant scaling limit. In
this part we show that the dimer model on the Temperleyan superposition of a
graph embedded on the surface and its dual is well posed, provided that we
remove an appropriate number of punctures. We further show that the resulting
dimer configuration is in bijection with an object which we call Temperleyan
forest, whose law is characterised in terms of a certain topological condition.
Finally we discuss the relation between height differences and Temperleyan
forest, and give a criterion guaranteeing the convergence of the height
fluctuations in terms of the Temperleyan forest.