The research of the first author was supported by NSA grant H98230-04-1-0039. The research of the second author was supported by NSF grant DMS-0111298.
<ph f="cmbx">On Chordal and Bilateral SLE in multiply connected domains</ph>

### Roland M. Friedrich

Department of Mathematics, University of Illinois at Urbana-Champaign, 1409 West Green Street, Urbana, IL 61801, USA E-mail address : rbauer@math.uiuc.edu Institute for Advanced Study, Princeton, NJ 08540, USA E-mail address : rolandf@ias.edu
• Abstract. We discuss the possible candidates for conformally invariant random non-self-crossing curves which begin and end on the boundary of a multiply connected planar domain, and which satisfy a Markovian-type property. We consider both, the case when the curve connects a boundary component to itself (chordal), and the case when the curve connects two different boundary components (bilateral). We establish appropriate extensions of Loewner's equation to multiply connected domains for the two cases. We show that a curve in the domain induces a motion on the boundary and that this motion is enough to first recover the motion of the moduli of the domain and then, second, the curve in the interior. For random curves in the interior we show that the induced random motion on the boundary is not Markov if the domain is multiply connected, but that the random motion on the boundary together with the random motion of the moduli forms a Markov process. In the chordal case, we show that this Markov process satisfies Brownian scaling and discuss how this limits the possible conformally invariant random non-self-crossing curves. We show that the possible candidates are labeled by a real constant and a function homogeneous of degree minus one which describes the interaction of the random curve with the boundary. We show that the random curve has the locality property if the interaction term vanishes and the real parameter equals six.

1 Introduction

