The Full Scaling Limit of Two-Dimensional Critical Percolation

Federico Camia Research partially supported by a Marie Curie Intra-European Fellowship under contract MEIF-CT-2003-500740. E-mail: camia@eurandom.tue.nl EURANDOM, P.O. Box 513, 5600 MB Eindhoven, The Netherlands

Charles M. Newman Research partially supported by the U.S. NSF under grant DMS-01-04278. E-mail: newman@courant.nyu.edu Courant Inst. of Mathematical Sciences, New York University, New York, NY 10012, USA

, which itself relies on Smirnov's result that crossing probabilities converge to Cardy's formula.

Keywords: continuum scaling limit, percolation, SLE, critical behavior, triangular lattice, conformal invariance.

AMS 2000 Subject Classification: 82B27, 60K35, 82B43, 60D05, 30C35.

1 Introduction and Motivation

In the theory of critical phenomena it is usually assumed that a physical system near a continuous phase transition is characterized by a single length scale (the “correlation length”) in terms of which all other lengths should be measured. When combined with the experimental observation that the correlation length diverges at the phase transition, this simple but strong assumption, known as the scaling hypothesis, leads to the belief that at criticality the system has no characteristic length, and is therefore invariant under scale transformations. This suggests that all thermodynamic functions at criticality are homogeneous functions, and predicts the appearance of power laws. It also means that it should be possible to rescale a critical system appropriately and obtain a continuum model (the “continuum scaling limit”) which may have more symmetries and be easier to study than the original discrete model defined on a lattice.

Indeed, thanks to the work of Polyakov [23] and others [4, 5] , it was understood by physicists since the early seventies that critical statistical mechanical models should possess continuum scaling limits with a global conformal invariance that goes beyond simple scale invariance, as long as the discrete models have “enough” rotation invariance.

This property gives important information, enabling the determination of twoand three-point functions at criticality, when they are nonvanishing. Because the conformal group is in general a finite dimensional Lie group, the resulting constraints are limited in number; however, the situation becomes particularly interesting in two dimensions, since there every analytic function $\omega =f\left(z\right)$  defines a conformal transformation, at least at points where $f^{\prime }\left(z\right)\ne 0$  . As a consequence, the conformal group in two dimensions is infinite-dimensional.

After this observation was made, a large number of critical problems in two dimensions were analyzed using conformal methods, which were applied, among others, to Ising and Potts models, Brownian motion, Self-Avoiding Walk (SAW), percolation, and Diffusion Limited Aggregation (DLA). The large body of knowledge and techniques that resulted, starting with the work of Belavin, Polyakov and Zamolodchikov [4, 5] in the early eighties, goes under the name of Conformal Field Theory (CFT). In two dimensions, one of the main goals of CFT and its most important application to statistical mechanics is a complete classification of all universality classes via irreducible representations of the infinite-dimensional Virasoro algebra.

Partly because of the success of CFT, work in recent years on critical phenomena seemed to slow down somewhat, probably due to the feeling that most of the leading problems had been resolved. Nonetheless, however powerful and successful it may be, CFT has some limitations and leaves various open problems. First of all, the theory deals primarily with correlation functions of local (or quasi-local) operators, and is therefore not always the best tool to investigate other quantities. Secondly, given some critical lattice model, there is no way, within the theory itself, of deciding to which CFT it corresponds. A third limitation, of a different nature, is due to the fact that the methods of CFT, although very powerful, are generally speaking not completely rigorous from a mathematical point of view.

In a somewhat surprising twist, the most recent developments in the area of two-dimensional critical phenomena have emerged in the mathematics literature and have followed a newdirection, which has provided new tools and a way of coping with at least some of the limitations of CFT. The new approach may even provide a reinterpretation of CFT, and seems to be complementary to the traditional one in the sense that questions that are difficult to pose and/or answer within CFT are easy and natural in this new approach and vice versa.

The main tool of this radically new approach is the Stochastic Loewner Evolution ( $SLE$  ), or Schramm Loewner Evolution, as it is also known, introduced by Schramm [28] .

The new approach, which is probabilistic in nature, focuses directly on non-local structures that characterize a given system, such as cluster boundaries in Ising, Potts and percolation models, or loops in the $O\left(n\right)$  model. At criticality, these non-local objects become, in the continuum limit, random curves whose distributions can be uniquely identified thanks to their conformal invariance and a certain “Markovian” property. There is a one-parameter family of $SLE$  s, indexed by a positive real number $\kappa$  , and they appear to be the only possible candidates for the scaling limits of interfaces of two-dimensional critical systems that are believed to be conformally invariant.

In particular, substantial progress has been made in recent years, thanks to $SLE$  , in understanding the fractal and conformally invariant nature of (the scaling limit of ) large percolation clusters, which has attracted much attention and is of interest both for intrinsic reasons, given the many applications of percolation, and as a paradigm for the behavior of other systems. The work of Schramm [28] and Smirnov [30] has identified the scaling limit of a certain percolation interface with $SL{E}_6$  , providing, along with the work of Lawler-Schramm-Werner [19, 20] and Smirnov-Werner [34] , a confirmation of many results in the physics literature, as well as some new results.

However, $SL{E}_6$  describes a single interface, which can be obtained by imposing special boundary conditions, and is not in itself sufficient to immediately describe the full scaling limit of the system. In fact, not only the nature and properties, but the very existence of the full scaling limit remained an open question. This is true of all models, such as Ising and Potts models, that are represented in terms of clusters. Werner [36] considered this problem in the context of $SL{E}_{\kappa }$  for values of $\kappa$  between $8/3$  and $4$  . For percolation (corresponding to $\kappa =6$  ), the same problem was addressed in [7] , where $SL{E}_6$  was used to construct a random process of continuous loops in the plane, which was identified with the full scaling limit of critical two-dimensional percolation, but without detailed proofs.

In this paper, we complete the analysis of [7] , making rigorous the connection between the construction given there and the full scaling limit of percolation, and we prove some properties of the full scaling limit, the Continuum Nonsimple Loop process, including (one version of ) conformal invariance. We do this in two parts. First, we give proofs in which we assume the validity of what we will call statement (S) (see Section  5 ), which is a specific version of the results of Schramm and of Smirnov [28, 30, 31, 32, 33] concerning convergence of percolation exploration paths to $SL{E}_6$  (see the discussion towards the end of Section  4.1 ). Since no detailed proof of statement (S) (or indeed, any version of convergence to $SL{E}_6$  ) has been available, in Appendix  A we give a proof based only on that part of Smirnov's results about the convergence of crossing probabilities to Cardy's formula [30] (see Theorem  4 in Appendix  A ). We note that statement (S) is restricted toJordan domains while no such restriction is indicated in [30, 31] .

The rest of the paper is organized as follows. In Section  2 , we give necessary definitions and introduce $SL{E}_6$  . Section  3 is devoted to the construction of the Continuum Nonsimple Loop process. In Section  4 , we introduce the discrete model and a discrete construction analogous to the continuum one presented in Section  3 . Most of the main results of this paper are stated in Section  5 , while Section  6 contains the proofs of those results, which use (S). The long Appendix  A contains the proof of statement (S) (it is a consequence of Corollary  A.1 there) and the short Appendix  B contains convergence results for sequences of conformal maps which are used throughout the paper.

We remark that although our proof in Appendix  A of convergence of exploration paths to $SL{E}_6$  roughly follows Smirnov's outline [30, 31] , based on his proof [30, 31] of convergence of crossing probabilities to Cardy's formula and on the Markovian properties of hulls and tips, there are at least two technically significant modifications. The first is that we use a different sequence of stopping times to obtain a Markov chain approximation to $SL{E}_6$  , which results in a different geometry for the approximation (see Remark  A.2 ). The second is that the control of “close encounters” by the exploration path to the domain boundary is not handled by general results for “three-arms” events at the boundary of a half-plane, but rather by an argument based on continuity of crossing probabilities with respect to domain boundaries (see Lemmas  A.2 ,  A.3 ,  A.4 and  A.5 ). Moreover, we cannot use directly Smirnov's result on convergence of crossing probabilities (see Theorem  4 ), but need an extended version which is given in Theorem  6 of Appendix  A .

We conclude by noting that the convergence results of Appendix  A are sufficient not only for our purposes of obtaining the full scaling limit, but also for obtaining the critical exponents (see [34] ).

2 Preliminary Definitions

We will find it convenient to identify the real plane ${\mathbb{R}}^2$  and the complex plane $\mathbb{C}$  . We will also refer to the Riemann sphere $\mathbb{C}\cup \infty$  and the open upper half-plane $\mathbb{H}=\{x+iy:y>0\}$  (and its closure $\overline{\mathbb{H}}$  ), where chordal $SLE$  will be defined (see Section  2.3 ). $\mathbb{D}$  will denote the open unit disc $\mathbb{D}=\{z\in \mathbb{C}:|z|<1\}$  .

A domain $D$  of the complex plane $\mathbb{C}$  is a nonempty, connected, open subset of $\mathbb{C}$  ; a simply connected domain $D$  is said to be a Jordan domain if its (topological) boundary $\partial D$  is a Jordan curve (i.e., a simple continuous loop).

We will make repeated use of Riemann's mapping theorem, which states that if $D$  is any simply connected domain other than the entire plane $\mathbb{C}$  and $z_0\in D$  , then there is a unique conformal map $f$  of $D$  onto $\mathbb{D}$  such that $f\left(z_0\right)=0$  and $f^{\prime }\left(z_0\right)>0$  .

2.1 Compactification of ${\mathbb{R}}^2$ 

When taking the scaling limit $\delta \to 0$  one can focus on fixed finite regions, $\Lambda \subset {\mathbb{R}}^2$  , or consider the whole ${\mathbb{R}}^2$  at once. The second option avoids dealing with boundary conditions, but requires an appropriate choice of metric.

A convenient way of dealing with the whole ${\mathbb{R}}^2$  is to replace the Euclidean metric with a distance function $\Delta (\cdot ,\cdot )$  defined on ${\mathbb{R}}^2\times {\mathbb{R}}^2$  by

$$\begin{array}{c}\Delta (u,v)={inf}_{\phi }\int (1+|\phi {|}^2{)}^{-1}ds,\end{array}$$

(1)

where the infimum is over all smooth curves $\phi \left(s\right)$  joining $u$  with $v$  , parametrized by arclength $s$  , and where $|\cdot |$  denotes the Euclidean norm. This metric is equivalent to the Euclidean metric in bounded regions, but it has the advantage of making ${\mathbb{R}}^2$  precompact.

Adding a single point at infinity yields the compact space ${\dot{\mathbb{R}}}^2$  which is isometric, via stereographic projection, to the two-dimensional sphere.

2.2 The Space of Curves

In dealing with the scaling limit we use the approach of Aizenman-Burchard [2] . Denote by ${\mathcal{S}}_R$  the complete separable metric space of continuous curves in the closure ${\overline{\mathbb{D}}}_R$  of the disc ${\mathbb{D}}_R$  of radius $R$  with the metric ( 2 ) defined below. Curves are regarded as equivalence classes of continuous functions from the unit interval to ${\overline{\mathbb{D}}}_R$  , modulo monotonic reparametrizations. $\gamma$  will represent a particular curve and $\gamma \left(t\right)$  a parametrization of $\gamma$  ; $\mathcal{ℱ}$  will represent a set of curves (more precisely, a closed subset of ${\mathcal{S}}_R$  ). $\text{d}(\cdot ,\cdot )$  will denote the uniform metric on curves, defined by

$$\begin{array}{c}\text{d}({\gamma }_1,{\gamma }_2)\equiv inf{sup}_{t\in [0,1]}\left|{\gamma }_1\right(t)-{\gamma }_2(t\left)\right|,\end{array}$$

(2)

where the infimum is over all choices of parametrizations of ${\gamma }_1$  and ${\gamma }_2$  from the interval $[0,1]$  . The distance between two closed sets of curves is defined by the induced Hausdorff metric as follows:

$$\begin{array}{c}\text{dist}(\mathcal{ℱ},{\mathcal{ℱ}}^{\prime })\le \varepsilon \iff \forall \gamma \in \mathcal{ℱ},\exists {\gamma }^{\prime }\in {\mathcal{ℱ}}^{\prime }\text{with}\text{d}(\gamma ,{\gamma }^{\prime })\le \varepsilon ,\text{and vice versa.}\end{array}$$

(3)

The space ${\Omega }_R$  of closed subsets of ${\mathcal{S}}_R$  (i.e., collections of curves in ${\overline{\mathbb{D}}}_R$  ) with the metric ( 3 ) is also a complete separable metric space. We denote by ${\mathcal{ℬ}}_R$  its Borel $\sigma$  -algebra.

For each fixed $\delta >0$  , the random curves that we consider are polygonal paths on the edges of the hexagonal lattice $\delta \mathcal{\mathscr{H}}$  , dual to the triangular lattice $\delta \mathcal{T}$  . A superscript $\delta$  is added to indicate that the curves correspond to a model with a “short distance cutoff ” of magnitude $\delta$  .

We will also consider the complete separable metric space $\mathcal{S}$  of continuous curves in ${\dot{\mathbb{R}}}^2$  with the distance

$$\begin{array}{c}\text{D}({\gamma }_1,{\gamma }_2)\equiv inf{sup}_{t\in [0,1]}\Delta \left({\gamma }_1\right(t),{\gamma }_2(t\left)\right),\end{array}$$

(4)

where the infimum is again over all choices of parametrizations of ${\gamma }_1$  and ${\gamma }_2$  from the interval $[0,1]$  . The distance between two closed sets of curves is again defined by the induced Hausdorff metric as follows:

$$\begin{array}{c}\text{Dist}(\mathcal{ℱ},{\mathcal{ℱ}}^{\prime })\le \varepsilon \iff \forall \gamma \in \mathcal{ℱ},\exists {\gamma }^{\prime }\in {\mathcal{ℱ}}^{\prime }\text{with}\text{D}(\gamma ,{\gamma }^{\prime })\le \varepsilon ,\text{and vice versa.}\end{array}$$

(5)

The space $\Omega$  of closed sets of $\mathcal{S}$  (i.e., collections of curves in ${\dot{\mathbb{R}}}^2$  ) with the metric ( 5 ) is also a complete separable metric space. We denote by $\mathcal{ℬ}$  its Borel $\sigma$  -algebra.

When we talk about convergence in distribution of random curves, we always mean with respect to the uniform metric ( 2 ), while when we deal with closed collections of curves, we always refer to the metric ( 3 ) or ( 5 ).

Remark 2.1. In this paper, the space $\Omega$  of closed sets of $\mathcal{S}$  is generally used for collections of exploration paths and cluster boundary loops and their scaling limits, $SL{E}_6$  paths and continuum nonsimple loops. There is one place however, in the statements and proofs of Lemmas  A.2 ,  A.4 and  A.5 , where we also apply $\Omega$  in essentially the original setting of Aizenman and Burchard [1, 2] , i.e., for collections of blue and yellow simple $\mathcal{T}$  -paths (see Section  4 for precise definitions) and their scaling limits. The slight modification needed to keep track of both the paths and their colors is easily managed.

2.3 Chordal $SLE$  in the Upper Half-Plane

The Stochastic Loewner Evolution ( $SLE$  ) was introduced by Schramm [28] as a tool for studying the scaling limit of two-dimensional discrete (defined on a lattice) probabilistic models whose scaling limits are expected to be conformally invariant. In this section we define the chordal version of $SLE$  ; for more on the subject, the interested reader can consult the original paper [28] as well as the fine reviews by Lawler [17] , Kager and Nienhuis [14] , and Werner [37] , and Lawler's book [18] .

Let $\mathbb{H}$  denote the upper half-plane. For a given continuous real function $U_t$  with $U_0=0$  , define, for each $z\in \overline{\mathbb{H}}$  , the function $g_t\left(z\right)$  as the solution to the ODE

$$\begin{array}{c}{\partial }_t{g}_t\left(z\right)=\frac{2}{g_t\left(z\right)-U_t},\end{array}$$

(6)

with $g_0\left(z\right)=z$  . This is well defined as long as $g_t\left(z\right)-U_t\ne 0$  , i.e., for all $t<T\left(z\right)$  , where

$$\begin{array}{c}T\left(z\right)\equiv sup\{t\ge 0:{min}_{s\in [0,t]}|g_s\left(z\right)-U_s|>0\}.\end{array}$$

(7)

Let $K_t\equiv \{z\in \overline{\mathbb{H}}:T(z)\le t\}$  and let ${\mathbb{H}}_t$  be the unbounded component of $\mathbb{H}\backslash K_t$  ; it can be shown that $K_t$  is bounded and that $g_t$  is a conformal map from ${\mathbb{H}}_t$  onto $\mathbb{H}$  . For each $t$  , it is possible to write $g_t\left(z\right)$  as

$$\begin{array}{c}{g}_t\left(z\right)=z+\frac{2t}{z}+o\left(\frac{1}{z}\right),\end{array}$$

(8)

when $z\to \infty$  . The family $(K_t,t\ge 0)$  is called the Loewner chain associated to the driving function $(U_t,t\ge 0)$  .

Definition 2.1. Chordal $SL{E}_{\kappa }$  is the Loewner chain $(K_t,t\ge 0)$  that is obtained when the driving function $U_t=\sqrt{\kappa }{B}_t$  is $\sqrt{\kappa }$  times a standard real-valued Brownian motion $(B_t,t\ge 0)$  with $B_0=0$  .

For all $\kappa \ge 0$  , chordal $SL{E}_{\kappa }$  is almost surely generated by a continuous random curve $\gamma$  in the sense that, for all $t\ge 0$  , ${\mathbb{H}}_t\equiv \mathbb{H}\backslash K_t$  is the unbounded connected component of $\mathbb{H}\backslash \gamma [0,t]$  ; $\gamma$  is called the trace of chordal $SL{E}_{\kappa }$  .

2.4 Chordal $SLE$  in an Arbitrary Simply Connected Domain

Let $D\subset \mathbb{C}$  ( $D\ne \mathbb{C}$  ) be a simply connected domain whose boundary is a continuous curve. By Riemann's mapping theorem, there are (many) conformal maps from the upper half-plane $\mathbb{H}$  onto $D$  . In particular, given two distinct points $a,b\in \partial D$  (or more accurately, two distinct prime ends), there exists a conformal map $f$  from $\mathbb{H}$  onto $D$  such that $f\left(0\right)=a$  and $f(\infty )\equiv {lim}_{\left|z\right|\to \infty }f\left(z\right)=b$  . In fact, the choice of the points $a$  and $b$  on the boundary of $D$  only characterizes $f(\cdot )$  up to a multiplicative factor, since $f(\lambda \cdot )$  would also do.

Suppose that $(K_t,t\ge 0)$  is a chordal $SL{E}_{\kappa }$  in $\mathbb{H}$  as defined above; we define chordal $SL{E}_{\kappa }$  $({\tilde{K}}_t,t\ge 0)$  in $D$  from $a$  to $b$  as the image of the Loewner chain $(K_t,t\ge 0)$  under $f$  . It is possible to show, using scaling properties of $SL{E}_{\kappa }$  , that the law of $({\tilde{K}}_t,t\ge 0)$  is unchanged, up to a linear time-change, if we replace $f(\cdot )$  by $f(\lambda \cdot )$  . This makes it natural to consider $({\tilde{K}}_t,t\ge 0)$  as a process from $a$  to $b$  in $D$  , ignoring the role of $f$  .

We are interested in the case $\kappa =6$  , for which $(K_t,t\ge 0)$  is generated by a continuous, nonsimple, non-self-crossing curve $\gamma$  with Hausdorff dimension $7/4$  . We will denote by ${\gamma }_{D,a,b}$  the image of $\gamma$  under $f$  and call it the trace of chordal $SL{E}_6$  in $D$  from $a$  to $b$  ; ${\gamma }_{D,a,b}$  is a continuous nonsimple curve inside $D$  from $a$  to $b$  , and it can be given a parametrization ${\gamma }_{D,a,b}\left(t\right)$  such that ${\gamma }_{D,a,b}\left(0\right)=a$  and ${\gamma }_{D,a,b}\left(1\right)=b$  , so that we are in the metric framework described in Section  2.2 . It will be convenient to think of ${\gamma }_{D,a,b}$  as an oriented path, with orientation from $a$  to $b$  .

3 Construction of the Continuum Nonsimple Loops

3.1 Construction of a Single Loop

As a preview to the full construction, we explain how to construct a single loop using two $SL{E}_6$  paths inside a domain $D$  whose boundary is assumed to have a given orientation (clockwise or counterclockwise). This is done in three steps (see Figure  1 ), of which the first consists in choosing two points $a$  and $b$  on the boundary $\partial D$  of $D$  and “running” a chordal $SL{E}_6$  , $\gamma ={\gamma }_{D,a,b}$  , from $a$  to $b$  inside $D$  . As explained in Section  2.4 , we consider $\gamma$  as an oriented path, with orientation from $a$  to $b$  . The set $D\backslash {\gamma }_{D,a,b}[0,1]$  is a countable union of its connected components, which are open and simply connected. If $z$  is a deterministic point in $D$  , then with probability one, $z$  is not touched by $\gamma$  [26] and so it belongs to a unique domain in $D\backslash {\gamma }_{D,a,b}[0,1]$  that we denote $D_{a,b}\left(z\right)$  .

The elements of $D\backslash {\gamma }_{D,a,b}[0,1]$  can be conveniently thought of in terms of how a point $z$  in the interior of the component was first “trapped” at some time $t_1$  by $\gamma [0,t_1]$  , perhaps together with either ${\partial }_{a,b}D$  or ${\partial }_{b,a}D$  (the portions of the boundary $\partial D$  from $a$  to $b$  counterclockwise or clockwise respectively): (1) those components whose boundary contains a segment of ${\partial }_{b,a}D$  between two successive visits at ${\gamma }_0\left(z\right)=\gamma \left(t_0\right)$  and ${\gamma }_1\left(z\right)=\gamma \left(t_1\right)$  to ${\partial }_{b,a}D$  (where here and below $t_0<t_1$  ), (2) the analogous components with ${\partial }_{b,a}D$  replaced by the other part of the boundary ${\partial }_{a,b}D$  , (3) those components formed when ${\gamma }_0\left(z\right)=\gamma \left(t_0\right)=\gamma \left(t_1\right)={\gamma }_1\left(z\right)\in D$  with $\gamma$  winding about $z$  in a counterclockwise direction between $t_0$  and $t_1$  , and finally (4) the analogous clockwise components.

We give to the boundary of a domain of type 3 or 4 the orientation induced by how the curve $\gamma$  winds around the points inside that domain. For a domain $D^{\prime }\ni z$  of type 1 or 2 which is produced by an “excursion” $\mathcal{ℰ}$  from ${\gamma }_0\left(z\right)\in \partial D$  to ${\gamma }_1\left(z\right)\in \partial D$  , the part of the boundary that corresponds to the inner perimeter of the excursion $\mathcal{ℰ}$  (i.e., the perimeter of $\gamma$  seen from $z$  ) is oriented according to the direction of $\gamma$  , i.e., from ${\gamma }_0\left(z\right)$  to ${\gamma }_1\left(z\right)$  .

If we assume that $\partial D$  is oriented from $a$  to $b$  clockwise, then the boundaries of domains of type 2 have a well defined orientation, while the boundaries of domains of type 1 do not, since they are composed of two parts which are both oriented from the beginning to the end of the excursion that produced the domain.

Now, let $D^{\prime }$  be a domain of type 1 and let $A$  and $B$  be respectively the starting and ending point of the excursion that generated $D^{\prime }$  . The second step to construct a loop is to run a chordal $SL{E}_6$  , ${\gamma }^{\prime }={\gamma }_{D^{\prime },B,A}$  , inside $D^{\prime }$  from $B$  to $A$  ; the third and final step consists in pasting together $\mathcal{ℰ}$  and ${\gamma }^{\prime }$  .

Running ${\gamma }^{\prime }$  inside $D^{\prime }$  from $B$  to $A$  partitions $D^{\prime }\backslash {\gamma }^{\prime }$  into new domains. Notice that if we assign an orientation to the boundaries of these domains according to the same rules used above, all of those boundaries have a well defined orientation, so that the construction of loops just presented can be iterated inside each one of these domains (as well as inside each of the domains of type 2, 3 and 4 generated by ${\gamma }_{D,a,b}$  in the first step). This will be done in the next section.

3.2 The Full Construction Inside The Unit Disc

In this section we define the Continuum Nonsimple Loop process inside the unit disc $\mathbb{D}={\mathbb{D}}_1$  via an inductive procedure. Later, in order to define the continuum nonsimple loops in the whole plane, the unit disc will be replaced by a growing sequence of large discs, ${\mathbb{D}}_R$  , with $R\to \infty$  (see Theorem  2 ). The basic ingredient in the algorithmic construction, given in the previous section, consists of a chordal $SL{E}_6$  path ${\gamma }_{D,a,b}$  between two points $a$  and $b$  of the boundary $\partial D$  of a given simply connected domain $D\subset \mathbb{C}$  .

We will organize the inductive procedure in steps, each one corresponding to one $SL{E}_6$  inside a certain domain generated by the previous steps. To do that, we need to order the domains present at the end of each step, so as to choose the one to use in the next step.

