Schramm–Loewner evolution

From HandWiki
Schramm–Loewner evolution on the upper half plane with hue indicating [math]\displaystyle{ \operatorname{Im}(\log(g_t(z))) }[/math]

In probability theory, the Schramm–Loewner evolution with parameter κ, also known as stochastic Loewner evolution (SLEκ), is a family of random planar curves that have been proven to be the scaling limit of a variety of two-dimensional lattice models in statistical mechanics. Given a parameter κ and a domain in the complex plane U, it gives a family of random curves in U, with κ controlling how much the curve turns. There are two main variants of SLE, chordal SLE which gives a family of random curves from two fixed boundary points, and radial SLE, which gives a family of random curves from a fixed boundary point to a fixed interior point. These curves are defined to satisfy conformal invariance and a domain Markov property.

It was discovered by Oded Schramm (2000) as a conjectured scaling limit of the planar uniform spanning tree (UST) and the planar loop-erased random walk (LERW) probabilistic processes, and developed by him together with Greg Lawler and Wendelin Werner in a series of joint papers.

Besides UST and LERW, the Schramm–Loewner evolution is conjectured or proven to describe the scaling limit of various stochastic processes in the plane, such as critical percolation, the critical Ising model, the double-dimer model, self-avoiding walks, and other critical statistical mechanics models that exhibit conformal invariance. The SLE curves are the scaling limits of interfaces and other non-self-intersecting random curves in these models. The main idea is that the conformal invariance and a certain Markov property inherent in such stochastic processes together make it possible to encode these planar curves into a one-dimensional Brownian motion running on the boundary of the domain (the driving function in Loewner's differential equation). This way, many important questions about the planar models can be translated into exercises in Itô calculus. Indeed, several mathematically non-rigorous predictions made by physicists using conformal field theory have been proven using this strategy.

The Loewner equation

Main page: Loewner differential equation

If [math]\displaystyle{ D }[/math] is a simply connected, open complex domain not equal to [math]\displaystyle{ \mathbb{C} }[/math], and [math]\displaystyle{ \gamma }[/math] is a simple curve in [math]\displaystyle{ D }[/math] starting on the boundary (a continuous function with [math]\displaystyle{ \gamma(0) }[/math] on the boundary of [math]\displaystyle{ D }[/math] and [math]\displaystyle{ \gamma((0,\infty)) }[/math] a subset of [math]\displaystyle{ D }[/math]), then for each [math]\displaystyle{ t\geq 0 }[/math], the complement [math]\displaystyle{ D_t = D\smallsetminus\gamma([0,t]) }[/math] of [math]\displaystyle{ \gamma([0,t]) }[/math] is simply connected and therefore conformally isomorphic to [math]\displaystyle{ D }[/math] by the Riemann mapping theorem. If [math]\displaystyle{ f_t }[/math] is a suitable normalized isomorphism from [math]\displaystyle{ D }[/math] to [math]\displaystyle{ D_t }[/math], then it satisfies a differential equation found by (Loewner 1923) in his work on the Bieberbach conjecture. Sometimes it is more convenient to use the inverse function [math]\displaystyle{ g_t }[/math] of [math]\displaystyle{ f_t }[/math], which is a conformal mapping from [math]\displaystyle{ D_t }[/math] to [math]\displaystyle{ D }[/math].

In Loewner's equation, [math]\displaystyle{ z\in D }[/math], [math]\displaystyle{ t\geq 0 }[/math], and the boundary values at time [math]\displaystyle{ t=0 }[/math] are [math]\displaystyle{ f_0(z)=z }[/math] or [math]\displaystyle{ g_0(z)=z }[/math]. The equation depends on a driving function [math]\displaystyle{ \zeta(t) }[/math] taking values in the boundary of [math]\displaystyle{ D }[/math]. If [math]\displaystyle{ D }[/math] is the unit disk and the curve [math]\displaystyle{ \gamma }[/math] is parameterized by "capacity", then Loewner's equation is

[math]\displaystyle{ \frac{\partial f_t(z)}{\partial t} = -z f^\prime_t(z)\frac{\zeta(t)+z}{\zeta(t)-z} }[/math]   or   [math]\displaystyle{ \dfrac{\partial g_t(z)}{\partial t} = g_t(z)\dfrac{\zeta(t)+g_t(z)}{\zeta(t)-g_t(z)}. }[/math]

When [math]\displaystyle{ D }[/math] is the upper half plane the Loewner equation differs from this by changes of variable and is

[math]\displaystyle{ \frac{\partial f_t(z)}{\partial t} = \frac{ 2f_t^\prime(z)}{\zeta(t)-z} }[/math]   or   [math]\displaystyle{ \dfrac{\partial g_t(z)}{\partial t} = \dfrac{2}{g_t(z)-\zeta(t)}. }[/math]

The driving function [math]\displaystyle{ \zeta }[/math] and the curve [math]\displaystyle{ \gamma }[/math] are related by

[math]\displaystyle{ f_t(\zeta(t)) = \gamma(t) \text{ or } \zeta(t) = g_t(\gamma(t)) }[/math]

where [math]\displaystyle{ f_t }[/math] and [math]\displaystyle{ g_t }[/math] are extended by continuity.

Example

Let [math]\displaystyle{ D }[/math] be the upper half plane and consider an SLE0, so the driving function [math]\displaystyle{ \zeta }[/math] is a Brownian motion of diffusivity zero. The function [math]\displaystyle{ \zeta }[/math] is thus identically zero almost surely and

[math]\displaystyle{ f_t(z) = \sqrt{z^2-4t} }[/math]
[math]\displaystyle{ g_t(z) = \sqrt{z^2+4t} }[/math]
[math]\displaystyle{ \gamma(t) = 2i\sqrt{t} }[/math]
[math]\displaystyle{ D_t }[/math] is the upper half-plane with the line from 0 to [math]\displaystyle{ 2i\sqrt{t} }[/math] removed.

Schramm–Loewner evolution

Schramm–Loewner evolution is the random curve γ given by the Loewner equation as in the previous section, for the driving function

[math]\displaystyle{ \zeta(t)=\sqrt{\kappa}B(t) }[/math]

where B(t) is Brownian motion on the boundary of D, scaled by some real κ. In other words, Schramm–Loewner evolution is a probability measure on planar curves, given as the image of Wiener measure under this map.

In general the curve γ need not be simple, and the domain Dt is not the complement of γ([0,t]) in D, but is instead the unbounded component of the complement.

There are two versions of SLE, using two families of curves, each depending on a non-negative real parameter κ:

  • Chordal SLEκ, which is related to curves connecting two points on the boundary of a domain (usually the upper half plane, with the points being 0 and infinity).
  • Radial SLEκ, which is related to curves joining a point on the boundary of a domain to a point in the interior (often curves joining 1 and 0 in the unit disk).

SLE depends on a choice of Brownian motion on the boundary of the domain, and there are several variations depending on what sort of Brownian motion is used: for example it might start at a fixed point, or start at a uniformly distributed point on the unit circle, or might have a built in drift, and so on. The parameter κ controls the rate of diffusion of the Brownian motion, and the behavior of SLE depends critically on its value.

The two domains most commonly used in Schramm–Loewner evolution are the upper half plane and the unit disk. Although the Loewner differential equation in these two cases look different, they are equivalent up to changes of variables as the unit disk and the upper half plane are conformally equivalent. However a conformal equivalence between them does not preserve the Brownian motion on their boundaries used to drive Schramm–Loewner evolution.

Special values of κ

  • For 0 ≤ κ < 4 the curve γ(t) is simple (with probability 1).
  • For 4 < κ < 8 the curve γ(t) intersects itself and every point is contained in a loop but the curve is not space-filling (with probability 1).
  • For κ ≥ 8 the curve γ(t) is space-filling (with probability 1).
  • κ = 2 corresponds to the loop-erased random walk, or equivalently, branches of the uniform spanning tree.
  • For κ = 8/3, SLEκ has the restriction property and is conjectured to be the scaling limit of self-avoiding random walks. A version of it is the outer boundary of Brownian motion.
  • κ = 3 is the limit of interfaces for the Ising model.
  • κ = 4 corresponds to the path of the harmonic explorer and contour lines of the Gaussian free field.
  • Percolation interface: Make a lozenge of equally sized hexagons in the plane. Color its upper and left side with black, and the lower and right side with white. Then color the other hexagons “white” or “black” independently with equal probability 1/2. There is a boundary between the black and the white, running from bottom left to top right. The scaling limit of the boundary is κ = 6.
  • For κ = 6, SLEκ has the locality property. This arises in the scaling limit of critical percolation on the triangular lattice and conjecturally on other lattices.
  • κ = 8 corresponds to the path separating the uniform spanning tree from its dual tree.

When SLE corresponds to some conformal field theory, the parameter κ is related to the central charge c of the conformal field theory by

[math]\displaystyle{ c = \frac{(8-3\kappa)(\kappa-6)}{2\kappa}. }[/math]

Each value of c < 1 corresponds to two values of κ, one value κ between 0 and 4, and a "dual" value 16/κ greater than 4. (see (Bauer Bernard) (Bauer Bernard))

(Beffara 2008) showed that the Hausdorff dimension of the paths (with probability 1) is equal to min(2, 1 + κ/8).

Left passage probability formulas for SLEκ

The probability of chordal SLEκ γ being on the left of fixed point [math]\displaystyle{ x_{0}+iy_{0}=z_{0}\in \mathbb{H} }[/math] was computed by (Schramm 2001a)[1]

[math]\displaystyle{ \mathbb{P}[\gamma \text{ passes to the left } z_0]=\frac{1}{2}+\frac{\Gamma(\frac{4}{\kappa})}{\sqrt{\pi} \, \Gamma(\frac{8-\kappa}{2\kappa})}\frac{x_0}{y_0} \, _2F_1 \left(\frac{1}{2},\frac{4}{\kappa}, \frac{3}{2}, - \left(\frac{x_0}{y_0}\right)^2 \right) }[/math]

where [math]\displaystyle{ \Gamma }[/math] is the Gamma function and [math]\displaystyle{ _2F_{1}(a,b,c,d) }[/math] is the hypergeometric function. This was derived by using the martingale property of

[math]\displaystyle{ h(x,y):=\mathbb{P}[\gamma \text{ passes to the left } x+iy] }[/math]

and Itô's lemma to obtain the following partial differential equation for [math]\displaystyle{ w:=\tfrac{x}{y} }[/math]

[math]\displaystyle{ \frac{\kappa}{2}\partial_{ww}h(w)+\frac{4w}{w^2+1}\partial_w h=0. }[/math]

For κ = 4, the RHS is [math]\displaystyle{ 1-\tfrac{1}{\pi}\arg(z_0) }[/math], which was used in the construction of the harmonic explorer,[2] and for κ = 6, we obtain Cardy's formula, which was used by Smirnov to prove conformal invariance in percolation.[3]

Applications

(Lawler Schramm) used SLE6 to prove the conjecture of (Mandelbrot 1982) that the boundary of planar Brownian motion has fractal dimension 4/3.

Critical percolation on the triangular lattice was proved to be related to SLE6 by Stanislav Smirnov.[4] Combined with earlier work of Harry Kesten,[5] this led to the determination of many of the critical exponents for percolation.[6] This breakthrough, in turn, allowed further analysis of many aspects of this model.[7][8]

Loop-erased random walk was shown to converge to SLE2 by Lawler, Schramm and Werner.[9] This allowed derivation of many quantitative properties of loop-erased random walk (some of which were derived earlier by Richard Kenyon[10]). The related random Peano curve outlining the uniform spanning tree was shown to converge to SLE8.[9]

Rohde and Schramm showed that κ is related to the fractal dimension of a curve by the following relation

[math]\displaystyle{ d = 1 + \frac{\kappa}{8}. }[/math]

Simulation

Computer programs (Matlab) are presented in this GitHub repository to simulate Schramm Loewner Evolution planar curves.

References

  1. Schramm, Oded (2001a), "Percolation formula.", Electron. Comm. 33 (6): 115–120, Bibcode2001math......7096S 
  2. Schramm, Oded; Sheffield, Scott (2005), "Harmonic explorer and its convergence to SLE4.", Annals of Probability 33 (6): 2127–2148, doi:10.1214/009117905000000477 
  3. Smirnov, Stanislav (2001). "Critical percolation in the plane: conformal invariance, Cardy's formula, scaling limits". Comptes Rendus de l'Académie des Sciences, Série I 333 (3): 239–244. doi:10.1016/S0764-4442(01)01991-7. ISSN 0764-4442. Bibcode2001CRASM.333..239S. 
  4. Smirnov, Stanislav (2001). "Critical percolation in the plane". Comptes Rendus de l'Académie des Sciences 333 (3): 239–244. doi:10.1016/S0764-4442(01)01991-7. Bibcode2001CRASM.333..239S. 
  5. Kesten, Harry (1987). "Scaling relations for 2D-percolation". Comm. Math. Phys. 109 (1): 109–156. doi:10.1007/BF01205674. Bibcode1987CMaPh.109..109K. http://conservancy.umn.edu/bitstream/11299/4435/1/241.pdf. 
  6. Smirnov, Stanislav; Werner, Wendelin (2001). "Critical exponents for two-dimensional percolation". Math. Res. Lett. 8 (6): 729–744. doi:10.4310/mrl.2001.v8.n6.a4. 
  7. Schramm, Oded; Steif, Jeffrey E. (2010). "Quantitative noise sensitivity and exceptional times for percolation". Ann. of Math. 171 (2): 619–672. doi:10.4007/annals.2010.171.619. 
  8. Garban, Christophe; Pete, Gábor; Schramm, Oded (2013). "Pivotal, cluster and interface measures for critical planar percolation". J. Amer. Math. Soc. 26 (4): 939–1024. doi:10.1090/S0894-0347-2013-00772-9. 
  9. 9.0 9.1 Lawler, Gregory F.; Schramm, Oded; Werner, Wendelin (2004). "Conformal invariance of planar loop-erased random walks and uniform spanning trees". Ann. Probab. 32 (1B): 939–995. doi:10.1214/aop/1079021469. 
  10. Kenyon, Richard (2000). "Long range properties of spanning trees". J. Math. Phys. 41 (3): 1338–1363. doi:10.1063/1.533190. Bibcode2000JMP....41.1338K. 

Further reading

External links