In this paper we discuss the possible candidates for a mathematically rigorous notion of conformally invariant random non-self-crossing curves which begin and end on the boundary of a multiply connected planar domain, and which satisfy a Markovian-type property. The Markovian type property means that the random curves can be developed dynamically as a (locally) growing family of random compacts. We aim to proceed in the spirit of Schramm, who deduced that, under an additional reflection symmetry, there is only a one parameter family of such random curves in simply connected domains, which he termed Stochastic Loewner Evolutions, see [24. As such conformally invariant random growing compacts are conjectured to arise as scaling limits of interfaces of 2-dimensional statistical mechanical systems at criticality, Schramm had with one stroke identified what those limits can be. This has many consequences and applications, see [14, [15, [16, [25, and references therein.
Statistical mechanical systems have been studied in discrete approximations of multiply connected domains and Riemann surfaces, see [13, and [2, and the connections with conformal field theory (CFT) indicate that the stochastic Loewner evolution should also extend to multiply connected domains and Riemann surfaces.
For multiply connected domains the situation is already more subtle when compared to the simply connected case, because moduli spaces enter the picture and, as we will show, one has to consider interactions with these moduli.
Families of random compacts from the boundary to the boundary now come in two flavors, as the random compact may grow to either connect a boundary component to itself (the chordal case) or it may grow to connect two different boundary components. We call the latter the bilateral case.
The radial case, treated in [4, where the random compact grows from the boundary to an interior point, can be considered as a limit of the bilateral case, when the boundary component the random compact grows towards shrinks to a point. This can be made precise, see [10.
Our procedure rests on an appropriate extension of Loewner's equation to the multiply connected case. In the simply connected case, Loewner's equation allows to encode a simple curve in a domain $D$  which has one endpoint on the boundary $\partial D$  by a continuous motion on the boundary, see [17. In the multiply connected case, we show in Theorem  3.1 and Theorem  3.2 that a simple curve induces a motion on the boundary of the domain. To recover the curve inside the domain requires also the knowledge of the moduli $\mathbf{M}$  (which describe the conformal equivalence class), as the curve grows. We show in Theorem  4.1 and Theorem  4.2 that these moduli can be recovered from the boundary motion and thus, once these moduli have been obtained, the curve in the interior itself.
A growing random non-self-crossing curve in a multiply connected domain can then also be encoded into a random motion $\xi \left(t\right)$  on the boundary. However, if the connectivity is greater than one, then $\xi$  cannot be Markov. We show in Section  5 that in the chordal case the boundary motion $\xi$  together with the motion of moduli $\mathbf{M}$  is a Markov process, and that it satisfies Brownian scaling.
These facts dramatically reduce the number of possible diffusions.
Indeed, in addition to a real parameter $\kappa$  one is only free to choose a function $A$  which is homogeneous of degree minus one in the variables $\xi$  and $\mathbf{M}$  . The term $A$  measures the interaction of the random growing compact with the boundary (for example if it is desired that the random set avoids the interior boundary components).
$\text{SLE}\left(\kappa ,\rho \right)$  , see [15, and [8, also fits naturally into this framework.
There, the random compact grows into the upper half-plane, the boundary is the real axis, and the interaction is with a finite number of points on the real axis and given in terms of the simplest homogeneous function of degree minus one, $1/x$  . Even though the upper half-plane is simply connected, the marked points on the boundary can serve as moduli and then $\text{SLE}\left(\kappa ,\rho \right)$  is given by a particular moduli diffusion.
For multiply connected domains it is natural to look for an interaction $A$  which is expressed in terms of domain functionals such as the Green function. Appropriate combinations of derivatives of the Green function are homogeneous of degree minus one in $\xi$  and the moduli. The `harmonic random Loewner chains' studied in [26are a particular example of this.
In our opinion the only further reduction in possible diffusions $\left(\xi ,\mathbf{M}\right)$  are regularity requirements on the homogenous function $A$  . In particular we may wish to allow only functions which are analytic. We propose to call the growing family of random compacts obtained by solving the chordal Loewner equation  18 for a diffusion $\left(\xi ,\mathbf{M}\right)$  associated to an analytic function $A$  homogenous of degree minus one by  43 , chordal stochastic Loewner evolution.
Finally, we would like to point out very briefly several physical aspects, which shall be discussed in more detail in a separate publication.
The classical Green function is basically the two-point function of the bosonic free field with central charge $1$  . This elucidates the appearance of the number $-1$  , the degree of homogeneity. Further, scale invariance is related to the operator ${L}_{0}$  , the infinitesimal generator of scale transformations. It is an element, as is well known [5, of the Virasoro algebra, and in particular an element of the sub-algebra corresponding to the group of global conformal transformations.
Interactions with boundary components such as those mentioned above, are conveniently modelled in the Coulomb gas formalism, and correspond to insertions of different operators (i.e. currents or vertex operators) into the correlator.
As we are interested in describing the measure on random paths in a dynamical way, the moduli have to change, due to the deformation of the domain and/or the presence of marked points. This is in agreement with the global framework introduced in [9, 12, where correlators are modelled as sections of a twisted version of a determinant line bundle over the appropriate moduli space.
In the case of the upper half plane, the CFT content of $\text{SLE}\left(\kappa ,\rho \right)$  has been studied and explained in [6.

2 Bilateral and Chordal standard domains

2.1 Harmonic Measures

Denote $D$  a region of connectivity $n>1$  in the complex plane.
The components of the complement in the extended complex plane are denoted by ${E}_{1},{E}_{2},...,{E}_{n}$  . We assume that no ${E}_{k}$  reduces to a point and that there is a unique unbounded component ${E}_{n}$  . By applying preliminary conformal maps, we may assume that $D$  is bounded by an outer contour ${C}_{n}$  and $n-1$  inner contours ${C}_{1},...,{C}_{n-1}$  , where the contours are oriented such that $D$  lies to the left in the direction of the contour. Denote ${\omega }_{k}\left(z\right)$  the solution to the Dirichlet problem in $D$  with the boundary values $1$  on ${C}_{k}$  and $0$  on the other contours. We have $0<{\omega }_{k}\left(z\right)<1$  in $D$  and
 $\begin{array}{c}{\omega }_{1}\left(z\right)+{\omega }_{2}\left(z\right)+\cdots +{\omega }_{n}\left(z\right)=1.\end{array}$ (1)
${\omega }_{k}\left(z\right)$  is called the harmonic measure of ${C}_{k}$  in $z$  . The conjugate harmonic differential of ${\omega }_{k}$  has periods
 $\begin{array}{c}{\alpha }_{kj}={\int }_{{C}_{j}}*d{\omega }_{k}={\int }_{{C}_{j}}\frac{\partial {\omega }_{k}}{\partial n}ds\end{array}$ (2)
along ${C}_{j}$  . Here, $\partial /\partial n$  denotes the normal derivative to the right of the direction of the contour, and $ds$  stands for arc-length measure. It is well known, [18, that the $\left(n-1\right)×\left(n-1\right)$  matrix $\mathbit{\alpha }$  with entries ${\alpha }_{kj}$  , $1\le k,j\le n-1$  , is positive definite and symmetric. In particular, the linear system
 $\begin{array}{cc}{\lambda }_{1}{\alpha }_{11}+{\lambda }_{2}{\alpha }_{21}+\cdots +{\lambda }_{n-1}{\alpha }_{n-1,1}& =2\pi \end{array}$
 $\begin{array}{cc}{\lambda }_{1}{\alpha }_{12}+{\lambda }_{2}{\alpha }_{22}+\cdots +{\lambda }_{n-1}{\alpha }_{n-1,2}& =0\end{array}$
 $\begin{array}{cc}& ⋮\end{array}$ (3)
 $\begin{array}{cc}{\lambda }_{1}{\alpha }_{1,n-1}+{\lambda }_{2}{\alpha }_{2,n-1}+\cdots +{\lambda }_{n-1}{\alpha }_{n-1,n-1}& =0\end{array}$
 $\begin{array}{}\end{array}$
has a unique solution. It follows from  1 that any solution of  3 also solves ${\lambda }_{1}{\alpha }_{1n}+{\lambda }_{2}{\alpha }_{2n}+\cdots +{\lambda }_{n-1}{\alpha }_{n-1,n}=-2\pi .$  Thus there is a multiple-valued integral $F\left(z\right)$  with periods $±2\pi i$  along ${C}_{1}$  and ${C}_{n}$  and all other periods equal to zero, the real part being constant equal to ${\lambda }_{k}$  on ${C}_{k}$  (we set ${\lambda }_{n}=0$  ). The function $f\left(z\right)={e}^{-F\left(z\right)}$  is then single-valued and one can show, [1, that $f$  maps $D$  conformally onto the annulus ${e}^{-{\lambda }_{1}}<|w|<1$  minus $n-2$  concentric arcs situated on the circles $|w|={e}^{-{\lambda }_{k}}$  , $k=2,...,n-1$  . We call such a circularly slit annulus a bilateral standard domain. By adding an imaginary constant to $F\left(z\right)$  we obtain another map onto a bilateral standard domain and we may normalize the map $f$  by requiring $f\left({z}_{0}\right)={e}^{-{\lambda }_{1}}$  for some ${z}_{0}\in {C}_{1}$  .
With this normalization we call $f$  the canonical map for $\left(D,{z}_{0},{C}_{n}\right)$  .

2.2 Green Function

Denote $D$  again a region of finite connectivity which is bounded by contours ${C}_{1},...,{C}_{n}$  ; this time the case $n=1$  is included.
We consider a point ${z}_{0}\in D$  and solve the Dirichlet problem in $D$  with the boundary values $ln|\zeta -{z}_{0}|$  . The solution is denoted by $h\left(z\right)$  .
The function $G\left(z\right)={G}_{D}\left(z,{z}_{0}\right)=h\left(z\right)-ln|z-{z}_{0}|$  is the Green function in $D$  with pole at ${z}_{0}$  . It is the unique function which is harmonic in $D$  except at ${z}_{0}$  , where it differs from $ln|z-{z}_{0}|$  by a harmonic function, and which vanishes on the boundary of $D$  . The Green function is conformally invariant in the sense that if $f:D\to {D}^{\prime }$  is conformal, then
 $\begin{array}{c}{G}_{D}\left(z,{z}_{0}\right)={G}_{{D}^{\prime }}\left(f\left(z\right),f\left({z}_{0}\right)\right).\end{array}$ (4)
The conjugate harmonic function of $G\left(z,{z}_{0}\right)$  is multiple-valued. It has the period $2\pi$  along a small circle about ${z}_{0}$  , and the periods ${p}_{k}\left({z}_{0}\right)={\int }_{{C}_{k}}*dG\left(z,{z}_{0}\right),k=1,...,n.$  It can be shown that ${p}_{k}\left({z}_{0}\right)=-2\pi {\omega }_{k}\left({z}_{0}\right)$  , [1. Let now ${z}_{0}\in {C}_{n}$  .
By linearity, $u\left(z\right)=\partial G\left(z,{z}_{0}\right)/\partial {n}_{{z}_{0}}$  is a harmonic function in $z$  . Its conjugate differential has periods
 $\begin{array}{cc}{A}_{k}\left({z}_{0}\right)& ={\int }_{{C}_{k}}*du=\frac{\partial }{\partial {n}_{{z}_{0}}}{\int }_{{C}_{k}}*{d}_{z}G\left(z,{z}_{0}\right)\end{array}$
 $\begin{array}{cc}& =-2\pi \frac{\partial }{\partial {n}_{{z}_{0}}}{\omega }_{k}\left({z}_{0}\right).\end{array}$ (5)
 $\begin{array}{}\end{array}$
Thus the linear combination $u+{\lambda }_{1}{\omega }_{1}+\cdots +{\lambda }_{n-1}{\omega }_{n-1}$  is free from periods provided that
 $\begin{array}{c}{\lambda }_{1}{\alpha }_{1k}+{\lambda }_{2}{\alpha }_{2k}+\cdots +{\lambda }_{n-1}{\alpha }_{n-1,k}=-{A}_{k},k=1,...,n-1.\end{array}$ (6)
If we write $\mathbf{P}$  for the matrix $\mathbit{\alpha }/2\pi$  , ${\mathbit{\lambda }}^{T}=\left({\lambda }_{1},...,{\lambda }_{n-1}\right)$  , and $\partial \mathbit{\omega }\left({z}_{0}{\right)}^{T}/\partial n=\left(\partial {\omega }_{1}\left({z}_{0}\right)/\partial n,...,\partial {\omega }_{n-1}\left({z}_{0}\right)/\partial n\right),$  then the solution to  6 is given by $\mathbit{\lambda }={\mathbf{P}}^{-1}\frac{\partial \mathbit{\omega }\left({z}_{0}\right)}{\partial n}.$  Hence $-i\left(\frac{\partial G\left(z,{z}_{0}\right)}{\partial {n}_{{z}_{0}}}+\mathbit{\omega }\left(z{\right)}^{T}{\mathbf{P}}^{-1}\frac{\partial \mathbit{\omega }\left({z}_{0}\right)}{\partial n}\right)$  is the imaginary part of a single-valued analytic function $\Psi \left(z\right)$  . It can be shown that $\Psi$  maps $D$  conformally onto the upper half-plane $\Im \left(w\right)>0$  minus $n-1$  horizontal slits with imaginary parts $\Im \left(w\right)=-\left[{\mathbf{P}}^{-1}\partial \mathbit{\omega }\left({z}_{0}\right)/\partial n{\right]}_{j},j=1,...,n-1.$  Under this map, $\Psi \left({C}_{n}\right)=\mathbb{R}$  , and $\Psi \left({z}_{0}\right)=\infty$  . We call the upper half-plane minus a finite number of horizontal slits a chordal standard domain. If $D$  is contained in the upper half-plane and for some $x>0$  we have $\left(\mathbb{R}\\left[-x,x\right]\right)\subset {C}_{n}$  and $\zeta =\infty$  , then, by adding an appropriate real constant, we may assume that $g\equiv \Psi /2$  satisfies the hydrodynamic normalization at infinity,
 $\begin{array}{c}{lim}_{z\to \infty }\left(g\left(z\right)-z\right)=0.\end{array}$ (7)
With this normalization, we call $g$  the canonical mapping for $D$  .

3 Evolution of slit mappings

3.1 Chordal Loewner equation

Consider a chordal standard domain $D$  . Let $\gamma :\left[0,{t}_{\gamma }\right]\to \overline{D}$  be a Jordan arc such that $\gamma \left(0\right)\in \mathbb{R}$  , and $\gamma \left(0,{t}_{\gamma }\right]\subset D$  . Let ${g}_{t}$  be the canonical mapping from $D\\gamma \left[0,t\right]$  with the normalization  7 , and denote ${D}_{t}$  the chordal standard domain ${g}_{t}\left(D\\gamma \left[0,t\right]\right)$  . It is well known, see [18, that ${g}_{t}$  solves the extremal problem ${a}_{1}=max$  among all univalent functions on $D\\gamma \left[0,t\right]$  with expansion $z+\frac{{a}_{1}}{z}+\frac{{a}_{2}}{{z}^{2}}+\cdots ,{a}_{k}\in \mathbb{R},$  near infinity. In particular, if ${g}_{t}\left(z\right)=z+{a}_{t}/z+o\left(1/|z|\right)$  , then ${a}_{{t}^{*}}\le {a}_{t}$  whenever $0<{t}^{*}  . In fact, a simple argument shows that
 $\begin{array}{c}{a}_{{t}^{*}}<{a}_{t}\text{if}{t}^{*} (8)
Thus we may assume that ${a}_{t}=2t$  . We wish to find a differential equation for the family $\left\{{g}_{t}:t\in \left[0,{t}_{\gamma }\right]\right\}$  .
Denote ${C}_{j}\left(t\right),j=1,...,n$  , the boundary components of ${D}_{t}$  . We always have ${C}_{n}\left(t\right)=\mathbb{R}$  . For $j=1,...,n-1$  , let ${y}_{j}\left(t\right)$  be the imaginary part of (points on) the slit ${C}_{j}\left(t\right)$  . Denote $\xi \left(t\right)$  the starting point on $\mathbb{R}$  of the Jordan arc ${g}_{t}\left(\gamma \left[t,{t}_{\gamma }\right]\right)$  in ${D}_{t}$  , i.e. ${g}_{t}\left({\gamma }_{t}\right)$  . For $0<{t}^{*}  , set ${g}_{t,{t}^{*}}={g}_{{t}^{*}}\circ {g}_{t}^{-1}.$  Then ${g}_{t,{t}^{*}}$  is a conformal map from ${D}_{t}$  onto ${D}_{{t}^{*}}\{g}_{{t}^{*}}\left(\gamma \left[{t}^{*},t\right]\right)$  . The point $\xi \left({t}^{*}\right)={g}_{{t}^{*}}\left({\gamma }_{{t}^{*}}\right)$  corresponds to two prime ends in ${D}_{{t}^{*}}\{g}_{{t}^{*}}\left(\gamma \left[{t}^{*},t\right]\right)$  .
Denote ${\beta }_{0}\left(t,{t}^{*}\right)$  and ${\beta }_{1}\left(t,{t}^{*}\right)$  , with ${\beta }_{0}\left(t,{t}^{*}\right)<{\beta }_{1}\left(t,{t}^{*}\right)$  , the pre-images of these prime ends under ${g}_{t,{t}^{*}}$  , i.e. ${g}_{t,{t}^{*}}\left({\beta }_{0}\left(t,{t}^{*}\right)\right)={g}_{t,{t}^{*}}\left({\beta }_{1}\left(t,{t}^{*}\right)\right)={g}_{{t}^{*}}\left({\gamma }_{{t}^{*}}\right).$  Then, if $x\in \mathbb{R}\\left[{\beta }_{0}\left(t,{t}^{*}\right),{\beta }_{1}\left(t,{t}^{*}\right)\right]$  , ${g}_{t,{t}^{*}}\left(x\right)\in \mathbb{R}.$  Consider the analytic function $z↦{g}_{t,{t}^{*}}\left(z\right)-z,$  which satisfies
 $\begin{array}{c}{g}_{t,{t}^{*}}\left(z\right)-z=\frac{2\left({t}^{*}-t\right)}{z}+o\left(1/|z|\right),\end{array}$ (9)
and note that $z↦\Im \left({g}_{t,{t}^{*}}\left(z\right)-z\right)$  is harmonic and constant on each boundary component. By Poisson's formula
 $\begin{array}{c}\Im \left({g}_{t,{t}^{*}}\left(z\right)-z\right)=-\frac{1}{2\pi }{\int }_{\partial {D}_{t}}\Im \left({g}_{t,{t}^{*}}\left(\zeta \right)-\zeta \right)\frac{\partial {G}_{t}\left(\zeta ,z\right)}{\partial {n}_{1}}ds,\end{array}$ (10)
where ${G}_{t}\left(\zeta ,z\right)$  is the Green function for ${D}_{t}$  with pole at $z$  . Note that there is no problem with integrability in  10 because $\Im \left({g}_{t,{t}^{*}}\left(\zeta \right)-\zeta \right)=\frac{-y}{{x}^{2}+{y}^{2}}+O\left(1/|\zeta {|}^{2}\right),\zeta =x+iy,$  and
 $\begin{array}{c}sup\left\{\Im \left(\zeta \right):\zeta \in \partial {D}_{t}\right\}<\infty .\end{array}$ (11)
Since $\Im \left({g}_{t,{t}^{*}}\left(z\right)-z\right)$  has a single-valued harmonic conjugate, it is orthogonal to the real part of any Abelian differential of the first kind, see [3, and we have
 $\begin{array}{cc}\Im & \left({g}_{t,{t}^{*}}\left(z\right)-z\right)\end{array}$
 $\begin{array}{cc}& =-\frac{1}{2\pi }{\int }_{\partial {D}_{t}}\Im \left({g}_{t,{t}^{*}}\left(\zeta \right)-\zeta \right)\left(\frac{\partial {G}_{t}\left(\zeta ,z\right)}{\partial {n}_{1}}+{\mathbit{\omega }}_{t}\left(z{\right)}^{T}{\mathbit{P}}_{t}^{-1}\frac{\partial {\mathbit{\omega }}_{t}\left(\zeta \right)}{\partial n}\right)ds.\end{array}$ (12)
 $\begin{array}{}\end{array}$
It follows from Section  2 that ${\int }_{{C}_{k}\left(t\right)}\left(\frac{\partial {G}_{t}\left(\zeta ,z\right)}{\partial {n}_{1}}+{\mathbit{\omega }}_{t}\left(z{\right)}^{T}{\mathbit{P}}_{t}^{-1}\frac{\partial {\mathbit{\omega }}_{t}\left(\zeta \right)}{\partial n}\right)ds=0,k=1,...,n-1,$  and also that $z↦-i\left(\frac{\partial {G}_{t}\left(\zeta ,z\right)}{\partial {n}_{1}}+{\mathbit{\omega }}_{t}\left(z{\right)}^{T}{\mathbit{P}}_{t}^{-1}\frac{\partial {\mathbit{\omega }}_{t}\left(\zeta \right)}{\partial n}\right)$  is the imaginary part of a single-valued analytic function ${\Psi }_{t}\left(z\right)={\Psi }_{t}\left(z,\zeta \right)$  . Thus, since $\Im \left({g}_{t,{t}^{*}}\left(\zeta \right)-\zeta \right)$  is constant on each ${C}_{k}\left(t\right)$  , $k=1,...,n-1$  , and identically zero on $\mathbb{R}\\left[{\beta }_{0}\left(t,{t}^{*}\right),{\beta }_{1}\left(t,{t}^{*}\right)\right]$  ,
 $\begin{array}{c}{g}_{t,{t}^{*}}\left(z\right)-z=\frac{1}{2\pi }{\int }_{{\beta }_{0}\left(t,{t}^{*}\right)}^{{\beta }_{1}\left(t,{t}^{*}\right)}\Im \left({g}_{t,{t}^{*}}\left(\zeta \right)-\zeta \right){\Psi }_{t}\left(z,\zeta \right)d\zeta +ic,\end{array}$ (13)
where $c$  is a real constant. Note that if $z↦{\stackrel{~}{\Psi }}_{t}\left(z,\zeta \right)$  is another analytic function with the same imaginary part as ${\Psi }_{t}$  , then ${\Psi }_{t}\left(z,\zeta \right)-{\stackrel{~}{\Psi }}_{t}\left(z,\zeta \right)=a\left(\zeta \right),$  where $a$  is real and depends only on $\zeta$  . We fix a normalization by requiring that
 $\begin{array}{c}{lim}_{z\to \infty }{\Psi }_{t}\left(z,\zeta \right)=0.\end{array}$ (14)
If we let $z\to \infty$  , then ${g}_{t,{t}^{*}}\left(z\right)-z\to 0$  . By bounded convergence, the integral in  13 converges to zero as well and it follows that $c=0$  . Next, $2\left({t}^{*}-t\right)={lim}_{z\to \infty }z\left({g}_{t,{t}^{*}}\left(z\right)-z\right)=f\left(0\right),$  where $w↦f\left(w\right)\equiv \frac{1}{w}\left[{g}_{t,{t}^{*}}\left(1/w\right)-\left(1/w\right)\right]$  is regular near zero. By the Schwarz reflection principle ${g}_{t,{t}^{*}}$  extends to the entire complex plane minus the slits ${C}_{1}\left(t\right),...,{C}_{n-1}\left(t\right)$  , their conjugates, and the real interval $\left[{\beta }_{0}\left(t,{t}^{*}\right),{\beta }_{1}\left(t,{t}^{*}\right)\right]$  . Denote $C$  the collection of these $2n-1$  finite slits. Then $f$  also extends to a corresponding domain with boundary $\stackrel{~}{C}$  . From Cauchy's integral formula we have
 $\begin{array}{cc}2\left({t}^{*}-t\right)& =\frac{1}{2\pi i}{\int }_{\stackrel{~}{C}}\frac{f\left(\zeta \right)}{\zeta }d\zeta =\frac{1}{2\pi i}{\int }_{\stackrel{~}{C}}\frac{{g}_{t,{t}^{*}}\left(1/\zeta \right)-1/\zeta }{{\zeta }^{2}}d\zeta \end{array}$
 $\begin{array}{cc}& =-\frac{1}{2\pi i}{\int }_{C}\left({g}_{t,{t}^{*}}\left(\eta \right)-\eta \right)d\eta =-\frac{1}{2\pi }{\int }_{C}\Im \left({g}_{t,{t}^{*}}\left(\eta \right)-\eta \right)d\eta ,\end{array}$ (15)
 $\begin{array}{}\end{array}$
where the final equality uses the fact that $d\eta$  is real for horizontal slits.
The slits ${C}_{1}\left(t\right),...,{C}_{n-1}\left(t\right)$  and their conjugates do not contribute to the last integral since $\Im \left({g}_{t,{t}^{*}}\left(\eta \right)-\eta \right)$  takes the same value on both “sides” of a given slit. For the slit $\left[{\beta }_{0}\left(t,{t}^{*}\right),{\beta }_{1}\left(t,{t}^{*}\right)\right]$  , $\Im \left({g}_{t,{t}^{*}}\left(\eta \right)-\eta \right)$  takes opposite values on the upper and lower “side” of the slit and, since the direction of integration is reversed, we finally get
 $\begin{array}{c}{t}^{*}-t=-\frac{1}{2\pi }{\int }_{{\beta }_{0}\left(t,{t}^{*}\right)}^{{\beta }_{1}\left(t,{t}^{*}\right)}\Im \left({g}_{t,{t}^{*}}\left(\eta \right)\right)d\eta .\end{array}$ (16)
Setting $z={g}_{t}\left(w\right)$  in  13 we have ${g}_{{t}^{*}}\left(w\right)-{g}_{t}\left(w\right)=\frac{1}{2\pi }{\int }_{{\beta }_{0}\left(t,{t}^{*}\right)}^{{\beta }_{1}\left(t,{t}^{*}\right)}\Im \left({g}_{t,{t}^{*}}\left(\eta \right)\right){\Psi }_{t}\left(z,\eta \right)d\eta .$  We are now ready to let ${t}^{*}↗t$  in  17 . Note first that, for $\eta \in \left[{\beta }_{0}\left(t,{t}^{*}\right),{\beta }_{1}\left(t,{t}^{*}\right)\right]$  , $\eta ↦\Im \left({g}_{t,{t}^{*}}\left(\eta \right)\right)$  is continuous and non-negative and that also $\eta ↦A\left(\eta \right):={\Psi }_{t}\left(z,\eta \right)$  is continuous. Thus it follows from the mean-value theorem of integration and  16 that
 $\begin{array}{cc}\frac{1}{2\pi }& {\int }_{{\beta }_{0}\left(t,{t}^{*}\right)}^{{\beta }_{1}\left(t,{t}^{*}\right)}\Im \left({g}_{t,{t}^{*}}\left(\eta \right)\right)A\left(\eta \right)d\eta \end{array}$
 $\begin{array}{cc}& =\frac{\Re \left(A\left({\eta }^{\prime }\right)\right)+i\Im \left(A\left({\eta }^{\prime \prime }\right)\right)}{2\pi }{\int }_{{\beta }_{0}\left(t,{t}^{*}\right)}^{{\beta }_{1}\left(t,{t}^{*}\right)}\Im \left({g}_{t,{t}^{*}}\left(\eta \right)d\eta \end{array}$
 $\begin{array}{cc}& =-\left[\Re \left(A\left({\eta }^{\prime }\right)\right)+i\Im \left(A\left({\eta }^{\prime \prime }\right)\right)\right]\left({t}^{*}-t\right),\end{array}$ (17)
 $\begin{array}{}\end{array}$
for some ${\eta }^{\prime },{\eta }^{\prime \prime }\in \left[{\beta }_{0}\left(t,{t}^{*}\right),{\beta }_{1}\left(t,{t}^{*}\right)\right]$  . Hence ${lim}_{{t}^{*}↗t}\frac{{g}_{{t}^{*}}\left(w\right)-{g}_{t}\left(w\right)}{{t}^{*}-t}=-{\Psi }_{t}\left(z,\xi \left(t\right)\right).$  By the same argument we may let $t↘{t}^{*}$  . On the right-hand side above we then only need to change $t$  to ${t}^{*}$  and introduce an overall minus sign. Thus we have established the following
Theorem 3.1 (Chordal Loewner equation). If $\gamma$  is a Jordan arc in a standard domain $D$  starting on $\mathbb{R}$  with the parametrization from above, and if ${g}_{t}$  is the canonical map for $D\\gamma \left[0,t\right]$  , then, using the notation from above, the family $\left\{{g}_{t}:t\in \left[0,{t}_{\gamma }\right]\right\}$  satisfies the equation
 $\begin{array}{c}{\partial }_{t}{g}_{t}\left(z\right)=-{\Psi }_{t}\left({g}_{t}\left(z\right),\xi \left(t\right)\right),\end{array}$ (18)
with initial condition ${g}_{0}\left(z\right)=z$  , and where ${\Psi }_{t}\left(z,\zeta \right)$  is the analytic function in $z$  with imaginary part $-\frac{\partial {G}_{t}\left(z,\zeta \right)}{\partial {n}_{\zeta }}-{\mathbit{\omega }}_{t}\left(z{\right)}^{T}{\mathbit{P}}_{t}^{-1}\frac{\partial {\mathbit{\omega }}_{t}\left(\zeta \right)}{\partial n},$  normalized by ${lim}_{z\to \infty }{\Psi }_{t}\left(z,\zeta \right)=0$  .
Remark 3.1. In the simply connected case, when $D=\mathbb{H}$  is the upper half-plane, the Green function is given by $G\left(z,w\right)=\Re \left(ln\frac{z-\overline{w}}{z-w}\right).$  Thus, if $w=x+iy$  , $-\Psi \left(z,w\right)=-i\frac{\partial }{\partial y}{|}_{y=0}ln\frac{z-x+iy}{z-x-iy}=\frac{2}{z-w},$  and  18 reduces to the well known chordal Loewner equation.

3.2 Bilateral Komatu-Loewner equation

The evolution of slit mappings in multiply connected domains was first studied by Komatu in [10for the doubly connected case, and in [11for general finite connectivity. Komatu treated this case by considering circular slit annuli.
Consider a bilateral standard domain $D$  with inner radius $Q$  . Let $\gamma :\left[0,{t}_{\gamma }\right]\to \overline{D}$  be a Jordan arc such that $\gamma \left(0\right)\in {S}^{1}$  , and $\gamma \left(0,{t}_{\gamma }\right]\subset D$  .
Let ${f}_{t}$  be the canonical mapping from $D\\gamma \left[0,t\right]$  with the normalization ${f}_{t}\left(Q\right)>0$  , and denote ${D}_{t}$  the chordal standard domain ${f}_{t}\left(D\\gamma \left[0,t\right]\right)$  .
If ${Q}_{t}={f}_{t}\left(Q\right)$  , then it can be shown that $t\in \left[0,{t}_{\gamma }\right]↦{Q}_{t}\in \left[Q,1\right]$  is continuous and strictly increasing, [11. Thus we may assume that $\gamma$  is parametrized such that $t=ln{Q}_{t}$  . For this parameter it is shown in [11that $t↦{f}_{t}\left(z\right)$  is differentiable. An expression for the derivative is also given. However, the expression given there is not explicit enough for the purposes we have in mind. In particular, we will need to know that the vector field is itself a Lipschitz function in the moduli of the domain.
We sketch a proof of what we call the bilateral Komatu-Loewner equation, leading to an expression of the derivative ${\partial }_{t}{f}_{t}$  in terms of the Green function, harmonic measures, their derivatives and harmonic conjugates.
The argument is similar to the radial case, [4. In fact, the radial case can be obtained as a limiting case from the bilateral case when $Q\to 0$  , [10.
Denote ${C}_{j}\left(t\right),j=1,...,n$  , the boundary components of ${D}_{t}$  . We always have ${C}_{n}\left(t\right)={S}^{1}$  , and ${C}_{1}\left(t\right)=\left\{|z|={e}^{t}\right\}$  . For $j=2,...,n-1$  , let ${m}_{j}\left(t\right)$  be the radial distance of the circular slit ${C}_{j}\left(t\right)$  from the origin.
Denote $\xi \left(t\right)$  the starting point on ${S}^{1}$  of the Jordan arc ${g}_{t}\left(\gamma \left[t,{t}_{\gamma }\right]\right)$  in ${D}_{t}$  , i.e. ${g}_{t}\left({\gamma }_{t}\right)$  . For $lnQ<{t}^{*}  , set ${g}_{t,{t}^{*}}={g}_{{t}^{*}}\circ {g}_{t}^{-1}.$  Then ${g}_{t,{t}^{*}}$  is a conformal map from ${D}_{t}$  onto ${D}_{{t}^{*}}\{g}_{{t}^{*}}\left(\gamma \left[{t}^{*},t\right]\right)$  . The point $\xi \left({t}^{*}\right)={g}_{{t}^{*}}\left({\gamma }_{{t}^{*}}\right)$  corresponds to two prime ends in ${D}_{{t}^{*}}\{g}_{{t}^{*}}\left(\gamma \left[{t}^{*},t\right]\right)$  .
Denote $exp\left(i{\beta }_{0}\left(t,{t}^{*}\right)\right)$  and $exp\left(i{\beta }_{1}\left(t,{t}^{*}\right)\right)$  , with ${\beta }_{0}\left(t,{t}^{*}\right)<{\beta }_{1}\left(t,{t}^{*}\right)$  , the pre-images of these prime ends under ${g}_{t,{t}^{*}}$  , i.e. ${g}_{t,{t}^{*}}\left(exp\left(i{\beta }_{0}\left(t,{t}^{*}\right)\right)\right)={g}_{t,{t}^{*}}\left(exp\left(i{\beta }_{1}\left(t,{t}^{*}\right)\right)\right)={g}_{{t}^{*}}\left({\gamma }_{{t}^{*}}\right).$  Then, if $|z|=1$  and ${\beta }_{1}\left(t,{t}^{*}\right)\le argz\le {\beta }_{0}\left(t,{t}^{*}\right)+2\pi$  , $|{g}_{t,{t}^{*}}\left(z\right)|=1.$  The function $z↦ln\frac{{g}_{t,{t}^{*}}\left(z\right)}{z}$  is analytic and single-valued throughout ${D}_{t}$  . By Poisson's formula
 $\begin{array}{c}ln|\frac{{g}_{t,{t}^{*}}\left(z\right)}{z}|=-\frac{1}{2\pi }{\int }_{\partial {D}_{t}}ln|\frac{{g}_{t,{t}^{*}}\left(\zeta \right)}{\zeta }|\frac{\partial {G}_{t}\left(\zeta ,z\right)}{\partial {n}_{1}}ds,\end{array}$ (19)
where ${G}_{t}\left(\zeta ,z\right)$  is the Green function for ${D}_{t}$  with pole at $z$  . Using orthogonality and the period relations as we did in the chordal case, it follows that
 $\begin{array}{c}ln\frac{{g}_{t,{t}^{*}}\left(z\right)}{z}=-\frac{i}{2\pi }{\int }_{{\beta }_{0}\left(t,{t}^{*}\right)}^{{\beta }_{1}\left(t,{t}^{*}\right)}ln|\frac{{g}_{t,{t}^{*}}\left(\zeta \right)}{\zeta }|{\Psi }_{t}\left(z,\zeta \right)ds+ic,\end{array}$ (20)
for some real constant $c$  . To eliminate $c$  , we evaluate the identity  20 at $z=q={e}^{t}$  and then take the difference:
 $\begin{array}{c}ln\frac{{g}_{t,{t}^{*}}\left(z\right)}{z}-ln\frac{{q}^{*}}{q}=-\frac{i}{2\pi }{\int }_{{\beta }_{0}\left(t,{t}^{*}\right)}^{{\beta }_{1}\left(t,{t}^{*}\right)}ln|\frac{{g}_{t,{t}^{*}}\left(\zeta \right)}{\zeta }|\left[{\Psi }_{t}\left(z,\zeta \right)-{\Psi }_{t}\left(q,\zeta \right)\right]ds.\end{array}$ (21)
By Cauchy's integral formula,
 $\begin{array}{c}0=\frac{1}{2\pi i}{\int }_{\partial {D}_{t}}ln\left(\frac{{g}_{t,{t}^{*}}\left(\zeta \right)}{\zeta }\right)\frac{d\zeta }{\zeta }.\end{array}$ (22)
In particular, the right-hand side of  22 is real. Since all boundary components are concentric circular arcs, $d\zeta /\zeta$  is purely imaginary along $\partial {D}_{t}$  , i.e. $\frac{d\zeta }{\zeta }=idarg\zeta ,\zeta \in \partial {D}_{t}.$  Hence
 $\begin{array}{cc}0& =\frac{1}{2\pi }{\int }_{\partial {D}_{t}}ln|\frac{{g}_{t,{t}^{*}}\left(\zeta \right)}{\zeta }|darg\zeta \end{array}$
 $\begin{array}{cc}& =\frac{1}{2\pi }{\int }_{{\beta }_{0}\left(t,{t}^{*}\right)}^{{\beta }_{1}\left(t,{t}^{*}\right)}ln|{g}_{t,{t}^{*}}\left({e}^{i\phi }\right)|d\phi -\frac{1}{2\pi }{\int }_{0}^{2\pi }ln\frac{{q}^{*}}{q}d\phi \end{array}$
 $\begin{array}{cc}& +\frac{1}{2\pi }{\sum }_{j=2}^{n-1}{\int }_{{C}_{j}\left(t\right)}ln\frac{{m}_{j}\left({t}^{*}\right)}{{m}_{j}\left(t\right)}darg\zeta .\end{array}$ (23)
 $\begin{array}{}\end{array}$
Since the two “sides” of ${C}_{j}\left(t\right)$  make opposite contributions, ${\int }_{{C}_{j}\left(t\right)}darg\zeta =0,j=2,...,n-1,$  and we finally get
 $\begin{array}{c}{t}^{*}-t=\frac{1}{2\pi }{\int }_{{\beta }_{0}\left(t,{t}^{*}\right)}^{{\beta }_{1}\left(t,{t}^{*}\right)}ln|{g}_{t,{t}^{*}}\left({e}^{i\phi }\right)|d\phi .\end{array}$ (24)
Letting $z={g}_{t}\left(w\right)$  in  21 , we have
 $\begin{array}{cc}ln\frac{{g}_{{t}^{*}}\left(w\right)}{{g}_{t}\left(w\right)}& -\left({t}^{*}-t\right)\end{array}$
 $\begin{array}{cc}& =-\frac{i}{2\pi }{\int }_{{\beta }_{0}\left(t,{t}^{*}\right)}^{{\beta }_{1}\left(t,{t}^{*}\right)}ln|{g}_{t,{t}^{*}}\left({e}^{i\phi }\right)|\left[{\Psi }_{t}\left(z,{e}^{i\phi }\right)-{\Psi }_{t}\left(q,{e}^{i\phi }\right)\right]ds.\end{array}$ (25)
 $\begin{array}{}\end{array}$
We now wish to let ${t}^{*}↗t$  in  25 . Note first that, for $\phi \in \left[0,2\pi \right]$  , $\phi ↦ln|{g}_{t,{t}^{*}}\left({e}^{i\phi }\right)|$  is continuous and non-positive and that also $\phi ↦A\left(\phi \right):={\Psi }_{t}\left(z,{e}^{i\phi }\right)-{\Psi }_{t}\left(q,{e}^{i\phi }\right)$  is continuous. Thus it follows from the mean-value theorem of integration that
 $\begin{array}{cc}\frac{1}{2\pi \left({t}^{*}-t\right)}& {\int }_{{\beta }_{0}\left(t,{t}^{*}\right)}^{{\beta }_{1}\left(t,{t}^{*}\right)}ln|{g}_{t,{t}^{*}}\left({e}^{i\phi }\right)|A\left(\phi \right)d\phi \end{array}$
 $\begin{array}{cc}& =\frac{\Re \left(A\left({\phi }^{\prime }\right)\right)+i\Im \left(A\left({\phi }^{\prime \prime }\right)\right)}{2\pi \left({t}^{*}-t\right)}{\int }_{{\beta }_{0}\left(t,{t}^{*}\right)}^{{\beta }_{1}\left(t,{t}^{*}\right)}ln|{g}_{t,{t}^{*}}\left({e}^{i\phi }\right)|d\phi \end{array}$
 $\begin{array}{cc}& =\Re \left(A\left({\phi }^{\prime }\right)\right)+i\Im \left(A\left({\phi }^{\prime \prime }\right)\right),\end{array}$ (26)
 $\begin{array}{}\end{array}$
for some ${\phi }^{\prime },{\phi }^{\prime \prime }\in \left[{\beta }_{0}\left(t,{t}^{*}\right),{\beta }_{1}\left(t,{t}^{*}\right)\right]$  . Hence
 $\begin{array}{c}{lim}_{{t}^{*}↗t}\frac{ln{g}_{{t}^{*}}\left(w\right)-ln{g}_{t}\left(w\right)}{{t}^{*}-t}=1+i\left[{\Psi }_{t}\left(z,{\xi }_{t}\right)-{\Psi }_{t}\left({e}^{t},\xi \left(t\right)\right)\right].\end{array}$ (27)
By the same argument we may let $t↘{t}^{*}$  . On the right-hand side above we then only need to change $t$  to ${t}^{*}$  and introduce an overall minus sign. Thus we have established the following
Theorem 3.2 (Bilateral Komatu-Loewner equation). If $\gamma$  is a Jordan arc in a standard domain $D$  starting on ${S}^{1}$  with the parametrization from above, and if ${g}_{t}$  is the canonical map for $D\\gamma \left[0,t\right]$  , then, using the notation from above, the family $\left\{{g}_{t}:t\in \left[lnQ,{t}_{\gamma }\right]\right\}$  satisfies the equation
 $\begin{array}{c}{\partial }_{t}ln{g}_{t}\left(z\right)=1+i\left[{\Psi }_{t}\left({g}_{t}\left(z\right),{\xi }_{t}\right)-{\Psi }_{t}\left({e}^{t},\xi \left(t\right)\right)\right],\end{array}$ (28)
with initial condition ${g}_{lnQ}\left(z\right)=z$  .

4 Motion of moduli

4.1 Chordal case

The right-hand side of the chordal Loewner equation, at time $t$  , involves the Green function of the domain ${D}_{t}$  , and also various functions derived from the Green function. Consequently, it does not make sense to ask for the solution of  18 for a given continuous curve $t↦\xi \left(t\right)$  , since the vector-field on the right-hand side of  18 is not specified by giving that information alone. To specify the Green function of ${D}_{t}$  we also need the moduli of the domain ${D}_{t}$  . We will now consider what the appropriate moduli space is for our purposes and find a system of equations these moduli satisfy. Once this system is found, we can solve it for a given input $t↦\xi \left(t\right)$  , and then, in a second step, solve the radial Komatu-Loewner equation using $\xi$  and the moduli.
The geometric description of ${D}_{t}$  requires $3n-3$  real parameters, three for each (interior) slit, given, for example, by the imaginary components of the slits, i.e ${y}_{j}\left(t\right)$  , $j=1,...,n-1$  , and the real components ${x}_{j}\left(t\right)<{x}_{j}^{\prime }\left(t\right),j=1,...,n-1,$  determining the endpoints of the slit ${C}_{j}\left(t\right)$  , $j=1...,n-1$  . On the other hand, it is well known that two $n$  -connected domains with non-degenerate boundary continua are conformally equivalent if $3n-6$  real parameters agree for $n>2$  . If $n=2$  then there is only one real parameter describing the conformal class, and if $n=1$  , then all such domains are conformally equivalent.
The slits we wish to grow mark two points on one of the boundary continua, the beginning ( $t=0$  ) and end point $\left(t=\infty \right)$  of the slit. Any $n$  -connected planar domain with two marked boundary points on one boundary component is conformally equivalent to the upper half-plane with $n-1$  horizontal slits and such that the marked boundary points are mapped to $0$  and $\infty$  . However, there is a one-parameter group of automorphisms, namely multiplication by $a>0$  , which maps the slit upper half-plane onto a slit upper half-plane, while fixing 0 and $\infty$  . It is now easy to see that the moduli space of $n$  -connected planar domains with two marked boundary points on one of the boundary components is $3n-4$  dimensional for all $n\ge 2$  , and zero dimensional if $n=1$  . Nonetheless, we will take $\mathbf{y}\left(t\right)=\left({y}_{1}\left(t\right),...,{y}_{n-1}\left(t\right)\right)$  , $\mathbf{x}\left(t\right)=\left({x}_{1}\left(t\right),...,{x}_{n-1}\left(t\right)\right)$  , and ${\mathbf{x}}^{\prime }\left(t\right)=\left({x}_{1}^{\prime }\left(t\right),...,{x}_{n-1}^{\prime }\left(t\right)\right)$  as the moduli of the domain ${D}_{t}$  and write $\mathbf{M}\left(t\right):=\left(\mathbf{y}\left(t\right),\mathbf{x}\left(t\right),{\mathbf{x}}^{\prime }\left(t\right)\right)$  . To obtain the conformal equivalence classes from this $3n-3$  dimensional parameter space, we need to identify $\left(\mathbf{y}\left(t\right),\mathbf{x}\left(t\right),{\mathbf{x}}^{\prime }\left(t\right)\right)$  and $\left(\stackrel{\mathbf{~}}{\mathbf{y}}\left(t\right),\stackrel{\mathbf{~}}{\mathbf{x}}\left(t\right),{\stackrel{\mathbf{~}}{\mathbf{x}}}^{\prime }\left(t\right)\right)$  , whenever there exists an $a>0$  such that $\mathbf{y}=a\stackrel{\mathbf{~}}{\mathbf{y}}$  , $\mathbf{x}=\mathbf{a}\stackrel{\mathbf{~}}{\mathbf{x}}$  , and ${\mathbf{x}}^{\prime }=a{\stackrel{\mathbf{~}}{\mathbf{x}}}^{\prime }$  .
The extra parameter $\mathbf{M}$  keeps track of will be reflected in a symmetry (invariance) of the moduli diffusion. For a standard domain the marked points are 0 and $\infty$  . For a point $\mathbf{M}$  in the “moduli space” we denote by $D=D\left(\mathbit{M}\right)$  the corresponding standard domain.
By boundary correspondence, if $z\in {C}_{j}$  , then ${g}_{t}\left(z\right)\in {C}_{j}\left(t\right)$  and $\Im \left({g}_{t}\left(z\right)\right)={y}_{j}\left(t\right).$  Thus, by considering the imaginary part of the chordal Loewner equation,
 $\begin{array}{c}{\partial }_{t}{y}_{j}\left(t\right)=-\Im \left({\Psi }_{t}\left({g}_{t}\left(z\right),\xi \left(t\right)\right)\right).\end{array}$ (29)
Further, if ${z}_{j}\left(t\right)={x}_{j}\left(t\right)+i{y}_{j}\left(t\right),{z}_{j}^{\prime }\left(t\right)={x}_{j}^{\prime }\left(t\right)+i{y}_{j}\left(t\right)$  are the endpoints of the slit ${C}_{j}\left(t\right)$  , then ${z}_{j}\left(t\right)={g}_{t}\left({\eta }_{j}\left(t\right)+i{y}_{j}\left(0\right)\right),{z}_{j}^{\prime }\left(t\right)={g}_{t}\left({\eta }_{j}^{\prime }\left(t\right)+i{y}_{j}\left(0\right)\right),$  where ${x}_{j}\left(0\right)<{\eta }_{j}\left(t\right),{\eta }_{j}^{\prime }\left(t\right)<{x}_{j}^{\prime }\left(0\right)$  . Indeed, the pre-images of the tips of ${C}_{j}\left(t\right)$  , that is ${\eta }_{j}\left(t\right)+i{y}_{j}\left(0\right)$  and ${\eta }_{j}^{\prime }\left(t\right)+i{y}_{j}\left(0\right)$  , are the solutions to the equation $\frac{\partial }{\partial z}{g}_{t}\left(z\right)=0,$  on the set of prime-ends corresponding to ${C}_{j}\\left\{{z}_{j}\left(0\right),{z}_{j}^{\prime }\left(0\right)\right\}$  . A tip of ${C}_{j}\left(t\right)$  cannot be the image of a tip of ${C}_{j}$  because then the analytic function $\partial {g}_{t}/\partial z$  would not have the required number of zeroes, $2n-2$  .
Lemma 4.1 (Motion of moduli—chordal case). The moduli $\mathbf{M}\left(t\right)=\left(\mathbf{y}\left(t\right),\mathbf{x}\left(t\right),{\mathbf{x}}^{\prime }\left(t\right)\right)$  satisfy the system of equations
 $\begin{array}{cc}{\partial }_{t}{y}_{j}\left(t\right)& ={\left[{\mathbit{P}}_{t}^{-1}\frac{\partial {\mathbit{\omega }}_{t}\left(\xi \left(t\right)\right)}{\partial n}\right]}_{j},\end{array}$
 $\begin{array}{cc}{\partial }_{t}{x}_{j}\left(t\right)& =-\Re \left({\Psi }_{t}\left({x}_{j}\left(t\right)+i{y}_{j}\left(t\right),\xi \left(t\right)\right)\right),\end{array}$
 $\begin{array}{cc}{\partial }_{t}{x}_{j}^{\prime }\left(t\right)& =-\Re \left({\Psi }_{t}\left({x}_{j}^{\prime }\left(t\right)+i{y}_{j}\left(t\right),\xi \left(t\right)\right)\right),\end{array}$ (30)
 $\begin{array}{}\end{array}$
for $j=1,...,n-1$  .
• Proof. We note that $\partial {g}_{t}/\partial z$  and ${\partial }^{2}{g}_{t}/\left(\partial z{\right)}^{2}$  are analytic functions that extend analytically to the prime-ends corresponding to ${C}_{1},...,{C}_{n-1}$  with the endpoints of the slits removed. By the implicit function theorem, $t↦{\eta }_{j}\left(t\right)+i{y}_{j}\left(0\right)$  is differentiable with derivative $DE{R}_{t}:={\left[\frac{{\partial }^{2}{g}_{t}}{\left(\partial z{\right)}^{2}}\left({\eta }_{j}\left(t\right)+i{y}_{j}\left(0\right)\right)\right]}^{-1}\frac{{\partial }^{2}{g}_{t}}{\partial t\partial z}\left({\eta }_{j}\left(t\right)+i{y}_{j}\left(0\right)\right).$  By counting zeroes we find that $\frac{{\partial }^{2}{g}_{t}}{\left(\partial z{\right)}^{2}}\left({\eta }_{j}\left(t\right)+i{y}_{j}\left(t\right)\right)\ne 0$  and so $DE{R}_{t}$  is finite. Hence
 $\begin{array}{cc}{\partial }_{t}{x}_{j}\left(t\right)& ={\partial }_{t}\Re \left({g}_{t}\left({\eta }_{j}\left(t\right)+i{y}_{j}\left(0\right)\right)\right)\end{array}$
 $\begin{array}{cc}& =-\Re \left({\Psi }_{t}\left({x}_{j}\left(t\right)+i{y}_{j}\left(t\right),\xi \left(t\right)\right)\right)+\Re \left(\left({\partial }_{z}{g}_{t}\right)\left({z}_{j}\left(t\right)\right)×DE{R}_{t}\right)\end{array}$
 $\begin{array}{cc}& =-\Re \left({\Psi }_{t}\left({x}_{j}\left(t\right)+i{y}_{j}\left(t\right),\xi \left(t\right)\right)\right).\end{array}$ (31)
 $\begin{array}{}\end{array}$
In a similar way we obtain the derivative of ${x}_{j}^{\prime }\left(t\right)$  . It remains to check that  29 agrees with the first equation in  30 . To this end we note that $\Im \left({\Psi }_{t}\left(z,\zeta \right)\right)=\frac{\partial {G}_{t}\left(\zeta ,z\right)}{\partial {n}_{1}}+{\mathbit{\omega }}_{t}\left(z{\right)}^{T}{\mathbit{P}}_{t}^{-1}\frac{\partial {\mathbit{\omega }}_{t}\left(\zeta \right)}{\partial n}.$  From the boundary behavior of the Green function and the harmonic measures, it follows that for $z\in {C}_{j}\left(t\right)$  $\frac{\partial G\left(\zeta ,z;t\right)}{\partial {n}_{1}}=0,\text{and}{\omega }_{k}\left(z\right)={\delta }_{jk}.$  The lemma follows.
We now have our main existence statement.
Theorem 4.1. Given a continuous function $t\in \left[0,\infty \right)↦\xi \left(t\right)\in \mathbb{R}$  and the moduli $\mathbf{M}$  of a standard domain $D$  , there exists a unique solution $\mathbf{M}\left(t\right)$  to the system  30 on an interval $\left[0,{t}_{\xi }\right)$  with $\mathbf{M}\left(0\right)=\mathbf{M}$  , and where ${t}_{\xi }$  is characterized by ${t}_{\xi }=inf\left\{\tau :{lim}_{t↗\tau }{y}_{j}\left(t\right)=0\text{for some}j\in \left\{1,...,n-1\right\}\right\}.$  Further, if ${D}_{t}$  is the standard domain determined by $\mathbf{M}\left(t\right)$  , and if ${\Psi }_{t}\left(z,\zeta \right)$  is the holomorphic vector field associated to ${D}_{t}$  (cf. Section  2.2 ), then, for any $z\in D$  , the equation ${\partial }_{t}{g}_{t}^{D}\left(z\right)={\Psi }_{t}\left({g}_{t}^{D}\left(z\right),\xi \left(t\right)\right),{g}_{0}^{D}\left(z\right)=z,$  has a unique solution on $\left[0,{t}_{z}\right)$  , where ${t}_{z}=sup\left\{t\le {t}_{\xi }:{inf}_{s\in \left[0,t\right]}|{g}_{s}^{D}\left(z\right)-\xi \left(s\right)|>0\right\}.$  Finally, for $t<{t}_{\xi }$  set ${K}_{t}=\left\{z\in D:{t}_{z}\le t\right\}$  . Then ${g}_{t}^{D}$  is the canonical conformal map from $D\{K}_{t}$  onto ${D}_{t}$  with hydrodynamic normalization at infinity.
• Proof. For the existence of the solution to the moduli equations  30 on $\left[0,{t}_{\xi }\right)$  we need to know that the vector field in  30 is Lipschitz as a function of $\mathbf{M}$  , with a Lipschitz constant that only depends on distance to $\xi \left(t\right)$  of the slit (or slits) nearest to $\xi \left(t\right)$  . Let $\mathbf{M}$  and $\stackrel{\mathbf{~}}{\mathbf{M}}$  be two points in moduli space with corresponding standard domains $D$  and $\stackrel{~}{D}$  , such that $|{y}_{j}-{\stackrel{~}{y}}_{j}|,|{x}_{j}-{\stackrel{~}{x}}_{j}|,|{x}_{j}^{\prime }-{\stackrel{~}{x}}_{j}^{\prime }|<\epsilon .$  We assume that $\epsilon$  is so small that ${C}_{j}\cap {\stackrel{~}{C}}_{k}=\varnothing ,\text{whenever}j\ne k.$  Denote ${z}_{j},{z}_{j}^{\prime }$  the endpoints of the slit ${C}_{j}$  and ${\stackrel{~}{z}}_{j},{\stackrel{~}{z}}_{j}^{\prime }$  the corresponding endpoints of ${\stackrel{~}{C}}_{j}$  . Denote $\Psi$  the canonical map for $D$  and $\stackrel{~}{\Psi }$  the canonical map for $\stackrel{~}{D}$  . Then we need to show that  $\begin{array}{c}\stackrel{~}{\Psi }\left({\stackrel{~}{z}}_{j}\right)-\Psi \left({z}_{j}\right),\stackrel{~}{\Psi }\left({{\stackrel{~}{z}}_{j}}^{\prime }\right)-\Psi \left({z}_{j}^{\prime }\right)=O\left(\epsilon \right),j=1,...,n-1.\end{array}$ (32)
This can be shown as in the radial case by the use of an interior variation that induces a smooth mapping $z↦\stackrel{~}{z}$  from $D$  to $\stackrel{~}{D}$  which maps slit-endpoints to corresponding slit-endpoints, see [3. The non-compactness of the upper half-plane is of no concern as the mapping from $D$  to $\stackrel{~}{D}$  may be assumed to be the identity outside of a compact.
The second part of the theorem now follows from general results about ordinary differential equations, exactly as in the simply connected case.

4.2 Bilateral case

As we mentioned before, the bilateral case is similar to the radial case. The geometric description of a bilateral standard domain with $n$  boundary components requires $1+3\left(n-2\right)$  real parameters: one for the radius $Q$  of the inner circle, and three for each concentric circular slit. If ${C}_{j}$  is one of the interior slits, then ${C}_{j}=\left\{{r}_{j}{e}^{i\theta },{\theta }_{j}\le \theta \le {\theta }_{j}^{\prime }\right\}$  , and we will take ${m}_{j}=ln{r}_{j}$  , and ${\theta }_{j},{\theta }_{j}^{\prime }$  as parameters to identify ${C}_{j}$  .
If, in an arbitrary $n$  -connected domain $D$  , where $n\ge 2$  , we choose a boundary point $w$  and a boundary component that does not contain $w$  , then there is a unique conformal map from $D$  onto a bilateral standard domain, which sends $w$  to 1, and the other distinguished boundary component to the inner boundary circle of the standard domain. Thus the conformal equivalence classes of $n$  -connected domains with one marked boundary point and one distinguished boundary component which does not include the marked point are given by $1+3\left(n-2\right)=3n-5$  parameters. We call the parameters $\left(lnQ,{m}_{2},...,{m}_{n-1},{\theta }_{2},...,{\theta }_{n-1},{\theta }_{2}^{\prime },...,{\theta }_{n-1}^{\prime }\right)$  the moduli of the domain. Note that, unlike in the chordal case, these are true moduli, in the sense that different sets of parameters correspond to different conformal equivalence classes.
In the bilateral case it was natural to choose the parameter $t=lnQ$  as time. For a bilateral standard domain ${D}_{t}$  , where $t=lnQ$  , we let $\mathbf{M}\left(t\right)=\left({m}_{2}\left(t\right),...,{m}_{n-1}\left(t\right),{\theta }_{2}\left(t\right),...,{\theta }_{n-1}\left(t\right),{\theta }_{2}^{\prime }\left(t\right),...,{\theta }_{n-1}^{\prime }\left(t\right)\right).$  We then can obtain the following results in the same way as in the chordal case.
Lemma 4.2 (Motion of moduli—bilateral case). The moduli $\mathbf{M}\left(t\right)$  satisfy the system
 $\begin{array}{cc}{\partial }_{t}{m}_{j}\left(t\right)& =1-\Im \left[{\Psi }_{t}\left({m}_{j}\left(t\right){e}^{i{\theta }_{j}\left(t\right)},\xi \left(t\right)\right)-{\Psi }_{t}\left({e}^{t},\xi \left(t\right)\right)\right],\end{array}$
 $\begin{array}{cc}{\partial }_{t}{\theta }_{j}\left(t\right)& =\Re \left[{\Psi }_{t}\left({m}_{j}\left(t\right){e}^{i{\theta }_{j}\left(t\right)},\xi \left(t\right)\right)-{\Psi }_{t}\left({e}^{t},\xi \left(t\right)\right)\right],\end{array}$
 $\begin{array}{cc}{\partial }_{t}{\theta }_{j}^{\prime }\left(t\right)& =\Re \left[{\Psi }_{t}\left({m}_{j}\left(t\right){e}^{i{\theta }_{j}^{\prime }\left(t\right)},\xi \left(t\right)\right)-{\Psi }_{t}\left({e}^{t},\xi \left(t\right)\right)\right],\end{array}$ (33)
 $\begin{array}{}\end{array}$
where $j=2,...,n-1$  .
As in the radial case, it can be shown that the vector field appearing on the right above is Lipschitz in the moduli and we obtain
Theorem 4.2. Given a continuous function $t\in \left[0,\infty \right)↦\xi \left(t\right)\in {S}^{1}$  and the moduli $\mathbf{M}$  of a bilateral standard domain $D$  with interior boundary circle of radius $Q$  , there exists a unique solution $\mathbf{M}\left(t\right)$  to the system  30 on an interval $\left[lnQ,{t}_{\xi }\right)$  with $\mathbf{M}\left(0\right)=\mathbf{M}$  , and where ${t}_{\xi }$  is characterized by ${t}_{\xi }=inf\left\{\tau :{lim}_{t↗\tau }{m}_{j}\left(t\right)=0\text{for some}j\in \left\{2,...,n-1\right\}\right\}.$  Further, if ${D}_{t}$  is the bilateral standard domain determined by $\mathbf{M}\left(t\right)$  , and if ${\Psi }_{t}\left(z,\zeta \right)$  is the holomorphic vector field associated to ${D}_{t}$  (cf. Section  2.2 ), then, for any $z\in D$  , the equation ${\partial }_{t}ln{g}_{t}^{D}\left(z\right)=1+\left[{\Psi }_{t}\left({g}_{t}^{D}\left(z\right),\xi \left(t\right)\right)-{\Psi }_{t}\left({e}^{t},\xi \left(t\right)\right)\right],{g}_{lnQ}^{D}\left(z\right)=z,$  has a unique solution on $\left[lnQ,{t}_{z}\right)$  , where ${t}_{z}=sup\left\{t\le {t}_{\xi }:{inf}_{s\in \left[lnQ,t\right]}|{g}_{s}^{D}\left(z\right)-\xi \left(s\right)|>0\right\}.$  Finally, for $t<{t}_{\xi }$  set ${K}_{t}=\left\{z\in D:{t}_{z}\le t\right\}$  . Then ${g}_{t}^{D}$  is the canonical conformal map from $D\{K}_{t}$  onto ${D}_{t}$  with ${g}_{t}^{D}\left(Q\right)={e}^{t}$  .

5 Chordal SLE in multiply connected domains

5.1 Conformal Invariance and Markovian-type Property

The purpose of this paper is 1) to give a “natural” construction of conformally invariant measures on “simple curves” in multiply connected domains, and 2) to study some of the properties of these random curves. We will now motivate, using informal arguments, our particular construction of conformally invariant measures on simple curves. The arguments lead to a small class of processes which contains chordal $SL{E}_{\kappa }$  in multiply connected domains.
For a domain $D$  with $n$  non-degenerate boundary continua and two boundary points (or, more generally, prime ends) $z$  and $w$  lying on the same boundary continuum, let $W\left(D,z,w\right)$  be the set of Jordan arcs in $D$  with endpoints $z$  and $w$  . Denote $\left\{{\mathcal{ℒ}}_{D,z,w}^{\mathbf{M}}{\right\}}_{D,z,w}$  a family of probability measures on Jordan arcs in the complex plane such that ${\mathcal{ℒ}}_{D,z,w}^{\mathbf{M}}\left(W\left(D,z,w\right)\right)=1,$  and where $\mathbf{M}=M\left(D\right)$  . Such families arise, or are conjectured to arise, as distributions of interfaces of statistical mechanical systems at criticality. Based on these models, e.g. percolation, one expects that the distributions describing the interfaces in different domains with different marked points are related by a Markovian-type property and conformal invariance. Denote $\gamma$  a random Jordan arc with law ${\mathcal{ℒ}}_{D,z,w}^{\mathbf{M}}$  .
The Markovian-type property says that if ${\gamma }^{\prime }$  is a sub-arc of $\gamma$  which has $z$  as one endpoint and whose other endpoint we denote by ${z}^{\prime }$  , and if ${\mathbf{M}}^{\prime }=M\left(D\{\gamma }^{\prime }\right)$  , then the conditional law of $\gamma$  given ${\gamma }^{\prime }$  is
 $\begin{array}{c}\text{law}\left(\gamma |{\gamma }^{\prime }\right)={\mathcal{ℒ}}_{D\{\gamma }^{\prime },{z}^{\prime },w}^{{\mathbf{M}}^{\prime }}.\end{array}$ (34)
Conformal invariance means that if $f:D\to {D}^{\prime }$  is conformal, ${z}^{\prime }=f\left(z\right)$  , ${w}^{\prime }=f\left(w\right)$  , then
 $\begin{array}{c}{\mathcal{ℒ}}_{{D}^{\prime },{z}^{\prime },{w}^{\prime }}^{\mathbf{M}}={f}_{*}{\mathcal{ℒ}}_{D,z,w}^{\mathbf{M}}.\end{array}$ (35)
If  35 holds, then to understand the family $\left\{{\mathcal{ℒ}}_{D,z,w}^{\mathbf{M}}\right\}$  it is enough to consider standard domains $D$  , take $w=\infty$  , $z=0$  , and, by the identification of standard domains with their moduli, we may write ${\mathcal{ℒ}}_{D,0,\infty }^{\mathbf{M}}={\mathcal{ℒ}}^{\mathbf{M}}.$  In this case there is a natural parametrization of the Jordan arcs we consider. Let $s\in \left[0,\infty \right)↦\gamma \left(s\right)\in \overline{D}$  be a Jordan arc in a standard domain $D$  such that $\gamma \left(0\right)\in \mathbb{R},\gamma \left(0,\infty \right)\subset D,\text{and}{lim}_{t\to \infty }\gamma \left(t\right)=\infty .$  Denote $\mathbf{M}=M\left(D\right)$  the point in the moduli space corresponding to $D$  and let ${g}_{t}^{\mathbf{M}}$  be the canonical mapping from $D\\gamma \left[0,t\right]$  onto a standard domain ${D}_{t}:={g}_{t}^{\mathbf{M}}\left(D\\gamma \left[0,t\right]\right)$  . Then
 $\begin{array}{c}{g}_{t}^{\mathbf{M}}\left(z\right)=z+\frac{{a}_{t}}{z}+o\left(1/|z|\right),z\to \infty ,\end{array}$ (36)
where ${a}_{t}$  is called the half-plane capacity. The function $t↦{a}_{t}$  is continuous, strictly increasing, starts at zero and satisfies ${a}_{t}\to \infty$  as $t\to \infty$  (this final statement is not true if the curve creeps along to infinity very close to the real axis and we exclude this case for the purpose of this argument). Thus we may and always will assume that $\gamma$  is parametrized by half-plane capacity, i.e. so that ${a}_{t}=2t$  . This parametrization is natural in the following sense. If $t\ge 0$  , $\mathbf{M}\left(t\right)=M\left({D}_{t}\right)$  , and $\stackrel{~}{\gamma }$  is the curve defined by $s\in \left[0,\infty \right)↦\stackrel{~}{\gamma }\left(s\right)={g}_{t}^{\mathbf{M}}\left(\gamma \left(t+s\right)\right),$  then the canonical mapping ${g}_{s}^{\mathbf{M}\left(t\right)}$  from ${D}_{t}\\stackrel{~}{\gamma }\left[0,s\right]$  is given by ${g}_{s}^{\mathbf{M}\left(t\right)}={g}_{t+s}^{\mathbf{M}}\circ {\left({g}_{t}^{\mathbf{M}}\right)}^{-1},$  and so ${g}_{s}^{\mathbf{M}\left(t\right)}\left({D}_{t}\\stackrel{~}{\gamma }\left[0,s\right]\right)={D}_{t+s}$  . In particular, it is easy to see that ${g}_{s}^{\mathbf{M}\left(t\right)}\left(z\right)=z+\frac{2s}{z}+o\left(1/|z|\right),z\to \infty ,$  i.e. $\stackrel{~}{\gamma }$  is also parametrized by half-plane capacity.
Let now $\left\{{g}_{s}^{\mathbf{M}}:s\ge 0\right\}$  be the random family of canonical maps corresponding to the random Jordan arcs $\left\{\gamma \left[0,s\right]:s\ge 0\right\}$  in a standard domain $D$  , and denote ${\mathcal{ℒ}}^{\mathbf{M}}=\text{law}\left(\left\{{g}_{s}^{\mathbf{M}}:s\ge 0\right\}\right).$  Then, applying first the Markovian-type property and then conformal invariance,  34 ,  35 , we find $\text{law}\left(\left\{{g}_{t+s}^{\mathbf{M}}:s\ge 0\right\}|{g}_{t}^{\mathbf{M}}\right)={\left({g}_{t}^{\mathbf{M}}\right)}_{*}^{-1}{\mathcal{ℒ}}^{\mathbf{M}\left(t\right)}.$  Equivalently,
 $\begin{array}{c}\text{law}\left(\left\{{g}_{t+s}^{\mathbf{M}}\circ {\left({g}_{t}^{\mathbf{M}}\right)}^{-1}:s\ge 0\right\}|{g}_{t}^{\mathbf{M}}\right)=\text{law}\left(\left\{{g}_{s}^{\mathbf{M}\left(t\right)}:s\ge 0\right\}\right).\end{array}$ (37)
By the chordal Loewner equation,  18 , for each $t\ge 0$  , the $\sigma$  -field generated by ${g}_{t}^{\mathbf{M}}$  is equal to $\sigma \left(\left(\xi \left(r\right),\mathbf{M}\left(r\right)\right):r\in \left[0,t\right]\right)$  , where $\xi \left(0\right)=0$  .
Similarly, it is easy to see that we can reconstruct ${g}_{t+s}^{\mathbf{M}}\circ {\left({g}_{t}^{\mathbf{M}}\right)}^{-1}$  from $\left\{\left(\xi \left(t+r\right)-\xi \left(t\right),\mathbf{M}\left(t+r\right)\right):r\in \left[0,s\right]\right\}$  . Thus  37 implies
 $\begin{array}{cc}\text{law}& \left(\left\{\left(\xi \left(t+s\right)-\xi \left(t\right),\mathbf{M}\left(t+s\right)\right):s\ge 0\right\}|\left\{\left(\xi \left(r\right),\mathbf{M}\left(r\right)\right):r\in \left[0,t\right]\right\}\right)\end{array}$
 $\begin{array}{cc}& =\text{law}\left(\left\{\left(\stackrel{~}{\xi }\left(s\right),\stackrel{~}{\mathbf{M}}\left(s\right)\right):s\ge 0\right\}\right),\end{array}$ (38)
 $\begin{array}{}\end{array}$
where $\stackrel{~}{\mathbf{M}}\left(s\right)=M\left({D}_{t}\\stackrel{~}{\gamma }\left[0,s\right]\right)$  , for a random Jordan arc $\stackrel{~}{\gamma }$  with law ${\mathcal{ℒ}}^{\mathbf{M}\left(t\right)}$  . The equality  38 is precisely the statement that $\left\{\left(\xi \left(t\right),\mathbf{M}\left(t\right)\right):t\ge 0\right\}$  is a Markov process. We note that in the simply connected case ( $n=1$  ),  38 reduces to $\text{law}\left(\left\{\xi \left(t+s\right)-\xi \left(t\right):s\ge 0\right\}|\left\{\xi \left(r\right):r\in \left[0,t\right]\right\}\right)=\text{law}\left(\left\{\stackrel{~}{\xi }\left(s\right):s\ge 0\right\}\right),$  from which it follows that $\xi$  is a process with independent, and identically distributed increments. From this, continuity, and the symmetry $\text{law}\left(\xi \right)=\text{law}\left(-\xi \right)$  , Schramm derived in [24that $\xi \left(t\right)=\sqrt{\kappa }{B}_{t}$  for a standard one-dimensional Brownian motion and a positive constant $\kappa$  . The continuity follows from the continuity of the Jordan arcs, and the symmetry is actually observed in various discrete models, such as the percolation exploration process.

5.2 Scaling

For chordal SLE in the upper half-plane $\mathbb{H}$  the scaling property is usually arrived at as a consequence of the scaling property of the driving function, Brownian motion. Indeed, denote ${\partial }_{t}{g}_{t}\left(z\right)=\frac{2}{{g}_{t}\left(z\right)-\sqrt{\kappa }{B}_{t}},{g}_{0}\left(z\right)=z,$  chordal SLE in $\mathbb{H}$  and let ${K}_{t}$  be its hull at time $t$  , i.e. ${g}_{t}$  maps $\mathbb{H}\{K}_{t}$  conformally onto $\mathbb{H}$  . If $c>0$  , then ${h}_{t}$  defined by ${h}_{t}\left(z\right)=\frac{1}{c}{g}_{{c}^{2}t}\left(cz\right)$  is the normalized conformal map from $\mathbb{H}\\frac{1}{c}{K}_{{c}^{2}t}$  onto $\mathbb{H}$  and ${\partial }_{t}{h}_{t}\left(z\right)=\frac{2}{{h}_{t}\left(z\right)-\sqrt{\kappa }\frac{1}{c}{B}_{{c}^{2}t}},{h}_{0}\left(z\right)=z.$  Since $\frac{1}{c}{B}_{{c}^{2}t}$  is also a standard Brownian motion, it follows that
 $\begin{array}{c}\text{law}\left(\frac{1}{c}{K}_{{c}^{2}t}:t\ge 0\right)=\text{law}\left({K}_{t}:t\ge 0\right).\end{array}$ (39)
However, we can also turn the argument around and ask for a law on growing compacts ${K}_{t}$  in the upper half-plane which is conformally invariant, the parameter $t$  being the half-plane capacity as above. For the conformal map $z↦cz$  , this implies  39 , as the half-plane capacity scales quadratically. Denote ${\gamma }_{t}$  the tip of the curve generating ${K}_{t}$  . Then the driving function for the Loewner equation is given by ${w}_{t}={g}_{t}\left({\gamma }_{t}\right)$  , and  39 implies $\text{law}\left(\frac{1}{c}{w}_{{c}^{2}t}:t\ge 0\right)=\text{law}\left({w}_{t}:t\ge 0\right),$  i.e. the driving function has Brownian scaling. Examples of diffusion processes with Brownian scaling are multiples of Brownian motion but also Bessel processes. More generally, if $w$  satisfies the stochastic differential equation $d{w}_{t}=\sigma \left({w}_{t}\right)d{B}_{t}+b\left({w}_{t}\right)dt,$  then $w$  has Brownian scaling if $\sigma \left(cx\right)=\sigma \left(x\right),cb\left(cx\right)=b\left(x\right),$  see [20. If we assume that the coefficients $\sigma$  and $b$  are continuous, then this is saying that $\sigma$  is constant, and $b$  homogeneous of degree minus one.
In the multiply connected case we can argue similarly. Denote ${g}_{t}^{\mathbf{M}}$  the normalized conformal map from $D\left(\mathbf{M}\right)\{K}_{t}^{\mathbf{M}}$  onto ${D}_{t}$  . The superscript $\mathbf{M}$  indicates that the random compact set is a hull in the domain $D\left(\mathbf{M}\right)$  .
Conformal invariance of the growing random compacts ${K}_{t}^{\mathbf{M}}$  requires that
 $\begin{array}{c}\text{law}\left(\frac{1}{c}{K}_{{c}^{2}t}^{c\mathbf{M}}:t\ge 0\right)=\text{law}\left({K}_{t}^{\mathbf{M}}:t\ge 0\right).\end{array}$ (40)
Let ${w}_{t}^{\mathbf{M}}={g}_{t}^{\mathbf{M}}\left({\gamma }_{t}\right)$  , where ${\gamma }_{t}$  is the tip of the curve generating ${K}_{t}^{\mathbf{M}}$  .
Then  40 implies $\text{law}\left(\frac{1}{c}{w}_{{c}^{2}t}^{c\mathbf{M}},\frac{1}{c}{\mathbf{M}}_{{c}^{2}t}^{c\mathbf{M}}:t\ge 0\right)=\text{law}\left({w}_{t}^{\mathbf{M}},{\mathbf{M}}_{t}^{\mathbf{M}}:t\ge 0\right),$  where the superscript $\mathbf{M}$  indicates that ${\mathbf{M}}_{0}=\mathbf{M}$  . Thus, the moduli diffusion $\left({w}_{t},{\mathbf{M}}_{t}\right)$  also satisfies Brownian scaling. As in the one dimensional (simply connected) case, this implies under mild regularity assumptions that the coefficients of the martingale parts of the stochastic differential equation are constant, and the drift coefficients all homogeneous of degree minus one. The drift coefficients of $d{\mathbf{M}}_{t}$  are given in  30 and we check immediately that they are indeed homogeneous of degree minus one.

5.3 Moduli diffusion and interactions with the boundary

For the purposes of this subsection a different normalization of mappings on standard domains is useful. We will change the normalization of the maps ${g}_{t}$  by changing the vector field in the chordal Loewner equation  18 . For a chordal standard domain $D$  and $w\in \mathbb{R}$  , define the real function $k\left(w\right)$  by
 $\begin{array}{c}k\left(w\right)={lim}_{z\to w}\left(\Psi \left(z,w\right)+\frac{2}{z-w}\right),\end{array}$ (41)
and the conformal map ${\Psi }^{0}\left(z\right)={\Psi }^{0}\left(z,w\right)$  by $\Psi \left(z,w\right)={\Psi }^{0}\left(z,w\right)+k\left(w\right).$  Then ${\Psi }^{0}\left(z,w\right)={\Psi }_{D}^{0}\left(z,w\right)$  is the unique conformal map from $D$  onto the upper half-plane with a finite number of horizontal slits which sends $w$  to $\infty$  and satisfies ${lim}_{z\to w}\left({\Psi }^{0}\left(z,w\right)+\frac{2}{z-w}\right)=0.$  Consider the modified chordal Loewner equation
 $\begin{array}{c}{\partial }_{t}{g}_{t}^{0}\left(z\right)=-{\Psi }^{0}\left({g}_{t}^{0}\left(z\right),{\xi }^{0}\left(t\right)\right),{g}_{0}^{0}\left(z\right)=z.\end{array}$ (42)
This is the normalization used in [26. Geometrically, this normalization means that if ${g}^{0}$  removes a small vertical slit from the boundary of the upper half-plane, then the images of the two sides of this slit under ${g}^{0}$  have the same length up to first order, see [4.
Let $\kappa$  be a positive real number and $A={A}_{\kappa }\left(w,\mathbf{M}\right)$  a function homogeneous of degree minus one in the variables $w\in \mathbb{R}$  , and $\mathbf{M}$  in an open subset of ${\mathbb{R}}^{3n-3}$  . Consider the system of stochastic differential equations
 $\begin{array}{cc}d\xi \left(t\right)& =\sqrt{\kappa }d{B}_{t}+{A}_{\kappa }\left(\xi \left(t\right),{\mathbf{M}}_{t}\right)dt,\end{array}$
 $\begin{array}{cc}d{y}_{j}\left(t\right)& =\Im \left({\Psi }_{t}^{0}\left({x}_{j}\left(t\right)+i{y}_{j}\left(t\right),\xi \left(t\right)\right)\right),\end{array}$
 $\begin{array}{cc}d{x}_{j}\left(t\right)& =\Re \left({\Psi }_{t}^{0}\left({x}_{j}\left(t\right)+i{y}_{j}\left(t\right),\xi \left(t\right)\right)\right),\end{array}$
 $\begin{array}{cc}d{x}_{j}^{\prime }\left(t\right)& =\Re \left({\Psi }_{t}^{0}\left({x}_{j}^{\prime }\left(t\right)+i{y}_{j}\left(t\right),\xi \left(t\right)\right)\right),j=1,...,n-1,\end{array}$ (43)
 $\begin{array}{}\end{array}$
where ${\mathbf{M}}_{t}=\left({y}_{1}\left(t\right),...,{y}_{n-1}\left(t\right),{x}_{1}\left(t\right),...,{x}_{n-1}\left(t\right),{x}_{1}^{\prime }\left(t\right),...,{x}_{n-1}^{\prime }\left(t\right)\right)$  .
If $A$  is Lipschitz, this system has a unique solution. Then we can solve the modified chordal Loewner equation  42 for $\left(\xi \left(t\right),{\mathbf{M}}_{t}\right)$  . Denote ${K}_{t}$  the random compact such that ${g}_{t}^{0}$  maps the complement of ${K}_{t}$  in $D$  conformally onto the standard domain ${D}_{t}$  . We can interpret the term $A$  as an interaction of the random growing compact set ${K}_{t}$  with the boundary components, and it may be possible to choose $A$  so that the set ${K}_{t}$  will avoid these interior boundary components. A similar situation arises for ${\text{SLE}}_{\kappa ,\rho }$  , see [8. In that case, a random growing compact set in a simply connected domain interacts with a finite number $n$  of boundary points, the interaction strength at point $j$  being given by a real constant ${\rho }_{j}$  . Then the driving function for the chordal Loewner equation is given by the diffusion
 $\begin{array}{cc}d\nu \left(t\right)& =\sqrt{\kappa }d{B}_{t}+{\sum }_{j=1}^{n}\frac{{\rho }_{j}}{\nu \left(t\right)-{Z}^{j}\left(t\right)}dt\end{array}$
 $\begin{array}{cc}d{Z}^{j}\left(t\right)& =\frac{2}{{Z}^{j}\left(t\right)-\nu \left(t\right)}dt,j=1,...,n,\end{array}$ (44)
 $\begin{array}{}\end{array}$
a system with drift coefficients homogeneous of degree minus one similar to  43 .
There are many possible candidates for the homogeneous function $A\left(w,\mathbf{M}\right)$  . If it is to be a domain functional of the domain $D=D\left(\mathbf{M}\right)$  , then natural candidates arise from derivatives of the Green function.
Indeed, if $G\left(z,w,\mathbf{M}\right)$  is the Green function for the domain $D=D\left(\mathbf{M}\right)$  and $c>0$  , then $G\left(z,w,\mathbf{M}\right)=G\left(cz,cw,c\mathbf{M}\right)$  by conformal invariance and so $\frac{{\partial }_{z}^{k}{\partial }_{w}^{l}G\left(z,w,\mathbf{M}\right)}{{\partial }_{z}^{m}{\partial }_{w}^{n}G\left(z,w,\mathbf{M}\right)}$  is homogeneous of degree minus one whenever $k+l=m+n+1,k,l,m,n\in \mathbb{N}.$  The “harmonic random Loewner chains” Zhan studies in his thesis, see [26, correspond to the choice $k=l=m=1$  , $n=0$  . Via integration, or directly by conformal invariance, we also see that $\frac{{\partial }_{z}^{k+1}{\omega }_{j}\left(z,\mathbf{M}\right)}{{\partial }_{z}^{k}{\omega }_{j}\left(z,\mathbf{M}\right)}$  is homogeneous of degree minus one.

5.4 Chordal SLE, percolation, and locality

The case of percolation is an example where there is no interaction, that is $A\equiv 0$  . For the following calculation we return to the original chordal Loewner equation  18 . Then $\xi$  in  43 has a nonzero drift coming from changing back the normalization.Thus, to model cluster-boundaries of percolation in a multiply connected domain $D$  we make the ansatz
 $\begin{array}{c}d\xi \left(t\right)=-{k}_{t}\left(\xi \left(t\right)\right)+\sqrt{\kappa }d{B}_{t},\end{array}$ (45)
where the subscript $t$  refers to the domain ${D}_{t}$  , ${k}_{t}$  to  41 , and where $\mathbf{M}\left(t\right)$  satisfies  30 .
This choice of drift reflects that the exploration process for percolation is as likely to turn right as it is to turn left. Other discrete models lead to different drifts. In this section we show that the ansatz  45 leads to random growing compacts satisfying the locality property if $\kappa =6$  .
Denote $\left\{{g}_{t}^{E},t\ge 0\right\}$  the solution of the chordal Loewner equation in a standard domain $E$  starting at $z=0$  for the diffusion  45 . Denote $\left\{{K}_{t},t\ge 0\right\}$  the associated growing compacts. Let $A$  be a hull in $E$  that does not contain zero. For the following calculations we restrict to the event $\left\{t<\tau \right\}$  , where $\tau :=inf\left\{t:{K}_{t}\cap A\ne \varnothing \right\}$  . Let ${\Phi }_{A}$  be the canonical mapping from $E\A$  , ${g}_{t}^{*}$  the canonical mapping from ${\Phi }_{A}\left(E\\left(A\cup {K}_{t}\right)\right)$  , and ${h}_{t}$  the canonical mapping from ${g}_{t}\left(E\\left(A\cup {K}_{t}\right)\right)$  . Since the canonical mapping for $E\\left(A\cup {K}_{t}\right)$  is unique, we have
 $\begin{array}{c}{h}_{t}\circ {g}_{t}={g}_{t}^{*}\circ {\Phi }_{A}.\end{array}$ (46)
Furthermore, up to a time change, the family $\left\{{g}_{t}^{*}\right\}$  also satisfies a chordal Loewner equation beginning with the standard domain ${E}^{*}:={\Phi }_{A}\left(E\A\right)$  . In fact, reasoning as in [4, it follows that
 $\begin{array}{c}{\partial }_{t}{g}_{t}^{*}\left(z\right)=-|{h}_{t}^{\prime }\left(\xi \left(t\right)\right){|}^{2}{\Psi }_{t}^{*}\left({\xi }^{*}\left(t\right),{w}_{t}^{*}\right),\end{array}$ (47)
where ${w}_{t}^{*}={g}_{t}^{*}\left(z\right)$  , and ${\xi }^{*}\left(t\right)={h}_{t}\left(\xi \left(t\right)\right)$  . The question we are interested in is whether $\left({\xi }^{*},{\mathbf{M}}^{*}\right)$  is a time change of $\left(\xi ,\mathbf{M}\right)$  . Since ${h}_{t}={g}_{t}^{*}\circ {\Phi }_{A}\circ {g}_{t}^{-1}$  , we have
 $\begin{array}{c}{\partial }_{t}{h}_{t}\left(z\right)=\left[{\partial }_{t}{g}_{t}^{*}\right]\left({\Phi }_{A}\left({g}_{t}^{-1}\left(z\right)\right)\right)+\left({g}_{t}^{*}\circ {\Phi }_{A}{\right)}^{\prime }\left({g}_{t}^{-1}\left(z\right)\right)\left({\partial }_{t}{g}_{t}^{-1}\left(z\right)\right),\end{array}$ (48)
and we note that
 $\begin{array}{c}{\partial }_{t}{g}_{t}^{-1}\left(z\right)=\left({g}_{t}^{-1}{\right)}^{\prime }\left(z\right){\Psi }_{t}\left(\xi \left(t\right),z\right).\end{array}$ (49)
Then  48 , 47 , and  49 imply
 $\begin{array}{c}{\partial }_{t}{h}_{t}\left(z\right)=-{h}_{t}^{\prime }\left(\xi \left(t\right){\right)}^{2}{\Psi }_{t}^{*}\left({\xi }^{*}\left(t\right),{h}_{t}\left(z\right)\right)+{h}_{t}^{\prime }\left(z\right){\Psi }_{t}\left(\xi \left(t\right),z\right).\end{array}$ (50)
Hence the stochastic differential ${\partial }_{t}{h}_{t}\left(z\right)dt+{h}_{t}^{\prime }\left(\xi \left(t\right)\right)d\xi \left(t\right)$  has martingale part ${h}_{t}^{\prime }\left(\xi \left(t\right)\right)\sqrt{\kappa }d{B}_{t}$  and its drift part can be grouped into the three components
 $\begin{array}{cc}I:& =-{h}_{t}^{\prime }\left(\xi \left(t\right){\right)}^{2}\left[{\Psi }_{t}^{*}\left({\xi }^{*}\left(t\right),{h}_{t}\left(z\right)\right)-{k}_{t}^{*}\left({\xi }^{*}\left(t\right)\right)\right]dt\end{array}$
 $\begin{array}{cc}& +{h}_{t}^{\prime }\left(z\right)\left[{\Psi }_{t}\left(\xi \left(t\right),z\right)-{k}_{t}\left(\xi \left(t\right)\right)\right]dt,\end{array}$
 $\begin{array}{cc}II:& =-{h}_{t}^{\prime }\left(\xi \left(t\right){\right)}^{2}{k}_{t}^{*}\left({\xi }^{*}\left(t\right)\right)dt,\end{array}$
 $\begin{array}{cc}III:& =-\left[{h}_{t}^{\prime }\left(\xi \left(t\right)\right)-{h}_{t}^{\prime }\left(z\right)\right]{k}_{t}\left(\xi \left(t\right)\right)dt.\end{array}$ (51)
 $\begin{array}{}\end{array}$
When $z\to \xi \left(t\right)$  , then part $III$  converges to zero, and part $II$  , together with the martingale part, converges to a time-change of  45 starting at ${E}^{*}$  . Finally, for part $I$  , by the definition of $k\left(\xi ;t\right)$  a double application of l'Hôpital's rule gives
 $\begin{array}{c}{lim}_{z\to \xi }\left(\frac{2{h}^{\prime }\left(\xi {\right)}^{2}}{h\left(z\right)-h\left(\xi \right)}-\frac{2{h}^{\prime }\left(z\right)}{z-\xi }\right)=-3{h}^{\prime \prime }\left(\xi \right).\end{array}$ (52)
Thus, by Itô's formula,
 $\begin{array}{c}d{h}_{t}\left(\xi \left(t\right)\right)=-{h}_{t}^{\prime }\left(\xi \left(t\right){\right)}^{2}{k}_{t}^{*}\left({\xi }^{*}\left(t\right)\right)dt+\frac{\kappa -6}{2}{h}_{t}^{\prime \prime }\left(\xi \left(t\right)\right)dt+{h}_{t}^{\prime }\left(\xi \left(t\right)\right)\sqrt{\kappa }d{B}_{t},\end{array}$ (53)
which is indeed a time-change of  45 if and only if $\kappa =6$  . From  47 it follows immediately that the equations for ${\mathbf{M}}^{*}$  are given by the same time change of the equations for $\mathbf{M}$  .
Theorem 5.1 (Chordal ${\text{SLE}}_{6}$  ). The solution to the chordal Loewner equation based on the diffusion  45 satisfies the locality property if and only if $\kappa =6$  .
References

1. L. Ahlfors, Complex Analysis, 2nd. ed., McGraw-Hill, New York, 1966.
2. M. Aizenman, The geometry of critical percolation and conformal invariance, Stat. Phys. 19 (1996), 104–120.
3. R. Bauer, R. Friedrich, Stochastic Loewner evolution in multiply connected domain, C. R. Acad. Sci. Paris, Ser. I 339, 579-584 (2004).
4. R. Bauer, R. Friedrich, On radial stochastic Loewner evolution in multiply connected domains, arXiv.
5. A.A. Belavin, A.M. Polyakov, A.B. Zamolodchikov, Infinite conformal symmetry in two-dimensional quantum field theory, Nuclear Phys. B 241 (1984), 333–380.
6. J. Cardy, SLE(kappa,rho) and Conformal Field Theory, arXiv math-ph/0412033.
7. R. Courant, Dirichlet's Principle, with an appendix by M. Schiffer, Interscience, New York, 1950.
8. J. Dubedat, Some remarks on commutation relations for SLE, arXiv, math.PR/0411299.
9. R. Friedrich, J. Kalkkinen, On conformal field theory and stochastic Loewner evolution, Prépublications de l'IHÉS P/03/28, (2003), which appeared in Nuclear Phys. B 687 (2004), no. 3, 279–302.
10. Y. Komatu, Untersuchungen über konforme Abbildung von zweifach zusammenhängenden Gebieten, Proc. Phys.-Math. Soc. Japan (3) 25 (1943), 1–42.
11. Y. Komatu, On conformal slit mapping of multiply-connected domains, Proc. Japan Acad. 26 (1950), no. 7, 26–31.
12. M. Kontsevich, Arbeitstagung 2003 “CFT, SLE and phase boundaries” MPIM (2003).
13. R. Langlands, Y. Pouillot, Y. Saint-Aubin, Conformal invariance in two-dimensional percolation, Bull. A.M.S. 30 (1994), 1–61.
14. G. F. Lawler, O. Schramm and W. Werner, Values of Brownian intersection exponents. I. Half-plane exponents, Acta Math. 187 (2001), no. 2, 275–308.
15. G. F. Lawler, O. Schramm and W. Werner, Conformal restriction: the chordal case, J. Amer. Math. Soc. 16 (2003), 917–955.
16. G. F. Lawler, O. Schramm and W. Werner, Conformal invariance of planar loop-erased random walks and uniform spanning trees, Ann. Probab. 32 (2004), no. 1B, 939–995.
17. K. Löwner, Untersuchungen über schlichte konforme Abbildungen des Einheitskreises I., Math. Ann. 89 (1923), 103–121.
18. Z. Nehari, Conformal Mapping, McGraw-Hill, New York, 1952.
19. C.Pommerenke, Boundary behaviour of conformal maps, Grundlehren der Mathematischen Wissenschaften, 299, Springer-Verlag, Berlin, 1992.
20. D. Revuz, M. Yor, Continuous Martingales and Brownian Motion, Grundlehren der mathematischen Wissenschaften, Vol. 293, 3rd edition, Springer, Heidelberg, 1999.
21. S. Rohde, O. Schramm, Basic properties of SLE, preprint, arXiv:math.PR/0106036 v2.
22. M. Schiffer, Hadamard's formula and variation of domain-functions Amer. J. Math. 68 (1946), 417–448.
23. M. Schiffer, D. Spencer, Functionals of finite Riemann surfaces, Princeton University Press, Princeton, New Jersey, 1954.
24. O. Schramm, Scaling limits of loop-erased random walks and uniform spanning trees, Israel J. Math. 118 (2000), 221–288.
25. W. Werner, Random planar curves and Schramm-Loewner evolutions, lecture notes from the 2002 St. Flour summer school, Springer, Berlin, 2003.
26. D. Zhan, Random Loewner chains in Riemann surfaces, thesis, California Institute of Technology, 2004.

Department of Mathematics, University of Illinois at Urbana-Champaign, 1409 West Green Street, Urbana, IL 61801, USA E-mail address : rbauer@math.uiuc.edu Institute for Advanced Study, Princeton, NJ 08540, USA E-mail address : rolandf@ias.edu