For this purpose, we introduce a deterministic countable set of points $\mathcal{P}$  that are dense in $\mathbb{C}$  and are endowed with a deterministic order (here and below by deterministic we mean that they are assigned before the beginning of the construction and are independent of the $SL{E}_6$  's).

The first step consists of an $SL{E}_6$  path, ${\gamma }_1={\gamma }_{\mathbb{D},-i,i}$  , inside $\mathbb{D}$  from $-i$  to $i$  , which produces many domains that are the connected components of the set $\mathbb{D}\backslash {\gamma }_1[0,1]$  . These domains can be priority-ordered according to the maximal $x$  or $y$  coordinate distances between points on their boundaries and using the rank of the points in $\mathcal{P}$  (contained in the domains) to break ties, as follows. For a domain $D$  , let ${\text{d}}_m\left(D\right)$  be the maximal $x$  or $y$  -distance between points on its boundary, whichever is greater. Domains with larger ${\text{d}}_m$  have higher priority, and if two domains have the same ${\text{d}}_m$  , the one containing the highest ranking point of $\mathcal{P}$  from those two domains has higher priority. The priority order of domains of course changes as the construction proceeds and new domains are formed.

The second step of the construction consists of an $SL{E}_6$  path, ${\gamma }_2$  , that is produced in the domain with highest priority (after the first step). Since all the domains that are produced in the construction are Jordan domains, as explained in the discussion following Corollary  5.1 , for all steps we can use the definition of chordal $SLE$  given in Section  2.4 .

As a result of the construction, the $SL{E}_6$  paths are naturally ordered: $\{{\gamma }_j{\}}_{j\in \mathbb{N}}$  . It will be shown (see especially the proof of Theorem  1 below) that every domain that is formed during the construction is eventually used (this is in fact one important requirement in deciding how to order the domains and therefore how to organize the construction).

So far we have not explained how to choose the starting and ending points of the $SL{E}_6$  paths on the boundaries of the domains. In order to do this, we give an orientation to the boundaries of the domains produced by the construction according to the rules explained in Section  3.1 . We call monochromatic a boundary which gets, as a consequence of those rules, a well defined (clockwise or counterclockwise) orientation; the choice of this term will be clarified when we discuss the lattice version of the loop construction below.

We will generally take our initial domain ${\mathbb{D}}_1$  (or ${\mathbb{D}}_R$  ) to have a monochromatic boundary (either clockwise or counterclockwise orientation).

It is easy to see by induction that the boundaries that are not monochromatic are composed of two “pieces” joined at two special points (call them A and B, as in theexample of Section  3.1 ), such that one piece is a portion of the boundary of a previous domain, and the other is the inner perimeter of an excursion (see again Section  3.1 ). Both pieces are oriented in the same direction, say from A to B (see Figure  1 ).

For a domain whose boundary is not monochromatic, we make the “natural” choice of starting and ending points, corresponding to the end and beginning of the excursion that produced the domain (the points B and A respectively, in the example above). As explained in Section  3.1 , when such a domain is used with this choice of points on the boundary, a loop is produced, together with other domains, whose boundaries are all monochromatic.

For a domain whose boundary is monochromatic, and therefore has a well defined orientation, there are various procedures which would yield the “correct” distribution for the resulting Continuum Nonsimple Loop process; one possibility is as follows.

Given a domain $D$  , $a$  and $b$  are chosen so that, of all pairs $(u,v)$  of points in $\partial D$  , they maximize $\left|\text{Re}\right(u-v\left)\right|$  if $\left|\text{Re}\right(u-v\left)\right|\ge \left|\text{Im}\right(u-v\left)\right|$  , or else they maximize $\left|\text{Im}\right(u-v\left)\right|$  . If the choice is not unique, to restrict the number of pairs one looks at those pairs, among the ones already obtained, that maximize the other of $\left\{\right|\text{Re}(u-v)|,|\text{Im}(u-v)\left|\right\}$  . Notice that this leaves at most two pairs of points; if that's the case, the pair that contains the point with minimal real (and, if necessary, imaginary) part is chosen. The iterative procedure produces a loop every time a domain whose boundary is not monochromatic is used. Our basic loop process consists of the collection of all loops generated by this inductive procedure (i.e., the limiting object obtained from the construction by letting the number of steps $k\to \infty$  ), to which we add a “trivial” loop for each $z$  in $\mathbb{D}$  , so that the collection of loops is closed in the appropriate sense [2] . The Continuum Nonsimple Loop process in the whole plane is introduced in Theorem  2 , Section  5 . There, a “trivial” loop for each $z\in \mathbb{C}\cup \infty$  has to be added to make the space of loops closed.

4 Lattices and Paths

We will denote by $\mathcal{T}$  the two-dimensional triangular lattice, whose sites we think of as the elementary cells of a regular hexagonal lattice $\mathcal{\mathscr{H}}$  embedded in the plane as in Figure  2 . A sequence $({\xi }_0,\dots ,{\xi }_n)$  of sites of $\mathcal{T}$  such that ${\xi }_{i-1}$  and ${\xi }_i$  are neighbors in $\mathcal{T}$  for all $i=1,\dots ,n$  and ${\xi }_i\ne {\xi }_j$  whenever $i\ne j$  will be called a $\mathcal{T}$  -path and denoted by $\pi$  .

If the first and last sites of the path are neighbors in $\mathcal{T}$  , the path will be called a $\mathcal{T}$  -loop.

We say that a finite subset $D$  of $\mathcal{T}$  is simply connected if both $D$  and $\mathcal{T}\backslash D$  are connected (by the edges of $\mathcal{T}$  ). For a simply connected set $D$  of hexagons, we denote by $\Delta D$  its external site boundary, or s-boundary (i.e., the set of hexagons that do not belong to $D$  but are adjacent to hexagons in $D$  ), and by $\partial D$  the topological boundary of $D$  when $D$  is considered as a domain of $\mathbb{C}$  . We will call a bounded, simply connected subset $D$  of $\mathcal{T}$  a Jordan set if its s-boundary $\Delta D$  is a $\mathcal{T}$  -loop.

For a Jordan set $D\subset \mathcal{T}$  , a vertex $x\in \mathcal{\mathscr{H}}$  that belongs to $\partial D$  can be either of two types, according to whether the edge incident on $x$  that is not in $\partial D$  belongs to a hexagon in $D$  or not. We call a vertex of the second type an e-vertex (e for “external” or “exposed”).

Given a Jordan set $D$  and two e-vertices $x,y$  in $\partial D$  , we denote by ${\partial }_{x,y}D$  the portion of $\partial D$  traversed counterclockwise from $x$  to $y$  , and call it the right boundary; the remaining part of the boundary is denote by ${\partial }_{y,x}D$  and is called the left boundary. Analogously, the portion of ${\Delta }_{x,y}D$  of $\Delta D$  whose hexagons are adjacent to ${\partial }_{x,y}D$  is called the right s-boundary and the remaining part the left s-boundary.

A percolation configuration $\sigma =\left\{\sigma \right(\xi ){\}}_{\xi \in \mathcal{T}}\in \{-1,+1{\}}^{\mathcal{T}}$  on $\mathcal{T}$  is an assignment of $-1$  (equivalently, yellow) or $+1$  (blue) to each site of $\mathcal{T}$  . For a domain $D$  of the plane, the restriction to the subset $D\cap \mathcal{T}$  of $\mathcal{T}$  of the percolation configuration $\sigma$  is denoted by ${\sigma }_D$  .

On the space of configurations $\Sigma =\{-1,+1{\}}^{\mathcal{T}}$  , we consider the usual product topology and denote by $\mathbb{P}$  the uniform measure, corresponding to Bernoulli percolation with equal density of yellow (minus) and blue (plus) hexagons, which is critical percolation in the case of the triangular lattice.

A (percolation) cluster is a maximal, connected, monochromatic subset of $\mathcal{T}$  ; we will distinguish between blue (plus) and yellow (minus) clusters. The boundary of a cluster $D$  is the set of edges of $\mathcal{\mathscr{H}}$  that surround the cluster (i.e., its Peierls contour); it coincides with the topological boundary of $D$  considered as a domain of $\mathbb{C}$  . The set of all boundaries is a collection of “nested” simple loops along the edges of $\mathcal{\mathscr{H}}$  .

Given a percolation configuration $\sigma$  , we associate an arrow to each edge of $\mathcal{\mathscr{H}}$  belonging to the boundary of a cluster in such a way that the hexagon to the right of the edge with respect to the direction of the arrow is blue (plus). The set of all boundaries then becomes a collection of nested, oriented, simple loops. A boundary path (or b-path) $\gamma$  is a sequence $(e_0,\dots ,e_n)$  of distinct edges of $\mathcal{\mathscr{H}}$  belonging to the boundary of a cluster and such that $e_{i-1}$  and $e_i$  meet at a vertex of $\mathcal{\mathscr{H}}$  for all $i=1,\dots ,n$  . To each b-path, we can associate a direction according to the direction of the edges in the path. Given a b-path $\gamma$  , we denote by ${\Gamma }_B\left(\gamma \right)$  (respectively, ${\Gamma }_Y\left(\gamma \right)$  ) the set of blue (resp., yellow) hexagons (i.e., sites of $\mathcal{T}$  ) adjacent to $\gamma$  ; we also let $\Gamma \left(\gamma \right)\equiv {\Gamma }_B\left(\gamma \right)\cup {\Gamma }_Y\left(\gamma \right)$  .

4.1 The Percolation Exploration Process and Path

For a Jordan set $D\subset \mathcal{T}$  and two e-vertices $x,y$  in $\partial D$  , imagine coloring blue all the hexagons in ${\Delta }_{x,y}D$  and yellow all those in ${\Delta }_{y,x}D$  . Then, for any percolation configuration ${\sigma }_D$  inside $D$  , there is a unique b-path $\gamma$  from $x$  to $y$  which separates the blue cluster adjacent to ${\Delta }_{x,y}D$  from the yellow cluster adjacent to ${\Delta }_{y,x}D$  . We call $\gamma ={\gamma }_{D,x,y}\left({\sigma }_D\right)$  a percolation exploration path (see Figure  3 ).

An exploration path $\gamma$  can be decomposed into left excursions $\mathcal{ℰ}$  , i.e., maximal b-subpaths of $\gamma$  that do not use edges of the left boundary ${\partial }_{y,x}D$  . Successive left excursions are separated by portions of $\gamma$  that contain only edges of the left boundary ${\partial }_{y,x}D$  .

Analogously, $\gamma$  can be decomposed into right excursions, i.e., maximal b-subpaths of $\gamma$  that do not use edges of the right boundary ${\partial }_{x,y}D$  . Successive right excursions are separated by portions of $\gamma$  that contain only edges of the right boundary ${\partial }_{x,y}D$  .

Notice that the exploration path $\gamma ={\gamma }_{D,x,y}\left({\sigma }_D\right)$  only depends on the percolation configuration ${\sigma }_D$  inside $D$  and the positions of the e-vertices $x$  and $y$  ; in particular, it does not depend on the color of the hexagons in $\Delta D$  , since it is defined by imposing fictitious $\pm$  boundary conditions on $D$  . To see this more clearly, we next show how to construct the percolation exploration path dynamically, via the percolation exploration process defined below.

Given a Jordan set $D\subset \mathcal{T}$  and two e-vertices $x,y$  in $\partial D$  , assign to ${\partial }_{x,y}D$  a counterclockwise orientation (i.e., from $x$  to $y$  ) and to ${\partial }_{y,x}D$  a clockwise orientation. Call $e_x$  the edge incident on $x$  that does not belong to $\partial D$  and orient it in the direction of $x$  ; this is the “starting edge” of an exploration procedure that will produce an oriented path inside $D$  along the edges of $\mathcal{\mathscr{H}}$  , together with two nonsimple monochromatic paths on $\mathcal{T}$  . From $e_x$  , the process moves along the edges of hexagons in $D$  according to the rules below. At each step there are two possible edges (left or right edge with respect to the current direction of exploration) to choose from, both belonging to the same hexagon $\xi$  contained in $D$  or $\Delta D$  .

• ∙ If $\xi$  belongs to $D$  and has not been previously “explored,” its color is determined by flipping a fair coin and then the edge to the left (with respect to the direction in which the exploration is moving) is chosen if $\xi$  is blue (plus), or the edge to the right is chosen if $\xi$  is yellow (minus).
• ∙ If $\xi$  belongs to $D$  and has been previously explored, the color already assigned to it is used to choose an edge according to the rule above.
• ∙ If $\xi$  belongs to the right external boundary ${\Delta }_{x,y}D$  , the left edge is chosen.
• ∙ If $\xi$  belongs to the left external boundary ${\Delta }_{y,x}D$  , the right edge is chosen.
• ∙ The exploration process stops when it reaches $b$  .

We can assign an arrow to each edge in the path in such a way that the hexagon to the right of the edge with respect to the arrow is blue; for edges in $\partial D$  , we assign the arrows according to the direction assigned to the boundary. In this way, we get an oriented path,whose shape and orientation depend solely on the color of the hexagons explored during the construction of the path.

When we present the discrete construction, we will encounter Jordan sets $D$  with two e-vertices $x,y\in \partial D$  assigned in some way to be discussed later. Such domains will have either monochromatic (plus or minus) boundaries or $\pm$  boundary conditions, corresponding to having both ${\Delta }_{x,y}D$  and ${\Delta }_{y,x}D$  monochromatic, but of different colors.

As explained, the exploration path ${\gamma }_{D,x,y}$  does not depend on the color of $\Delta D$  , but the interpretation of ${\gamma }_{D,x,y}$  does. For domains with $\pm$  boundary conditions, the exploration path represents the interface between the yellow cluster containing the yellow portion of the s-boundary of $D$  and the blue cluster containing its blue portion.

For domains with monochromatic blue (resp., yellow) boundary conditions, the exploration path represents portions of the boundaries of yellow (resp., blue) clusters touching ${\partial }_{y,x}D$  and adjacent to blue (resp., yellow) hexagons that are the starting point of a blue (resp., yellow) path (possibly an empty path) that reaches ${\partial }_{x,y}D$  , pasted together using portions of ${\partial }_{y,x}D$  .

In order to study the continuum scaling limit of an exploration path, we introduce the following definitions.

Definition 4.1. Given a bounded, simply connected domain $D$  of the plane, we denote by $D^{\delta }$  the largest Jordan set of hexagons of the scaled hexagonal lattice $\delta \mathcal{\mathscr{H}}$  that is contained in $D$  , and call it the $\delta$  -approximation of $D$  .

It is clear that if $D$  is a Jordan domain, then as $\delta \to 0$  , $\partial D^{\delta }$  converges to $\partial D$  in the metric ( 2 ).

Definition 4.2. Let $D$  be a bounded domain of the plane and $D^{\delta }$  its $\delta$  -approximation. For $a,b\in \partial D$  , choose the pair $(x_a,x_b)$  of e-vertices in $\partial D^{\delta }$  closest to, respectively, $a$  and $b$  (if there are two such vertices closest to $a$  , we choose, say, the first one encountered going clockwise along $\partial D^{\delta }$  , and analogously for $b$  ). Given a percolation configuration $\sigma$  , we define the exploration path ${\gamma }_{D,a,b}^{\delta }\left(\sigma \right)\equiv {\gamma }_{D^{\delta },x_a,x_b}\left(\sigma \right)$  .

For a fixed $\delta >0$  , the measure $\mathbb{P}$  on percolation configurations $\sigma$  induces a measure ${\mu }_{D,a,b}^{\delta }$  on exploration paths ${\gamma }_{D,a,b}^{\delta }\left(\sigma \right)$  . In the continuum scaling limit, $\delta \to 0$  , one is interested in the weak convergence of ${\mu }_{D,a,b}^{\delta }$  to a measure ${\mu }_{D,a,b}$  supported on continuous curves, with respect to the uniform metric ( 2 ) on continuous curves.

One of the main tools in this paper is the result on convergence to $SL{E}_6$  announced by Smirnov [30] (see also [31] ), whose detailed proof is to appear [32] : The distribution of ${\gamma }_{D,a,b}^{\delta }$  converges, as $\delta \to 0$  , to that of the trace of chordal $SL{E}_6$  inside $D$  from $a$  to $b$  , with respect to the uniform metric ( 2 ) on continuous curves.

Actually, we will rather use a slightly stronger conclusion, given as statement (S) at the beginning of Section  5 below, a version of which, according to [34] (see p. 734 there), and [33] , will be contained in [32] . This stronger statement is that the convergence of the percolation process to $SL{E}_6$  takes place locally uniformly with respect to the shape of the domain $D$  and the positions of the starting and ending points $a$  and $b$  on its boundary $\partial D$  .

We will use this version of convergence to $SL{E}_6$  to identify the Continuum Nonsimple Loop process with the scaling limit of all critical percolation clusters. Statement (S) is a direct consequence of Corollary  A.1 , which is proved in Appendix  A . Although the convergence statements in Corollary  A.1 and in (S) are stronger than those in [30, 31] , we note that they are restricted to Jordan domains, a restriction not present in [30, 31] .

Before concluding this section, we give one more definition. Consider the exploration path $\gamma ={\gamma }_{D,x,y}^{\delta }$  and the set $\Gamma \left(\gamma \right)={\Gamma }_Y\left(\gamma \right)\cup {\Gamma }_B\left(\gamma \right)$  . The set $D^{\delta }\backslash \Gamma \left(\gamma \right)$  is the union of its connected components (in the lattice sense), which are simply connected. If the domain $D$  is large and the e-vertices $x_a,y_a\in \partial D^{\delta }$  are not too close to each other, then with high probability the exploration process inside $D^{\delta }$  will make large excursions into $D^{\delta }$  , so that $D^{\delta }\backslash \Gamma \left(\gamma \right)$  will have more than one component. Given a point $z\in \mathbb{C}$  contained in $D^{\delta }\backslash \Gamma \left(\gamma \right)$  , we will denote by $D_{a,b}^{\delta }\left(z\right)$  the domain corresponding to the unique element of $D^{\delta }\backslash \Gamma \left(\gamma \right)$  that contains $z$  (notice that for a deterministic $z\in D$  , $D_{a,b}^{\delta }\left(z\right)$  is well defined with high probability for $\delta$  small, i.e., when $z\in D^{\delta }$  and $z/\in \Gamma \left(\gamma \right)$  ).

4.2 Discrete Loop Construction

Next, we show how to construct, by twice using the exploration process described in Section  4.1 , a loop $\Lambda$  along the edges of $\mathcal{\mathscr{H}}$  corresponding to the external boundary of a monochromatic cluster contained in a large, simply connected, Jordan set $D$  with monochromatic blue (say) boundary conditions (see Figures  4 and  5 ).

Consider the exploration path $\gamma ={\gamma }_{D,x,y}$  and the sets ${\Gamma }_Y\left(\gamma \right)$  and ${\Gamma }_B\left(\gamma \right)$  (see Figure  4 ).

The set $D\backslash \left\{{\Gamma }_Y\right(\gamma )\cup {\Gamma }_B(\gamma \left)\right\}$  is the union of its connected components (in the lattice sense), which are simply connected. If the domain $D$  is large and the e-vertices $x,y\in \partial D$  are chosen not too close to each other, with large probability the exploration process inside $D$  will make large excursions into $D$  , so that $D\backslash \left\{{\Gamma }_Y\right(\gamma )\cup {\Gamma }_B(\gamma \left)\right\}$  will have many components.

There are four types of components which may be usefully thought of in terms of their external site boundaries: (1) those components whose site boundary contains both sites in ${\Gamma }_Y\left(\gamma \right)$  and ${\Delta }_{y,x}D$  , (2) the analogous components with ${\Delta }_{y,x}D$  replaced by ${\Delta }_{x,y}D$  and ${\Gamma }_Y\left(\gamma \right)$  by ${\Gamma }_B\left(\gamma \right)$  , (3) those components whose site boundary only contains sites in ${\Gamma }_Y\left(\delta \right)$  , and finally (4) the analogous components with ${\Gamma }_Y\left(\gamma \right)$  replaced by ${\Gamma }_B\left(\gamma \right)$  .

Notice that the components of type 1 are the only ones with $\pm$  boundary conditions, while all other components have monochromatic s-boundaries. For a given component $D^{\prime }$  of type 1, we can identify the two edges that separate the yellow and blue portions of its s-boundary. The vertices $x^{\prime }$  and $y^{\prime }$  of $\mathcal{\mathscr{H}}$  where those two edges intersect $\partial D^{\prime }$  are e-vertices and are chosen to be the starting and ending points of the exploration path ${\gamma }_{D^{\prime },x^{\prime },y^{\prime }}$  inside $D^{\prime }$  .

If $x^{\prime \prime },y^{\prime \prime }\in \partial D$  are respectively the ending and starting points of the left excursion $\mathcal{ℰ}$  of ${\gamma }_{D,x,y}$  that “created” $D^{\prime }$  , by pasting together $\mathcal{ℰ}$  and ${\gamma }_{D^{\prime },x^{\prime },y^{\prime }}$  with the help of the edges of $\partial D$  contained between $x^{\prime }$  and $x^{\prime \prime }$  and between $y^{\prime }$  and $y^{\prime \prime }$  , we get a loop $\Lambda$  which corresponds to the boundary of a yellow cluster adjacent to ${\partial }_{y,x}D$  (see Figure  5 ). Notice that the path ${\gamma }_{D^{\prime },x^{\prime },y^{\prime }}$  in general splits $D^{\prime }$  into various other domains, all of which have monochromatic boundary conditions.

4.3 Full Discrete Construction

We now give the algorithmic construction for discrete percolation which is the analogue of the continuum one. Each step of the construction is a single percolation exploration process; the order of successive steps is organized as in the continuum construction detailed in Section  3.2 . We start with the smallest Jordan set $D_0^{\delta }={\mathbb{D}}^{\delta }$  of hexagons that covers the unit disc $\mathbb{D}$  . We will also make use of the countable set $\mathcal{P}$  of points dense in $\mathbb{C}$  that was introduced earlier.

The first step consists of an exploration process inside $D_0^{\delta }$  . For this, we need to select two points $x$  and $y$  in $\partial D_0^{\delta }$  (which identify the starting and ending edges). We choose for $x$  the e-vertex closest to $-i$  , and for $y$  the e-vertex closest to $i$  (if there are two such vertices closest to $-i$  , we can choose, say, the one with smallest real part, and analogously for $i$  ). The first exploration produces a path ${\gamma }_1^{\delta }$  and, for $\delta$  small, many new domains of all four types. These domains are ordered according to the maximal $x$  or $y$  distance ${\text{d}}_m$  between points on their boundaries and, if necessary, with the help of points in $\mathcal{P}$  , as in the continuum case, and that order is used, at each step of the construction, to determine the next exploration process. With this choice, the exploration processes and paths are naturally ordered: ${\gamma }_1^{\delta },{\gamma }_2^{\delta },\dots$  .

Each exploration process of course requires choosing a starting and ending vertex and edge. For domains of type 1, with a $\pm$  or $\mp$  boundary condition, the choice is the natural one, explained before.

For a domain $D_k^{\delta }$  (used at the $k$  th step) of type other than 1, and therefore with a monochromatic boundary, the starting and ending edges are chosen with a procedure that mimics what is done in the continuum case. Once again, the exact procedure used to choose the pair of points is not important, as long as they are not chosen too close to each other. This is clear in the discrete case because the procedure that we are presenting is only “discovering” the cluster boundaries. In more precise terms, it is clear that one could couple the processes obtained with different rules by means of the same percolation configuration, thus obtaining exactly the same cluster boundaries.

As in the continuum case, we can choose the following procedure. (In Theorem  1 we will slightly reorganize the procedure by using a coupling to the continuum construction to guarantee that the order of exploration of domains of the discrete and continuum procedures match despite the rules for breaking ties.) Given a domain $D$  , $x$  and $y$  are chosen so that, of all pairs $(u,v)$  of points in $\partial D$  , they maximize $\left|\text{Re}\right(u-v\left)\right|$  if $\left|\text{Re}\right(u-v\left)\right|\ge \left|\text{Im}\right(u-v\left)\right|$  , or else they maximize $\left|\text{Im}\right(u-v\left)\right|$  . If the choice is not unique, to restrict the number of pairs one looks at those pairs, among the ones already obtained, that maximize the other of $\left\{\right|\text{Re}(u-v)|,|\text{Im}(u-v)\left|\right\}$  . Notice that this leaves at most two pairs of points; if that's the case, the pair that contains the point with minimal real (and, if necessary, imaginary) part is chosen.

The procedure continues iteratively, with regions that have monochromatic boundaries playing the role played in the first step by the unit disc. Every time a region with $\pm$  boundary conditions is used, a new loop, corresponding to the outer boundary contour ofa cluster, is formed by pasting together, as explained in Section  3.1 , the new exploration path and the excursion containing the region where the last exploration was carried out.

All the new regions created at a step when a loop is formed have monochromatic boundary conditions.

5 Main Results

In this section we collect our main results about the Continuum Nonsimple Loop process. Before doing that, we state a precise version, called statement (S), on convergence of exploration paths to $SL{E}_6$  that we will use in the proofs of these results, presented in Section  6 . (A proof of statement (S) is given in Appendix  A ; it is an immediate consequence of Corollary  A.1 there. The proof relies, among other things, on the result of Smirnov [30] concerning convergence of crossing probabilities to Cardy's formula [10, 11] – see Theorem  4 .) We note that (S) or (Corollary  A.1 ) is both more general and more special than the convergence statements in [30, 31] — more general in that the domain can vary with $\delta$  as $\delta \to 0$  , but more special in the restriction to Jordan domains.

Given a Jordan domain $D$  with two distinct points $a,b\in \partial D$  on its boundary, let ${\mu }_{D,a,b}$  denote the law of ${\gamma }_{D,a,b}$  , the trace of chordal $SL{E}_6$  , and let ${\mu }_{D,a,b}^{\delta }$  denote the law of the percolation exploration path ${\gamma }_{D,a,b}^{\delta }$  . Let $X$  be the space of continuous curves inside $D$  from $a$  to $b$  . We define $\rho ({\mu }_{D,a,b},{\mu }_{D,a,b}^{\delta })\equiv inf\{\varepsilon >0:{\mu }_{D,a,b}(U)\le {\mu }_{D,a,b}^{\delta }({\cup }_{x\in U}{B}_{\text{d}}(x,\varepsilon ))+\varepsilon \text{for all Borel}U\subset X\}$  (where $B_{\text{d}}(x,\varepsilon )$  denotes the open ball of radius $\varepsilon$  centered at $x$  in the metric ( 2 )) and denote by ${\text{d}}_{\text{P}}({\mu }_{D,a,b},{\mu }_{D,a,b}^{\delta })\equiv max\left\{\rho \right({\mu }_{D,a,b},{\mu }_{D,a,b}^{\delta }),\rho ({\mu }_{D,a,b}^{\delta },{\mu }_{D,a,b}\left)\right\}$  the Prohorov distance; weak convergence is equivalent to convergence in the Prohorov metric. Statement (S) is the following; it is used in the proofs of all the results of this section except for Lemmas  5.1 - 5.2 .

• (S) For Jordan domains, there is convergence in distribution of the percolation exploration path to the trace of chordal $SL{E}_6$  that is locally uniform in the shape of the boundary with respect to the uniform metric on continuous curves ( 2 ), and in the location of the starting and ending points with respect to the Euclidean metric; i.e., for $(D,a,b)$  a Jordan domain with $a,b\in \partial D$  , $\forall \varepsilon >0$  , $\exists {\alpha }_0={\alpha }_0\left(\varepsilon \right)$  and ${\delta }_0={\delta }_0\left(\varepsilon \right)$  such that for all $(D^{\prime },a^{\prime },b^{\prime })$  with $D^{\prime }$  Jordan and with $max\left(\text{d}\right(\partial D,\partial D^{\prime }),|a-a^{\prime }|,|b-b^{\prime }\left|\right)\le {\alpha }_0$  and $\delta \le {\delta }_0$  , ${\text{d}}_{\text{P}}({\mu }_{D^{\prime },a^{\prime },b^{\prime }},{\mu }_{D^{\prime },a^{\prime },b^{\prime }}^{\delta })\le \varepsilon$  .

5.1 Preliminary Results

We first give some important results which are needed in the proofs of the main theorems. We start with two lemmas which are consequences of [2] , of standard bounds on the probability of events corresponding to having a certain number of monochromatic crossings of an annulus (see Lemma 5 of [16] and Appendix A of [20] ), but which do not depend on statement (S).

Lemma 5.1. Let ${\gamma }_{\mathbb{D},-i,i}^{\delta }$  be the percolation exploration path on the edges of $\delta \mathcal{\mathscr{H}}$  inside (the $\delta$  -approximation of ) $\mathbb{D}$  between (the e-vertices closest to) $-i$  and $i$  . For any fixed point $z\in \mathbb{D}$  , chosen independently of ${\gamma }_{\mathbb{D},-i,i}^{\delta }$  , as $\delta \to 0$  , ${\gamma }_{\mathbb{D},-i,i}^{\delta }$  and the boundary $\partial {\mathbb{D}}_{-i,i}^{\delta }\left(z\right)$  of the domain ${\mathbb{D}}_{-i,i}^{\delta }\left(z\right)$  that contains $z$  jointly have limits in distribution along subsequences of $\delta$  with respect to the uniform metric ( 2 ) on continuous curves. Moreover, any subsequence limit of $\partial {\mathbb{D}}_{-i,i}^{\delta }\left(z\right)$  is almost surely a simple loop [3] .

Lemma 5.2. Using the notation of Lemma  5.1 , let ${\gamma }_{\mathbb{D},-i,i}$  be the limit in distribution of ${\gamma }_{\mathbb{D},-i,i}^{\delta }$  as $\delta \to 0$  along some convergent subsequence $\left\{{\delta }_k\right\}$  and $\partial {\mathbb{D}}_{-i,i}\left(z\right)$  the boundary of the domain ${\mathbb{D}}_{-i,i}\left(z\right)$  of $\mathbb{D}\backslash {\gamma }_{D,-i,i}[0,1]$  that contains $z$  .

Then, as $k\to \infty$  , $({\gamma }_{\mathbb{D},-i,i}^{{\delta }_k},\partial {\mathbb{D}}_{-i,i}^{{\delta }_k}(z\left)\right)$  converges in distribution to $({\gamma }_{\mathbb{D},-i,i},\partial {\mathbb{D}}_{-i,i}(z\left)\right)$  .

The two lemmas above are important ingredients in the proof of Theorem  1 below. The second one says that, for every subsequence limit, the discrete boundaries converge to the boundaries of the domains generated by the limiting continuous curve. If we use statement (S), then the limit ${\gamma }_{\mathbb{D},-i,i}$  of ${\gamma }_{\mathbb{D},-i,i}^{{\delta }_k}$  is the trace of chordal $SL{E}_6$  for every subsequence ${\delta }_k\downarrow 0$  , and we can use Lemmas  5.2 and  5.1 to deduce that all the domains produced in the continuum construction are Jordan domains. The key step in that direction is represented by the following result, our proof of which relies on (S).

Corollary 5.1. For any deterministic $z\in \mathbb{D}$  , the boundary $\partial {\mathbb{D}}_{-i,i}\left(z\right)$  of a domain ${\mathbb{D}}_{-i,i}\left(z\right)$  of the continuum construction is almost surely a Jordan curve.

The corollary says that the domains that appear after the first step of the continuum construction are Jordan domains. The steps in the second stage of the continuum construction consist of $SL{E}_6$  paths inside Jordan domains, and therefore Corollary  5.1 , combined with Riemann's mapping theorem and the conformal invariance of $SL{E}_6$  , implies that the domains produced during the second stage are also Jordan. By induction, we deduce that all the domains produced in the continuum construction are Jordan domains.

We end this section with one more lemma which is another key ingredient in the proof of Theorem  1 ; we remark that its proof requires (S) in a fundamental way.

Lemma 5.3. Let $(D,a,b)$  denote a random Jordan domain, with $a,b$  two points on $\partial D$  . Let $\left\{\right(D_k,a_k,b_k){\}}_{k\in \mathbb{N}},a_k,b_k\in \partial D_k$  , be a sequence of random Jordan domains with points on their boundaries such that, as $k\to \infty$  , $(\partial D_k,a_k,b_k)$  converges in distribution to $(\partial D,a,b)$  with respect to the uniform metric ( 2 ) on continuous curves, and the Euclidean metric on $(a,b)$  . For any sequence $\{{\delta }_k{\}}_{k\in \mathbb{N}}$  with ${\delta }_k\downarrow 0$  as $k\to \infty$  , ${\gamma }_{D_k,a_k,b_k}^{{\delta }_k}$  converges in distribution to ${\gamma }_{D,a,b}$  with respect to the uniform metric ( 2 ) on continuous curves.

5.2 The Main Theorems

In this section we state the main theorems of this paper and a corollary, our most important result, that the Continuum Nonsimple Loop process is the scaling limit of the set of all cluster boundaries for critical site percolation on the triangular lattice. The corollary is obtained by combining the first two theorems. The proofs of all these results rely on statement (S). As noted before, statement (S) is proved in Appendix  A .

Theorem 1. For any $k\in \mathbb{N}$  , the first $k$  steps of (a suitably reorganized version of ) the full discrete construction inside the unit disc (of Section  4.3 ) converge, jointly in distribution, to the first $k$  steps of the full continuum construction inside the unit disc (of Section  3.2 ). Furthermore, the scaling limit of the full (original or reorganized) discrete construction is the full continuum construction.

Moreover, if for any fixed $\varepsilon >0$  we let $K_{\delta }\left(\varepsilon \right)$  denote the number of steps needed to find all the cluster boundaries of Euclidean diameter larger than $\varepsilon$  in the discrete construction, then $K_{\delta }\left(\varepsilon \right)$  is bounded in probability as $\delta \to 0$  ; i.e., ${lim}_{C\to \infty }{limsup}_{\delta \to 0}\mathbb{P}\left(K_{\delta }\right(\varepsilon )>C)=0$  . This is so in both the original and reorganized versions of the discrete construction.

The second part of Theorem  1 means that both versions of the discrete construction used in the theorem find all large contours in a number of steps which does not diverge as $\delta \to 0$  . This, together with the first part of the same theorem, implies that the continuum construction does indeed describe all macroscopic contours contained inside the unit disc (with blue boundary conditions) as $\delta \to 0$  .

The construction presented in Section  3.2 can of course be repeated for the disc ${\mathbb{D}}_R$  of radius $R$  , for any $R$  , so we should take a “thermodynamic limit” by letting $R\to \infty$  . In this way, we would eliminate the boundary (and the boundary conditions) and obtain a process on the whole plane. Such an extension from the unit disc to the plane is contained in the next theorem.

Let $P_R$  be the (limiting) distribution of the set of curves (all continuum nonsimple loops) generated by the continuum construction inside ${\mathbb{D}}_R$  (i.e., the limiting measure, defined by the inductive construction, on the complete separable metric space ${\Omega }_R$  of collections of continuous curves in ${\mathbb{D}}_R$  ).

For a domain $D$  , we denote by $I_D$  the mapping (on $\Omega$  or ${\Omega }_R$  ) in which all portions of curves that exit $D$  are removed. When applied to a configuration of loops in the plane, $I_D$  gives a set of curves which either start and end at points on $\partial D$  or form closed loops completely contained in $D$  . Let ${\hat{I}}_D$  be the same mapping lifted to the space of probability measures on $\Omega$  or ${\Omega }_R$  .

Theorem 2. Theorem  1 implies that there exists a unique probability measure $P$  on the space $\Omega$  of collections of continuous curves in ${\dot{\mathbb{R}}}^2$  such that $P_R\to P$  as $R\to \infty$  in the sense that for every bounded domain $D$  , as $R\to \infty$  , ${\hat{I}}_D{P}_R\to {\hat{I}}_{D}P$  .

Remark 5.1. We remark that we will generally take monochromatic blue boundary conditions on the disc ${\mathbb{D}}_R$  of radius $R$  . But one could also takemonochromatic boundary conditions with color depending on $R$  or even non-monochromatic boundary conditions without any essential change in the results or the proofs.

Corollary 5.2. The Continuum Nonsimple Loop process $P$  in the plane defined in Theorem  2 is the scaling limit of the collection of all boundary contours for critical site percolation on the triangular lattice.

The next theorem states some properties of the Continuum Nonsimple Loop process in the plane.

Theorem 3. The Continuum Nonsimple Loop process in the plane has the following properties, the first three of which are valid with probability one:

• 1. The Continuum Nonsimple Loop process is a random collection of noncrossing continuous loops in the plane. The loops can and do touch themselves and each other many times, but there are no triple points; i.e. no three or more loops can come together at the same point, and a single loop cannot touch the same point more than twice, nor can a loop touch a point where another loop touches itself.
• 2. Any deterministic point (i.e., chosen independently of the loop process) of the plane is surrounded by an infinite family of nested loops with diameters going to both zero and infinity; any annulus about that point with inner radius $r_1>0$  and outer radius $r_2<\infty$  contains only a finite number of those loops. Consequently, any two distinct deterministic points of the plane are separated by loops winding around each of them.
• 3. Any two loops are connected by a finite “path” of touching loops.
• 4. The Continuum Nonsimple Loop process is conformally invariant in the sense that, given a Jordan domain $D$  and a conformal homeomorphism $f:D\to D^{\prime }$  onto $D^{\prime }$  , the scaling limits, $P_D$  and $P_{D^{\prime }}$  , of the loops inside $D$  and $D^{\prime }$  taken with, say, blue boundary conditions are related by $f*P_D=P_{D^{\prime }}$  . (Here $f*P_D$  denotes the probability distribution of the loop process $f\left(X\right)$  when $X$  is distributed by $P_D$  .)

To conclude this section, we show how to recover chordal $SL{E}_6$  from the Continuum Nonsimple Loop process, i.e., given a (deterministic) Jordan domain $D$  with two boundary points $a$  and $b$  , we give a construction that uses the continuum nonsimple loops of $P$  to generate a process distributed like chordal $SL{E}_6$  inside $D$  from $a$  to $b$  .

Remember, first of all, that each continuum nonsimple loop has either a clockwise or counterclockwise direction, with the set of all loops surrounding any deterministic point alternating in direction. For convenience, let us suppose that $a$  is at the “bottom” and $b$  is at the “top” of $D$  so that the boundary is divided into a left and right part by thesetwo points. Fix $\varepsilon >0$  and call $LR\left(\varepsilon \right)$  the set of all the directed segments of loops that connect from the left to the right part of the boundary touching $\partial D$  at a distance larger than $\varepsilon$  from both $a$  and $b$  , and $RL\left(\varepsilon \right)$  the analogous set of directed segments from the right to the left portion of $\partial D$  . For a fixed $\varepsilon >0$  , there is only a finite number of such segments, and, if they are ordered moving along the left boundary of $D$  from $a$  to $b$  , they alternate in direction (i.e., a segment in $LR\left(\varepsilon \right)$  is followed by one in $RL\left(\varepsilon \right)$  and so on).

Between a segment in $RL\left(\varepsilon \right)$  and the next segment in $LR\left(\varepsilon \right)$  , there are countably many portions of loops intersecting $D$  which start and end on $\partial D$  and are maximal in the sense that they are not contained inside any other portion of loop of the same type; they all have counterclockwise direction and can be used to make a “bridge” between the right-to-left segment and the next one (in $LR\left(\varepsilon \right)$  ). This is done by pasting the portions of loops together with the help of points in $\partial D$  and a limit procedure to produce a connected (nonsimple) path.

If we do this for each pair of successive segments on both sides of the boundary of $D$  , we get a path that connects two points on $\partial D$  . By letting $\varepsilon \to 0$  and taking the limit of this procedure, since almost surely $a$  and $b$  are surrounded by an infinite family of nested loops with diameters going to zero, we obtain a path that connects $a$  with $b$  ; this path is distributed as chordal $SL{E}_6$  inside $D$  from $a$  to $b$  . The last claim follows from considering the analogous procedure for percolation on the discrete lattice $\delta \mathcal{\mathscr{H}}$  , using segments of boundaries. It is easy to see that in the discrete case this procedure produces exactly the same path as the percolation exploration process. By Corollary  5.2 , the scaling limit of this discrete procedure is the continuum one described above, therefore the claim follows from (S).

6 Proofs

In this section we present the proofs of the results stated in Section  5 .

Proof of Lemma  5.1 . The first part of the lemma is a direct consequence of [2] ; it is enough to notice that the (random) polygonal curves ${\gamma }_{\mathbb{D},-i,i}^{\delta }$  and $\partial {\mathbb{D}}_{-i,i}^{\delta }\left(z\right)$  satisfy the conditions in [2] and thus have a scaling limit in terms of continuous curves, at least along subsequences of $\delta$  .

To prove the second part, we use standard percolation bounds (see Lemma 5 of [16] and Appendix A of [20] ) to show that, in the limit $\delta \to 0$  , the loop $\partial {\mathbb{D}}_{-i,i}^{\delta }\left(z\right)$  does not collapse on itself but remains a simple loop.

Let us assume that this is not the case and that the limit $\tilde{\gamma }$  of $\partial {\mathbb{D}}_{-i,i}^{{\delta }_k}\left(z\right)$  along some subsequence $\{{\delta }_k{\}}_{k\in \mathbb{N}}$  touches itself, i.e., $\tilde{\gamma }\left(t_0\right)=\tilde{\gamma }\left(t_1\right)$  for $t_0\ne t_1$  with positive probability.

If that happens, we can take $\varepsilon >{\varepsilon }^{\prime }>0$  small enough so that the annulus $B\left(\tilde{\gamma }\right(t_1),\varepsilon )\backslash B\left(\tilde{\gamma }\right(t_1),{\varepsilon }^{\prime })$  is crossed at least four times by $\tilde{\gamma }$  (here $B(u,r)$  is the ball of radius $r$  centered at $u$  ).

Because of the choice of topology, the convergence in distribution of $\partial {\mathbb{D}}_{-i,i}^{{\delta }_k}\left(z\right)$  to $\tilde{\gamma }$  implies that we can find coupled versions of $\partial {\mathbb{D}}_{-i,i}^{{\delta }_k}\left(z\right)$  and $\tilde{\gamma }$  on some probability space $({\Omega }^{\prime },{\mathcal{ℬ}}^{\prime },{\mathbb{P}}^{\prime })$  such that $\text{d}(\partial {\mathbb{D}}_{-i,i}^{\delta }(z),\tilde{\gamma })\to 0$  , for all ${\omega }^{\prime }\in {\Omega }^{\prime }$  as $k\to \infty$  (see, for example, Corollary 1 of [6] ).

Using this coupling, we can choose $k$  large enough (depending on ${\omega }^{\prime }$  ) so that $\partial {\mathbb{D}}_{-i,i}^{{\delta }_k}\left(z\right)$  stays in an ${\varepsilon }^{\prime }/2$  -neighborhood $\mathcal{N}(\tilde{\gamma },{\varepsilon }^{\prime }/2)\equiv {\cup }_{u\in \tilde{\gamma }}B(u,{\varepsilon }^{\prime }/2)$  of $\tilde{\gamma }$  . This event however would correspond to (at least) four paths of one color (corresponding to the four crossings by $\partial {\mathbb{D}}_{-i,i}^{{\delta }_k}\left(z\right)$  ) and two of the other color of the annulus $B\left(\tilde{\gamma }\right(t_1),\varepsilon -{\varepsilon }^{\prime }/2)\backslash B\left(\tilde{\gamma }\right(t_1),3{\varepsilon }^{\prime }/2)$  .

As ${\delta }_k\to 0$  , we can let ${\varepsilon }^{\prime }\to 0$  , in which case the probability of seeing the event just described somewhere inside $\mathbb{D}$  goes to zero [16, 20] , leading to a contradiction.

Proof of Lemma  5.2 . Let $\{{\delta }_k{\}}_{k\in \mathbb{N}}$  be a convergent subsequence for ${\gamma }_{\mathbb{D},-i,i}^{\delta }$  and $\gamma \equiv {\gamma }_{\mathbb{D},-i,i}$  the limit in distribution of ${\gamma }_{\mathbb{D},-i,i}^{{\delta }_k}$  as $k\to \infty$  . For simplicity of notation, in the rest of the proof we will drop the $k$  and write $\delta$  instead of ${\delta }_k$  . Because of the choice of topology, the convergence in distribution of ${\gamma }^{\delta }\equiv {\gamma }_{\mathbb{D},-i,i}^{\delta }$  to $\gamma$  implies that we can find coupled versions of ${\gamma }^{\delta }$  and $\gamma$  on some probability space $({\Omega }^{\prime },{\mathcal{ℬ}}^{\prime },{\mathbb{P}}^{\prime })$  such that $\text{d}\left({\gamma }^{\delta }\right({\omega }^{\prime }),\gamma ({\omega }^{\prime }\left)\right)\to 0$  , for all ${\omega }^{\prime }$  as $k\to \infty$  (see, for example, Corollary 1 of [6] ). Using this coupling, our first task will be to prove the following claim:

• (C) For two (deterministic) points $u,v\in \mathbb{D}$  , the probability that ${\mathbb{D}}_{-i,i}\left(u\right)={\mathbb{D}}_{-i,i}\left(v\right)$  but ${\mathbb{D}}_{-i,i}^{\delta }\left(u\right)\ne {\mathbb{D}}_{-i,i}^{\delta }\left(v\right)$  or vice versa goes to zero as $\delta \to 0$  .

Let us consider first the case of $u,v$  such that ${\mathbb{D}}_{-i,i}\left(u\right)={\mathbb{D}}_{-i,i}\left(v\right)$  but ${\mathbb{D}}_{-i,i}^{\delta }\left(u\right)\ne {\mathbb{D}}_{-i,i}^{\delta }\left(v\right)$  . Since ${\mathbb{D}}_{-i,i}\left(u\right)$  is an open subset of $\mathbb{C}$  , there exists a continuous curve ${\gamma }_{u,v}$  joining $u$  and $v$  and a constant $\varepsilon >0$  such that the $\varepsilon$  -neighborhood $\mathcal{N}({\gamma }_{u,v},\varepsilon )$  of the curve is contained in ${\mathbb{D}}_{-i,i}\left(u\right)$  , which implies that $\gamma$  does not intersect $\mathcal{N}({\gamma }_{u,v},\varepsilon )$  . Now, if ${\gamma }^{\delta }$  does not intersect $\mathcal{N}({\gamma }_{u,v},\varepsilon /2)$  , for $\delta$  small enough, then there is a $\mathcal{T}$  -path $\pi$  of unexplored hexagons connecting the hexagon that contains $u$  with the hexagon that contains $v$  , and we conclude that ${\mathbb{D}}_{-i,i}^{\delta }\left(u\right)={\mathbb{D}}_{-i,i}^{\delta }\left(v\right)$  . This shows that the event that ${\mathbb{D}}_{-i,i}\left(u\right)={\mathbb{D}}_{-i,i}\left(v\right)$  but ${\mathbb{D}}_{-i,i}^{\delta }\left(u\right)\ne {\mathbb{D}}_{-i,i}^{\delta }\left(v\right)$  implies the existence of a curve ${\gamma }_{u,v}$  whose $\varepsilon$  -neighborhood $\mathcal{N}({\gamma }_{u,v},\varepsilon )$  is not intersected by $\gamma$  but whose $\varepsilon /2$  -neighborhood $\mathcal{N}({\gamma }_{u,v},\varepsilon /2)$  is intersected by ${\gamma }^{\delta }$  . This implies that $\forall u,v\in \mathbb{D}$  , $\exists \varepsilon >0$  such that ${\mathbb{P}}^{\prime }\left({\mathbb{D}}_{-i,i}\right(u)={\mathbb{D}}_{-i,i}(v\left)\text{but}{\mathbb{D}}_{-i,i}^{\delta }\right(u)\ne {\mathbb{D}}_{-i,i}^{\delta }(v\left)\right)\le {\mathbb{P}}^{\prime }\left(\text{d}\right({\gamma }^{\delta },\gamma )\ge \varepsilon /2)$  . But the right hand side goes to zero for every $\varepsilon >0$  as $\delta \to 0$  , which concludes the proof of one direction of the claim.

To prove the other direction, we consider two points $u,v\in D$  such that $D_{-i,i}\left(u\right)\ne D_{-i,i}\left(v\right)$  but $D_{-i,i}^{\delta }\left(u\right)=D_{-i,i}^{\delta }\left(v\right)$  . Assume that $u$  is trapped before $v$  by $\gamma$  and suppose for the moment that ${\mathbb{D}}_{-i,i}\left(u\right)$  is a domain of type 3 or 4; the case of a domain of type 1 or 2 is analogous and will be treated later. Let $t_1$  be the first time $u$  is trapped by $\gamma$  with $\gamma \left(t_0\right)=\gamma \left(t_1\right)$  the double point of $\gamma$  where the domain ${\mathbb{D}}_{-i,i}\left(u\right)$  containing $u$  is “sealed off.” At time $t_1$  , a new domain containing $u$  is created and $v$  is disconnected from $u$  .

Choose $\varepsilon >0$  small enough so that neither $u$  nor $v$  is contained in the ball $B\left(\gamma \right(t_1),\varepsilon )$  of radius $\varepsilon$  centered at $\gamma \left(t_1\right)$  , nor in the $\varepsilon$  -neighborhood $\mathcal{N}\left(\gamma \right[t_0,t_1],\varepsilon )$  of the portion of $\gamma$  which surrounds $u$  . Then it follows from the coupling that, for $\delta$  small enough, there are appropriate parameterizations of $\gamma$  and ${\gamma }^{\delta }$  such that the portion ${\gamma }^{\delta }[t_0,t_1]$  of ${\gamma }^{\delta }\left(t\right)$  is inside $\mathcal{N}\left(\gamma \right[t_0,t_1],\varepsilon )$  , and ${\gamma }^{\delta }\left(t_0\right)$  and ${\gamma }^{\delta }\left(t_1\right)$  are contained in $B\left(\gamma \right(t_1),\varepsilon )$  .

For $u$  and $v$  to be contained in the same domain in the discrete construction, there must be a $\mathcal{T}$  -path $\pi$  of unexplored hexagons connecting the hexagon that contains $u$  to the hexagon that contains $v$  . From what we said in the previous paragraph, any such $\mathcal{T}$  -path connecting $u$  and $v$  would have to go though a “bottleneck” in $B\left(\gamma \right(t_1),\varepsilon )$  .

Assume now, for concreteness but without loss of generality, that ${\mathbb{D}}_{-i,i}\left(u\right)$  is a domain of type 3, which means that $\gamma$  winds around $u$  counterclockwise, and consider the hexagons to the “left” of ${\gamma }^{\delta }[t_0,t_1]$  . Those hexagons form a “quasi-loop” around $u$  since they wind around it (counterclockwise) and the first and last hexagons are both contained in $B\left(\gamma \right(t_1),\varepsilon )$  . The hexagons to the left of ${\gamma }^{\delta }[t_0,t_1]$  belong to the set ${\Gamma }_Y\left({\gamma }^{\delta }\right)$  , which can be seen as a (nonsimple) path by connecting the centers of the hexagons in ${\Gamma }_Y\left({\gamma }^{\delta }\right)$  by straight segments. Such a path shadows ${\gamma }^{\delta }$  , with the difference that it can have double (or even triple) points, since the same hexagon can be visited more than once. Consider ${\Gamma }_Y\left({\gamma }^{\delta }\right)$  as a path ${\hat{\gamma }}^{\delta }$  with a given parametrization ${\hat{\gamma }}^{\delta }\left(t\right)$  , chosen so that ${\hat{\gamma }}^{\delta }\left(t\right)$  is inside $B\left(\gamma \right(t_1),\varepsilon )$  when ${\gamma }^{\delta }\left(t\right)$  is, and it winds around $u$  together with ${\gamma }^{\delta }\left(t\right)$  .

Now suppose that there were two times, ${\hat{t}}_0$  and ${\hat{t}}_1$  , such that ${\hat{\gamma }}^{\delta }\left({\hat{t}}_1\right)={\hat{\gamma }}^{\delta }\left({\hat{t}}_0\right)\in B\left(\gamma \right(t_1),\varepsilon )$  and ${\hat{\gamma }}^{\delta }[{\hat{t}}_0,{\hat{t}}_1]$  winds around $u$  . This would imply that the “quasi-loop” of explored yellow hexagons around $u$  is actually completed, and that $D_{a,b}^{\delta }\left(v\right)\ne D_{a,b}^{\delta }\left(u\right)$  .

Thus, for $u$  and $v$  to belong to the same discrete domain, this cannot happen.

For any $0<{\varepsilon }^{\prime }<\varepsilon$  , if we take $\delta$  small enough, ${\hat{\gamma }}^{\delta }$  will be contained inside $\mathcal{N}(\gamma ,{\varepsilon }^{\prime })$  , due to the coupling. Following the considerations above, the fact that $u$  and $v$  belong to the same domain in the discrete construction but to different domains in the continuum construction implies, for $\delta$  small enough, that there are four disjoint yellow $\mathcal{T}$  -paths crossing the annulus $B\left(\gamma \right(t_1),\varepsilon )\backslash B\left(\gamma \right(t_1),{\varepsilon }^{\prime })$  (the paths have to be disjoint because, as we said, ${\hat{\gamma }}^{\delta }$  cannot, when coming back to $B\left(\gamma \right(t_1),\varepsilon )$  after winding around $u$  , touch itself inside $B\left(\gamma \right(t_1),\varepsilon )$  ). Since $B\left(\gamma \right(t_1),\varepsilon )\backslash B\left(\gamma \right(t_1),{\varepsilon }^{\prime })$  is also crossed by at least two blue $\mathcal{T}$  -paths from ${\Gamma }_B\left({\gamma }^{\delta }\right)$  , there is a total of at least six $\mathcal{T}$  -paths, not all of the same color,crossing the annulus $B\left(\gamma \right(t_1),\varepsilon )\backslash B\left(\gamma \right(t_1),{\varepsilon }^{\prime })$  .

Let us call ${\mathcal{A}}_w(\varepsilon ,{\varepsilon }^{\prime })$  the event described above, where $\gamma \left(t_1\right)=w$  ; a standard bound [16] on the probability of six disjoint crossings (not all of the same color) of an annulus gives that the probability of ${\mathcal{A}}_w(\varepsilon ,{\varepsilon }^{\prime })$  scales as $(\frac{{\varepsilon }^{\prime }}{\varepsilon }{)}^{2+\alpha }$  with $\alpha >0$  . As $\delta \to 0$  , we can let ${\varepsilon }^{\prime }$  go to zero (keeping $\varepsilon$  fixed); when we do this, the probability of ${\mathcal{A}}_w(\varepsilon ,{\varepsilon }^{\prime })$  goes to zero sufficiently rapidly with ${\varepsilon }^{\prime }$  to conclude, like in the proof of Lemma  5.1 , that the probability to see such an event anywhere in $\mathbb{D}$  goes to zero.

In the case in which $u$  belongs to a domain of type 1 or 2, let $\mathcal{ℰ}$  be the excursion that traps $u$  and $\gamma \left(t_0\right)\in \partial \mathbb{D}$  be the point on the boundary of $\mathbb{D}$  where $\mathcal{ℰ}$  starts and $\gamma \left(t_1\right)\in \partial \mathbb{D}$  the point where it ends. Choose $\varepsilon >0$  small enough so that neither $u$  nor $v$  is contained in the balls $B\left(\gamma \right(t_0),\varepsilon )$  and $B\left(\gamma \right(t_1),\varepsilon )$  of radius $\varepsilon$  centered at $\gamma \left(t_0\right)$  and $\gamma \left(t_1\right)$  , nor in the $\varepsilon$  -neighborhood $\mathcal{N}(\mathcal{ℰ},\varepsilon )$  of the excursion $\mathcal{ℰ}$  . Because of the coupling, for $\delta$  small enough (depending on $\varepsilon$  ), ${\gamma }^{\delta }$  shadows $\gamma$  along $\mathcal{ℰ}$  , staying within $\mathcal{N}(\mathcal{ℰ},\varepsilon )$  . If this is the case, any $\mathcal{T}$  -path of unexplored hexagons connecting the hexagon that contains $u$  with the hexagon that contains $v$  would have to go through one of two “bottlenecks,” one contained in $B\left(\gamma \right(t_0),\varepsilon )$  and the other in $B\left(\gamma \right(t_1),\varepsilon )$  .

Assume for concreteness (but without loss of generality) that $u$  is in a domain of type 1, which means that $\gamma$  winds around $u$  counterclockwise. If we parameterize $\gamma$  and ${\gamma }^{\delta }$  so that ${\gamma }^{\delta }\left(t_0\right)\in B\left(\gamma \right(t_0),\varepsilon )$  and ${\gamma }^{\delta }\left(t_1\right)\in B\left(\gamma \right(t_1),\varepsilon )$  , ${\gamma }^{\delta }[t_0,t_1]$  forms a “quasi-excursion” around $u$  since it winds around it (counterclockwise) and it starts inside $B_{\varepsilon }\left(\gamma \right(t_0\left)\right)$  and ends inside $B_{\varepsilon }\left(\gamma \right(t_1\left)\right)$  . Notice that if ${\gamma }^{\delta }$  touched $\partial {\mathbb{D}}^{\delta }$  , inside both $B_{\varepsilon }\left(\gamma \right(t_0\left)\right)$  and $B_{\varepsilon }\left(\gamma \right(t_1\left)\right)$  , this would imply that the “quasi-excursion” is a real excursion and that $D_{a,b}^{\delta }\left(v\right)\ne D_{a,b}^{\delta }\left(u\right)$  .

For any $0<{\varepsilon }^{\prime }<\varepsilon$  , if we take $\delta$  small enough, ${\gamma }^{\delta }$  will be contained inside $\mathcal{N}(\gamma ,{\varepsilon }^{\prime })$  , due to the coupling. Therefore, the fact that ${\mathbb{D}}_{a,b}^{\delta }\left(v\right)={\mathbb{D}}_{a,b}^{\delta }\left(u\right)$  implies, with probability going to one as $\delta \to 0$  , that for $\varepsilon >0$  fixed and any $0<{\varepsilon }^{\prime }<\varepsilon$  , ${\gamma }^{\delta }$  enters the ball $B\left(\gamma \right(t_i),{\varepsilon }^{\prime })$  and does not touch $\partial {\mathbb{D}}^{\delta }$  inside the larger ball $B\left(\gamma \right(t_i),\varepsilon )$  , for $i=0$  or $1$  . This is equivalent to having at least two yellow and one blue $\mathcal{T}$  -paths (contained in ${\mathbb{D}}^{\delta }$  ) crossing the annulus $B\left(\gamma \right(t_i),\varepsilon )\backslash B\left(\gamma \right(t_i),{\varepsilon }^{\prime })$  . Let us call ${\mathcal{ℬ}}_w(\varepsilon ,{\varepsilon }^{\prime })$  the event described above, where $\gamma \left(t_i\right)=w$  ; a standard bound [20] (this bound can also be derived from the one obtained in [16] ) on the probability of disjoint crossings (not all of the same color) of a semi-annulus in the upper half-plane gives that the probability of ${\mathcal{ℬ}}_w(\varepsilon ,{\varepsilon }^{\prime })$  scales as $(\frac{{\varepsilon }^{\prime }}{\varepsilon }{)}^{1+\beta }$  with $\beta >0$  . (We can apply the bound to our case because the unit disc is a convex subset of the half-plane $\{x+iy:y>-1\}$  and therefore the intersection of an annulus centered at say $-i$  with the unit disc is a subset of the intersection of the same annulus with the half-plane $\{x+iy:y>-1\}$  .) As $\delta \to 0$  , we can let ${\varepsilon }^{\prime }$  go to zero (keeping $\varepsilon$  fixed), concluding that the probability that such an event occurs anywhere on the boundary of the disc goes to zero.

We have shown that, for two fixed points $u,v\in \mathbb{D}$  , having ${\mathbb{D}}_{-i,i}\left(u\right)\ne {\mathbb{D}}_{-i,i}\left(v\right)$  but ${\mathbb{D}}_{-i,i}^{\delta }\left(u\right)={\mathbb{D}}_{-i,i}^{\delta }\left(v\right)$  or vice versa implies the occurrence of an event whose probability goes to zero as $\delta \to 0$  , and the proof of the claim is concluded.

We now introduce the Hausdorff distance ${\text{d}}_{\text{H}}(A,B)$  between two closed nonemptysubsets of $\overline{\mathbb{D}}$  :

$$\begin{array}{c}{\text{d}}_{\text{H}}(A,B)\equiv inf\{\ell \ge 0:B\subset {\cup }_{a\in A}B(a,\ell ),A\subset {\cup }_{b\in B}B(b,\ell \left)\right\}.\end{array}$$

(9)

With this metric, the collection of closed subsets of $\overline{\mathbb{D}}$  is a compact space. We will next prove that $\partial {\mathbb{D}}_{-i,i}^{\delta }\left(z\right)$  converges in distribution to $\partial {\mathbb{D}}_{-i,i}\left(z\right)$  as $\delta \to 0$  , in the topology induced by ( 9 ). (Notice that the coupling between ${\gamma }^{\delta }$  and $\gamma$  provides a coupling between $\partial {\mathbb{D}}_{-i,i}^{\delta }\left(z\right)$  and $\partial {\mathbb{D}}_{-i,i}\left(z\right)$  , seen as boundaries of domains produced by the two paths.) We will now use Lemma  5.1 and take a further subsequence $k_n$  of the $\delta$  's that for simplicity of notation we denote by $\{{\delta }_n{\}}_{n\in \mathbb{N}}$  such that, as $n\to \infty$  , $\{{\gamma }^{{\delta }_n},\partial {\mathbb{D}}_{-i,i}^{{\delta }_n}(z\left)\right\}$  converge jointly in distribution to $\{\gamma ,\tilde{\gamma }\}$  , where $\tilde{\gamma }$  is a simple loop. For any $\varepsilon >0$  , since $\tilde{\gamma }$  is a compact set, we can find a covering of $\tilde{\gamma }$  by a finite number of balls of radius $\varepsilon /2$  centered at points on $\tilde{\gamma }$  . Each ball contains both points in the interior $\text{int}\left(\tilde{\gamma }\right)$  of $\tilde{\gamma }$  and in the exterior $\text{ext}\left(\tilde{\gamma }\right)$  of $\tilde{\gamma }$  , and we can choose (independently of $n$  ) one point from $\text{int}\left(\tilde{\gamma }\right)$  and one from $\text{ext}\left(\tilde{\gamma }\right)$  inside each ball.

Once again, the convergence in distribution of $\partial {\mathbb{D}}_{-i,i}^{{\delta }_n}\left(z\right)$  to $\tilde{\gamma }$  implies the existence of a coupling such that, for $n$  large enough, the selected points that are in $\text{int}\left(\tilde{\gamma }\right)$  are contained in ${\mathbb{D}}_{-i,i}^{{\delta }_n}\left(z\right)$  , and those that are in $\text{ext}\left(\tilde{\gamma }\right)$  are contained in the complement of $\overline{{\mathbb{D}}_{-i,i}^{{\delta }_n}\left(z\right)}$  . But by claim (C), each one of the selected points that is contained in ${\mathbb{D}}_{-i,i}^{{\delta }_n}\left(z\right)$  is also contained in ${\mathbb{D}}_{-i,i}\left(z\right)$  with probability going to $1$  as $n\to \infty$  ; analogously, each one of the selected points contained in the complement of $\overline{{\mathbb{D}}_{-i,i}^{{\delta }_n}\left(z\right)}$  is also contained in the complement of $\overline{{\mathbb{D}}_{-i,i}\left(z\right)}$  with probability going to $1$  as $n\to \infty$  . This implies that $\partial {\mathbb{D}}_{-i,i}\left(z\right)$  crosses each one of the balls in the covering of $\tilde{\gamma }$  , and therefore $\tilde{\gamma }\subset {\cup }_{u\in \partial {\mathbb{D}}_{-i,i}\left(z\right)}B(u,\varepsilon )$  . From this and the coupling between $\partial {\mathbb{D}}_{-i,i}^{{\delta }_n}\left(z\right)$  and $\tilde{\gamma }$  , it follows immediately that, for $n$  large enough, $\partial {\mathbb{D}}_{-i,i}^{{\delta }_n}\left(z\right)\subset {\cup }_{u\in \partial {\mathbb{D}}_{-i,i}\left(z\right)}B(u,\varepsilon )$  with probability close to one.

A similar argument (analogous to the previous one but simpler, since it does not require the use of $\tilde{\gamma }$  ), with the roles of ${\mathbb{D}}_{-i,i}^{{\delta }_n}\left(z\right)$  and ${\mathbb{D}}_{-i,i}\left(z\right)$  inverted, shows that $\partial {\mathbb{D}}_{-i,i}\left(z\right)\subset {\cup }_{u\in \partial {\mathbb{D}}_{-i,i}^{{\delta }_n}\left(z\right)}B(u,\varepsilon )$  with probability going to $1$  as $n\to \infty$  . Therefore, for all $\varepsilon >0$  , $\mathbb{P}\left({\text{d}}_{\text{H}}\right(\partial {\mathbb{D}}_{-i,i}^{{\delta }_n}\left(z\right),\partial {\mathbb{D}}_{-i,i}\left(z\right))>\varepsilon )\to 0$  as $n\to \infty$  , which implies convergence in distribution of $\partial {\mathbb{D}}_{-i,i}^{{\delta }_n}\left(z\right)$  to $\partial {\mathbb{D}}_{-i,i}\left(z\right)$  , as ${\delta }_n\to 0$  , in the topology induced by ( 9 ). But Lemma  5.1 implies that $\partial {\mathbb{D}}_{-i,i}^{{\delta }_n}\left(z\right)$  converges in distribution (using ( 2 )) to a simple loop, therefore $\partial {\mathbb{D}}_{-i,i}\left(z\right)$  must also be a simple loop; and we have convergence in the topology induced by ( 2 ).

It is also clear that the argument above is independent of the subsequence $\left\{{\delta }_n\right\}$  , so the limit of $\partial {\mathbb{D}}_{-i,i}^{\delta }\left(z\right)$  is unique and coincides with $\partial {\mathbb{D}}_{-i,i}\left(z\right)$  . Hence, we have convergence in distribution of $\partial {\mathbb{D}}_{-i,i}^{\delta }\left(z\right)$  to $\partial {\mathbb{D}}_{-i,i}\left(z\right)$  , as $\delta \to 0$  , in the topology induced by ( 2 ), and indeed joint convergence of $({\gamma }^{\delta },\partial {\mathbb{D}}_{-i,i}^{\delta }(z\left)\right)$  to $(\gamma ,\partial {\mathbb{D}}_{-i,i}(z\left)\right)$  .

Proof of Corollary  5.1 . The corollary follows immediately from Lemma  5.1 and Lemma  5.2 , as already seen in the proof of Lemma  5.2 .

Proof of Lemma  5.3 . First of all recall that the convergence of $(\partial D_k,a_k,b_k)$  to $(\partial D,a,b)$  in distribution implies the existence of coupled versions of $(\partial D_k,a_k,b_k)$  and $(\partial D,a,b)$  onsome probability space $({\Omega }^{\prime },{\mathcal{ℬ}}^{\prime },{\mathbb{P}}^{\prime })$  such that $\text{d}(\partial D({\omega }^{\prime }),\partial D_k({\omega }^{\prime }\left)\right)\to 0$  , $a_k\left({\omega }^{\prime }\right)\to a\left({\omega }^{\prime }\right)$  , $b_k\left({\omega }^{\prime }\right)\to b\left({\omega }^{\prime }\right)$  for all ${\omega }^{\prime }$  as $k\to \infty$  (see, for example, Corollary 1 of [6] ). This immediately implies that the conditions to apply Radó's theorem (see Theorem  9 of Appendix  B ) are satisfied. Let $f_k$  be the conformal map that takes the unit disc $\mathbb{D}$  onto $D_k$  with $f_k\left(0\right)=0$  and $f_k^{\prime }\left(0\right)>0$  , and let $f$  be the conformal map from $\mathbb{D}$  onto $D$  with $f\left(0\right)=0$  and $f^{\prime }\left(0\right)>0$  . Then, by Theorem  9 , $f_k$  converges to $f$  uniformly in $\overline{\mathbb{D}}$  , as $k\to \infty$  .

Let $\gamma$  (resp., ${\gamma }_k$  ) be the chordal $SL{E}_6$  inside $D$  (resp., $D_k$  ) from $a$  to $b$  (resp., from $a_k$  to $b_k$  ), $\tilde{\gamma }=f^{-1}\left(\gamma \right)$  , $\tilde{a}=f^{-1}\left(a\right)$  , $\tilde{b}=f^{-1}\left(b\right)$  , and ${\tilde{\gamma }}_k=f_k^{-1}\left({\gamma }_k\right)$  , ${\tilde{a}}_k=f_k^{-1}\left(a_k\right)$  , ${\tilde{b}}_k=f_k^{-1}\left(b_k\right)$  .

We note that, because of the conformal invariance of chordal $SL{E}_6$  , $\tilde{\gamma }$  (resp., ${\tilde{\gamma }}_k$  ) is distributed as chordal $SL{E}_6$  in $\mathbb{D}$  from $\tilde{a}$  to $\tilde{b}$  (resp., from ${\tilde{a}}_k$  to ${\tilde{b}}_k$  ). Since $|a-a_k|\to 0$  and $|b-b_k|\to 0$  for all ${\omega }^{\prime }$  , and $f_k\to f$  uniformly in $\overline{\mathbb{D}}$  , we conclude that $|\tilde{a}-{\tilde{a}}_k|\to 0$  and $|\tilde{b}-{\tilde{b}}_k|\to 0$  for all ${\omega }^{\prime }$  .

Later we will prove a “continuity” property of $SL{E}_6$  (Lemma  6.1 ) that allows us to conclude that, under these conditions, ${\tilde{\gamma }}_k$  converges in distribution to $\tilde{\gamma }$  in the uniform metric ( 2 ) on continuous curves. Once again, this implies the existence of coupled versions of ${\tilde{\gamma }}_k$  and $\tilde{\gamma }$  on some probability space $({\Omega }^{\prime },{\mathcal{ℬ}}^{\prime },{\mathbb{P}}^{\prime })$  such that $\text{d}\left(\tilde{\gamma }\right({\omega }^{\prime }),{\tilde{\gamma }}_k({\omega }^{\prime }\left)\right)\to 0$  , for all ${\omega }^{\prime }$  as $k\to \infty$  . Therefore, thanks to the convergence of $f_k$  to $f$  uniformly in $\overline{\mathbb{D}}$  , $\text{d}\left(f\right(\tilde{\gamma }\left({\omega }^{\prime }\right)),f_k({\tilde{\gamma }}_k\left({\omega }^{\prime }\right)\left)\right)\to 0$  , for all ${\omega }^{\prime }$  as $k\to \infty$  . But since $f\left({\tilde{\gamma }}_k\right)$  is distributed as ${\gamma }_{D_k,a_k,b_k}$  and $f\left(\tilde{\gamma }\right)$  is distributed as ${\gamma }_{D,a,b}$  , we conclude that, as $k\to \infty$  , ${\gamma }_{D_k,a_k,b_k}$  converges in distribution to ${\gamma }_{D,a,b}$  in the uniform metric ( 2 ) on continuous curves.

We now note that (S) implies that, as $\delta \to 0$  , ${\gamma }_{D_k,a_k,b_k}^{\delta }$  converges in distribution to ${\gamma }_{D_k,a_k,b_k}$  uniformly in $k$  , for $k$  large enough. Therefore, as $k\to \infty$  , ${\gamma }_{D_k,a_k,b_k}^{{\delta }_k}$  converges in distribution to ${\gamma }_{D,a,b}$  , and the proof is concluded.

Lemma 6.1. Let $\mathbb{D}\subset \mathbb{C}$  be the unit disc, $a$  and $b$  two distinct points on its boundary, and $\gamma$  the trace of chordal $SL{E}_6$  inside $\mathbb{D}$  from $a$  to $b$  . Let $\left\{a_k\right\}$  and $\left\{b_k\right\}$  be two sequences of points in $\partial \mathbb{D}$  such that $a_k\to a$  and $b_k\to b$  . Then, as $k\to \infty$  , the trace ${\gamma }_k$  of chordal $SL{E}_6$  inside $\mathbb{D}$  from $a_k$  to $b_k$  converges in distribution to $\gamma$  in the uniform topology ( 2 ) on continuous curves.

Proof. Let $f_k\left(z\right)=e^{i{\alpha }_k}\frac{z-z_k}{1-{\overline{z}}_{k}z}$  be the (unique) linear fractional transformation that takes the unit disc $\mathbb{D}$  onto itself, mapping $a$  to $a_k$  , $b$  to $b_k$  , and a third point $c\in \partial \mathbb{D}$  distinct from $a$  and $b$  to itself. ${\alpha }_k$  and $z_k$  depend continuously on $a_k$  and $b_k$  . As $k\to \infty$  , since $a_k\to a$  and $b_k\to b$  , $f_k$  converges uniformly to the identity in $\overline{\mathbb{D}}$  .

Using the conformal invariance of chordal $SL{E}_6$  , we couple ${\gamma }_k$  and $\gamma$  by writing ${\gamma }_k=f_k\left(\gamma \right)$  . The uniform convergence of $f_k$  to the identity implies that $\text{d}(\gamma ,{\gamma }_k)\to 0$  as $k\to \infty$  , which is enough to conclude that ${\gamma }_k$  converges to $\gamma$  in distribution.

Proof of Theorem  1 . Let us prove the second part of the theorem first. We will do this for the original version of the discrete construction, but essentially the same proof works for the reorganized version we will describe below, as we will explain later. Suppose that at step $k$  of this discrete construction an exploration process ${\gamma }_k^{\delta }$  is run inside a domain $D_{k-1}^{\delta }$  , and write $D_{k-1}^{\delta }\backslash \Gamma \left({\gamma }_k^{\delta }\right)={\cup }_j{D}_{k,j}^{\delta }$  , where $\left\{D_{k,j}^{\delta }\right\}$  are the maximal connected domains of unexplored hexagons into which $D_{k-1}^{\delta }$  is split by removing the set $\Gamma \left({\gamma }_k^{\delta }\right)$  of hexagonsexplored by ${\gamma }_k^{\delta }$  .

Let ${\text{d}}_x\left(D_{k-1}^{\delta }\right)$  and ${\text{d}}_y\left(D_{k-1}^{\delta }\right)$  be respectively the maximal $x$  and $y$  -distances between pairs of points in $\partial D_{k-1}^{\delta }$  . Suppose, without loss of generality, that ${\text{d}}_x\left(D_{k-1}^{\delta }\right)\ge {\text{d}}_y\left(D_{k-1}^{\delta }\right)$  , and consider the rectangle $\mathcal{ℛ}$  (see Figure  6 ) whose vertical sides are aligned to the $y$  -axis, have length ${\text{d}}_x\left(D_{k-1}^{\delta }\right)$  , and are each placed at $x$  -distance $\frac{1}{3}{\text{d}}_x\left(D_{k-1}^{\delta }\right)$  from points of $\partial D_{k-1}^{\delta }$  with minimal or maximal $x$  -coordinate in such a way that the horizontal sides of $\mathcal{ℛ}$  have length $\frac{1}{3}{\text{d}}_x\left(D_{k-1}^{\delta }\right)$  ; the bottom and top sides of $\mathcal{ℛ}$  are placed in such a way that they are at equal $y$  -distance from the points of $\partial D_{k-1}^{\delta }$  with minimal or maximal $y$  -coordinate, respectively.

It follows from the Russo-Seymour-Welsh lemma [27, 29] (see also [15, 13] ) that the probability to have two vertical $\mathcal{T}$  -crossings of $\mathcal{ℛ}$  of different colors is bounded away from zero by a positive constant $p_0$  that does not depend on $\delta$  (for $\delta$  small enough). If that happens, then ${max}_j{\text{d}}_x\left(D_{k,j}^{\delta }\right)\le \frac{2}{3}{\text{d}}_x\left(D_{k-1}^{\delta }\right)$  . The same argument of course applies to the maximal $y$  -distance when ${\text{d}}_y\left(D_{k-1}^{\delta }\right)\ge {\text{d}}_x\left(D_{k-1}^{\delta }\right)$  . We can summarize the above observation in the following lemma.

Lemma 6.2. Suppose that at step $k$  of the full discrete construction an exploration process ${\gamma }_k^{\delta }$  is run inside a domain $D_{k-1}^{\delta }$  . If ${\text{d}}_x\left(D_{k-1}^{\delta }\right)\ge {\text{d}}_y\left(D_{k-1}^{\delta }\right)$  , then for $\delta$  small enough (i.e., $\delta \le C{\text{d}}_x\left(D_{k-1}^{\delta }\right)$  for some constant $C$  ), ${max}_j{\text{d}}_x\left(D_{k,j}^{\delta }\right)\le \frac{2}{3}{\text{d}}_x\left(D_{k-1}^{\delta }\right)$  with probability at least $p_0$  independent of $\delta$  . The same holds for the maximal $y$  -distances when ${\text{d}}_y\left(D_{k-1}^{\delta }\right)\ge {\text{d}}_x\left(D_{k-1}^{\delta }\right)$  .

Here is another lemma that will be useful later on.

Lemma 6.3. Two “daughter” subdomains, $D_{k,j}^{\delta }$  and $D_{k,j^{\prime }}^{\delta }$  , either have disjoint s-boundaries, or else their common s-boundary consists of exactly two adjacent hexagons (of the same color) where the exploration path ${\gamma }_k^{\delta }$  came within $2$  hexagons of touching itself just when completing the s-boundary of one of the two subdomains.

Proof. Suppose that the two daughter subdomains have s-boundaries $\Delta D_{k,j}^{\delta }$  and $\Delta D_{k,j^{\prime }}^{\delta }$  that are not disjoint and let $S=\{{\xi }_1,\dots ,{\xi }_i\}$  be the set of (sites of $\mathcal{T}$  that are the centers of the) hexagons that belong to both s-boundaries. $S$  can be partitioned into subsets consisting of single hexagons that are not adjacent to any another hexagon in $S$  and groups of hexagons that form simple $\mathcal{T}$  -paths (because the s-boundaries of the two subdomains are simple $\mathcal{T}$  -loops). Let $\{{\xi }_l,\dots ,{\xi }_m\}$  be such a subset of hexagons of $S$  that form a simple $\mathcal{T}$  -path ${\pi }_0=({\xi }_l,\dots ,{\xi }_m)$  . Then there is a $\mathcal{T}$  -path ${\pi }_1$  of hexagons in $\Delta D_{k,j}^{\delta }$  that goes from ${\xi }_l$  to ${\xi }_m$  without using any other hexagon of ${\pi }_0$  and a different $\mathcal{T}$  -path ${\pi }_2$  in $\Delta D_{k,j^{\prime }}^{\delta }$  that goes from ${\xi }_m$  to ${\xi }_l$  without using any other hexagon of ${\pi }_0$  . But then, all the hexagons in ${\pi }_0$  other than ${\xi }_l$  and ${\xi }_m$  are “surrounded” by ${\pi }_1\cup {\pi }_2$  and therefore cannot have been explored by the exploration process that produced $D_{k,j}^{\delta }$  and $D_{k,j^{\prime }}^{\delta }$  , andcannot belong to $\Delta D_{k,j}^{\delta }$  or $\Delta D_{k,j^{\prime }}^{\delta }$  , leading to a contradiction, unless ${\pi }_0=({\xi }_l,{\xi }_m)$  . Similar arguments lead to a contradiction if $S$  is partitioned into more than one subset.

If ${\xi }_i\in S$  is not adjacent to any other hexagon in $S$  , then it is adjacent to two other hexagons of $\Delta D_{k,j}^{\delta }$  and two hexagons of $\Delta D_{k,j^{\prime }}^{\delta }$  . Since ${\xi }_i$  has only six neighbors and neither the two hexagons of $\Delta D_{k,j}^{\delta }$  adjacent to ${\xi }_i$  nor those of $\Delta D_{k,j^{\prime }}^{\delta }$  can be adjacent to each other, each hexagon of $\Delta D_{k,j}^{\delta }$  is adjacent to one of $\Delta D_{k,j^{\prime }}^{\delta }$  . But then, as before, ${\xi }_i$  is “surrounded” by $\{\Delta D_{k,j}^{\delta }\cup \Delta D_{k,j^{\prime }}^{\delta }\}\backslash {\xi }_i$  and therefore cannot have been explored by the exploration process that produced $D_{k,j}^{\delta }$  and $D_{k,j^{\prime }}^{\delta }$  , and cannot belong to $\Delta D_{k,j}^{\delta }$  or $\Delta D_{k,j^{\prime }}^{\delta }$  , leading once again to a contradiction. The proof is now complete, since the only case remaining is the one where $S$  consists of a single pair of adjacent hexagons as stated in the lemma.

With these lemmas, we can now proceed with the proof of the second part of the theorem. Lemma  6.2 tells us that large domains are “chopped” with bounded away from zero probability ( $\ge p_0>0$  ), but we need to keep track of domains of diameter larger than $\varepsilon$  in such a way as to avoid “double counting” as the lattice construction proceeds.

More accurately, we will keep track of domains ${\tilde{D}}^{\delta }$  having ${\text{d}}_m\left({\tilde{D}}^{\delta }\right)\ge \frac{1}{\sqrt{2}}\varepsilon$  , since only these can have diameter larger than $\varepsilon$  . To do so, we will associate with each domain ${\tilde{D}}^{\delta }$  having ${\text{d}}_m\left({\tilde{D}}^{\delta }\right)\ge \frac{1}{\sqrt{2}}\varepsilon$  that we encounter as we do the lattice construction a non-negative integer label. The first domain is $D_0^{\delta }={\mathbb{D}}^{\delta }$  (see the beginning of Section  4.3 ) and this gets label $1$  . After each exploration process in a domain ${\tilde{D}}^{\delta }$  with ${\text{d}}_m\left({\tilde{D}}^{\delta }\right)\ge \frac{1}{\sqrt{2}}\varepsilon$  , if the number $\tilde{m}$  of “daughter” subdomains ${\tilde{D}}_j^{\delta }$  with ${\text{d}}_m\left({\tilde{D}}_j^{\delta }\right)\ge \frac{1}{\sqrt{2}}\varepsilon$  is $0$  , then the label of ${\tilde{D}}^{\delta }$  is no longer used, if instead $\tilde{m}\ge 1$  , then one of these $\tilde{m}$  subdomains (chosen by any procedure – e.g., the one with the highest priority for further exploration) is assigned the same label as ${\tilde{D}}^{\delta }$  and the rest are assigned the next $\tilde{m}-1$  integers that have never before been used as labels. Note that once all domains have ${\text{d}}_m<\frac{1}{\sqrt{2}}\varepsilon$  , there are no more labelled domains.

Lemma 6.4. Let $M_{\varepsilon }^{\delta }$  denote the total number of labels used in the above procedure; then for any fixed $\varepsilon >0$  , $M_{\varepsilon }^{\delta }$  is bounded in probability as $\delta \to 0$  ; i.e., ${lim}_{M\to \infty }{limsup}_{\delta \to 0}\mathbb{P}(M_{\varepsilon }^{\delta }>M)=0$  .

Proof. Except for $D_0^{\delta }$  , every domain comes with (at least) a “physically correct” monochromatic “half-boundary” (notice that we are considering s-boundaries and that a half-boundary coming from the “artificially colored” boundary of $D_0^{\delta }$  is not considered a physically correct monochromatic half-boundary). Let us assume, without loss of generality, that $M_{\varepsilon }^{\delta }>1$  .

If we associate with each label the “last” (in terms of steps of the discrete construction) domain which used that label (its daughter subdomains all had ${\text{d}}_m<\frac{1}{\sqrt{2}}\varepsilon$  ), then we claim that it follows from Lemma  6.3 that (with high probability) any two such last domains that are labelled have disjoint s-boundaries. This is a consequence of the fact that the two domains are subdomains of two “ancestors” that are distinct daughter subdomains of the same domain (possibly $D_0^{\delta }$  ) and whose s-boundaries are therefore (by Lemma  6.3 ) either disjoint or else overlap at a pair of hexagons where an exploration path had a close encounter of distance two hexagons with itself. But since we are dealing only with macroscopic domains (of diameter at least order $\varepsilon$  ), such a close encounter would imply,like in Lemmas  5.1 and  5.2 , the existence of six crossings, not all of the same color, of an annulus whose outer radius can be kept fixed while the inner radius is sent to zero together with $\delta$  . The probability of such an event goes to zero as $\delta \to 0$  and hence the unit disc $\mathbb{D}$  contains, with high probability, at least $M_{\varepsilon }^{\delta }$  disjoint monochromatic $\mathcal{T}$  -paths of diameter at least $\frac{1}{\sqrt{2}}\varepsilon$  , corresponding to the physically correct half-boundaries of the $M_{\varepsilon }^{\delta }$  labelled domains.

Now take the collection of squares $s_j$  of side length ${\varepsilon }^{\prime }>0$  centered at the sites $c_j$  of a scaled square lattice ${\varepsilon }^{\prime }{\mathbb{Z}}^2$  of mesh size ${\varepsilon }^{\prime }$  , and let $N\left({\varepsilon }^{\prime }\right)$  be the number of squares of side ${\varepsilon }^{\prime }$  needed to cover the unit disc. Let ${\varepsilon }^{\prime }<\varepsilon /2$  and consider the event $\{M_{\varepsilon }^{\delta }\ge 6N({\varepsilon }^{\prime }\left)\right\}$  , which implies that, with high probability, the unit disc contains at least $6N\left({\varepsilon }^{\prime }\right)$  disjoint monochromatic $\mathcal{T}$  -paths of diameter at least $\frac{1}{\sqrt{2}}\varepsilon$  and that, for at least one $j=j_0$  , the square $s_{j_0}$  intersects at least six disjoint monochromatic $\mathcal{T}$  -paths of diameter larger that $\frac{1}{\sqrt{2}}\varepsilon$  , so that the “annulus” $B(c_{j_0},\frac{1}{2\sqrt{2}}\varepsilon )\backslash s_{j_0}$  is crossed by at least six disjoint monochromatic $\mathcal{T}$  -paths contained inside the unit disc.

If all these $\mathcal{T}$  -paths crossing $B(c_{j_0},\frac{1}{2\sqrt{2}}\varepsilon )\backslash s_{j_0}$  have the same color, say blue, then since they are portions of boundaries of domains discovered by exploration processes, they are “shadowed” by exploration paths and therefore between at least one pair of blue $\mathcal{T}$  -paths, there is at least one yellow $\mathcal{T}$  -path crossing $B(c_{j_0},\frac{1}{2\sqrt{2}}\varepsilon )\backslash s_{j_0}$  . Therefore, whether the original monochromatic $\mathcal{T}$  -paths are all of the same color or not, $B(c_{j_0},\frac{1}{2\sqrt{2}}\varepsilon )\backslash s_{j_0}$  is crossed by at least six disjoint monochromatic $\mathcal{T}$  -paths not all of the same color contained in the unit disc. Let $g(\varepsilon ,{\varepsilon }^{\prime })$  denote the $limsup$  as $\delta \to 0$  of the probability that such an event happens anywhere inside the unit disc. We have shown that the event $\{M_{\varepsilon }^{\delta }\ge 6N({\varepsilon }^{\prime }\left)\right\}$  implies a “six-arms” event unless not all labelled domains have disjoint s-boundaries. But the latter also implies a “six-arms” event, as discussed before; therefore

$$\begin{array}{c}{limsup}_{\delta \to 0}\mathbb{P}(M_{\varepsilon }^{\delta }\ge 6N({\varepsilon }^{\prime }\left)\right)\le 2g(\varepsilon ,{\varepsilon }^{\prime }).\end{array}$$

(10)

Since $B(c_{j_0},\frac{1}{2\sqrt{2}}\varepsilon )\backslash B(c_{j_0},\frac{1}{\sqrt{2}}{\varepsilon }^{\prime })\subset B(c_{j_0},\frac{1}{2\sqrt{2}}\varepsilon )\backslash s_{j_0}$  , bounds in [16] imply that, for $\varepsilon$  fixed, $g(\varepsilon ,{\varepsilon }^{\prime })\to 0$  as ${\varepsilon }^{\prime }\to 0$  , which shows that

$$\begin{array}{c}{lim}_{M\to \infty }{limsup}_{\delta \to 0}\mathbb{P}(M_{\varepsilon }^{\delta }>M)=0\end{array}$$

(11)

and concludes the proof of the lemma.

Now, let $N_i^{\delta }$  denote the number of distinct domains that had label $i$  (this is equal to the number of steps that label $i$  survived). Let us also define $H\left(\varepsilon \right)$  to be the smallest integer $h\ge 1$  such that $(\frac{2}{3}{)}^h<\frac{1}{\sqrt{2}}\varepsilon$  and $G_h$  to be the random variable corresponding to how many Bernoulli trials (with probability $p_0$  of success) it takes to have $h$  successes.

Then, we may apply (sequentially) Lemma  6.2 to conclude that for any $i$ 

$$\begin{array}{c}\mathbb{P}(N_i^{\delta }\ge k+1)\le \mathbb{P}(G_{H\left(\varepsilon \right)}+G_{H\left(\varepsilon \right)}^{\prime }\ge k),\end{array}$$

(12)

where $G_h^{\prime }$  is an independent copy of $G_h$  . Now let ${\tilde{N}}_1\left(\varepsilon \right),{\tilde{N}}_2\left(\varepsilon \right),\dots$  be i.i.d. random variables equidistributed with $G_{H\left(\varepsilon \right)}+G_{H\left(\varepsilon \right)}^{\prime }$  .

Let ${\tilde{K}}_{\delta }\left(\varepsilon \right)$  be the number of steps needed so that all domains left to explore have ${\text{d}}_m<\frac{1}{\sqrt{2}}\varepsilon$  .

Then, for any positive integer $M$  ,

$$\begin{array}{c}\mathbb{P}\left({\tilde{K}}_{\delta }\right(\varepsilon )>C)\le \mathbb{P}(M_{\varepsilon }^{\delta }\ge M+1)+\mathbb{P}\left({\tilde{N}}_1\right(\varepsilon )+\dots +{\tilde{N}}_M(\varepsilon )\ge C).\end{array}$$

(13)

Notice that, for fixed $M$  , $\mathbb{P}\left({\tilde{N}}_1\right(\varepsilon )+\dots +{\tilde{N}}_M(\varepsilon )\ge C)\to 0$  as $C\to \infty$  . Moreover, for any $\hat{\varepsilon }>0$  , by Lemma  6.4 , we can choose $M_0=M_0\left(\hat{\varepsilon }\right)$  large enough so that ${limsup}_{\delta \to 0}\mathbb{P}(M_{\varepsilon }^{\delta }>M_0)<\hat{\varepsilon }$  . So, for any $\hat{\varepsilon }>0$  , it follows that

$$\begin{array}{c}{limsup}_{C\to \infty }{limsup}_{\delta \to 0}\mathbb{P}\left({\tilde{K}}_{\delta }\right(\varepsilon )>C)<\hat{\varepsilon },\end{array}$$

(14)

which implies that

$$\begin{array}{c}{lim}_{C\to \infty }{limsup}_{\delta \to 0}\mathbb{P}\left({\tilde{K}}_{\delta }\right(\varepsilon )>C)=0.\end{array}$$

(15)

To conclude this part of the proof, notice that the discrete construction cannot “skip” a contour and move on to explore its interior, so that all the contours with diameter larger than $\varepsilon$  must have been found by step $k$  if all the domains present at that step have diameter smaller than $\varepsilon$  . Therefore, $K_{\delta }\left(\varepsilon \right)\le {\tilde{K}}_{\delta }\left(\varepsilon \right)$  , which shows that $K_{\delta }\left(\varepsilon \right)$  is bounded in probability as $\delta \to 0$  .

For the first part of the theorem, we need to prove, for any fixed $k\in \mathbb{N}$  , joint convergence in distribution of the first $k$  steps of a suitably reorganized discrete construction to the first $k$  steps of the continuum one. Later we will explain why this reorganized construction has the same scaling limit as the one defined in Section  4.3 . For each $k$  , the first $k$  steps of the reorganized discrete construction will be coupled to the first $k$  steps of the continuum one with suitable couplings in order to obtain the convergence in distribution of those steps of the discrete construction to the analogous steps of the continuum one; the proof will proceed by induction in $k$  . We will explain how to reorganize the discrete construction as we go along; in order to explain the idea of the proof, we will consider first the cases $k=1$  , $2$  and $3$  , and then extend to all $k>3$  .

$k=1$  . The first step of the continuum construction consists of an $SL{E}_6$  ${\gamma }_1$  from $-i$  to $i$  inside $\mathbb{D}$  . Correspondingly, the first step of the discrete construction consists of an exploration path ${\gamma }_1^{\delta }$  inside ${\mathbb{D}}^{\delta }$  from the e-vertex closest to $-i$  to the e-vertex closest to $i$  .

The convergence in distribution of ${\gamma }_1^{\delta }$  to ${\gamma }_1$  is covered by statement (S).

$k=2$  . The convergence in distribution of the percolation exploration path to chordal $SL{E}_6$  implies that we can couple ${\gamma }_1^{\delta }$  and ${\gamma }_1$  generating them as random variables on some probability space $({\Omega }^{\prime },{\mathcal{ℬ}}^{\prime },{\mathbb{P}}^{\prime })$  such that $\text{d}\left({\gamma }_1\right({\omega }^{\prime }),{\gamma }_1^{\delta }({\omega }^{\prime }\left)\right)\to 0$  for all ${\omega }^{\prime }$  as $k\to \infty$  (see, for example, Corollary 1 of [6] ).

Now, let $D_1$  be the domain generated by ${\gamma }_1$  that is chosen for the second step of the continuum construction, and let $c_1\in \mathcal{P}$  be the highest ranking point of $\mathcal{P}$  contained in $D_1$  .For $\delta$  small enough, $c_1$  is also contained in ${\mathbb{D}}^{\delta }$  ; let $D_1^{\delta }=D_1^{\delta }\left(c_1\right)$  be the unique connected component of the set ${\mathbb{D}}^{\delta }\backslash \Gamma \left({\gamma }_1^{\delta }\right)$  containing $c_1$  (this is well-defined with probability close to $1$  for small $\delta$  ); $D_1^{\delta }$  is the domain where the second exploration process is to be carried out.

From the proof of Lemma  5.2 , we know that the boundaries $\partial D_1^{\delta }$  and $\partial D_1$  of the domains $D_1^{\delta }$  and $D_1$  produced respectively by the path ${\gamma }_1^{\delta }$  and ${\gamma }_1$  are close with probability close to one for $\delta$  small enough.

For the next step of the discrete construction, we choose the two e-vertices $x_1$  and $y_1$  in $\partial D_1^{\delta }$  that are closest to the points $a_1$  and $b_1$  of $\partial D_1$  selected for the coupled continuum construction (if the choice is not unique, we can select the e-vertices with any rule to break the tie) and call ${\gamma }_2^{\delta }$  the percolation exploration path inside $D_1^{\delta }$  from $x_1$  to $y_1$  . It follows from [2] that $\{{\gamma }_1^{\delta },\partial D_1^{\delta },{\gamma }_2^{\delta }\}$  converge jointly in distribution along some subsequence to some limit $\{{\tilde{\gamma }}_1,\partial {\tilde{D}}_1,{\tilde{\gamma }}_2\}$  . We already know that ${\tilde{\gamma }}_1$  is distributed like ${\gamma }_1$  and we can deduce from the joint convergence in distribution of $({\gamma }_1^{\delta },\partial D_1^{\delta })$  to $({\gamma }_1,\partial D_1)$  (Lemma  5.2 ), that $\partial {\tilde{D}}_1$  is distributed like $\partial D_1$  . Therefore, if we call ${\gamma }_2$  the $SL{E}_6$  path inside $D_1$  from $a_1$  to $b_1$  , Lemma  5.3 implies that ${\tilde{\gamma }}_2$  is distributed like ${\gamma }_2$  and indeed that, as $\delta \to 0$  , $\{{\gamma }_1^{\delta },\partial D_1^{\delta },{\gamma }_2^{\delta }\}$  converge jointly in distribution to $\{{\gamma }_1,\partial D_1,{\gamma }_2\}$  .

$k=3$  . So far, we have proved the convergence in distribution of the (paths and boundaries produced in the) first two steps of the discrete construction to the (paths and boundaries produced in the) first two steps of the discrete construction. The third step of the continuum construction consists of an $SL{E}_6$  path ${\gamma }_3$  from $a_2\in \partial D_2$  to $b_2\in \partial D_2$  , inside the domain $D_2$  with highest priority after the second step has been completed. Let $c_2\in \mathcal{P}$  be the highest ranking point of $\mathcal{P}$  contained in $D_2$  , $D_2^{\delta }$  the domain of the discrete construction containing $c_2$  after the second step of the discrete construction has been completed (this is well defined with probability close to $1$  for small $\delta$  ), and choose the two e-vertices $x_2$  and $y_2$  in $\partial D_2^{\delta }$  that are closest to the points $a_2$  and $b_2$  of $\partial D_2$  selected for the coupled continuum construction (if the choice is not unique, we can select the e-vertices with any rule to break the tie). The third step of the discrete construction consists of an exploration path ${\gamma }_3^{\delta }$  from $x_2$  to $y_2$  inside $D_2^{\delta }$  .

It follows from [2] that $\{{\gamma }_1^{\delta },\partial D_1^{\delta },{\gamma }_2^{\delta },\partial D_2^{\delta },{\gamma }_3^{\delta }\}$  converge jointly in distribution along some subsequence to some limit $\{{\tilde{\gamma }}_1,\partial {\tilde{D}}_1,{\tilde{\gamma }}_2,\partial {\tilde{D}}_2,{\tilde{\gamma }}_3\}$  . We already know that ${\tilde{\gamma }}_1$  is distributed like ${\gamma }_1$  , $\partial {\tilde{D}}_1$  like $\partial D_1$  and ${\tilde{\gamma }}_2$  like ${\gamma }_2$  , and we would like to apply Lemma  5.3 to conclude that ${\tilde{\gamma }}_3$  is distributed like ${\gamma }_3$  and indeed that, as $\delta \to 0$  , $({\gamma }_1^{\delta },\partial D_1^{\delta },{\gamma }_2^{\delta },\partial D_2^{\delta },{\gamma }_3^{\delta })$  converges in distribution to $({\gamma }_1,\partial D_1,{\gamma }_2,\partial D_2,{\gamma }_3)$  . In order to do so, we have to first show that $\partial {\tilde{D}}_2$  is distributed like $\partial D_2$  . If $D_2^{\delta }$  is a subset of ${\mathbb{D}}^{\delta }\backslash \Gamma \left({\gamma }_1^{\delta }\right)$  , this follows from Lemma  5.2 , as in the previous case, but if the s-boundary of $D_2^{\delta }$  contains hexagons of $\Gamma \left({\gamma }_2^{\delta }\right)$  , then we cannot use Lemma  5.2 directly, although the proof of the lemma can be easily adapted to the present case, as we now explain.

Indeed, the only difference is in the proof of claim (C) and is due to the fact that, when dealing with a domain of type 1 or 2, we cannot use the bound on the probability of three disjoint crossings of a semi-annulus because the domains we are dealing with may not be convex (like the unit disc). On the other hand, the discrete domains like $D_1^{\delta }$  and $D_2^{\delta }$  where we have to run exploration processes at various steps of the discreteconstruction are themselves generated by previous exploration processes, so that any hexagon of the s-boundary of such a domain has three adjacent hexagons which are the starting points of three disjoint $\mathcal{T}$  -paths (two of one color and one of the other). Two of these $\mathcal{T}$  -paths belong to the s-boundary of the domain, while the third belongs to the adjacent percolation cluster (see Figure  7 ). This allows us to use the bound on the probability of six disjoint crossings of an annulus.

To see this, let ${\pi }_1,{\pi }_2$  be the $\mathcal{T}$  -paths contained in the s-boundary of the discrete domain (i.e., $D_1^{\delta }$  in the present context) and ${\pi }_3$  the $\mathcal{T}$  -path belonging to the adjacent cluster, all starting from hexagons adjacent to some hexagon $\xi$  (centered at $u$  ) in the s-boundary of $D_1^{\delta }$  . For $0<{\varepsilon }^{\prime }<\varepsilon$  and $\delta$  small enough, let ${\mathcal{A}}_u(\varepsilon ,{\varepsilon }^{\prime })$  be the event that the exploration path ${\gamma }_2^{\delta }$  enters the ball $B(u,{\varepsilon }^{\prime })$  without touching $\partial D_1^{\delta }$  inside the larger ball $B(u,\varepsilon )$  . ${\mathcal{A}}_u(\varepsilon ,{\varepsilon }^{\prime })$  implies having (at least) three disjoint $\mathcal{T}$  -paths (two of one color and one of the other), ${\pi }_4,{\pi }_5$  and ${\pi }_6$  , contained in $D_1^{\delta }$  and crossing the annulus $B(u,\varepsilon )\backslash B(u,{\varepsilon }^{\prime })$  , with ${\pi }_4,{\pi }_5$  and ${\pi }_6$  disjoint from ${\pi }_1,{\pi }_2$  and ${\pi }_3$  . Hence, ${\mathcal{A}}_u(\varepsilon ,{\varepsilon }^{\prime })$  implies the event that there are (at least) six disjoint crossings (not all of the same color) of the annulus $B(u,\varepsilon )\backslash B(u,{\varepsilon }^{\prime })$  .

Once claim (C) is proved, the rest of the proof of Lemma  5.2 applies to the present case. Therefore, we have convergence in distribution of $\partial D_2^{\delta }$  to $\partial D_2$  , which allows us to use Lemma  5.3 and conclude that $({\gamma }_1^{\delta },{\gamma }_2^{\delta },{\gamma }_3^{\delta })$  converges in distribution to $({\gamma }_1,{\gamma }_2,{\gamma }_3)$  .

$k>3$  . We proceed by induction in $k$  , iterating the steps explained above; there are no new difficulties; all steps for $k\ge 4$  are analogous to the case $k=3$  .

To conclude the proof of the theorem, we need to show that the scaling limit of the original full discrete construction defined in Section  4.3 is the same as that of the reorganized one just used in the proof of the first part of the theorem. In order to do so, we can couple the two constructions by using the same percolation configuration for both, so that the two constructions have at their disposal the same set of loops to discover. We proved above that the original discrete construction finds all the “macroscopic” loops, so we have to show that this is true also for the reorganized version of the discrete construction. This is what we will do next, using essentially the same arguments as those employed for the original discrete construction; we present these arguments for the sake of completeness since there are some changes.

Consider the reorganized discrete construction described above, where the starting and ending points of the exploration processes at each step are chosen to be close to those of the corresponding (coupled) continuum construction. Suppose that at step $k$  of this discrete construction an exploration process ${\gamma }_k^{\delta }$  is run inside a domain $D_{k-1}^{\delta }$  , and write $D_{k-1}^{\delta }\backslash \Gamma \left({\gamma }_k^{\delta }\right)={\cup }_j{D}_{k,j}^{\delta }$  , where $\left\{D_{k,j}^{\delta }\right\}$  are the connected domains into which $D_{k-1}^{\delta }$  is split by the set $\Gamma \left({\gamma }_k^{\delta }\right)$  of hexagons explored by ${\gamma }_k^{\delta }$  .

Let ${\text{d}}_x\left(D_{k-1}\right)$  (resp., ${\text{d}}_x\left(D_{k-1}^{\delta }\right)$  ) and ${\text{d}}_y\left(D_{k-1}\right)$  (resp., ${\text{d}}_y\left(D_{k-1}^{\delta }\right)$  ) be respectively the maximal $x$  and $y$  -distance between pairs of points in $\partial D_{k-1}$  (resp., $\partial D_{k-1}^{\delta }$  ). If ${\text{d}}_x\left(D_{k-1}^{\delta }\right)\ge {\text{d}}_y\left(D_{k-1}^{\delta }\right)$  and the e-vertices on $\partial D_{k-1}^{\delta }$  are chosen to be closest to two points of $\partial D_{k-1}$  with maximal $x$  -distance, then the same construction and argument spelled out earlier in the first part of the proof (corresponding to the second part of the theorem) show that ${max}_j{\text{d}}_x\left(D_{k,j}^{\delta }\right)\le \frac{2}{3}{\text{d}}_x\left(D_{k-1}^{\delta }\right)$  with bounded away from zero probability.

If the e-vertices on $\partial D_{k-1}^{\delta }$  are chosen to be closest to two points of $\partial D_{k-1}$  with maximal $x$  -distance but ${\text{d}}_x\left(D_{k-1}^{\delta }\right)\le {\text{d}}_y\left(D_{k-1}^{\delta }\right)$  , then consider the rectangle ${\mathcal{ℛ}}^{\prime }$  whose vertical sides are aligned to the $y$  -axis, have length ${\text{d}}_y\left(D_{k-1}^{\delta }\right)$  , and are each placed at the same $x$  -distance from the points of $\partial D_{k-1}^{\delta }$  with minimal or maximal $x$  -coordinate in such a way that the horizontal sides of ${\mathcal{ℛ}}^{\prime }$  have length $\frac{1}{3}{\text{d}}_y\left(D_{k-1}^{\delta }\right)$  ; the bottom and top sides of ${\mathcal{ℛ}}^{\prime }$  are placed in such a way that they touch the points of $\partial D_{k-1}^{\delta }$  with minimal or maximal $y$  -coordinate, respectively. Notice that, because of the coupling between the continuum and discrete constructions, for any $\tilde{\varepsilon }>0$  , for $k$  large enough, $\left|{\text{d}}_x\right(D_{k-1}^{\delta })-{\text{d}}_x(D_{k-1}\left)\right|\le \tilde{\varepsilon }$  and $\left|{\text{d}}_y\right(D_{k-1}^{\delta })-{\text{d}}_y(D_{k-1}\left)\right|\le \tilde{\varepsilon }$  . Since in the case under consideration we have ${\text{d}}_y\left(D_{k-1}^{\delta }\right)\ge {\text{d}}_x\left(D_{k-1}^{\delta }\right)$  and ${\text{d}}_x\left(D_{k-1}\right)\ge {\text{d}}_y\left(D_{k-1}\right)$  , for $\delta$  large enough, we must also have $\left|{\text{d}}_y\right(D_{k-1}^{\delta })-{\text{d}}_x(D_{k-1}^{\delta }\left)\right|\le 2\tilde{\varepsilon }$  .

Once again, it follows from the Russo-Seymour-Welsh lemma that the probability to have two vertical $\mathcal{T}$  -crossings of ${\mathcal{ℛ}}^{\prime }$  of different colors is bounded away from zero by a positive constant that does not depend on $\delta$  (for $\delta$  small enough). If that happens, then ${max}_j{\text{d}}_x\left(D_{k,j}^{\delta }\right)\le \frac{2}{3}{\text{d}}_x\left(D_{k-1}^{\delta }\right)+\frac{1}{3}\tilde{\varepsilon }$  .

All other cases are handled in the same way, implying that the maximal $x$  and $y$  -distances of domains that appear in the discrete construction have a positive probability (bounded away from zero) to decrease by (approximately) a factor $2/3$  at each step of the discrete construction in which an exploration process is run in that domain.

With this result at our disposal, the rest of the proof, that for any $\varepsilon >0$  the number of steps needed to find all the loops of diameter larger than $\varepsilon$  is bounded in probability as $\delta \to 0$  (which implies that all the “macroscopic” loops are discovered), proceeds exactly like for the original discrete construction.

Proof of Theorem  2 . First of all, we want to show that $P^D\equiv {\hat{I}}_D{P}_R$  does not depend on $R$  , provided $D$  is strictly contained in ${\mathbb{D}}_R$  and $\partial D\cap \partial {\mathbb{D}}_R=\varnothing$  . In order to do this, we assume that the above conditions are satisfied for the pair $D,R$  and show that ${\hat{I}}_D{P}_R={\hat{I}}_D{P}_{R^{\prime }}$  for all $R^{\prime }>R$  .

Take two copies of the scaled hexagonal lattice, $\delta \mathcal{\mathscr{H}}$  and $\delta {\mathcal{\mathscr{H}}}^{\prime }$  , their dual lattices $\delta \mathcal{T}$  and $\delta {\mathcal{T}}^{\prime }$  , and two percolation configurations, ${\sigma }_{{\mathbb{D}}_R}$  and ${\sigma }_{{\mathbb{D}}_{R^{\prime }}}^{\prime }$  , both with blue boundary conditions and coupled in such a way that ${\sigma }_{{\mathbb{D}}_R}={\sigma }_{{\mathbb{D}}_R}^{\prime }$  . The laws of the boundaries of $\sigma$  and ${\sigma }^{\prime }$  are also coupled, in such a way that the boundaries or portions of boundaries contained inside $D$  are identical for all small enough $\delta$  . Therefore, letting $\delta \to 0$  and using the convergence of the percolation boundaries inside ${\mathbb{D}}_R$  and ${\mathbb{D}}_{R^{\prime }}$  to the continuum nonsimple loop processes $P_R$  and $P_{R^{\prime }}$  respectively, we conclude that ${\hat{I}}_D{P}_R={\hat{I}}_D{P}_{R^{\prime }}$  .

From what we have just proved, it follows that the probability measures $P^{{\mathbb{D}}_R}$  on $({\Omega }_R,{\mathcal{ℬ}}_R)$  , for $R\in {\mathbb{R}}_{+}$  , satisfy the consistency conditions $P^{{\mathbb{D}}_{R_1}}={\hat{I}}_{{\mathbb{D}}_{R_1}}{P}^{{\mathbb{D}}_{R_2}}$  for all $R_1\le R_2$  .

Since ${\Omega }_R$  , $\Omega$  are complete separable metric spaces, the measurable spaces $({\Omega }_R,{\mathcal{ℬ}}_R)$  , $(\Omega ,\mathcal{ℬ})$  are standard Borel spaces and so we can apply Kolmogorov's extension theorem (see, for example, [12] ) and conclude that there exists a unique probability measure on $(\Omega ,\mathcal{ℬ})$  with $P^{{\mathbb{D}}_R}={\hat{I}}_{{\mathbb{D}}_R}P$  for all $R\in {\mathbb{R}}_{+}$  . It follows that, for $R^{\prime }>R$  and all $D$  strictly contained in ${\mathbb{D}}_R$  and such that $\partial D\cap \partial {\mathbb{D}}_R=\varnothing$  , ${\hat{I}}_D{P}_R=P^D={\hat{I}}_D{P}_{R^{\prime }}={\hat{I}}_D{\hat{I}}_{{\mathbb{D}}_R}{P}_{R^{\prime }}={\hat{I}}_D{P}^{{\mathbb{D}}_R}={\hat{I}}_D{\hat{I}}_{{\mathbb{D}}_R}P={\hat{I}}_{D}P$  , which concludes the proof.

Proof of Corollary  5.2 . The corollary is an immediate consequence of Theorems  1 and  2 , where the full scaling limit is intended in the topology induced by ( 5 ).

Proof of Theorem  3 . 1. The fact that the Continuum Nonsimple Loop process is a random collection of noncrossing continuous loops is a direct consequence of its definition.

The fact that the loops touch themselves is a consequence of their being constructed out of $SL{E}_6$  , while the fact that they touch each other follows from the observation that a chordal $SL{E}_6$  path ${\gamma }_{D,a,b}$  touches $\partial D$  with probability one. Therefore, each new loop in the continuum construction touches one or more previous ones (many times).

The nonexistence of triple points follows directly from Lemma 5 of [16] on the number of crossings of an annulus, combined with Corollary  5.2 , which allows to transport discrete results to the continuum case. In fact, a triple point would imply, for discrete percolation, at least six crossings (not all of the same color) of an annulus whose ratio of inner to outer radius goes to zero in the scaling limit, leading to a contradiction.

2. This follows from straightforward Russo-Seymour-Welsh type arguments for percolation (for more details, see, for example, Lemma 3 of [16] ), combined with Corollary  5.2 .

3. Combining Russo-Seymour-Welsh type arguments for percolation (see, for example, Lemma 3 of [16] ) with Corollary  5.2 , we know that $P$  -a.s. there exists a (random) $R^{*}=R^{*}\left(R\right)$  , with $R^{*}<\infty$  , such that ${\mathbb{D}}_R$  is surrounded by a continuum nonsimple loop contained in ${\mathbb{D}}_{R^{*}}$  . From (the proof of ) Theorem  2 , we also know that ${\hat{I}}_{{\mathbb{D}}_{R^{\prime \prime }}}P=P^{{\mathbb{D}}_{R^{\prime \prime }}}=I_{{\mathbb{D}}_{R^{\prime \prime }}}{P}_{R^{\prime }}$  for all $R^{\prime }>R^{\prime \prime }$  . This implies that by taking $R^{\prime }$  large enough and performing the continuum construction inside ${\mathbb{D}}_{R^{\prime }}$  , we have a positive probability of generating a loop $\lambda$  contained in the annulus ${\mathbb{D}}_{R^{\prime \prime }}\backslash {\mathbb{D}}_R$  , with $R^{\prime }>R^{\prime \prime }>R$  . If that is the case, all the loops contained inside ${\mathbb{D}}_R$  are connected, by construction, to the loop $\lambda$  surrounding ${\mathbb{D}}_R$  by a finite sequence (a “path”) of loops (remember that in the continuum construction each loop is generated by pasting together portions of $SL{E}_6$  paths inside domains whose boundaries are determined by previously formed loops or excursions). Therefore, any two loops contained inside ${\mathbb{D}}_R$  are connected to each other by a “path” of loops.

Using again the fact that ${\hat{I}}_{{\mathbb{D}}_{R^{\prime \prime }}}P=P^{{\mathbb{D}}_{R^{\prime \prime }}}=I_{{\mathbb{D}}_{R^{\prime \prime }}}{P}_{R^{\prime }}$  for all $R^{\prime }>R^{\prime \prime }$  , and letting first $R^{\prime }$  and then $R^{\prime \prime }$  go to $\infty$  , we see from the discussion above (with $R\to \infty$  as well) that any two loops are connected by a finite “path” of intermediate loops, $P$  -a.s.

4. In order to prove the claim, we will define a discrete construction inside $D^{\prime }$  coupled to the continuum construction inside $D$  , by means of the conformal map $f$  from $D$  to $D^{\prime }$  .

Roughly speaking, this new discrete construction for $D^{\prime }$  is one in which the $(x,y)$  pairs at each step are chosen to be closest to the $\left(f\right(a),f(b\left)\right)$  points in $D^{\prime }$  mapped from $D$  via $f$  ,where the pairs $(a,b)$  are those that appear at the corresponding steps of the continuum construction inside $D$  .

More precisely, let ${\gamma }_1$  be the first $SL{E}_6$  path in $D$  from $a_1$  to $b_1$  . Because of the conformal invariance of $SL{E}_6$  , the image $f\left({\gamma }_1\right)$  of ${\gamma }_1$  under $f$  is a path distributed as the trace of chordal $SL{E}_6$  in $D^{\prime }$  from $f\left(a_1\right)$  to $f\left(b_1\right)$  . Therefore, the exploration path ${\gamma }_1^{\delta }$  inside $D^{\prime }$  from $x_1$  to $y_1$  , chosen to be closest to $f\left(a_1\right)$  and $f\left(b_1\right)$  respectively, converges in distribution to $f\left({\gamma }_1\right)$  , as $\delta \to 0$  , which means that there exist a coupling between ${\gamma }_1^{\delta }$  and $f\left({\gamma }_1\right)$  such that the paths stay close for $\delta$  small.

We see already that one can use the same strategy as in the proof of the first part of Theorem  1 , and obtain a discrete construction whose exploration paths are coupled to the $SL{E}_6$  paths $f\left({\gamma }_k\right)$  that are the images of the paths ${\gamma }_k$  in $D$  . Then, for this discrete construction, the scaling limit of the exploration paths will be distributed as the images of the $SL{E}_6$  paths in $D$  .

In order to conclude the proof, we just have to show that the discrete construction inside $D^{\prime }$  defined above finds all the boundaries in a number of steps that is bounded in probability as $\delta \to 0$  (this is equivalent to the second part of Theorem  1 ). To do that, we use the second part of Theorem  1 , which implies that, for any fixed $\varepsilon >0$  and $C<\infty$  , the number of steps of the discrete construction in $D$  that are necessary to ensure that only domains with diameter less than $\varepsilon /C$  are present is bounded in probability as $\delta \to 0$  .

Since $f$  (can be extended to a function that) is continuous in the compact set $\overline{D}$  , $f$  is uniformly continuous and so we can now choose $C=C\left(f\right)<\infty$  such that any subdomain of $D$  of diameter at most $\varepsilon /C$  is mapped by $f$  to a subdomain of $D^{\prime }$  of diameter at most $\varepsilon$  . This, combined with the coupling between $SL{E}_6$  paths and exploration paths inside $D^{\prime }$  , assures that the number of steps necessary for the new discrete construction inside $D^{\prime }$  to find all the loops of diameter at least $\varepsilon$  is bounded in probability as $\delta \to 0$  .

Therefore, the scaling limit, as $\delta \to 0$  , of this new discrete construction for $D^{\prime }$  gives the measure $P_{D^{\prime }}$  . It follows by construction that $f*P_D=P_{D^{\prime }}$  , which concludes the proof. Appendix A: Convergence of the Percolation Exploration Path In this first appendix, we provide a detailed proof of statement (S). The existence of subsequential limits for the percolation exploration path, which follows from the work of Aizenman and Burchard [2] , means that the proof can be divided into two parts: first we will give a characterization of chordal $SL{E}_6$  in terms of two properties that determine it uniquely; then we will show that any subsequential scaling limit of the percolation exploration path satisfies these two properties.

The characterization part will follow from known properties of hulls and of $SL{E}_6$  (see [22] and [35] ). The second part will follow from an extension of Smirnov's result about the convergence of crossing probabilities to Cardy's formula [10, 11] (see Theorem  6 below) for sequences of Jordan domains $D_k$  , with the domain $D_k$  changing together with the mesh ${\delta }_k$  of the lattice, combined with the proof of a certain spatial Markov property forsubsequential limits of percolation exploration hulls (Theorem  7 ). We note that although Theorem  6 represents only a slight extension to Smirnov's result on convergence of crossing probabilities, this extension and its proof play a major role in the technically important Lemmas  A.2 ,  A.4 and  A.5 , which control the “close encounters” of exploration paths to domain boundaries. The proof of Theorem  6 is modelled after a simpler geometric argument involving only rectangles used in [8] .

Let $D^{\prime }$  be a bounded simply connected domain containing the origin whose boundary $\partial D^{\prime }$  is a continuous curve. Let $f:\mathbb{D}\to D^{\prime }$  be the (unique) conformal map from the unit disc to $D^{\prime }$  with $f\left(0\right)=0$  and $f^{\prime }\left(0\right)>0$  ; note that by Theorem  10 of Appendix  B , $f$  has a continuous extension to $\overline{\mathbb{D}}$  . Let $z_1,z_2,z_3,z_4$  be four points of $\partial D^{\prime }$  (or more accurately, four prime ends) in counterclockwise order – i.e., such that $z_j=f\left(w_j\right),j=1,2,3,4$  , with $w_1,\dots ,w_4$  in counterclockwise order. Also, let $\eta =\frac{(w_1-w_2)(w_3-w_4)}{(w_1-w_3)(w_2-w_4)}$  . Cardy's formula [10, 11] for the probability ${\Phi }_{D^{\prime }}(z_1,z_2;z_3,z_4)$  of a “crossing” inside $D^{\prime }$  from the counterclockwise arc $\overline{z_1{z}_2}$  to the counterclockwise arc $\overline{z_3{z}_4}$  is

$$\begin{array}{c}{\Phi }_{D^{\prime }}(z_1,z_2;z_3,z_4)=\frac{\Gamma (2/3)}{\Gamma (4/3)\Gamma (1/3)}{\eta }_2^{1/3}{F}_1(1/3,2/3;4/3;\eta ),\end{array}$$

(16)

where ${}_2{F}_1$  is a hypergeometric function.

For a given mesh $\delta >0$  , the probability of a crossing inside $D^{\prime }$  from the counterclockwise arc $\overline{z_1{z}_2}$  to the counterclockwise arc $\overline{z_3{z}_4}$  is the probability of the existence of a blue $\mathcal{T}$  -path contained in $(D^{\prime }{)}^{\delta }$  , the $\delta$  -approximation of $D^{\prime }$  (see Definition  4.1 above), that starts at a hexagon adjacent to one intersecting $\overline{z_1{z}_2}$  and ends at a hexagon adjacent to one intersecting $\overline{z_3{z}_4}$  . Smirnov proved the following major theorem, concerning the conjectured behavior [10, 11] of crossing probabilities in the scaling limit.

Theorem 4. (Smirnov [30] ) In the limit $\delta \to 0$  , the crossing probability becomes conformally invariant and converges to Cardy's formula ( 16 ).

Remark A.1. We actually only need Theorem  4 for Jordan domains, as can be seen by a careful reading of the proof of Theorem  6 . We note that Smirnov does not restrict attention to that case.

Let us now specify the objects that we are interested in. Suppose $D$  is a simply connected domain whose boundary $\partial D$  is a continuous curve, and $a,b$  are two distinct points in $\partial D$  (or more accurately, two distinct prime ends), and let ${\tilde{\mu }}_{D,a,b}$  be a probability measure on continuous (non-self-crossing) curves $\tilde{\gamma }={\tilde{\gamma }}_{D,a,b}:[0,\infty )\to \overline{D}$  with $\tilde{\gamma }\left(0\right)=a$  and $\tilde{\gamma }(\infty )\equiv {lim}_{t\to \infty }\tilde{\gamma }\left(t\right)=b$  (we remark that the use of $[0,\infty )$  instead of $[0,1]$  for the time parametrization is purely for convenience). Let $D_t\equiv D\backslash {\tilde{K}}_t$  denote the (unique) connected component of $D\backslash \tilde{\gamma }[0,t]$  whose closure contains $b$  , where ${\tilde{K}}_t$  , the filling of $\tilde{\gamma }[0,t]$  , is a closed connected subset of $\overline{D}$  . ${\tilde{K}}_t$  is called a hull if it satisfies the condition

$$\begin{array}{c}\overline{\stackrel{t}{\stackrel{~}{K}}\cap D}={\tilde{K}}_t.\end{array}$$

(17)

We will generally be interested in curves $\tilde{\gamma }$  such that ${\tilde{K}}_t$  is a hull for each $t$  , although we normally only consider ${\tilde{K}}_T$  at certain stopping times $T$  . (An example of such a curve that we are particularly interested in is the trace of chordal $SL{E}_6$  .) Let $C^{\prime }\subset D$  be a closed subset of $\overline{D}$  such that $a/\in C^{\prime }$  , $b\in C^{\prime }$  , and $D^{\prime }=D\backslash C^{\prime }$  is a bounded simply connected domain whose boundary contains the counterclockwise arc $\overline{cd}$  that does not belong to $\partial D$  (except for its endpoints $c$  and $d$  – see Figure  8 ).

Let $T^{\prime }=inf\{t:{\tilde{K}}_t\cap C^{\prime }\ne \varnothing \}$  be the first time that $\tilde{\gamma }\left(t\right)$  hits $C^{\prime }$  and assume that the filling ${\tilde{K}}_{T^{\prime }}$  of $\tilde{\gamma }[0,T^{\prime }]$  is a hull; we denote by ${\tilde{\nu }}_{D^{\prime },a,c,d}$  the distribution of ${\tilde{K}}_{T^{\prime }}$  . To explain what we mean by the distribution of a hull, consider the set $\tilde{\mathcal{A}}$  of closed subsets $\tilde{A}$  of $\overline{D^{\prime }}$  that do not contain $a$  and such that $\partial \tilde{A}\backslash \partial D^{\prime }$  is a simple (continuous) curve contained in $D^{\prime }$  starting at a point on $\partial D^{\prime }\cap D$  and ending at a point on $\partial D$  (see Figure  9 ).

Let $\mathcal{A}$  be the set of closed subsets of $\overline{D^{\prime }}$  of the form ${\tilde{A}}_1\cup {\tilde{A}}_2$  , where ${\tilde{A}}_1,{\tilde{A}}_2\in \tilde{\mathcal{A}}$  and ${\tilde{A}}_1\cap {\tilde{A}}_2=\varnothing$  .

For a given $C^{\prime }$  and corresponding $T^{\prime }$  , let $\mathcal{K}$  be the set whose elements are possible hulls at time $T^{\prime }$  ; we claim that the events $\{K\in \mathcal{K}:K\cap A=\varnothing \}$  , for $A\in \mathcal{A}$  , form a $\pi$  -system $\Pi$  (i.e., they are closed under finite intersections; we also include the empty set in $\Pi$  ), and we consider the $\sigma$  -algebra $\Sigma =\sigma (\Pi )$  generated by these events. To see that $\Pi$  is closed under pairwise intersections, notice that, if $A_1,A_2\in \mathcal{A}$  , then $\{K\in \mathcal{K}:K\cap A_1=\varnothing \}\cap \{K\in \mathcal{K}:K\cap A_2=\varnothing \}=\{K\in \mathcal{K}:K\cap \{A_1\cup A_2\}=\varnothing \}$  and $A_1\cup A_2\in \mathcal{A}$  (or else $\{K\in \mathcal{K}:K\cap \{A_1\cup A_2\}=\varnothing \}$  is empty). We are interested in probability spaces of the form $(\mathcal{K},\Sigma ,{\mathbb{P}}^{*})$  .

We say that the exit distribution of $\tilde{\gamma }\left(t\right)$  is determined by Cardy's formula if, for any $C^{\prime }$  and any counterclockwise arc $\overline{xy}$  of $\overline{cd}$  , the probability that $\tilde{\gamma }$  hits $C^{\prime }$  at time $T^{\prime }$  on $\overline{xy}$  is given by

$$\begin{array}{c}{\mathbb{P}}^{*}\left(\tilde{\gamma }\right(T^{\prime })\in \overline{xy})={\Phi }_{D^{\prime }}(a,c;x,d)-{\Phi }_{D^{\prime }}(a,c;y,d).\end{array}$$

(18)

It is easy to see that if the hitting distribution of $\tilde{\gamma }\left(t\right)$  is determined by Cardy's formula, then the probabilities of events in $\Pi$  are also determined by Cardy's formula in the following way. Let $A\in \mathcal{A}$  be the union of ${\tilde{A}}_1,{\tilde{A}}_2\in \tilde{\mathcal{A}}$  , with $\partial {\tilde{A}}_1\backslash \partial D^{\prime }$  given by a curve from $u_1\in \partial D^{\prime }\cap D$  to $v_1\in \partial D$  and $\partial {\tilde{A}}_2\backslash \partial D^{\prime }$  given by a curve from $u_2\in \partial D^{\prime }\cap D$  to $v_2\in \partial D$  ; then, assuming that $a$  , $v_1$  , $u_1$  , $u_2$  , $v_2$  are ordered counterclockwise around $\partial D^{\prime }$  ,

$$\begin{array}{c}{\mathbb{P}}^{*}({\tilde{K}}_{T^{\prime }}\cap A=\varnothing )={\Phi }_{D^{\prime }\backslash A}(a,v_1;u_1,v_2,)-{\Phi }_{D^{\prime }\backslash A}(a,v_1;u_2,v_2).\end{array}$$

(19)

Since $\Pi$  is a $\pi$  -system, the probabilities of the events in $\Pi$  determine uniquely the distribution of the hull in the sense described above. Therefore, if we let ${\gamma }_{D,a,b}$  denote the trace of chordal $SL{E}_6$  inside $D$  from $a$  to $b$  , $K_t$  its hull up to time $t$  , and $\tau =inf\{t:K_t\cap C^{\prime }\ne \varnothing \}$  the first time that ${\gamma }_{D,a,b}$  hits $C^{\prime }$  , we have the following simple but useful lemma.

Lemma A.1. With the notation introduced above, if the hitting distribution of ${\tilde{\gamma }}_{D,a,b}$  is determined by Cardy's formula and ${\tilde{K}}_{T^{\prime }}$  is a hull, then ${\tilde{K}}_{T^{\prime }}$  is distributed like the hull $K_{\tau }$  of chordal $SL{E}_6$  .

Proof. It is enough to note that the hitting distribution for chordal $SL{E}_6$  is determined by Cardy's formula [19] .

Now let ${\tilde{f}}_0$  be a conformal map from $D$  to the upper half-plane $\mathbb{H}$  such that ${\tilde{f}}_0\left(a\right)=0$  and ${lim}_{z\to b}{\tilde{f}}_0\left(z\right)=\infty$  (these two conditions determine ${\tilde{f}}_0$  only up to a multiplicative constant). For $\varepsilon >0$  fixed, let $C(u,\varepsilon )=\{z:|u-z|<\varepsilon \}\cap \mathbb{H}$  denote the semi-ball of radius $\varepsilon$  centered at $u$  on the real line and let ${\tilde{T}}_1={\tilde{T}}_1\left(\varepsilon \right)$  denote the first time $\tilde{\gamma }\left(t\right)$  hits $D\backslash {\tilde{G}}_1$  , where ${\tilde{G}}_1\equiv {\tilde{f}}_0^{-1}\left(C\right(0,\varepsilon \left)\right)$  . Define recursively ${\tilde{T}}_{j+1}$  as the first time $\tilde{\gamma }[{\tilde{T}}_j,\infty )$  hits ${\tilde{D}}_{{\tilde{T}}_j}\backslash {\tilde{G}}_{j+1}$  , where ${\tilde{D}}_{{\tilde{T}}_j}\equiv D\backslash {\tilde{K}}_{{\tilde{T}}_j}$  , ${\tilde{G}}_{j+1}\equiv {\tilde{f}}_{{\tilde{T}}_j}^{-1}\left(C\right(0,\varepsilon \left)\right)$  , and ${\tilde{f}}_{{\tilde{T}}_j}$  is a conformal map from ${\tilde{D}}_{{\tilde{T}}_j}$  to $\mathbb{H}$  that maps $\tilde{\gamma }\left({\tilde{T}}_j\right)$  to $0$  and $b$  to $\infty$  . We also define ${\tilde{\tau }}_{j+1}\equiv {\tilde{T}}_{j+1}-{\tilde{T}}_j$  , so that ${\tilde{T}}_j={\tilde{\tau }}_1+\dots +{\tilde{\tau }}_j$  . We note that, like ${\tilde{f}}_0$  , the conformal maps ${\tilde{f}}_{{\tilde{T}}_j}$  are only defined up to a multiplicative factor.

Notice that ${\tilde{G}}_{j+1}$  is a bounded simply connected domain chosen so that the conformal transformation which maps ${\tilde{D}}_{{\tilde{T}}_j}$  to $\mathbb{H}$  maps ${\tilde{G}}_{j+1}$  to the semi-ball $C(0,\varepsilon )$  centered at the origin on the real line. With these definitions, consider the (discrete-time) stochastic process ${\tilde{X}}_j\equiv ({\tilde{K}}_{{\tilde{T}}_j},\tilde{\gamma }({\tilde{T}}_j\left)\right)$  for $j=1,2,\dots$  ; we say that ${\tilde{K}}_t$  satisfies the spatial Markov property if each ${\tilde{K}}_{{\tilde{T}}_j}$  is a hull and ${\tilde{X}}_j$  for $j=1,2,\dots$  is a Markov chain (for any choice of the multitplicative factors for ${\tilde{f}}_0,{\tilde{f}}_{{\tilde{T}}_1},{\tilde{f}}_{{\tilde{T}}_2},\dots$  ). Notice that the trace of chordal $SL{E}_6$  satisfies the spatial Markov property, due to the conformal invariance and Markovian properties [28] of $SL{E}_6$  .

Next, we give the main characterization theorem.

Theorem 5. If the filling process ${\tilde{K}}_t$  of ${\tilde{\gamma }}_{D,a,b}$  satisfies the spatial Markov property and its hitting distribution is determined by Cardy's formula, then ${\tilde{\gamma }}_{D,a,b}$  is distributed like the trace ${\gamma }_{D,a,b}$  of chordal $SL{E}_6$  inside $D$  started at $a$  and aimed at $b$  .

Proof. Since the trace ${\gamma }_{D,a,b}$  of chordal $SL{E}_6$  in a bounded Jordan domain $D$  is defined (up to a linear time change) as $f\left(\gamma \right)$  , where $\gamma ={\gamma }_{\mathbb{H},0,\infty }$  is the trace of chordal $SL{E}_6$  in the upper half-plane started at $0$  and $f$  is any conformal map from the upper half-plane $\mathbb{H}$  to $D$  such that $f\left(0\right)=a$  and $f(\infty )=b$  , it is enough to show that $\hat{\gamma }=f^{-1}\left({\tilde{\gamma }}_{D,a,b}\right)$  is distributed like the trace of chordal $SL{E}_6$  in the upper half-plane. Let ${\hat{K}}_t$  denote the filling of $\hat{\gamma }\left(t\right)$  at time $t$  and let ${\hat{g}}_t\left(z\right)$  be the unique conformal transformation that maps $\mathbb{H}\backslash {\hat{K}}_t$  onto $\mathbb{H}$  with the following expansion at infinity:

$$\begin{array}{c}{\hat{g}}_t\left(z\right)=z+\frac{\hat{a}\left(t\right)}{z}+o\left(\frac{1}{z}\right).\end{array}$$

(20)

We choose to parametrize $\hat{\gamma }\left(t\right)$  so that $t=\frac{\hat{a}\left(t\right)}{2}$  (this is often called parametrization by capacity, $\hat{a}\left(t\right)$  being the capacity of the filling up to time $t$  ). We want to compare $\hat{\gamma }\left(t\right)$  with the trace $\gamma \left(t\right)$  of chordal $SL{E}_6$  in the upper half-plane parametrized in the same way (i.e., with $a\left(t\right)=2t$  ), so that, if $K_t$  denotes the filling of $\gamma$  at time $t$  , $\mathbb{H}\backslash K_t$  is mapped onto $\mathbb{H}$  by a conformal transformation with the following expansion at infinity:

$$\begin{array}{c}{g}_t\left(z\right)=z+\frac{2t}{z}+o\left(\frac{1}{z}\right).\end{array}$$

(21)

Our strategy will be to construct suitable polygonal approximations ${\hat{\gamma }}_{\varepsilon }$  and ${\gamma }_{\varepsilon }$  of $\hat{\gamma }$  and $\gamma$  which converge, as $\varepsilon \to 0$  , to the original curves (in the uniform metric on continuous curves ( 2 )), and show that ${\hat{\gamma }}_{\varepsilon }$  and ${\gamma }_{\varepsilon }$  have the same distribution. This implies that the distributions of $\hat{\gamma }$  and $\gamma$  must coincide.

Let us describe first the construction for ${\gamma }_{\varepsilon }\left(t\right)$  ; we use exactly the same construction for ${\hat{\gamma }}_{\varepsilon }\left(t\right)$  . We remark that the important features in the construction of the polygonal approximations are the spatial Markov property of the fillings and Cardy's formula, which are valid for both $\gamma$  and $\hat{\gamma }$  .

For $\varepsilon >0$  fixed, as above let $C(u,\varepsilon )=\{z:|u-z|<\varepsilon \}\cap \mathbb{H}$  denote the semi-ball of radius $\varepsilon$  centered at $u$  on the real line. Let $T_1=T_1\left(\varepsilon \right)$  denote the first time $\gamma \left(t\right)$  hits $\mathbb{H}\backslash G_1$  , where $G_1\equiv C(0,\varepsilon )$  , and define recursively $T_{j+1}$  as the first time $\gamma [T_j,\infty )$  hits ${\mathbb{H}}_{T_j}\backslash G_{j+1}$  , where ${\mathbb{H}}_{T_j}=\mathbb{H}\backslash K_{T_j}$  and $G_{j+1}\equiv g_{T_j}^{-1}\left(C\right(g_{T_j}\left(\gamma \right(T_j\left)\right),\varepsilon \left)\right)$  . Notice that $G_{j+1}$  is a bounded simply connected domain chosen so that the conformal transformation which maps ${\mathbb{H}}_{T_j}$  to $\mathbb{H}$  maps $G_{j+1}$  to the semi-ball $C\left(g_{T_j}\right(\gamma \left(T_j\right)),\varepsilon )$  centered at the point of the real line where the “tip” $\gamma \left(T_j\right)$  of the hull $K_{T_j}$  is mapped. The spatial Markov property and the conformal invariance of the hull of $SL{E}_6$  imply that if we write $T_j={\tau }_1+\dots +{\tau }_j$  , with ${\tau }_{j+1}\equiv T_{j+1}-T_j$  , the ${\tau }_j$  's are i.i.d. random variables, and also that the distribution of $K_{T_{j+1}}$  is the same as the distribution of $K_{T_j}\cup g_{T_j}^{-1}(K_{T_1}^{\prime }+g_{T_j}(\gamma \left(T_j\right)\left)\right)$  , where $K_{T_1}^{\prime }$  is equidistributed with $K_{T_1}$  , but also is independent of $K_{T_1}$  . The polygonal approximation ${\gamma }_{\varepsilon }$  is obtained by joining, for all $j$  , $\gamma \left(T_j\right)$  to $\gamma \left(T_{j+1}\right)$  with a straight segment, where the speed ${\gamma }_{\varepsilon }^{\prime }\left(t\right)$  is chosen to be constant.

Now let ${\hat{T}}_1={\hat{T}}_1\left(\varepsilon \right)$  denote the first time $\hat{\gamma }\left(t\right)$  hits $\mathbb{H}\backslash {\hat{G}}_1$  , where ${\hat{G}}_1\equiv C(0,\varepsilon )$  , and define recursively ${\hat{T}}_{j+1}$  as the first time $\hat{\gamma }[{\hat{T}}_j,\infty )$  hits ${\hat{\mathbb{H}}}_{{\hat{T}}_j}\backslash {\hat{G}}_{j+1}$  , where ${\hat{\mathbb{H}}}_{{\hat{T}}_j}\equiv \mathbb{H}\backslash {\hat{K}}_{{\hat{T}}_j}$  and ${\hat{G}}_{j+1}\equiv {\hat{g}}_{T_j}^{-1}\left(C\right({\hat{g}}_{{\hat{T}}_j}\left(\hat{\gamma }\right({\hat{T}}_j\left)\right),\varepsilon \left)\right)$  . We also define ${\hat{\tau }}_{j+1}\equiv {\hat{T}}_{j+1}-{\hat{T}}_j$  , so that ${\hat{T}}_j={\hat{\tau }}_1+\dots +{\hat{\tau }}_j$  .

Once again, ${\hat{G}}_{j+1}$  is a bounded simply connected domain chosen so that the conformal transformation which maps ${\hat{\mathbb{H}}}_{{\hat{T}}_j}$  to $\mathbb{H}$  maps ${\hat{G}}_{j+1}$  to the semi-ball $C\left({\hat{g}}_{{\hat{T}}_j}\right(\hat{\gamma }\left({\hat{T}}_j\right)),\varepsilon )$  centered at the point on the real line where the “tip” $\hat{\gamma }\left({\hat{T}}_j\right)$  of the hull ${\hat{K}}_{{\hat{T}}_j}$  is mapped. The polygonal approximation ${\hat{\gamma }}_{\varepsilon }$  is obtained by joining, for all $j$  , $\hat{\gamma }\left({\hat{T}}_j\right)$  to $\hat{\gamma }\left({\hat{T}}_{j+1}\right)$  with a straight segment, where the speed ${\hat{\gamma }}_{\varepsilon }^{\prime }\left(t\right)$  is chosen to be constant.

Consider the sequence of times ${\tilde{T}}_j$  defined in the natural way so that $\tilde{\gamma }\left({\tilde{T}}_j\right)=f\left(\hat{\gamma }\right({\hat{T}}_j\left)\right)$  and the (discrete-time) stochastic processes ${\hat{X}}_j\equiv ({\hat{K}}_{{\hat{T}}_j},\hat{\gamma }({\hat{T}}_j\left)\right)$  and ${\tilde{X}}_j\equiv ({\tilde{K}}_{{\tilde{T}}_j},\tilde{\gamma }({\tilde{T}}_j\left)\right)$  related by ${\hat{X}}_j=f^{-1}\left({\tilde{X}}_j\right)$  . If for $x\in \mathbb{R}$  we let $\theta \left[x\right]$  denote the translation that maps $x$  to $0$  and define the family of conformal maps ${\tilde{f}}_{{\tilde{T}}_j}=\theta \left[g_{{\hat{T}}_j}\right(\hat{\gamma }\left({\hat{T}}_j\right)\left)\right]\circ g_{{\hat{T}}_j}\circ f^{-1}$  from $D\backslash {\tilde{K}}_{{\tilde{T}}_j}$  to $\mathbb{H}$  , then ${\tilde{f}}_{{\tilde{T}}_j}$  sends $\tilde{\gamma }\left({\tilde{T}}_j\right)$  to $0$  and $b$  to $\infty$  , and ${\tilde{T}}_{j+1}$  is the first time $\tilde{\gamma }[{\tilde{T}}_j,\infty )$  hits ${\tilde{\mathbb{H}}}_{{\hat{T}}_j}\backslash {\tilde{G}}_{j+1}$  ,where ${\tilde{\mathbb{H}}}_{{\tilde{T}}_j}=\mathbb{H}\backslash {\tilde{K}}_{{\hat{T}}_j}$  and ${\tilde{G}}_{j+1}={\tilde{f}}_{{\tilde{T}}_j}^{-1}\left(C\right(0,\varepsilon \left)\right)$  . Therefore, $\left\{{\tilde{T}}_j\right\}$  is a sequence of stopping times like those used in the definition of the spatial Markov property and, thanks to the relation ${\hat{X}}_j=f^{-1}\left({\tilde{X}}_j\right)$  , the fact that ${\tilde{K}}_t$  satisfies the spatial Markov property implies that ${\hat{X}}_j$  is a Markov chain. We also note that the fact that the hitting distribution of $\tilde{\gamma }\left(t\right)$  is determined by Cardy's formula implies the same for the hitting distribution of $\hat{\gamma }\left(t\right)$  , thanks to the conformal invariance of Cardy's formula. We next use these properties to show that ${\hat{\gamma }}_{\varepsilon }$  is distributed like ${\gamma }_{\varepsilon }$  .

To do so, we first note that the conformal transformations $g_{T_j}$  and ${\hat{g}}_{{\hat{T}}_j}$  are random and that their distributions are functionals of the distributions of the hulls $K_{T_j}$  and ${\hat{K}}_{{\hat{T}}_j}$  , since there is a one-to-one correspondence between hulls and conformal maps (with the normalization we have chosen in ( 20 )–( 21 )). Therefore, since ${\hat{K}}_{{\hat{T}}_1}$  is distributed like $K_{T_1}$  (see Lemma  A.1 ), $g_{T_1}$  and ${\hat{g}}_{{\hat{T}}_1}$  have the same distribution, which also implies that ${\hat{T}}_1$  is distributed like $T_1$  because, due to the parametrization by capacity of $\gamma$  and $\hat{\gamma }$  , $2T_1$  is exactly the coefficient of the term $1/z$  in the expansion at infinity of $g_{T_1}$  , and $2{\hat{T}}_1$  is exactly the coefficient of the term $1/z$  in the expansion at infinity of ${\hat{g}}_{{\hat{T}}_1}$  . Moreover, it is also clear that $\hat{\gamma }\left({\hat{T}}_1\right)$  is distributed like $\gamma \left(T_1\right)$  , because their distributions are both determined by Cardy's formula, and so ${\hat{g}}_{{\hat{T}}_1}\left(\hat{\gamma }\right({\hat{T}}_1\left)\right)$  is distributed like $g_{T_1}\left(\gamma \right(T_1\left)\right)$  . Notice that the law of the hull ${\hat{K}}_{{\hat{T}}_1}$  is conformally invariant because, by Lemma  A.1 , it coincides with the law of the $SL{E}_6$  hull $K_{T_1}$  .

Using now the Markovian character of ${\hat{X}}_j$  , which implies that, conditioned on ${\hat{X}}_1=({\hat{K}}_{{\hat{T}}_1},\hat{\gamma }({\hat{T}}_1\left)\right)$  , ${\hat{K}}_{{\hat{T}}_2}\backslash {\hat{K}}_{{\hat{T}}_1}$  and $\hat{\gamma }\left({\hat{T}}_2\right)$  are determined by Cardy's formula in ${\hat{G}}_2$  , from the fact that ${\hat{K}}_{{\hat{T}}_1}$  is equidistributed with $K_{T_1}$  and therefore ${\hat{G}}_2$  is equidistributed with $G_2$  , we obtain that the hull ${\hat{K}}_{{\hat{T}}_2}$  is distributed like $K_{T_2}$  and its “tip” $\hat{\gamma }\left({\hat{T}}_2\right)$  is distributed like the “tip” $\gamma \left(T_2\right)$  of the hull $K_{T_2}$  , and we can conclude that the joint distribution of $\left\{\hat{\gamma }\right({\hat{T}}_1),\hat{\gamma }({\hat{T}}_2\left)\right\}$  is the same as the joint distribution of $\left\{\gamma \right(T_1),\gamma (T_2\left)\right\}$  . It also follows immediately that ${\hat{g}}_{{\hat{T}}_2}$  is equidistributed with $g_{T_2}$  and ${\hat{\tau }}_2$  is equidistributed with ${\tau }_2$  or indeed with ${\tau }_1$  .

By repeating this argument recursively, using at each step the Markovian character of the hulls and tips, we obtain that, for all $j$  , the joint distribution of $\left\{\hat{\gamma }\right({\hat{T}}_1),\dots ,\hat{\gamma }({\hat{T}}_j\left)\right\}$  is the same as the joint distribution of $\left\{\gamma \right(T_1),\dots ,\gamma (T_j\left)\right\}$  . This immediately implies that ${\hat{\gamma }}_{\varepsilon }$  has the same distribution as ${\gamma }_{\varepsilon }$  .

In order to conclude the proof, we just have to show that, as $\varepsilon \to 0$  , ${\hat{\gamma }}_{\varepsilon }$  converges to $\hat{\gamma }$  and ${\gamma }_{\varepsilon }$  to $\gamma$  in the uniform metric on continuous curves ( 2 ). This, however, follows easily from properties of continuous curves, if we can show that the time intervals ${\hat{T}}_{j+1}-{\hat{T}}_j={\hat{\tau }}_{j+1}$  and $T_{j+1}-T_j={\tau }_{j+1}$  go to $0$  as $\varepsilon \to 0$  . To see this, we recall that ${\hat{\tau }}_{j+1}$  and ${\tau }_{j+1}$  are distributed like ${\tau }_1$  and use Lemma 2.1 of [21] , which implies the (deterministic) bound ${\tau }_1\left(\varepsilon \right)\le \frac{1}{2}{\varepsilon }^2$  , which follows from the well-know bound $a\left(t\right)\le {\varepsilon }^2$  for the capacity $a\left(t\right)=2t$  of ( 21 ).

Remark A.2. The procedure for constructing the polygonal approximations of $\hat{\gamma }$  and $\gamma$  and the recursive strategy for proving that they have the same distribution include significant modifications to the sketched argument forconvergence of the percolation exploration process to chordal $SL{E}_6$  given by Smirnov in [30] and [31] . One modification is that we use “conformal semi-balls” instead of balls (see [30, 31] ) to define the sequences of stopping times $\left\{{\hat{T}}_j\right\}$  and $\left\{T_j\right\}$  . Since the paths we are dealing with touch themselves (or almost do), if one were to use balls, some of them would intersect multiple disjoint pieces of the past hull, making it impossible to use Cardy's formula in the “triangular setting” proposed by (Carleson and) Smirnov and used here. The use of conformally mapped semi-balls ensures, thanks to the choice of the conformal maps, that the domains used to define the stopping times intersect a single piece of the past hull. This seems to be a natural choice (exploiting the conformal invariance) to obtain a good polygonal approximation of the paths while still being able to use Cardy's formula to determine hitting distributions.

We will next prove a version of Smirnov's result (Theorem  4 above) extended to cover the convergence of crossing probabilities to Cardy's formula for the case of sequences of domains. The statement of Theorem  6 below is certainly not optimal, but it is sufficient for our purposes. We remark that a weaker statement restricted, for instance, only to Jordan domains would not be sufficient – see Figure  19 and the discussion referring to it in the proof of Theorem  7 below. First, we introduce some definitions that will simplify the notation in the rest of the paper.

We will consider bounded simply connected domains $D$  whose boundaries $\partial D$  are continuous curves. Let $f:\mathbb{D}\to D$  be the (unique) conformal map from the unit disc to $D$  with $f\left(0\right)=0$  and $f^{\prime }\left(0\right)>0$  ; note that by Theorem  10 of Appendix  B , $f$  has a continuous extension to $\overline{\mathbb{D}}$  . Let $a,c,d$  be three points of $\partial D$  (or more accurately, three prime ends) in counterclockwise order – i.e., such that $a=f\left(a^{*}\right)$  , $c=f\left(c^{*}\right)$  and $d=f\left(d^{*}\right)$  , with $a^{*}$  , $c^{*}$  and $d^{*}$  in counterclockwise order. We will call $D$  admissible with respect to $(a,c,d)$  if the counterclockwise arcs $\overline{ac}$  , $\overline{cd}$  and $\overline{da}$  are simple curves, $\overline{cd}$  does not touch the interior of either $\overline{ac}$  or $\overline{da}$  , and from each point in $\overline{cd}$  there is a path to infinity that does not cross $\partial D$  .

Notice that, according to our definition, the interiors of the arcs $\overline{ac}$  and $\overline{da}$  can touch.

If that happens, the double-points of the boundary (belonging to both arcs) are counted twice and considered as two distinct points (and are two different prime ends). The significance of the notion of admissible is that certain domains arising naturally in the proof of Theorem  7 (and thus in the proof of convergence to $SL{E}_6$  ) are not Jordan but are admissible – see Figure  19 .

Consider a sequence of admissible domains $D_k$  with $j$  distinct selected points $u_k^1,\dots ,u_k^j$  on $\partial D_k$  on each of their boundaries. If $D$  is an admissible domain with $j$  selected distinct points $u^1,\dots ,u^j$  on its boundary such that, as $k\to \infty$  , $\text{d}(\partial D_k,\partial D)\to 0$  and $|u_k^1-u^1|,\dots ,|u_k^j-u^j|\to 0$  , we say that $(D_k,u_k^1,\dots ,u_k^j)$  converges to $(D,u^1,\dots ,u^j)$  and write $(D_k,u_k^1,\dots ,u_k^j)\to (D,u^1,\dots ,u^j)$  .

Theorem 6. Consider a sequence $\left\{\right(D_k,a_k,b_k,c_k,d_k\left)\right\}$  of domains containing the origin, admissible with respect to $(a_k,c_k,d_k)$  , and with $b_k$  belonging to the interior of the counterclockwise arc $\overline{c_k{d}_k}$  . Assume that $(D_k,a_k,b_k,c_k,d_k)$  converges, as $k\to \infty$  , to $(D,a,b,c,d)$  , where $D$  is a domain containing the origin, admissible with respect to $(a,c,d)$  , and $b$  belongs to the interior of the counterclockwise arc $\overline{cd}$  . Then, for any sequence ${\delta }_k\downarrow 0$  , the probability ${\Phi }_k^{{\delta }_k}(\equiv {\Phi }_{D_k}^{{\delta }_k})$  of a blue crossing inside $D_k$  from the counterclockwise segment $\overline{a_k{c}_k}$  of $\partial D_k$  to the counterclockwise segment $\overline{b_k{d}_k}$  of $\partial D_k$  converges, as $k\to \infty$  , to Cardy's formula ${\Phi }_D$  (see ( 16 )) for a crossing inside $D$  from the counterclockwise segment $\overline{ac}$  of $\partial D$  to the counterclockwise segment $\overline{cd}$  of $\partial D$  .

Proof. We will construct for each small $\varepsilon >0$  , two domains with boundary points, $(\tilde{D},\tilde{a},\tilde{b},\tilde{c},\tilde{d})$  and $(\hat{D},\hat{a},\hat{b},\hat{c},\hat{d})$  , approximating $(D,a,b,c,d)$  in such a way that not only does ${\tilde{\Phi }}_{\varepsilon }\equiv {\Phi }_{\tilde{D}}(\tilde{a},\tilde{c};\tilde{b},\tilde{d})\stackrel{\varepsilon \to 0}{⟶}\Phi \equiv {\Phi }_D(a,c;b,d)$  (by the continuity of Cardy's formula – see Lemma  B.2 ) and the same for ${\hat{\Phi }}_{\varepsilon }$  , but also so that

$$\begin{array}{c}{\tilde{\Phi }}_{\varepsilon }={liminf}_{k\to \infty }{\Phi }_{\tilde{D}\left(\varepsilon \right)}^{{\delta }_k}\le {liminf}_{k\to \infty }{\Phi }_k^{{\delta }_k}\le {limsup}_{k\to \infty }{\Phi }_k^{{\delta }_k}\le {limsup}_{k\to \infty }{\Phi }_{\hat{D}\left(\varepsilon \right)}^{{\delta }_k}={\hat{\Phi }}_{\varepsilon }.\end{array}$$

(22)

This yields the desired result. The construction of the approximating domains uses fairly straightforward conformal mapping arguments. We provide details for $\tilde{D}$  ; the construction of $\hat{D}$  is analogous. Before providing the construction details for $\tilde{D}$  , we give several paragraphs of overview.

To construct $\tilde{D}$  , we will first need to take an inner approximation $\underline{E}=\underline{E}(D,\varepsilon )$  of $\partial D$  , where $\underline{E}$  is a simple loop surrounded by $\partial D$  and with $\text{d}(\partial D,\underline{E})\le \varepsilon$  . We will then construct a Jordan domain $G\left(\varepsilon \right)$  whose boundary $\partial G$  is composed of pieces of $\partial D$  and of $\underline{E}$  , plus four segments joining $\partial D$  with $\underline{E}$  (see Figure  10 ), as we explain below.

Next, we will take an outer approximation $\overline{E}=\overline{E}(G,\varepsilon )$  of $\partial G$  , where $\overline{E}$  is a simple loop surrounding $\partial G$  and with $\text{d}(\partial G,\overline{E})\le \varepsilon$  . We will also need four simple curves $\{{\partial }_a,{\partial }_b,{\partial }_c,{\partial }_d\}$  in the interior of the (topological) annulus between $\underline{E}$  and $\overline{E}$  and connecting their endpoints $\left\{\right(\underline{a},\overline{a}),(\underline{b},\overline{b}),(\underline{c},\overline{c}),(\underline{d},\overline{d}\left)\right\}$  on $\underline{E}$  and $\overline{E}$  with each of the four curves touching $\partial D$  at exactly one point which is either in the interior of the counterclockwise segment $\overline{ac}$  (for ${\partial }_a$  and ${\partial }_c$  ) or else the counterclockwise segment $\overline{bd}$  (for ${\partial }_b$  and ${\partial }_d$  ). Furthermore each of these connecting curves is close to its corresponding point $a,b,c$  , or $d$  ; i.e., $\text{d}({\partial }_a,a)\le \varepsilon$  , etc. (see Figure  10 where each connecting segment represents “half ” of one of these connecting curves).

We will take $\tilde{a}=\overline{a}$  , $\tilde{b}=\overline{b}$  , $\tilde{c}=\overline{c}$  , and $\tilde{d}=\overline{d}$  with $\partial \tilde{D}$  the concatenation of: ${\partial }_a$  from $\underline{a}$  to $\overline{a}$  , the portion of $\overline{E}$  from $\overline{a}$  to $\overline{c}$  counterclockwise, ${\partial }_c$  from $\overline{c}$  to $\underline{c}$  , the portion of $\underline{E}$  from $\underline{c}$  to $\underline{b}$  counterclockwise, ${\partial }_b$  from $\underline{b}$  to $\overline{b}$  , the portion of $\overline{E}$  from $\overline{b}$  to $\overline{d}$  counterclockwise, ${\partial }_d$  from $\overline{d}$  to $\underline{d}$  , and the portion of $\underline{E}$  from $\underline{d}$  to $\underline{a}$  counterclockwise. It is important that (for fixed $\varepsilon$  and $\tilde{D}$  ) there is a strictly positive minimal distance between $\underline{E}$  and $\partial D$  , and between $\partial \tilde{D}$  and the union of the two counterclockwise segments $\overline{cb}$  and $\overline{da}$  of $\partial D$  . These features will guarantee (as we explain with more detail below) that for fixed $\varepsilon$  , once $k$  is large enough, a continuous curve within $\tilde{D}$  that corresponds to the crossing event whose probability is ${\Phi }_{\tilde{D}}(\tilde{a},\tilde{c};\tilde{b},\tilde{d})$  must have a subpath corresponding to the crossing event in $D_k$  whose probability is ${\Phi }_{D_k}(a_k,c_k;b_k,d_k)$  . This is the key feature of $\tilde{D}$  , which will yield the first inequality of ( 22 ).

We now give a more detailed explanation of the construction of $\tilde{D}$  . We first construct the parts of $\partial \tilde{D}$  that are inside $\partial D$  . Let $f$  be the conformal map from $\mathbb{D}$  onto $D$  with $f\left(0\right)=0$  and $f^{\prime }\left(0\right)>0$  , and consider the image $f(\partial {\mathbb{D}}_{1-{\varepsilon }^{\prime }})$  of the circle $\partial {\mathbb{D}}_{1-{\varepsilon }^{\prime }}=\{z:|z|=1-{\varepsilon }^{\prime }\}$  under $f$  and the inverse images, $a^{*},b^{*},c^{*},d^{*}$  , under $f^{-1}$  of $a,b,c,d$  . Let ${\partial }_a^{*}({\varepsilon }^{\prime },{\varepsilon }_a)$  be the straight segment between $e^{-i{\varepsilon }_a}{a}^{*}$  on the unit circle $\partial \mathbb{D}$  , and $(1-{\varepsilon }^{\prime })e^{-i{\varepsilon }_a}{a}^{*}$  on the circle $\partial {\mathbb{D}}_{1-{\varepsilon }^{\prime }}$  , and define ${\partial }_b^{*},{\partial }_c^{*}$  , and ${\partial }_d^{*}$  similarly, but using clockwise rotations by $e^{+i{\varepsilon }_c}$  and $e^{+i{\varepsilon }_d}$  for $c^{*}$  and $d^{*}$  (see Figure  11 ).

$f(\partial {\mathbb{D}}_{1-{\varepsilon }^{\prime }})$  is a candidate for $\underline{E}(D,\varepsilon )$  and $f\left({\partial }_{♯}^{*}\right({\varepsilon }^{\prime },{\varepsilon }_{♯}\left)\right)$  is a candidate for half of ${\partial }_{♯}$  (where $♯=a$  or $b$  or $c$  or $d$  ), so we must choose ${\varepsilon }^{\prime }$  and the ${\varepsilon }_{♯}$  's small enough so that $\text{d}(\partial D,f({\mathbb{D}}_{1-{\varepsilon }^{\prime }}\left)\right)\le \varepsilon$  , $\text{d}\left(f\right({\partial }_a^{*}({\varepsilon }^{\prime },{\varepsilon }_a)),a)\le \varepsilon$  , etc. We then define $\underline{a}=f\left(\right(1-{\varepsilon }^{\prime }\left)e^{-i{\varepsilon }_a}{a}^{*}\right)$  and similarly for $\underline{b}$  , $\underline{c}$  and $\underline{d}$  (see Figure  12 ).

Consider now the Jordan domain $G=G\left(\varepsilon \right)$  whose boundary $\partial G$  is given by the concatenation of: $f\left({\partial }_a^{*}\right)$  from $\underline{a}$  to $f\left(e^{-i{\varepsilon }_a}{a}^{*}\right)$  , the portion of $\partial D$  from $f\left(e^{-i{\varepsilon }_a}{a}^{*}\right)$  to $f\left(e^{i{\varepsilon }_c}{c}^{*}\right)$  counterclockwise, $f\left({\partial }_c^{*}\right)$  from $f\left(e^{i{\varepsilon }_c}{c}^{*}\right)$  to $\underline{c}$  , the portion of $\underline{E}$  from $\underline{c}$  to $\underline{b}$  counterclockwise, $f\left({\partial }_b^{*}\right)$  from $\underline{b}$  to $f\left(e^{-i{\varepsilon }_b}{b}^{*}\right)$  , the portion of $\partial D$  from $f\left(e^{-i{\varepsilon }_b}{b}^{*}\right)$  to $f\left(e^{i{\varepsilon }_d}{d}^{*}\right)$  counterclockwise, $f\left({\partial }_d^{*}\right)$  from $f\left(e^{i{\varepsilon }_d}{d}^{*}\right)$  to $\underline{d}$  , and the portion of $\underline{E}$  from $\underline{d}$  to $\underline{a}$  counterclockwise (see Figure  13 ).

The exterior of this new Jordan domain $G$  is a connected domain for which we can do a construction analogous to the one for the original domain $D$  using a conformal mapfrom $\mathbb{D}$  to obtain candidates for $\overline{E}(D,\varepsilon )$  and for the exterior halves of the ${\partial }_{♯}$  's. To do this, we use a conformal map from $\mathbb{D}$  to the exterior of $\partial G$  and use $\underline{a},\underline{b},\underline{c},\underline{d}$  as replacements for $a,b,c,d$  (see Figure  14 ).

Finally, we use the freedom to choose the exterior replacements for ${\varepsilon }^{\prime }$  and the ${\varepsilon }_{♯}$  's differently from the interior values to make sure that the interior and exterior halves of the ${\partial }_{♯}$  's match up. We also choose the exterior values for ${\varepsilon }^{\prime }$  small enough so that $\overline{E}$  stays between $a$  and $f\left(e^{-i{\varepsilon }_a}{a}^{*}\right)$  , between $b$  and $f\left(e^{-i{\varepsilon }_b}{b}^{*}\right)$  , between $c$  and $f\left(e^{i{\varepsilon }_c}{c}^{*}\right)$  , between $d$  and $f\left(e^{i{\varepsilon }_d}{d}^{*}\right)$  (see Figure  15 ).

Once all the pieces of $\partial \tilde{D}\left(\varepsilon \right)$  are available, they are put together as explained above, and as show in Figures  16 and  17 below.

It should be clear that for a given approximation $\partial \tilde{D}$  of $\partial D$  constructed as described above there is a strictly positive $\tilde{\varepsilon }$  such that the distance between $\partial D$  and the portions of $\partial \tilde{D}$  that belong to $\underline{E}$  and $\overline{E}$  is not smaller than $\tilde{\varepsilon }$  , and the distance between the union $\overline{cb}\cup \overline{da}$  , two counterclockwise segments of $\partial D$  , and ${\partial }_a\cup {\partial }_b\cup {\partial }_c\cup {\partial }_d$  is also not smaller than $\tilde{\varepsilon }$  . On the other hand, for any $\tilde{\varepsilon }>0$  , there exists $k_0=k_0\left(\tilde{\varepsilon }\right)$  such that for all $k\ge k_0$  , $\partial D_k$  is contained inside the $\tilde{\varepsilon }$  -neighborhood of $\partial D$  with the counterclockwise segment $\overline{c_k{b}_k}$  (resp., $\overline{d_k{a}_k}$  ) in the $\tilde{\varepsilon }$  -neighborhood of the counterclockwise segment $\overline{cb}$  (resp., $\overline{da}$  ). This implies that for $k$  large enough, any blue path crossing inside $\tilde{D}$  from the counterclockwise segment $\overline{\overline{a}\overline{c}}$  of $\partial \tilde{D}$  to the counterclockwise segment $\overline{\overline{b}\overline{d}}$  of $\partial \tilde{D}$  must have a subpath that stays inside $D_k$  and crosses between the counterclockwise segment $\overline{a_k{c}_k}$  of $\partial D_k$  and thecounterclockwise segment $\overline{b_k{d}_k}$  of $\partial D_k$  . Therefore the crossing probability ${\tilde{\Phi }}_{\tilde{D}\left(\varepsilon \right)}^{{\delta }_k}$  is a lower bound for ${\Phi }_k^{{\delta }_k}$  for all $k\ge k_0$  , so that

$$\begin{array}{c}{liminf}_{k\to \infty }{\Phi }_k^{{\delta }_k}\ge {lim}_{k\to \infty }{\tilde{\Phi }}_{\tilde{D}\left(\varepsilon \right)}^{{\delta }_k}={\tilde{\Phi }}_{\varepsilon },\end{array}$$

(23)

as desired (the equality uses Smirnov's result, Theorem  4 , for fixed $\tilde{D}\left(\varepsilon \right)$  ).

We now note that as $\varepsilon \to 0$  , $(\tilde{D},\tilde{a},\tilde{b},\tilde{c},\tilde{d})\to (D,a,b,c,d)$  . This allows us to use the continuity of Cardy's formula (Lemma  B.2 in Appendix  B ) to obtain

$$\begin{array}{c}{lim}_{\varepsilon \to 0}{\tilde{\Phi }}_{\varepsilon }=\Phi .\end{array}$$

(24)

From this and ( 23 ) it follows that

$$\begin{array}{c}{liminf}_{k\to \infty }{\Phi }_k^{{\delta }_k}\ge \Phi .\end{array}$$

(25)

The remaining part of the proof involves defining a domain $\hat{D}$  analogous to $\tilde{D}$  but with the property that the probability ${\hat{\Phi }}_{\varepsilon }^{{\delta }_k}$  of an appropriate crossing, such that

$$\begin{array}{c}{lim}_{\varepsilon \to 0}{lim}_{k\to \infty }{\hat{\Phi }}_{\hat{D}\left(\varepsilon \right)}^{{\delta }_k}=\Phi ,\end{array}$$

(26)

is an upper bound for ${\Phi }_k^{{\delta }_k}$  for all $k$  large enough. (The details of the construction of $\hat{D}$  are analogous to those of $\tilde{D}$  ; we leave them to the reader.) This shows that

$$\begin{array}{c}{limsup}_{k\to \infty }{\Phi }_k^{{\delta }_k}\le \Phi ,\end{array}$$

(27)

which, combined with ( 25 ), implies

$$\begin{array}{c}{lim}_{k\to \infty }{\Phi }_k^{{\delta }_k}=\Phi \end{array}$$

(28)

and concludes the proof.

We note that Theorem  6 , combined with the continuity of Cardy's formula in the shape of the domain (for admissible domains) and positions of the four points on the boundary (Lemma  B.2 ), implies the convergence of crossing probabilities to Cardy's formula locally uniformly in the shape of the domain with respect to the uniform metric on curves, and in the location of the four points on the boundary with respect to the Euclidean metric; i.e., for $(D,a,b,c,d)$  an admissible domain with $a,b,c,d\in \partial D$  (with the notation used in Theorem  6 ), $\forall \varepsilon >0$  , $\exists {\alpha }_0={\alpha }_0\left(\varepsilon \right)$  and ${\delta }_0={\delta }_0\left(\varepsilon \right)$  such that for all admissible domains $(D^{\prime },a^{\prime },b^{\prime },c^{\prime },d^{\prime })$  with $max\left(\text{d}\right(\partial D,\partial D^{\prime }),|a-a^{\prime }|,|b-b^{\prime }|,|c-c^{\prime }|,|d-d^{\prime }\left|\right)\le {\alpha }_0$  and $\delta \le {\delta }_0$  , $\left|{\Phi }_{D^{\prime }}\right(a^{\prime },c^{\prime };b^{\prime },d^{\prime })-{\Phi }_{D^{\prime }}^{\delta }(a^{\prime },c^{\prime };b^{\prime },d^{\prime }\left)\right|\le \varepsilon$  , where ${\Phi }_{D^{\prime }}(a^{\prime },c^{\prime };b^{\prime },d^{\prime })$  is Cardy's formula and ${\Phi }_{D^{\prime }}^{\delta }(a^{\prime },c^{\prime };b^{\prime },d^{\prime })$  is the corresponding crossing probability.

Our next task is to show that the filling of any subsequential scaling limit of the percolation exploration process satisfies the spatial Markov property. Let us start withsome notation. First of all, suppose that $D$  is a Jordan domain, and $a$  and $b$  are two points on $\partial D$  . Define the discrete filling (or simply filling) at time $t$  of a percolation exploration path ${\gamma }_{D,a,b}^{\delta }$  (with a given parametrization) inside (the $\delta$  -approximation $D^{\delta }$  of ) $D$  from (the e-vertex closest to) $a$  to (the e-vertex closest to) $b$  to be the union of the hexagons explored up to time $t$  and those unexplored hexagons from which it is not possible to reach $b$  without crossing an explored hexagon or $\partial D$  (in other words, this is the set of hexagons that at time $t$  have been explored or are disconnected from $b$  by the exploration path).

Consider a sequence $\left\{\right(D_k,a_k,b_k\left)\right\}$  of Jordan domains such that $(D_k,a_k,b_k)\to (D,a,b)$  , where $a_k,b_k\in \partial D_k$  are two distinct points on $\partial D_k$  . Denote by ${\gamma }_k^{\delta }\equiv {\gamma }_{D_k,a_k,b_k}^{\delta }$  the percolation exploration path inside (the $\delta$  -approximation $D_k^{\delta }$  of ) $D_k$  from (the e-vertex closest to) $a_k$  to (the e-vertex closest to) $b_k$  .

Notice that we can couple the paths ${\gamma }_k^{\delta }$  simultaneously for all values of $k$  and $\delta$  by using the same percolation configuration to generate all of them. We can then apply the results of [2] to conclude that there exists a subsequence ${\delta }_k\downarrow 0$  such that the law of ${\gamma }_k^{{\delta }_k}$  converges to some limiting law for a process $\tilde{\gamma }$  supported on (Hölder) continuous curves inside $D$  from $a$  to $b$  . The filling ${\tilde{K}}_t$  of $\tilde{\gamma }[0,t]$  , appearing in the next theorem, is defined just above Equation ( 17 ).

Theorem 7. For any subsequential limit process $\tilde{\gamma }$  of the percolation exploration path ${\gamma }_k^{\delta }$  defined above, the filling ${\tilde{K}}_t$  of $\tilde{\gamma }[0,t]$  , as a process, satisfies the spatial Markov property.

Proof. Let ${\delta }_k\downarrow 0$  be a subsequence such that the law of ${\gamma }_k^{{\delta }_k}$  converges to some limiting law supported on continuous curves $\tilde{\gamma }$  inside $D$  from $a$  to $b$  . We will prove the spatial Markov property by showing that $({\tilde{K}}_{{\tilde{T}}_j},\tilde{\gamma }({\tilde{T}}_j\left)\right)$  as defined in the proof of Theorem  5 are jointly distributed like the corresponding $SL{E}_6$  hull variables, which do have the spatial Markov property. For each $k$  , let $K_t^k$  denote the filling at time $t$  of ${\gamma }_k^{{\delta }_k}$  (with some parametrization – we do not need to worry about the choice of parametrization here). It follows from the Markovian character of the percolation exploration process that, for all $k$  , the filling $K_t^k$  of the percolation exploration path ${\gamma }_k^{{\delta }_k}$  satisfies a suitably adapted (to the discrete setting) spatial Markov property. (In fact, the percolation exploration path satisfies a stronger property – roughly speaking, that the future of the path given the filling of the past is distributed as a percolation exploration path in the original domain from which the filling of the past has been removed.) To be more precise, let $f_0^k$  be a conformal transformation that maps $D_k$  to $\mathbb{H}$  such that $f_0^k\left(a_k\right)=0$  and $f_0^k\left(b_k\right)=\infty$  and let $T_1^k=T_1^k\left(\varepsilon \right)$  denote the first exit time of ${\gamma }_k^{{\delta }_k}\left(t\right)$  from $G_1^k\equiv \left(f_0^k{)}^{-1}\right(C(0,\varepsilon ))$  defined as the first time an explored hexagon intersects the image under $(f_0^k{)}^{-1}$  of the semi-circle $\{z:|z|=\varepsilon \}\cap \mathbb{H}$  . Define recursively $T_{j+1}^k$  as the first exit time of ${\gamma }_k^{{\delta }_k}[T_j^k,\infty )$  from $G_{j+1}^k\equiv \left(f_{T_j^k}^k{)}^{-1}\right(C(0,\varepsilon ))$  , where $f_{T_j^k}^k$  is a conformal map from $D_k\backslash K_{T_j^k}^k$  to $\mathbb{H}$  that maps ${\gamma }_k^{{\delta }_k}\left(T_j^k\right)$  to $0$  and $b_k$  to $\infty$  . We also define ${\tau }_{j+1}^k\equiv T_{j+1}^k-T_j^k$  , so that $T_j^k={\tau }_1^k+\dots +{\tau }_j^k$  , and the (discrete-time) stochastic process $X_j^k\equiv (K_{T_j^k}^k,{\gamma }_k^{{\delta }_k}(T_j^k\left)\right)$  for $j=1,2,\dots$  . The Markovian character of the percolation exploration process impliesthat, for every $k$  , $X_j^k$  is a Markov process (in $j$  ). In order to study the limit as $k\to \infty$  of $X_1^k,X_2^k,\dots$  , we first need to analyze in more detail the mappings $f_{T_j^k}^k$  .

The conformal transformations $f_{T_j^k}^k$  are defined as follows, with an arbitrary multiplicative factor ${\lambda }_j$  . We choose $f_{T_j^k}^k$  to be the composition ${\lambda }_j{\psi }_j^k\circ {\phi }_j^k$  of two maps, where ${\lambda }_j>0$  , ${\phi }_j^k$  is the conformal transformation that maps $D_k\backslash K_{T_j^k}^k$  onto $\mathbb{D}$  with ${\phi }_j^k\left(0\right)=0$  and $\left({\phi }_j^k{)}^{\prime }\right(0)>0$  (we are assuming for simplicity that the domain $D_k\backslash K_{T_j^k}^k$  contains the origin; if that is not the case, one can think of a translated domain that does contain the origin), while ${\psi }_j^k$  is the inverse of the transformation

$$\begin{array}{c}w=e^{i{\theta }_j^k}\left(\frac{(z+1)-z_j^k}{(z+1)-\overline{z_j^k}}\right)\end{array}$$

(29)

that maps $\overline{\mathbb{H}}$  onto $\overline{\mathbb{D}}$  , where ${\theta }_j^k$  is chosen so that $e^{i{\theta }_j^k}={\phi }_j^k\left(b_k\right)$  and $z_j^k$  can be chosen so that $|1-z_j^k|=1$  , $\text{Im}\left(z_j^k\right)>0$  and ${\phi }_j^k\left(b_k\right)\left(\frac{1-z_j^k}{1-\overline{z_j^k}}\right)={\phi }_j^k\left({\gamma }_k^{{\delta }_k}\right(T_j^k\left)\right)$  , which means that ${\lambda }_j{\psi }_j^k$  maps ${\phi }_j^k\left({\gamma }_k^{{\delta }_k}\right(T_j^k\left)\right)$  to $0$  and ${\phi }_j^k\left(b_k\right)$  to $\infty$  , so that $f_{T_j^k}^k={\lambda }_j{\psi }_j^k\circ {\phi }_j^k$  indeed maps ${\gamma }_k^{{\delta }_k}\left(T_j^k\right)$  to $0$  and $b_k$  to $\infty$  .

Since ${\gamma }_k^{{\delta }_k}$  converges in distribution to $\tilde{\gamma }$  , we can find two coupled versions of ${\gamma }_k^{{\delta }_k}$  and $\tilde{\gamma }$  on some probability space $({\Omega }^{\prime },{\mathcal{ℬ}}^{\prime },{\mathbb{P}}^{\prime })$  such that ${\gamma }_k^{{\delta }_k}$  converges to $\tilde{\gamma }$  for all ${\omega }^{\prime }\in {\Omega }^{\prime }$  ; in the rest of the proof we work with these new versions which, with a slight abuse of notation, we denote with the same names as the original ones. Let ${\tilde{f}}_0$  be a conformal transformation that maps $D$  to $\mathbb{H}$  such that ${\tilde{f}}_0\left(a\right)=0$  and ${\tilde{f}}_0\left(b\right)=\infty$  and let ${\tilde{T}}_1={\tilde{T}}_1\left(\varepsilon \right)$  denote the first time $\tilde{\gamma }\left(t\right)$  hits $D\backslash {\tilde{G}}_1$  , with ${\tilde{G}}_1\equiv {\tilde{f}}_0^{-1}\left(C\right(0,\varepsilon \left)\right)$  . Define recursively ${\tilde{T}}_{j+1}$  as the first time $\tilde{\gamma }\left(t\right)$  hits $D\backslash {\tilde{G}}_{j+1}$  , with ${\tilde{G}}_{j+1}\equiv {\tilde{f}}_{{\tilde{T}}_j}^{-1}\left(C\right(0,\varepsilon \left)\right)$  , where ${\tilde{f}}_{{\tilde{T}}_j}$  is a conformal map from $D\backslash {\tilde{K}}_{{\tilde{T}}_j}$  to $\mathbb{H}$  that maps $\tilde{\gamma }\left({\tilde{T}}_j\right)$  to $0$  and $b$  to $\infty$  . We also define ${\tilde{\tau }}_j\equiv {\tilde{T}}_{j+1}-{\tilde{T}}_j$  , so that ${\tilde{T}}_j={\tilde{\tau }}_1+\dots +{\tilde{\tau }}_j$  , and the (discrete-time) stochastic process ${\tilde{X}}_j\equiv ({\tilde{K}}_{{\tilde{T}}_j},\tilde{\gamma }({\tilde{T}}_j\left)\right)$  . As above, we choose the conformal transformation ${\tilde{f}}_{{\tilde{T}}_j}$  to be the composition ${\lambda }_j{\tilde{\psi }}_j\circ {\tilde{\phi }}_j$  of two maps, where ${\lambda }_j>0$  , ${\tilde{\phi }}_j$  is the conformal transformation that maps $D\backslash {\tilde{K}}_{T_j}$  onto $\mathbb{D}$  with ${\tilde{\phi }}_j\left(0\right)=0$  and ${\tilde{\phi }}_j^{\prime }\left(0\right)>0$  (once again, we are assuming for simplicity that the domain $D\backslash {\tilde{K}}_{T_j}$  contains the origin; if that is not the case, one can think of a translated domain that does contain the origin), while ${\tilde{\psi }}_j$  is the inverse of the transformation

$$\begin{array}{c}w=e^{i{\tilde{\theta }}_j}\left(\frac{(z+1)-{\tilde{z}}_j}{(z+1)-\overline{\stackrel{j}{\stackrel{~}{z}}}}\right)\end{array}$$

(30)

that maps $\overline{\mathbb{H}}$  onto $\overline{\mathbb{D}}$  , where ${\tilde{\theta }}_j$  is chosen so that $e^{i{\tilde{\theta }}_j}={\tilde{\phi }}_j\left(b\right)$  and ${\tilde{z}}_j$  can be chosen so that $|1-{\tilde{z}}_j|=1$  , $\text{Im}\left({\tilde{z}}_j\right)>0$  and ${\tilde{\phi }}_j\left(b\right)\left(\frac{1-{\tilde{z}}_j}{1-\overline{\stackrel{j}{\stackrel{~}{z}}}}\right)={\tilde{\phi }}_j\left(\tilde{\gamma }\right({\tilde{T}}_j\left)\right)$  , which means that ${\lambda }_j{\tilde{\psi }}_j$  maps ${\tilde{\phi }}_j\left(\tilde{\gamma }\right({\tilde{T}}_j\left)\right)$  to $0$  and ${\tilde{\phi }}_j\left(b\right)$  to $\infty$  , so that ${\tilde{f}}_{{\tilde{T}}_j}={\lambda }_j{\tilde{\psi }}_j\circ {\tilde{\phi }}_j$  indeed maps $\tilde{\gamma }\left({\tilde{T}}_j\right)$  to $0$  and $b$  to $\infty$  .

Analogous quantities can be defined for the trace of chordal $SL{E}_6$  . For clarity, they will be indicated here by the superscript $SL{E}_6$  ; e.g., $f_j^{SL{E}_6}$  , $K_{T_j}^{SL{E}_6}$  , $G_j^{SL{E}_6}$  and $X_j^{SL{E}_6}$  .We want to show recursively that, for any $j$  , as $k\to \infty$  , $\{X_1^k,\dots ,X_j^k\}$  converge jointly in distribution to $\{{\tilde{X}}_1,\dots ,{\tilde{X}}_j\}$  . By recursively applying Theorem  6 and Lemma  A.1 , we can then conclude that $\{{\tilde{X}}_1,\dots ,{\tilde{X}}_j\}$  are jointly equidistributed with the corresponding $SL{E}_6$  hull variables (at the corresponding stopping times) $\{X_1^{SL{E}_6},\dots ,X_j^{SL{E}_6}\}$  . Since the latter do satisfy the spatial Markov property, so will the former, as desired.

The zeroth step consists in noticing that the convergence of $(D_k,a_k,b_k)$  to $(D,a,b)$  as $k\to \infty$  allows us to apply Radó's theorem (i.e., Theorem  9 of Appendix  B ) to show that $({\phi }_0^k{)}^{-1}$  converges to ${\tilde{\phi }}_0^{-1}$  uniformly in $\overline{\mathbb{D}}$  . This, together with the convergence of $a_k$  to $a$  and $b_k$  to $b$  , implies that ${\phi }_0^k\left(a_k\right)$  converges to ${\tilde{\phi }}_0\left(a\right)$  and ${\phi }_0^k\left(b_k\right)$  to ${\tilde{\phi }}_0\left(b\right)$  . Therefore, we also have the convergence of ${\lambda }_0{\psi }_0^k$  to ${\lambda }_0{\tilde{\psi }}_0$  and we can conclude that $(f_0^k{)}^{-1}$  converges to ${\tilde{f}}_0^{-1}$  uniformly on compact subsets of $\overline{\mathbb{H}}$  , which implies that the boundary $\partial G_1^k$  of $G_1^k=\left(f_0^k{)}^{-1}\right(C(0,\varepsilon ))$  converges to the boundary $\partial {\tilde{G}}_1$  of ${\tilde{G}}_1={\tilde{f}}_0^{-1}\left(C\right(0,\varepsilon \left)\right)$  in the uniform metric on continuous curves.

Starting from there, the first step of our recursion argument is organized as follows:

• (1) $K_{T_1^k}^k\to {\tilde{K}}_{{\tilde{T}}_1}$  by “number of arms” percolation bounds [16] and Lemma  A.2 below, but also $K_{T_1^k}^k\to K_{T_1}^{SL{E}_6}$  by Lemma  A.1 (Theorem  6 is used here).
• (2) $D_k\backslash K_{T_1^k}^k\to D\backslash {\tilde{K}}_{{\tilde{T}}_1}$  , but also $D_k\backslash K_{T_1^k}^k\to D\backslash K_{T_1}^{SL{E}_6}$  , by (1).
• (3) $f_{T_1^k}^k\to {\tilde{f}}_{{\tilde{T}}_1}$  , but also $f_{T_1^k}^k\to f_{T_1}^{SL{E}_6}$  , by Corollary  B.1 .
• (4) $G_2^k\to {\tilde{G}}_2$  , but also $G_2^k\to G_2^{SL{E}_6}$  , by (3).

At this point, we are in the same situation as at the zeroth step, but with $G_1^k$  , ${\tilde{G}}_1$  and $G_1^{SL{E}_6}$  replaced by $G_2^k$  , ${\tilde{G}}_2$  and $G_2^{SL{E}_6}$  respectively, and we can proceed by recursion. As explained above, the theorem then follows from the fact that the $SL{E}_6$  hull variables do posses the spatial Markov property.

In what follows we will show that the “number of arms” bounds [16] and Lemma  A.2 below imply the convergence of $K_{T_j^k}^k$  to ${\tilde{K}}_{{\tilde{T}}_j}$  , and that the conditions to apply Theorem  6 , Lemma  A.2 and Corollary  B.1 are always satisfied. (This last point boils down to showing that the domains $D\backslash {\tilde{K}}_{{\tilde{T}}_j}$  and ${\tilde{G}}_j$  are admissible for all $j$  .) We begin by showing first that, as $k\to \infty$  , $K_{T_1^k}^k$  converges in distribution to ${\tilde{K}}_{{\tilde{T}}_1}$  (which also implies that ${\gamma }_k^{{\delta }_k}\left(T_1^k\right)$  converges in distribution to $\tilde{\gamma }\left({\tilde{T}}_1\right)$  ), from which it follows that $X_1^k$  converges in distribution to ${\tilde{X}}_1$  .

Consider $G_1^k\backslash K_{T_1^k}^k$  and ${\tilde{G}}_1\backslash {\tilde{K}}_{{\tilde{T}}_1}$  ; they are both composed of two domains (which “meet” at ${\gamma }_k^{{\delta }_k}\left(T_1^k\right)$  and $\tilde{\gamma }\left({\tilde{T}}_1\right)$  respectively), which we denote by $A_{1,1}^k$  and $A_{1,2}^k$  and by ${\tilde{A}}_{1,1}$  and ${\tilde{A}}_{1,2}$  , respectively (see Figure  18 ).

It follows from [2] that for some further subsequence $k_n$  of the $k$  's (which we denote by simply replacing $k$  by $n$  ), $({\gamma }_n^{{\delta }_n},\partial A_{1,1}^n,\partial A_{1,2}^n)$  converge jointly in distribution to some limit;we already know that ${\gamma }_n^{{\delta }_n}$  must converge to $\tilde{\gamma }$  and want to use this fact and a suitably adapted triangular array version of Lemma  5.2 (whose validity does not rely on statement (S)) to conclude that the limit is unique and coincides with $(\tilde{\gamma },\partial {\tilde{A}}_{1,1},\partial {\tilde{A}}_{1,2})$  . (Notice that $A_{1,1}^n$  and $A_{1,2}^n$  are two of the (many) domains in $G_1^n$  produced by the exploration process started at $a_n$  and stopped when it first hits the image under $(f_0^n{)}^{-1}$  of the semi-circle $\{z:|z|=\varepsilon \}\cap \mathbb{H}$  , so we are in a context close to that of Lemma  5.2 ).

First of all, we need to show that the scaling limit of $K_{T_1^n}^{{\delta }_n}$  touches the image of the semi-circle $\{z:|z|=\varepsilon \}\cap \mathbb{H}$  under $(f_0{)}^{-1}$  at a single point. This follows immediately from the definition of the stopping time $T_1^n$  for every fixed $n$  (with the map $(f_0^n{)}^{-1}$  ), but it could fail to be true in the limit $n\to \infty$  . The fact that it holds true in the limit is a direct consequence of Lemma  A.2 below (in fact, a simpler version concerning Jordan domains would suffice here, but not when we iterate the argument – see below), which also implies that the single point at which the scaling limit of $K_{T_1^n}^{{\delta }_n}$  touches the image of the semi-circle $\{z:|z|=\varepsilon \}\cap \mathbb{H}$  under $(f_0{)}^{-1}$  coincides with the limit of ${\gamma }^{{\delta }_n}\left(T_1^n\right)$  and with $\tilde{\gamma }\left({\tilde{T}}_1\right)$  .

Therefore, if we remove the single point $\tilde{\gamma }\left({\tilde{T}}_1\right)$  , the scaling limit of the boundary of $K_{T_1^n}^{{\delta }_n}$  splits into a left and a right part (corresponding to the scaling limit of the leftmost yellow and the rightmost blue $\mathcal{T}$  -paths of hexagons explored by ${\gamma }_n^{{\delta }_n}$  , respectively) that do not touch the image of the semi-circle $\{z:|z|=\varepsilon \}\cap \mathbb{H}$  under $(f_0{)}^{-1}$  .

Moreover, Lemma  A.3 below implies that if ${\gamma }_n^{{\delta }_n}$  has a “close encounter” with $\partial D_n$  , then it touches $\partial D_n^{{\delta }_n}$  . Analogously, the standard bound on the probability of six crossings of an annulus [16] , used repeatedly before, implies that wherever ${\gamma }_n^{{\delta }_n}$  has a “close encounter” with itself, there is touching (see the proof of Lemma  5.1 ). These two observations assure that the scaling limit of $K_{T_1^n}^{{\delta }_n}$  is almost surely a filling (of $\tilde{\gamma }[0,{\tilde{T}}_1]$  ), i.e., a closed connected set whose complement in $D$  is simply connected. From the same bound on the probability of six crossings of an annulus, we can also conclude that the scaling limits of the left and right boundaries of $K_{T_1^n}^n$  are almost surely simple (continuous) curves, as in the proof of Lemma  5.1 .

It is also possible to conclude that the intersection of the scaling limit of the left and right boundaries of $K_{T_1^n}^n$  with the boundary of $D$  almost surely does not contain arcs of positive length. In fact, if that were the case, it would be possible to find a subdomain $D^{\prime }$  with three counterclockwise points $z_1,z_2,z_3$  on its boundary such that the probability that an exploration path started at $z_1$  and stopped when it first hits the arc $\overline{z_2{z}_3}$  of $\partial D^{\prime }$  has a positive probability, in the scaling limit, of hitting at $z_2$  or $z_3$  , contradicting Cardy's formula (which, by Theorem  6 , holds for all subsequential scaling limits). This means that the scaling limit of $K_{T_1^n}^n$  almost surely satisfies the condition in ( 17 ) and is therefore a hull. It says as well that, almost surely, the scaling limit $\tilde{\gamma }$  of ${\gamma }_n^{{\delta }_n}$  does not “stick” to the boundary of ${\tilde{G}}_1$  , which implies that also ${\tilde{K}}_{{\tilde{T}}_1}$  satisfies the condition in ( 17 ) and is therefore a hull.

It also implies that $D\backslash {\tilde{K}}_{{\tilde{T}}_1}$  and ${\tilde{G}}_2$  are admissible domains since the part of the boundary of either $D\backslash {\tilde{K}}_{{\tilde{T}}_1}$  or ${\tilde{G}}_2$  that belongs to the boundary of ${\tilde{K}}_{{\tilde{T}}_1}$  can be split up, by removing the single point $\tilde{\gamma }\left({\tilde{T}}_1\right)$  , into two pieces which are, by an application of theproof of Lemma  5.1 , simple continuous curves, while the remaining part of the boundary of either $D\backslash {\tilde{K}}_{{\tilde{T}}_1}$  or ${\tilde{G}}_2$  is a Jordan arc whose interior does not touch the hull ${\tilde{K}}_{{\tilde{T}}_1}$  .

(Notice however, that they need not be Jordan domains because ${\tilde{K}}_{{\tilde{T}}_1}$  has cut-points with positive probability – see Figure  19 ). This will be important later, when we need to apply Lemma  A.1 (and therefore Theorem  6 ), Corollary  B.1 and Lemmas  A.2 - A.3 again.

Then, since hulls are characterized by their “envelope” (see Lemma  A.1 and the discussion preceding it), the joint convergence in distribution of $\{\partial A_{1,1}^n,\partial A_{1,2}^n\}$  to $\{\partial {\tilde{A}}_{1,1},\partial {\tilde{A}}_{1,2}\}$  would be enough to conclude that $K_{T_1^n}^n$  converges to ${\tilde{K}}_{{\tilde{T}}_1}$  as $n\to \infty$  , and in fact that $({\gamma }_n^{{\delta }_n},K_{T_1^n}^n)$  converges in distribution to $(\tilde{\gamma },K_{{\tilde{T}}_1})$  (and this will be valid also for the original subsequence $k$  and not just for the further subsequence $k_n$  ). In order to get that, as explained before, we can use the convergence in distribution of ${\gamma }_n^{{\delta }_n}$  to $\tilde{\gamma }$  and apply almost the same arguments as used in the proof of Lemma  5.2 . The only difference is that, in proving claim (C), we cannot use the bound on the probability of three crossings of an annulus centered at a boundary point because we are not necessarily dealing with a convex domain. To replace that bound we use once again Lemmas  A.2 - A.3 below (a simpler version concerning Jordan domains would again suffice here, but not when we iterate the argument – see below).

We can then conclude that $K_{T_1^n}^n$  converges in distribution to ${\tilde{K}}_{{\tilde{T}}_1}$  , which in turn implies the joint convergence in distribution of $(K_{T_1^k}^k,{\gamma }_k^{{\delta }_k}(T_1^k\left)\right)$  to $({\tilde{K}}_{{\tilde{T}}_1},\tilde{\gamma }({\tilde{T}}_1\left)\right)$  and concludes the first step of the argument.

We next need to prove that $\left(\right(K_{T_1^k}^k,{\gamma }_k^{{\delta }_k}\left(T_1^k\right)),(K_{T_2^k}^k,{\gamma }_k^{{\delta }_k}\left(T_2^k\right)\left)\right)$  converges in distribution to $\left(\right({\tilde{K}}_{{\tilde{T}}_1},\tilde{\gamma }\left({\tilde{T}}_1\right)),({\tilde{K}}_{{\tilde{T}}_2},\tilde{\gamma }\left({\tilde{T}}_2\right)\left)\right)$  . Since we have already proved the convergence of $(K_{T_1^k}^k,{\gamma }_k^{{\delta }_k}(T_1^k\left)\right)$  to $({\tilde{K}}_{{\tilde{T}}_1},\tilde{\gamma }({\tilde{T}}_1\left)\right)$  , we claim that all we really need to prove is the convergence of $(K_{T_2^k}^k\backslash K_{T_1^k}^k,{\gamma }_k^{{\delta }_k}(T_2^k\left)\right)$  to $({\tilde{K}}_{{\tilde{T}}_2}\backslash {\tilde{K}}_{{\tilde{T}}_1},\tilde{\gamma }({\tilde{T}}_2\left)\right)$  . To see this, notice that $K_{T_2^k}^k\backslash K_{T_1^k}^k$  is distributed like the hull of a percolation exploration path inside $D_k\backslash K_{T_1^k}^k$  . Besides, the convergence in distribution of $(K_{T_1^k}^k,{\gamma }_k^{{\delta }_k}(T_1^k\left)\right)$  to $({\tilde{K}}_{{\tilde{T}}_1},\tilde{\gamma }({\tilde{T}}_1\left)\right)$  implies that we can find versions of $({\gamma }_k^{{\delta }_k},K_{T_1^k}^k)$  and $(\tilde{\gamma },{\tilde{K}}_{{\tilde{T}}_1})$  on some probability space $({\Omega }^{\prime },{\mathcal{ℬ}}^{\prime },{\mathbb{P}}^{\prime })$  such that ${\gamma }_k^{{\delta }_k}\left({\omega }^{\prime }\right)$  converges to $\tilde{\gamma }\left({\omega }^{\prime }\right)$  and $(K_{T_1^k}^k,{\gamma }_k^{{\delta }_k}(T_1^k\left)\right)$  converges to $({\tilde{K}}_{{\tilde{T}}_1},\tilde{\gamma }({\tilde{T}}_1\left)\right)$  for all ${\omega }^{\prime }\in {\Omega }^{\prime }$  . These two observations imply that, if we work with the coupled versions of $({\gamma }_k^{{\delta }_k},K_{T_1^k}^k)$  and $(\tilde{\gamma },{\tilde{K}}_{{\tilde{T}}_1})$  , we are in the same situation as before, but with $D_k$  (resp., $D$  ) replaced by $D_k\backslash K_{T_1^k}^k$  (resp., $D\backslash {\tilde{K}}_{{\tilde{T}}_1}$  ) and $a_k$  (resp., $a$  ) by ${\gamma }_k^{{\delta }_k}\left(T_1^k\right)$  (resp., $\tilde{\gamma }\left({\tilde{T}}_1\right)$  ). As already remarked, $D\backslash {\tilde{K}}_{{\tilde{T}}_1}$  and ${\tilde{G}}_2$  are admissible domains, which allows us to use Theorem  6 (and therefore Lemma  A.1 ), Corollary  B.1 and Lemmas  A.2 - A.3 .

Then the conclusion that $\left(\right(K_{T_1^k}^k,{\gamma }_k^{{\delta }_k}\left(T_1^k\right)),(K_{T_2^k}^k,{\gamma }_k^{{\delta }_k}\left(T_2^k\right)\left)\right)$  converges in distribution to $\left(\right({\tilde{K}}_{{\tilde{T}}_1},\tilde{\gamma }\left({\tilde{T}}_1\right)),({\tilde{K}}_{{\tilde{T}}_2},\tilde{\gamma }\left({\tilde{T}}_2\right)\left)\right)$  follows from the same arguments as before, again using Corollary  B.1 and Lemmas  A.2 - A.3 . In order to get claim (C), in places where the exploration path comes close to the boundary of the past hull we can use the bound on the probability of six crossings of an annulus in the plane (as already seen in the case $k=3$  of the proof of Theorem  1 ), while in places where it comes close to the remaining portionof the boundary (i.e., $\partial D$  or the Jordan arc $\overline{cd}$  in Figure  19 ) we can use Lemmas  A.2 - A.3 .

We can now iterate those same arguments $j$  times, for any $j>1$  . It is in fact easy to see by induction that the domains $D\backslash {\tilde{K}}_{{\tilde{T}}_j}$  and ${\tilde{G}}_j$  that appear in the successive steps are admissible for all $j$  . Therefore we can keep using Theorem  6 (and therefore Lemma  A.1 ), Corollary  B.1 and Lemmas  A.2 - A.3 . If we keep track at each step of the previous ones, in the spirit of Theorem  1 , this provides the joint convergence of all the curves and fillings involved at each step and concludes the proof of Theorem  7 .

Lemma A.2. Let $\left\{\right(D_k,a_k,c_k,d_k\left)\right\}$  be a sequence of domains admissible with respect to $(a_k,c_k,d_k)$  and let ${\gamma }_k^{\delta }$  be the percolation exploration path in $D_k$  started at $a_k$  and stopped when it first hits the counterclockwise arc $J_k^{\prime }=\overline{c_k^{\prime }{d}_k^{\prime }}\subset J_k=\overline{c_k{d}_k}$  of $\partial D_k$  . Assume that, as $k\to \infty$  , $(D_k,a_k,c_k,c_k^{\prime },d_k,d_k^{\prime })$  converges to $(D,a,c,c^{\prime },d,d^{\prime })$  , where $D$  is a domain admissible with respect to $(a,c,d)$  and $J=\overline{c^{\prime }{d}^{\prime }}\subset J=\overline{cd}$  . Let ${\mathcal{ℰ}}_k^{\delta }(J_k;\varepsilon ,{\varepsilon }^{\prime })=\left\{{\cup }_{v\in J_k\backslash J_k^{\prime }}{\mathcal{ℬ}}_k^{\delta }\right(v;\varepsilon ,{\varepsilon }^{\prime }\left)\right\}\cup \left\{{\cup }_{v\in J_k^{\prime }}{\mathcal{A}}_k^{\delta }\right(v;\varepsilon ,{\varepsilon }^{\prime }\left)\right\}$  , where ${\mathcal{A}}_k^{\delta }(v;\varepsilon ,{\varepsilon }^{\prime })$  is the event that ${\gamma }_k^{\delta }$  contains a segment that stays within $B(v,\varepsilon )$  and has a double crossing of the annulus $B(v,\varepsilon )\backslash B(v,{\varepsilon }^{\prime })$  without that segment touching $\partial D_k^{\delta }$  , and ${\mathcal{ℬ}}_k^{\delta }(v;\varepsilon ,{\varepsilon }^{\prime })$  is the event that ${\gamma }_k^{\delta }$  enters $B(v,{\varepsilon }^{\prime })$  , but is stopped outside $B(v,\varepsilon )$  and does not touch $\partial D_k\cap B(v,\varepsilon )$  . Then, for any $\varepsilon >0$  ,

$$\begin{array}{c}{lim}_{{\varepsilon }^{\prime }\to 0}{limsup}_{k\to \infty \delta \to 0}\mathbb{P}\left({\mathcal{ℰ}}_k^{\delta }\right(J_k;\varepsilon ,{\varepsilon }^{\prime }\left)\right)=0.\end{array}$$

(31)

Essentially, this means that as $\delta \to 0$  (and $k\to \infty$  ), it becomes increasingly unlikely that the exploration path ever comes close to $J_k$  without quickly touching $\partial D_k^{\delta }$  nearby.

Lemma  A.2 easily implies the following result, used in Theorem  7 to show that if $\tilde{\gamma }$  touches $\partial D$  , then ${\gamma }_n^{{\delta }_n}$  touches $\partial D_n^{{\delta }_n}$  nearby, for $n$  large enough.

Lemma A.3. With the notation and assumptions of Lemma  A.2 ,

$$\begin{array}{c}{lim}_{\varepsilon \to 0}{lim}_{{\varepsilon }^{\prime }\to 0}{limsup}_{k\to \infty \delta \to 0}\mathbb{P}\left({\bigcup }_{v\in J_k\backslash J_k^{\prime }}{\mathcal{A}}_k^{\delta }\right(v;\varepsilon ,{\varepsilon }^{\prime }\left)\right)=0.\end{array}$$

(32)

Proof. First of all notice that for $v$  in $J_k\backslash J_k^{\prime }$  but not in $B(c^{\prime },\varepsilon )$  and not in $B(d^{\prime },\varepsilon )$  , the events ${\mathcal{A}}_k^{\delta }(v;\varepsilon ,{\varepsilon }^{\prime })$  and ${\mathcal{ℬ}}_k^{\delta }(v;\varepsilon ,{\varepsilon }^{\prime })$  are exactly the same because the exploration path is, by definition, stopped on $J_k^{\prime }$  . Therefore, we only have to prove that the event corresponding to the union over $v\in \{J_k\backslash J_k^{\prime }\}\cap \left\{B\right(c^{\prime },\varepsilon )\cup B(d^{\prime },\varepsilon \left)\right\}$  of ${\mathcal{A}}_k^{\delta }(v;\varepsilon ,{\varepsilon }^{\prime })$  has probability going to zero as $\varepsilon \to 0$  . We already know from Lemma  A.2 that ${\mathcal{ℬ}}_k^{\delta }(v;\varepsilon ,{\varepsilon }^{\prime })$  happens with small probability for those points. This is, however, not sufficient because the exploration path could enter $B(v,{\varepsilon }^{\prime })$  , then exit $B(v,\varepsilon )$  , and then re-enter it and touch $J_k^{\prime }$  inside $B(v,\varepsilon )$  ,which is not an event in ${\mathcal{ℬ}}_k^{\delta }(v;\varepsilon ,{\varepsilon }^{\prime })$  . But such an event would imply that ${\gamma }_k^{\delta }$  first touches $J_k^{\prime }$  inside one of the two balls of radius $\varepsilon$  centered at $c_k^{\prime }$  and $d_k^{\prime }$  , and by an application of Cardy's formula the probability that the latter happens goes to zero as $\varepsilon \to 0$  .

The proof of Lemma  A.2 is partly based on relating the failure of ( 31 ) to the occurrence with strictly positive probability of certain continuum limit “mushroom” events (see Lemma  A.5 ) that we will show must have zero probability because otherwise there would be a contradiction to Lemma  A.4 , which itself is a consequence of the continuity of Cardy's formula with respect to the domain boundary. In both of the next two lemmas, we denote by $\mu$  any subsequence limit of the probability measures for the collection of all colored (blue and yellow) $\mathcal{T}$  -paths on all of ${\mathbb{R}}^2$  , in the Aizenman-Burchard sense (see Remark  2.1 ). We recall that in our notation, $D$  represents an open domain and $\overline{z_1{z}_2}$  , $\overline{z_3{z}_4}$  represent closed segments of its boundary. In Lemma  A.4 below, we restrict attention to a Jordan domain $D$  since that case suffices for the use of Lemma  A.4 in the proof of Lemma  A.2 .

Lemma A.4. For $(D,z_1,z_2,z_3,z_4)$  , with $D$  a Jordan domain, consider the following crossing events, ${\mathcal{C}}_i^{*}={\mathcal{C}}_i^{*}(D,z_1,z_2,z_3,z_4)$  , where $*$  denotes either blue or yellow and $i=1,2,3$  : ${\mathcal{C}}_1^{*}=\{\exists \text{a}*\text{path in the closure}\overline{D}\text{from}\overline{z_1{z}_2}\text{to}\overline{z_3{z}_4}\},$  ${\mathcal{C}}_2^{*}=\{\exists \text{a}*\text{path in}D\text{from the interior of}\overline{z_1{z}_2}\text{to the interior of}\overline{z_3{z}_4}\},$  ${\mathcal{C}}_3^{*}=\{\exists \text{a}*\text{path starting and ending outside}\overline{D}\text{whose restriction to}D\text{is as in}{\mathcal{C}}_2^{*}\}.$  Then $\mu \left({\mathcal{C}}_1^{*}\right)=\mu \left({\mathcal{C}}_2^{*}\right)=\mu \left({\mathcal{C}}_3^{*}\right)={\Phi }_D(z_1,z_2;z_3,z_4)$  .

Proof. The proof is similar to that of Theorem  6 , but easier because $D$  is here a Jordan domain. Indeed, it is enough to construct a new Jordan domain $\tilde{D}\left(\varepsilon \right)$  (with appropriately selected points ${\tilde{z}}_1\left(\varepsilon \right),{\tilde{z}}_2\left(\varepsilon \right),{\tilde{z}}_3\left(\varepsilon \right),{\tilde{z}}_4\left(\varepsilon \right)$  on the boundary and corresponding events ${\tilde{\mathcal{C}}}_i^{*}$  ) such that the occurrence of ${\tilde{\mathcal{C}}}_1^{*}$  in $\tilde{D}\left(\varepsilon \right)$  implies the occurrence of ${\mathcal{C}}_3^{*}$  in $D$  and with $\left(\tilde{D}\right(\varepsilon ),{\tilde{z}}_1(\varepsilon ),{\tilde{z}}_2(\varepsilon ),{\tilde{z}}_3(\varepsilon ),{\tilde{z}}_4(\varepsilon \left)\right)\to (D,z_1,z_2,z_3,z_4)$  as $\varepsilon \to 0$  . The continuity of Cardy's formula (Lemma  B.2 in Appendix  B ) does the rest.

Lemma A.5. For $(D,a,c,d)$  as in Lemma  A.2 , $v\in J\equiv \overline{cd}$  , and $\varepsilon >0$  , we define $U^{yellow}(D,\varepsilon ,v)$  , the yellow “mushroom” event (at $v$  ), to be the event that there is a yellow path in $\overline{D}$  from $v$  to $\partial B(v,\varepsilon )$  and a blue path in $\overline{D}$  , between some pair of distinct points $v_1,v_2$  in $\partial D\cap \left\{B\right(v,\varepsilon /3)\backslash B(v,\varepsilon /8\left)\right\}$  , that passes through $v$  and such that this blue path is between $\partial D$  and the yellow path (see Figure  20 ). We similarly define $U^{blue}(D,\varepsilon ,v)$  with the colors interchanged and $U^{*}(D,\varepsilon ,J)={\cup }_{v\in J}{U}^{*}(D,\varepsilon ,v)$  where $*$  denotes blue or yellow. Then for any deterministic domain $D$  and any $0<\varepsilon <min\left\{\right|a-c|,|a-d\left|\right\}$  , $\mu \left(U^{*}\right(D,\varepsilon ,J\left)\right)=0$  .

Proof. If $\mu \left(U^{*}\right(D,\varepsilon ,J\left)\right)>0$  for some $\varepsilon >0$  , then there is some segment $\overline{a^{\prime }{b}^{\prime }}\subset J$  of $\partial D$  of diameter not larger than $\varepsilon /10$  such that

$$\begin{array}{c}\mu \left({\cup }_{v\in \overline{a^{\prime }{b}^{\prime }}}{U}^{*}\right(D,\varepsilon ,v\left)\right)>0.\end{array}$$

(33)

Choose any point $v_0\in \overline{a^{\prime }{b}^{\prime }}$  and consider the new domain $D^{\prime }$  whose boundary consists of the correctly chosen (as we explain below) segment of the circle $\partial B(v_0,\varepsilon /2)$  between the two points $c^{\prime },d^{\prime }$  where $\partial D$  first hits $\partial B(v_0,\varepsilon /2)$  on either side of $v_0$  , together with the segment from $d^{\prime }$  to $c^{\prime }$  of $\partial D$  (see Figure  21 ).

The correct circle segment between $c^{\prime }$  and $d^{\prime }$  is the (counter) clockwise one if $v_0$  is between $c^{\prime }$  and $d^{\prime }$  along $\partial D$  when $\partial D$  is oriented (counter) clockwise. It is also not hard to see that since $\varepsilon <min\left\{\right|a-c|,|a-d\left|\right\}$  , $D^{\prime }$  is a Jordan domain, so that Lemma  A.4 can be applied. In the new domain $D^{\prime }$  , $\overline{a^{\prime }{b}^{\prime }}$  is the same curve segment as it was in the old domain $D$  , but $\overline{c^{\prime }{d}^{\prime }}$  is now a segment of the circle $\partial B(v_0,\varepsilon /2)$  . It should be clear that

$$\begin{array}{c}{\cup }_{v\in \overline{a^{\prime }{b}^{\prime }}}{U}^{*}(D,\varepsilon ,v)\subset {\mathcal{C}}_1^{*}(D^{\prime },a^{\prime },b^{\prime },c^{\prime },d^{\prime })\backslash {\mathcal{C}}_3^{*}(D^{\prime },a^{\prime },b^{\prime },c^{\prime },d^{\prime })\end{array}$$

(34)

which yields a contradiction of Lemma  A.4 if $\mu \left(U^{*}\right(D,\varepsilon ,J\left)\right)>0$  .

Proof of Lemma  A.2 . We first note that since the probability in ( 31 ) is nonincreasing in $\varepsilon$  , we may assume that $\varepsilon <min\left\{\right|a-c|,|a-d\left|\right\}$  , as requested by Lemma  A.5 .

Let us first consider the simpler case of ${\mathcal{ℬ}}^{{\delta }_k}(v;\varepsilon ,{\varepsilon }^{\prime })$  in which $v\in J_k^{\prime }$  . We follow the exploration process until time $T$  , when it first touches $\partial B(v,{\varepsilon }^{\prime })$  for some $v\in J_k^{\prime }$  , and consider the annulus $B(v,\varepsilon )\backslash B(v,{\varepsilon }^{\prime })$  . Let ${\pi }_Y$  be the leftmost yellow $\mathcal{T}$  -path and ${\pi }_B$  the rightmost blue $\mathcal{T}$  -path in $\Gamma \left({\gamma }_k^{\delta }\right)$  at time $T$  that cross $B(v,\varepsilon )\backslash B(v,{\varepsilon }^{\prime })$  . ${\pi }_Y$  and ${\pi }_B$  split the annulus $B(v,\varepsilon )\backslash B(v,{\varepsilon }^{\prime })$  into three sectors that, for simplicity, we will call the central sector, containing the crossing segment of the exploration path, the yellow (left) sector, with ${\pi }_Y$  as part of its boundary, and the blue (right) sector, the remaining one, with ${\pi }_B$  as part of its boundary.

We then look for a yellow “lateral” crossing within the yellow sector from ${\pi }_Y$  to $\partial D_k$  and a blue lateral crossing within the blue sector from ${\pi }_B$  to $\partial D_k$  . Notice that the yellow sector may contain “excursions” of the exploration path coming off $\partial B(v,\varepsilon )$  , producing nested yellow and blue excursions off $\partial B(v,\varepsilon )$  , and the same for the blue sector. But for topological reasons, those excursions are such that for every group of nested excursions, the outermost one is always yellow in the yellow sector and blue in the blue sector. Therefore, by standard percolation theory arguments, the conditional probability (conditioned on $\Gamma \left({\gamma }_k^{\delta }\right)$  at time $T$  ) to find a yellow lateral crossing of the yellow sector from ${\pi }_Y$  to $\partial D_k$  is bounded below by the probability to find a yellow circuit in an annulus with inner radius ${\varepsilon }^{\prime }$  and outer radius $\varepsilon$  . An analogous statement holds for the conditional probability (conditioned on $\Gamma \left({\gamma }_k^{\delta }\right)$  at time $T$  and also on the entire percolation configuration in theyellow sector) to find a blue lateral crossing of the blue sector from ${\pi }_B$  to $\partial D_k$  . Thus for any fixed $\varepsilon >0$  , by an application of the Russo-Seymour-Welsh lemma [27, 29] , the conditional probability to find both a yellow lateral crossing within the yellow sector from ${\pi }_Y$  to $\partial D_k$  and a blue lateral crossing within the blue sector from ${\pi }_B$  to $\partial D_k$  goes to one as ${\varepsilon }^{\prime }\to 0$  .

But if such yellow and blue crossings are present, the exploration path is forced to touch $J_k^{\prime }$  before exiting $B(v,\varepsilon )$  , and if that happens, the exploration process is stopped, so that it will never exit $B(v,\varepsilon )$  and the union over $v\in J_k^{\prime }$  of ${\mathcal{ℬ}}_k^{\delta }(v;\varepsilon ,{\varepsilon }^{\prime })$  cannot occur.

This concludes the proof of this case.

Let us now consider the remaining case in which $v/\in J_k^{\prime }$  . The basic idea of the proof is then that by straightforward weak convergence and related coupling arguments, the failure of ( 31 ) would imply that some subsequence limit $\mu$  would satisfy $\mu \left(U^{yellow}\right(D,\varepsilon ,J)\cup U^{blue}(D,\varepsilon ,J\left)\right)>0$  , which would contradict Lemma  A.5 . This is essentially because the close approach of an exploration path on the $\delta$  -lattice to $J_k\backslash J_k^{\prime }$  without quickly touching nearby yields one two-sided colored $\mathcal{T}$  -path (the “perimeter” of the portion of the hull of the exploration path seen from a boundary point of close approach) and a one-sided $\mathcal{T}$  -path of the other color belonging to the percolation cluster not seen from the boundary point (i.e., shielded by the two-sided path). Both the two-sided path and the one-sided one are subsets of $\Gamma \left({\gamma }_k^{\delta }\right)$  .

Assume by contradiction that ( 31 ) is false, so that close encounters without touching happen with bounded away from zero probability. Consider for concreteness an exploration path ${\gamma }_k^{\delta }$  that has a close approach to a point $v$  in the counterclockwise arc $\overline{d_k^{\prime }{d}_k}$  . The exploration path may have multiple close approaches to $v$  with differing colors of the perimeter as seen from $v$  , but for topological reasons, the last time the exploration path comes close to $v$  , it must do so in such a way as to produce a yellow $\mathcal{T}$  -path ${\pi }_Y$  (seen from $v$  ) that crosses $B(v,\varepsilon )\backslash B(v,{\varepsilon }^{\prime })$  twice, and a blue path ${\pi }_B$  that crosses it once (see Figure  22 ). This is so because the exploration process that produced ${\gamma }_k^{\delta }$  ended somewhere on $J_k^{\prime }$  (and outside $B(v,\varepsilon )$  ), which is to the right of (i.e., clockwise to) $v$  .

The presence of ${\pi }_Y$  implies that there are a yellow leftmost $\mathcal{T}$  -path ${\pi }_L$  and a yellow rightmost $\mathcal{T}$  -path ${\pi }_R$  (looking at $v$  from inside $D_k$  ) crossing the annulus $B(v,\varepsilon )\backslash B(v,{\varepsilon }^{\prime })$  .

The paths ${\pi }_L$  and ${\pi }_R$  split the annulus $B(v,\varepsilon )\backslash B(v,{\varepsilon }^{\prime })$  into three sectors, that we will call the central sector, containing ${\pi }_B$  , the left sector, with ${\pi }_L$  as part of its boundary, and the right sector, with ${\pi }_R$  as part of its boundary. Again for topological reasons, all other monochromatic crossings of the annulus are contained in the central sector, including at least one blue path ${\pi }_B$  . As in the previous case, the left and right sectors can contain nested monochromatic excursions off $\partial B(v,\varepsilon )$  , but this time for every group of excursions, the outermost one is yellow in both sectors. Now consider the annulus $B(v,\varepsilon /3)\backslash B(v,\varepsilon /8)$  . We look for a yellow lateral crossing within the left sector from ${\pi }_L$  to $\partial D_k$  and a yellow lateral crossing within the right sector from ${\pi }_R$  to $\partial D_k$  . Since the outermost excursions in both sectors are yellow, the conditional probability to find a yellow lateral crossing within the left sector from ${\pi }_L$  to $\partial D_k$  is bounded below by the probability to find a yellow circuit in an annulus with inner radius ${\varepsilon }^{\prime }$  and outer radius $\varepsilon$  , and an analogous statement holds for the conditional probability to find a yellow lateral crossing within the right sector from ${\pi }_R$  to $\partial D_k$  . Thus for any fixed $\varepsilon >0$  , by an application of the Russo-Seymour-Welsh lemma [27, 29] , the conditional probability to find both yellow lateral crossings goes to one as ${\varepsilon }^{\prime }\to 0$  . But the presence of such yellow crossings would produce a (blue) mushroom event, leading to a contradiction with Lemma  A.5 .

We are finally ready to prove the main result of this section which implies statement (S) at the beginning of Section  5 .

Corollary A.1. Consider a sequence $\left\{\right(D_k,a_k,b_k\left)\right\}$  of Jordan domains with two distinct selected points $a_k,b_k$  on their boundaries $\partial D_k$  . Assume that $(D_k,a_k,b_k)\to (D,a,b)$  , where $D$  is a Jordan domain with two distinct selected points on its boundary $\partial D$  . Denote by ${\gamma }_k^{\delta }\equiv {\gamma }_{D_k,a_k,b_k}^{\delta }$  the percolation exploration path inside (the $\delta$  -approximation $D_k^{\delta }$  of ) $D_k$  from (the e-vertex closest to) $a_k$  to (the e-vertex closest to) $b_k$  . Then for any sequence ${\delta }_k\downarrow 0$  , as $k\to \infty$  , ${\gamma }_k^{{\delta }_k}$  converges in distribution to the trace ${\gamma }_{D,a,b}$  of chordal $SL{E}_6$  inside $D$  from $a$  to $b$  .

Proof. It follows from [2] that ${\gamma }_k^{{\delta }_k}$  converges in distribution along subsequence limits $k_n$  .

Since we have proved that the filling of any such subsequence limit $\tilde{\gamma }$  satisfies the spatial Markov property (Theorem  7 ) and the exit distribution of $\tilde{\gamma }$  is determined by Cardy's formula (Theorem  6 ), we can deduce from Theorem  5 that the limit is unique and that the law of ${\gamma }_k^{{\delta }_k}$  converges, as $k\to \infty$  , to the law of the trace ${\gamma }_{D,a,b}$  of chordal $SL{E}_6$  inside $D$  from $a$  to $b$  .

Appendix B: Sequences of Conformal Maps In this appendix, we give some results about sequences of conformal maps that are used in various places throughout the paper. For more details, the interested reader should consult [24] .

Definition B.1. (see Section 1.4 of [24] ) Let $w_0\in \mathbb{C}$  be given and let $\left\{G_n\right\}$  be domains with $w_0\in G_n\subset \mathbb{C}$  . We say that $G_n\to G$  as $n\to \infty$  with respect to $w_0$  in the sense of kernel convergence if

• 1. either $G=\left\{w_0\right\}$  , or else $G$  is a domain $\ne \mathbb{C}$  with $w_0\in G$  such that some neighborhood of every $w\in G$  lies in $G_n$  for large $n$  ; and
• 2. for $w\in \partial G$  there exist $w_n\in \partial G_n$  such that $w_n\to w$  as $n\to \infty$  .

It is clear from the definition that every subsequence limit also converges to $G$  and it is also easy to see that the limit is uniquely determined. With this definition we can now state Carathéodory's kernel theorem [9] .

Theorem 8. (see Theorem 1.8 of [24] ) Let $f_n$  map $\mathbb{D}$  conformally onto $G_n$  with $f_n\left(0\right)=w_0$  and $f_n^{\prime }\left(0\right)>0$  . If $G=\left\{w_0\right\}$  , let $f\left(z\right)\equiv w_0$  ; otherwise let $f$  map $\mathbb{D}$  conformally onto $G$  with $f\left(0\right)=w_0$  and $f^{\prime }\left(0\right)>0$  . Then, as $n\to \infty$  , $f_n\to f$  locally uniformly in $\mathbb{D}$  if and only if $G_n\to G$  with respect to $w_0$  .

The next result, Radó's theorem [25] , deals with sequences of Jordan domains and is used in the main body of the paper. In this case the conformal maps have a continuous extension to $\mathbb{D}\cup \partial \mathbb{D}$  .

Theorem 9. (see Theorem 2.11 of [24] ) For $n=1,2,\dots$  , let $J_n$  and $J$  be Jordan curves parametrized respectively by ${\phi }_n\left(t\right)$  and $\phi \left(t\right)$  , $t\in [0,1]$  , and let $f_n$  and $f$  be conformal maps from $\mathbb{D}$  onto the inner domains of $J_n$  and $J$  such that $f_n\left(0\right)=f\left(0\right)$  and $f_n^{\prime }\left(0\right)>0$  , $f^{\prime }\left(0\right)>0$  for all $n$  . If ${\phi }_n\to \phi$  as $n\to \infty$  uniformly in $[0,1]$  then $f_n\to f$  as $n\to \infty$  uniformly in $\overline{\mathbb{D}}$  .

The type of convergence of sequences of Jordan domains $\left\{G_n\right\}$  to a Jordan domain $G$  encountered in the main body of the paper (i.e., in the sense that $\partial G_n$  converges, as $n\to \infty$  , to $\partial G$  in the uniform metric ( 2 ) on continuous curves) is clearly sufficient to apply Theorem  9 . In Appendix  A , however, we have to deal with domains that are not Jordan, and therefore we cannot use Radó's theorem. The tools needed to deal with those situations are described below.

Definition B.2. (see Section 2.2 of [24] ) The closed set $A\subset \mathbb{C}$  is called locally connected if for every $\varepsilon >0$  there is $\delta >0$  such that, for any two points $a,b\in A$  with $|a-b|<\delta$  , we can find a continuum $B$  with diameter smaller than $\varepsilon$  and with $a,b\in B\subset A$  .

In the definition above, a continuum denotes a compact connected set with more than one point. We remark that every continuous curve (with more than one point) is a locally connected continuum (the converse is also true: every locally connected continuum is a curve). The concept of local connectedness gives a topological answer to the problem of global extension of a conformal map to the domain boundary, as follows.

Theorem 10. (see Continuity Theorem in Section 2.1 of [24] ) Let $f$  map the unit disk $\mathbb{D}$  conformally onto $G\subset \mathbb{C}\cup \{\infty \}$  . Then the function $f$  has a continuous extension to $\mathbb{D}\cup \partial \mathbb{D}$  if and only if $\partial G$  is locally connected.

When $f$  has a continuous extension to $\mathbb{D}\cup \partial \mathbb{D}$  , we do not distinguish between $f$  and its extension. This is always the case for the conformal maps considered in this paper. The problem wether this extension is injective on $\overline{\mathbb{D}}$  has also a topological answer, as follows.

Theorem 11. (see Carathhéodory Theorem in Section 2.1 of [24] ) In the notation of Theorem  10 , the function $f$  has a continuous and injective extension if and only if $\partial G$  is a Jordan curve.

When considering sequences of domains whose boundaries are locally connected the following definition is useful.

Definition B.3. (see Section 2.2 of [24] ) The closed sets $A_n\subset \mathbb{C}$  are uniformly locally connected if, for every $\varepsilon >0$  , there exists $\delta >0$  independent of $n$  such that any two points $a_n,b_n\in A_n$  with $|a_n-b_n|<\delta$  can be joined by continua $B_n\subset A_n$  of diameter smaller than $\varepsilon$  .

The convergence of domains used in this paper (i.e., $G_n\to G$  if $\partial G_n\to \partial G$  in the uniform metric ( 2 ) on continuous curves) clearly implies kernel convergence, which immediately allows us to use Theorem  8 . However, we need uniform convergence in $\overline{\mathbb{D}}$  .

This is guaranteed by Radó's theorem in the case of Jordan domains; in the non-Jordan case, sufficient conditions to have uniform convergence are stated in the next theorem.

Theorem 12. (see Corollary 2.4 of [24] ) Let $\left\{G_n\right\}$  be a sequence of bounded domains such that, for some $0<r<R<\infty$  , $B(0,r)\subset G_n\subset B(0,R)$  for all $n$  and such that $\{\mathbb{C}\backslash G_n\}$  is uniformly locally connected. Let $f_n$  map $\mathbb{D}$  conformally onto $G_n$  with $f_n\left(0\right)=0$  . If $f_n\left(z\right)\to f\left(z\right)$  as $n\to \infty$  for each $z\in \mathbb{D}$  , then the convergence is uniform in $\overline{\mathbb{D}}$  .

In order to use Theorem  12 in Appendix  A we need the following lemma. The definitions of admissible domain and the related notion of convergence are given just before Theorem  6 in Appendix  A .

Lemma B.1. Let $\left\{\right(G_n,a_n,c_n,d_n\left)\right\}$  be a sequence of domains admissible with respect to $(a_n,c_n,d_n)$  and assume that, as $n\to \infty$  , $(G_n,a_n,c_n,d_n)\to (G,a,c,d)$  , where $G$  is a domain admissible with respect to $(a,c,d)$  . Then the sequence of closed sets $\{\mathbb{C}\backslash G_n\}$  is uniformly locally connected.

Proof. In order to prove the lemma, we claim that it suffices to focus on pairs of points on the boundaries, i.e., to show that: $(*)$  for every $\varepsilon >0$  , there exists $\delta =\delta \left(\varepsilon \right)>0$  independent of $n$  such that any two points $u_n,v_n\in \partial G_n$  with $|u_n-v_n|<\delta$  can be joined by a continuum of diameter $<\varepsilon$  contained in the complement $\mathbb{C}\backslash G_n$  of $G_n$  .

To verify our claim, let us assume $(*)$  for the moment, and consider two points $u_n,v_n\in \mathbb{C}\backslash G_n$  (but not necessarily in $\partial G_n$  ) with $|u_n-v_n|<{\delta }^{\prime }$  , where ${\delta }^{\prime }=min\left\{\frac{1}{3}\delta \right(\frac{\varepsilon }{3}),\frac{\varepsilon }{3}\}$  .

If (at least) one of the two points, say $u_n$  , is at distance greater than ${\delta }^{\prime }$  from $\partial G_n$  , then we can connect $u_n$  and $v_n$  using the closed ball of radius ${\delta }^{\prime }$  centered at $u_n$  , since $v_n\in \overline{B(u_n,{\delta }^{\prime })}\subset \mathbb{C}\backslash G_n$  . If both points are at distance smaller than ${\delta }^{\prime }$  from $\partial G_n$  , we can connect each point to a closest point on $\partial G_n$  by a straight segment of length smaller than ${\delta }^{\prime }$  . Those two points on $\partial G_n$  can then be connected to each other by a continuum $B_n$  of diameter $<\varepsilon /3$  contained in $\mathbb{C}\backslash G_n$  , and the union of $B_n$  with the two straight segments gives a continuum of diameter $<\varepsilon$  connecting $u_n$  with $v_n$  and contained in $\mathbb{C}\backslash G_n$  .

We now prove $(*)$  . Since $\partial G_n\to \partial G$  in the uniform metric ( 2 ) on continuous curves, for every $\varepsilon >0$  there exists $n_0=n_0\left(\varepsilon \right)$  such that for all $n\ge n_0$  , $\text{d}(\partial G_n,\partial G)<\varepsilon$  . The admissibility of $G$  implies that we can split its boundary into three Jordan arcs, $J_1=\overline{da}$  , $J_2=\overline{ac}$  , $J_3=\overline{cd}$  , such that $J_3$  does not touch the interior of either $J_1$  or $J_2$  . We can do the same with $\partial G_n$  , letting $J_{1,n}=\overline{d_n{a}_n}$  , $J_{2,n}=\overline{a_n{c}_n}$  and $J_{3,n}=\overline{c_n{d}_n}$  . Let ${\phi }_{i,n}\left(t\right)$  and ${\phi }_i\left(t\right)$  , $t\in [0,1]$  , $i=1,2,3$  be parametrizations of $J_{i,n}$  and $J_i$  respectively, with ${sup}_{t\in [0,1]}\left|{\phi }_{i,n}\right(t)-{\phi }_i(t\left)\right|<\varepsilon$  for $n\ge n_0$  and $i=1,2,3$  .

Let us assume, by contradiction, that $(*)$  is false. Then there are indices $k$  (actually $n_k$  , but we abuse notation a bit) and points $u_k,v_k\in \partial G_k$  with $|u_k-v_k|\to 0$  (as $k\to \infty$  ) that cannot be joined by a continuum of diameter $<\varepsilon$  contained in $\mathbb{C}\backslash G_k$  . By compactness considerations, we may assume that $u_k\to u$  and $v_k\to v$  as $k\to \infty$  , with $u=v$  . Suppose that $u_k$  and $v_k$  belong to the interior of the same Jordan arc $J_{i,k}$  for all $k$  large enough.

Let $u_k={\phi }_{i,k}\left({\tau }_k\right)$  , $v_k={\phi }_{i,k}\left({\tau }_k^{\prime }\right)$  , $u={\phi }_i\left(\tau \right)$  and $v={\phi }_i\left({\tau }^{\prime }\right)$  . It follows that ${\tau }_k\to \tau$  and ${\tau }_k^{\prime }\to {\tau }^{\prime }$  , and since $J_i$  is a Jordan arc, $\tau ={\tau }^{\prime }$  . For $k$  large enough, the function ${\phi }_{i,k}$  maps the closed segment of $[0,1]$  between $\tau$  and ${\tau }^{\prime }$  onto a continuum in $J_{i,k}$  containing $u_k$  and $v_k$  whose diameter tends to zero as $k\to \infty$  , leading to a contradiction with our assumption.

Similar reasoning gives a contradiction if $u_k$  and $v_k$  both belong to $J_{1,k}\cup J_{3,k}$  or both belong to $J_{2,k}\cup J_{3,k}$  for all $k$  large enough, since the concatenation of $J_{1,k}$  with $J_{3,k}$  or of $J_{2,k}$  with $J_{3,k}$  is still a Jordan arc. The above reasoning applies except when $u(=v)$  is on both $J_1$  and $J_2$  . When $u=v=a$  , one can paste together small Jordan arcs on $J_{1,k}$  and $J_{2,k}$  to get a suitable continuum leading to a contradiction. The sole remaining case is when for all $k$  large enough, $u_k$  belongs to the interior of $J_{1,k}$  and $v_k$  belongs to the interior of $J_{2,k}$  .

(Notice that we are ignoring the “degenerate” case in which $c=d$  coincides with the “last” [from $a$  ] double-point on $\partial G$  , and $J_3$  is a simple loop. In that case $u_k$  and $v_k$  could converge to $u=v=c=d\in J_1\cap J_2$  and $u_k$  or $v_k$  could still belong to $J_{3,k}$  for arbitrarily large $k$  's. However, in that case one can find two distinct points on $J_3$  , $c^{\prime }$  and $d^{\prime }$  , such that $D$  is admissible with respect to $(a,c^{\prime },d^{\prime })$  , and points $c_k^{\prime }$  and $d_k^{\prime }$  on $J_{3,k}$  converging to $c^{\prime }$  and $d^{\prime }$  respectively, and define accordingly new Jordan arcs, $J_1^{\prime },J_2^{\prime },J_3^{\prime }$  and $J_{1,k}^{\prime },J_{2,k}^{\prime },J_{3,k}^{\prime }$  , so that $u_k\in J_{1,k}^{\prime }$  and $v_k\in J_{2,k}^{\prime }$  for $k$  large enough. We assume that this has been done if necessary, and for simplicity of notation drop the primes.) In this case let $\left[u_k{v}_k\right]$  denote the closed straight line segment in the plane between $u_k$  and $v_k$  . Imagine that $\left[u_k{v}_k\right]$  is oriented from $u_k$  to $v_k$  and let $v_k^{\prime }$  be the first point of $J_{2,k}$  intersected by $\left[u_k{v}_k\right]$  and $u_k^{\prime }$  be the previous intersection of $\left[u_k{v}_k\right]$  with $\partial G_k$  . Clearly, $u_k^{\prime }/\in J_{2,k}$  . For $k$  large enough, $u_k^{\prime }$  cannot belong to $J_{3,k}$  either, or otherwise in the limit $k\to \infty$  , $J_3$  would touch the interior of $J_1$  and $J_2$  . We deduce that for all $k$  large enough, $u_k^{\prime }\in J_{1,k}$  . Since $J_{1,k}$  and $J_{2,k}$  are continuous curves and therefore locally connected, $u_k$  and $u_k^{\prime }$  belong to a continuum $B_{1,k}$  contained in $J_{1,k}$  whose diameter goes to zero as $k\to \infty$  , and the same for $v_k$  and $v_k^{\prime }$  (with $B_{1,k}$  and $J_{1,k}$  replaced by $B_{2,k}$  and $J_{2,k}$  ).

Since the interior of $\left[u_k^{\prime }{v}_k^{\prime }\right]$  does not intersect any portion of $\partial G_k$  , it is either contained in $G_k$  or in its complement $\mathbb{C}\backslash G_k$  . If $\left[u_k^{\prime }{v}_k^{\prime }\right]\subset \mathbb{C}\backslash G_k$  , we have a contradiction since the union of $\left[u_k^{\prime }{v}_k^{\prime }\right]$  with $B_{1,k}$  and $B_{2,k}$  is contained in $\mathbb{C}\backslash G_k$  and is a continuum containing $u_k$  and $v_k$  whose diameter goes to zero as $k\to \infty$  .

If the interior of $\left[u_k^{\prime }{v}_k^{\prime }\right]$  is contained in $G_k$  , let us consider a conformal map $f_k$  from $\mathbb{D}$  onto $G_k$  . Since $\partial G_k$  is locally connected, the conformal map $f_k$  extends continuously to the boundary of the unit disc. Let $u_k^{\prime }=f_k\left(u_k^{*}\right)$  , $v_k^{\prime }=f_k\left(v_k^{*}\right)$  , $a_k=f_k\left(a_k^{*}\right)$  , $c_k=f_k\left(c_k^{*}\right)$  and $d_k=f_k\left(d_k^{*}\right)$  . The points $c_k^{*},d_k^{*},u_k^{*},a_k^{*},v_k^{*}$  are in counterclockwise order on $\partial \mathbb{D}$  , so that any curve in $\mathbb{D}$  from $a_k^{*}$  to the counterclockwise arc $\overline{c_k^{*}{d}_k^{*}}$  must cross the curve from $u_k^{*}$  to $v_k^{*}$  whose image under $f_k$  is $\left[u_k^{\prime }{v}_k^{\prime }\right]$  . This implies that any curve in $G_k$  going from $a_k$  to the counterclockwise arc $\overline{c_k{d}_k}$  of $\partial G_k$  must cross the (interior of the) line segment $\left[u_k^{\prime }{v}_k^{\prime }\right]$  .

Then, in the limit $k\to \infty$  , any curve in $G$  from $a$  to the counterclockwise arc $\overline{cd}$  must contain the limit point $u={lim}_{k\to \infty }{u}_k^{\prime }={lim}_{k\to \infty }{v}_k^{\prime }=v$  . On the other hand, except for its starting and ending point, any such curve is completely contained in $G$  , which implies that either $u=v=a$  or else that (in the limit $k\to \infty$  ) the counterclockwise arc $\overline{cd}$  is the single point at $u=v=c=d$  . We have already dealt with the former case. In the latter case, one can paste together small Jordan arcs from $u_k^{\prime }$  to $d_k$  , from $d_k$  to $c_k$  , and from $c_k$  to $v_k^{\prime }$  , and take the union with $B_{1,k}$  and $B_{2,k}$  (defined above) to get a suitable continuum in $\mathbb{C}\backslash G_k$  containing $u_k$  and $v_k$  , leading to a contradiction. This concludes the proof.

Theorem  12 , together with Theorem  8 and Lemma  B.1 , implies the following result, which is used in Appendix  A .

Corollary B.1. With the notation and assumptions of Lemma  B.1 (and also assuming that $G_n$  and $G$  contain the origin), let $f_n$  map $\mathbb{D}$  conformally onto $G_n$  with $f_n\left(0\right)=0$  and $f_n^{\prime }\left(0\right)>0$  , and $f$  map $\mathbb{D}$  conformally onto $G$  with $f\left(0\right)=0$  and $f^{\prime }\left(0\right)>0$  . Then, as $n\to \infty$  , $f_n\to f$  uniformly in $\overline{\mathbb{D}}$  .

Proof. As already remarked, the convergence of $\partial G_n$  to $\partial G$  in the uniform metric ( 2 ) on continuous curves (which is part of the definition of $(G_n,a_n,c_n,d_n)\to (G,a,c,d)$  ) easily implies that the conditions in Carathéodory's kernel theorem (Theorem  8 ) are satisfied and therefore that $f_n$  converges to $f$  locally uniformly in $\mathbb{D}$  , as $n\to \infty$  . By an application of Lemma  B.1 , the sequence $\{\mathbb{C}\backslash D_n\}$  is uniformly locally connected, so that we can apply Theorem  12 to conclude that, as $n\to \infty$  , $f_n$  converges to $f$  uniformly in $\overline{\mathbb{D}}$  .

We conclude this appendix with a simple lemma, used in the proof of Theorem  6 , about the continuity of Cardy's formula with respect to the shape of the domain and the positions of the four points on the boundary.

Lemma B.2. For $\left\{\right(D_n,a_n,b_n,c_n,d_n\left)\right\}$  and $(D,a,b,c,d)$  as in Theorem  6 , let ${\Phi }_n$  denote Cardy's formula (see ( 16 )) for a crossing inside $D_n$  from the counterclockwise segment $\overline{a_n{c}_n}$  of $\partial D_n$  to the counterclockwise segment $\overline{b_n{d}_n}$  of $\partial D_n$  and $\Phi$  the corresponding Cardy's formula for the limiting domain $D$  .

Then, as $n\to \infty$  , ${\Phi }_n\to \Phi$  .

Proof. Let $f_n$  be the conformal map that takes $\mathbb{D}$  onto $D_n$  with $f_n\left(0\right)=0$  and $f_n^{\prime }\left(0\right)>0$  , and let $f$  denote the conformal map from $\mathbb{D}$  onto $D$  with $f\left(0\right)=0$  and $f^{\prime }\left(0\right)>0$  ; let $z_1=f^{-1}\left(a\right)$  , $z_2=f^{-1}\left(c\right)$  , $z_3=f^{-1}\left(b\right)$  , $z_4=f^{-1}\left(d\right)$  , $z_1^n=f_n^{-1}\left(a_n\right)$  , $z_2^n=f_n^{-1}\left(c_n\right)$  , $z_3^n=f_n^{-1}\left(b_n\right)$  , and $z_4^n=f_n^{-1}\left(d_n\right)$  . We can apply Corollary  B.1 to conclude that, as $n\to \infty$  , $f_n$  converges to $f$  uniformly in $\overline{\mathbb{D}}$  . This, in turn, implies that, as $n\to \infty$  , $z_1^n\to z_1$  , $z_2^n\to z_2$  , $z_3^n\to z_3$  , and $z_4^n\to z_4$  .

Cardy's formula for a crossing inside $D_n$  from the counterclockwise segment $\overline{a_n{c}_n}$  of $\partial D_n$  to the counterclockwise segment $\overline{b_n{d}_n}$  of $\partial D_n$  is given by

$$\begin{array}{c}{\Phi }_n=\frac{\Gamma (2/3)}{\Gamma (4/3)\Gamma (1/3)}{{\eta }_n^{1/3}}_2{F}_1(1/3,2/3;4/3;{\eta }_n),\end{array}$$

(35)

where

$$\begin{array}{c}{\eta }_n=\frac{(z_1^n-z_2^n)(z_3^n-z_4^n)}{(z_1^n-z_3^n)(z_2^n-z_4^n)}.\end{array}$$

(36)

Because of the continuity of ${\eta }_n$  in $z_1^n$  , $z_2^n$  , $z_3^n$  , $z_4^n$  , and the continuity of Cardy's formula ( 35 ) in ${\eta }_n$  , the convergence of $z_1^n\to z_1$  , $z_2^n\to z_2$  , $z_3^n\to z_3$  and $z_4^n\to z_4$  immediately implies the convergence of ${\Phi }_n$  to $\Phi$  .

Acknowledgements. We are grateful to Greg Lawler, Oded Schramm and Wendelin Werner for various interesting and useful conversations and to Stas Smirnov for communications about a paper in preparation. We note in particular that a discussion with Oded Schramm at the November 2004 Northeast Probability Seminar at the CUNY Graduate Center, about dependence of exploration paths with respect to small changes of domain boundaries, pointed us in a direction that eventually led to Lemmas  A.2 - A.5 . F. C. thanks Wendelin Werner for an invitation to Université Paris-Sud 11, and Vincent Beffara and Luiz Renato Fontes for many helpful discussions. We are especially grateful to Vincent Beffara for pointing out a gap in a preliminary version of Appendix  A . We thank Lai-Sang Young for comments about the presentation of our results. F. C. acknowledges the kind hospitality of the Courant Institute where part of this work was completed.

References

1. M. Aizenman, Scaling limit for the incipient spanning clusters, in Mathematics of Multiscale Materials; the IMA Volumes in Mathematics and its Applications (K. Golden, G. Grimmett, R. James, G. Milton and P. Sen, eds.), Springer (1998).
2. M. Aizenman and A. Burchard, Hölder regularity and dimension bounds for random curves, Duke Math. J. 99, 419-453 (1999).
3. M. Aizenman, B. Duplantier and A. Aharony, Connectivity exponents and the external perimeter in $2D$  independent percolation, Phys. Rev. Lett. 83, 1359-1362 (1999).
4. A. A. Belavin, A. M. Polyakov and A. B. Zamolodchikov, Infinite conformal symmetry of critical fluctuations in two dimensions, J. Stat. Phys. 34, 763-774 (1984).
5. A. A. Belavin, A. M. Polyakov and A. B. Zamolodchikov, Infinite conformal symmetry in two-dimensional quantum field theory, Nucl. Phys. B 241, 333-380 (1984).
6. P. Billingsley, Weak Convergence of Measures: Applications in Probability, Society for Industrial and Applied Mathematics, Philadelphia (1971).
7. F. Camia and C. M. Newman, Continuum Nonsimple Loops and 2D Critical Percolation, J. Stat. Phys. 116, 157-173 (2004).
8. F. Camia, C. M. Newman and V. Sidoravicius, Cardy's formula for some dependent percolation models, Bull. Brazilian Math. Soc. 33, 147-156 (2002).
9. C. Carathéodory, Untersuchungen über die konformen Abbildungen von festen und veränderlichen Gebieten, Math. Ann. 72, 107-144 (1912).
10. J. L. Cardy, Critical percolation in finite geometries, J. Phys. A 25, L201-L206 (1992).
11. J. Cardy, Lectures on Conformal Invariance and Percolation, available at arXiv:math-ph/0103018 (2001).
12. R. Durrett, Probability: Theory and Examples, Third Edition, Duxbury Advanced Series (2004).
13. G. R. Grimmett, Percolation, second edition, Springer, Berlin (1999).
14. W. Kager and B. Nienhuis, A Guide to Stochastic Löwner Evolution and Its Applications, J. Phys. A 115, 1149-1229 (2004).
15. H. Kesten, Percolation Theory for Mathematicians, Birkhäuser, Boston (1982).
16. H. Kesten, V. Sidoravicius and Y. Zhang, Almost all words are seen in critical site percolation on the triangular lattice, Electr. J. Probab. 3, paper no. 10 (1998).
17. G. Lawler, Conformally Invariant Processes, in Lecture notes for the 2002 ICTP School and Conference on Probability, ICTP Lecture Notes Series, Vol. XVII (2004).
18. G. F. Lawler, Conformally Invariant Processes in the Plane, to be published by the American Mathematical Society (2005).
19. G. Lawler, O. Schramm and W. Werner, Values of Brownian intersection exponents I: Half-plane exponents, Acta Math. 187, 237-273 (2001).
20. G. Lawler, O. Schramm and W. Werner, One arm exponent for critical 2D percolation, Electronic J. Probab. 7, paper no. 2 (2002).
21. G. Lawler, O. Schramm and W. Werner, Conformal invariance of planar loop-erased random walk and uniform spanning trees, Ann. Prob. 32, 939-995 (2004).
22. G. Lawler, O. Schramm and W. Werner, Conformal restriction: the chordal case, J. Amer. Math. Soc. 16, 917-955 (2003).
23. A. M. Polyakov, Conformal symmetry of critical fluctuations, JETP Letters 12, 381-383 (1970).
24. Ch. Pommerenke, Boundary Behaviour of Conformal Maps, Springer-Verlag, Berlin (1992).
25. T. Radó, Sur la représentation conforme de domaines variables, Acta Sci. Math. (Szeged) 1, 180-186 (1923).
26. S. Rohde and O. Schramm, Basic properties of SLE, Ann. Math., to appear, preprint arXiv:math.PR0106036 (2003).
27. L. Russo, A note on percolation, Z. Wahrsch. Ver. Geb. 43, 39-48 (1978).
28. O. Schramm, Scaling limits of loop-erased random walks and uniform spanning trees, Israel J. Math. 118, 221-288 (2000).
29. P. D. Seymour, D. J. A. Welsh, Percolation probabilities on the square lattice, in Advances in Graph Theory (B. Bollobás ed.), Annals of Discrete Mathematics 3, North-Holland, Amsterdam, pp. 227-245 (1978).
30. S. Smirnov, Critical percolation in the plane: Conformal invariance, Cardy's formula, scaling limits, C. R. Acad. Sci. Paris 333, 239-244 (2001).
31. S. Smirnov, Critical percolation in the plane. I. Conformal invariance and Cardy's formula. II. Continuum scaling limit. (long version of [30] , dated Nov. 15, 2001), available at http://www.math.kth.se/ $\sim$  stas/papers/index.html.
32. S. Smirnov, in preparation.
33. S. Smirnov, private communication.
34. S. Smirnov and W. Werner, Critical exponents for two-dimensional percolation, Math. Rev. Lett. 8, 729-744 (2001).
35. W. Werner, Critical exponents, conformal invariance and planar Brownian motion, in Proceedings of the 3rd Europ. Congress of Math., Prog. Math., Vol. 202, 87-103 (2001).
36. W. Werner, SLEs as boundaries of clusters of Brownian loops, C. R. Math. Acad. Sci. Paris 337, 481–486 (2003).
37. W. Werner, Random planar curves and Schramm-Loewner Evolutions, in Lectures on probability theory and statistics, Lecture Notes in Math., Vol. 1840, Springer, Berlin, pp. 107-195 (2004).