Loop-erased random walk


Loop-erased random walk

In mathematics, loop-erased random walk is a model for a random simple path with important applications in combinatorics and, in physics, quantum field theory. It is intimately connected to the uniform spanning tree, a model for a random tree. [See also "random walk" for more general treatment of this topic.]

Definition

Assume "G" is some graph and gamma is some path of length "n" on "G". In other words, gamma(1),dots,gamma(n) are vertices of "G" such that gamma(i) and gamma(i+1) are neighbors. Then the loop erasure of gamma is a new simple path created by erasing all the loops of gamma in chronological order. Formally, we define indices i_j inductively using

:i_1 = 1,:i_{j+1}=max{i:gamma(i)=gamma(i_j)}+1,

where "max" here means up to the length of the path gamma. The induction stops when for some i_j we have gamma(i_j)=gamma(n). Assume this happens at "J" i.e. i_J is the last i_j. Then the loop erasure of gamma, denoted by mathrm{LE}(gamma) is a simple path of length "J" defined by

:mathrm{LE}(gamma)(j)=gamma(i_j).,

Now let "G" be some graph, let "v" be a vertex of "G", and let "R" be a random walk on "G" starting from "v". Let "T" be some stopping time for "R". Then the loop-erased random walk until time "T" is LE("R"( [1,"T"] )). In other words, take "R" from its beginning until "T" — that's a (random) path — erase all the loops in chronological order as above — you get a random simple path.

The stopping time "T" may be fixed, i.e. one may perform "n" steps and then loop-erase. However, it is usually more natural to take "T" to be the hitting time in some set. For example, let "G" be the graph Z2 and let "R" be a random walk starting from the point (0,0). Let "T" be the time when "R" first hits the circle of radius 100 (we mean here of course a "discretized" circle). LE("R") is called the loop-erased random walk starting at (0,0) and stopped at the circle.

The uniform spanning tree

Let "G" again be a graph. A spanning tree of "G" is a subgraph of "G" containing all vertices and some of the edges, which is a tree, i.e. connected and with no cycles. The uniform spanning tree (UST for short) is a "random" spanning tree chosen among all the possible spanning trees of "G" with equal probability.

Let now "v" and "w" be two vertices in "G". Any spanning tree contains precisely one simple path between "v" and "w". Taking this path in the "uniform" spanning tree gives a random simple path. It turns out that the distribution of this path is identical to the distribution of the loop-erased random walk starting at "v" and stopped at "w".

An immediate corollary is that loop-erased random walk is symmetric in its start and end points. More precisely, the distribution of the loop-erased random walk starting at "v" and stopped at "w" is identical to the distribution of the reversal of loop-erased random walk starting at "w" and stopped at "v". This is not a trivial fact at all! Loop-erasing a path and the reverse path do not give the same result. It is only the "distributions" that are identical

A-priori sampling a UST seems difficult. Even a relatively modest graph (say a 100x100 grid) has way too many spanning trees to prepare a complete list. Therefore a different approach is needed. There are a number of algorithms for sampling a UST, but we will concentrate on Wilson's algorithm.

Take any two vertices and perform loop-erased random walk from one to the other. Now take a third vertex (not on the constructed path) and perform loop-erased random walk until hitting the already constructed path. This gives a tree with three leaves. Choose a fourth vertex and do loop-erased random walk until hitting this tree. Continue until the tree spans all the vertices. It turns out that "no matter which method you use to choose the starting vertices" you always end up with the same distribution on the spanning trees, namely the uniform one.

The Laplacian random walk

Another representation of loop-erased random walk stems from solutions of the discrete Laplace equation. Let again "G" be a graph and let "v" and "w" be two vertices in "G". Construct a random path from "v" to "w" inductively using the following procedure. Assume we have already defined gamma(1),...,gamma(n). Let "f" be a function from "G" to R satisfying

:f(gamma(i))=0 for all ileq n and f(w)=1:"f" is discretely harmonic everywhere else

Where a function "f" on a graph is discretely harmonic at a point "x" if "f"("x") equals the average of "f" on the neighbors of "x".

With "f" defined choose gamma(n+1) using "f" at the neighbors of gamma(n) as weights. In other words, if x_1,...,x_d are these neighbors, choose x_i with probability

:frac{f(x_i)}{sum_{j=1}^d f(x_j)}.

Continue this process, (notice that at each step you have to calculate "f" again) and you will end up with a random simple path from "v" to "w". It turns out that the distribution of this path is identical to that of a loop-erased random walk from "v" to "w".

An alternative view is that the distribution of a loop-erased random walk conditioned to start in some path β is identical to the loop-erasure of a random walk conditioned not to hit β. This property is often referred to as the Markov property of loop-erased random walk (the relation to the usual Markov property is somewhat vague).

It is important to notice that while the proof of the equivalence is quite easy, normally models which involve dynamically changing harmonic functions or measures are extremely difficult to analyze. Practically nothing is known about the p-Laplacian walk or diffusion-limited aggregation. Another somewhat related model is the harmonic explorer.

Finally there is another link that should be mentioned: Kirchhoff's theorem relates the number of spanning trees of a graph "G" to the eigenvalues of the discrete Laplacian. See spanning tree for details.

Grids

Let "d" be the dimension, which we will assume to be at least 2. Examine Z"d" i.e. all the points (a_1,...,a_d) with integer a_i. This is an infinite graph with degree 2"d" when you connect each point to its nearest neighbors. From now on we will consider loop-erased random walk on this graph or its subgraphs.

High dimensions

The easiest case to analyze is dimension 5 and above. In this case it turns out that there the intersections are only local. A calculation shows that if one takes a random walk of length "n", its loop-erasure has length of the same order of magnitude, i.e. "n". Scaling accordingly, it turns out that loop-erased random walk converges (in an appropriate sense) to Brownian motion as "n" goes to infinity. Dimension 4 is more complicated, but the general picture is still true. It turns out that the loop-erasure of a random walk of length "n" has approximately n/log^{1/3}n vertices, but again, after scaling (that takes into account the logarithmic factor) the loop-erased walk converges to Brownian motion.

Two dimensions

In two dimensions, arguments from conformal field theory and simulation results led to a number of exciting conjectures. Assume "D" is some simply connected domain in the plane and "x" is a point in "D". Take the graph "G" to be

:G:=Dcap varepsilon mathbb{Z}^2,

that is, a grid of side length ε restricted to "D". Let "v" be the vertex of "G" closest to "x". Examine now a loop-erased random walk starting from "v" and stopped when hitting the "boundary" of "G", i.e. the vertices of "G" which correspond to the boundary of "D". Then the conjectures are
* As ε goes to zero the distribution of the path converges to some distribution on simple paths from "x" to the boundary of "D" (different from Brownian motion, of course — in 2 dimensions paths of Brownian motion are not simple). This distribution (denote it by S_{D,x}) is called the scaling limit of loop-erased random walk.
* These distributions are conformally invariant. Namely, if φ is a Riemann map between "D" and a second domain "E" then

:phi(S_{D,x})=S_{E,phi(x)}.,

*The Hausdorff dimension of these paths is 5/4 almost surely.

The first attack at these conjectures came from the direction of domino tilings. Taking a spanning tree of "G" and adding to it its planar dual one gets a domino tiling of a special derived graph (call it "H"). Each vertex of "H" corresponds to a vertex, edge or face of "G", and the edges of "H" show which vertex lies on which edge and which edge on which face. It turns out that taking a uniform spanning tree of "G" leads to a uniformly distributed random domino tiling of "H". The number of domino tilings of a graph can be calculated using the determinant of special matrices, which allow to connect it to the discrete Green function which is approximately conformally invariant. These arguments allowed to show that certain measurables of loop-erased random walk are (in the limit) conformally invariant, and that the expected number of vertices in a loop-erased random walk stopped at a circle of radius "r" is of the order of r^{5/4}.

In 2002 these conjectures were resolved (positively) using Stochastic Löwner Evolution. Very roughly, it is a stochastic conformally invariant partial differential equation which allows to catch the Markov property of loop-erased random walk (and many other probabilistic processes).

Three dimensions

A [http://arxiv.org/abs/math.PR/0508344 paper] by Gady Kozma claimsto prove that the scaling limit exists and is invariant under rotations and dilations. If we denote by L(r) the expected number of vertices in the loop-erased random walk until it gets to a distance of "r", then it was proved that

:cr^{1+varepsilon}leq L(r)leq Cr^{5/3},

where ε, "c" and "C" are some positive numbers (the numbers can, in principle, be calculated from the proofs, but the authors did not do it). This suggests that the scaling limit should have Hausdorff dimension between 1+varepsilon and 5/3 almost surely. Numerical experiments show that it should be 1.62pm 0.01.

Notes

References

*Richard Kenyon, "The asymptotic determinant of the discrete Laplacian", Acta Math. 185:2 (2000), 239-286, [http://arxiv.org/abs/math-ph/0011042 online version] .
*Richard Kenyon, "Conformal invariance of domino tiling", Ann. Probab. 28:2 (2000), 759-795, [http://arxiv.org/abs/math-ph/9910002 online version] .
*Richard Kenyon, "Long-range properties of spanning trees", Probabilistic techniques in equilibrium and nonequilibrium statistical physics, J. Math. Phys. 41:3 (2000), 1338-1363, [http://www.math.ubc.ca/~kenyon/papers/long.ps.Z online version] .
*Gady Kozma, "The scaling limit of loop-erased random walk in three dimensions", [http://arxiv.org/abs/math.PR/0508344 online version] .
*Gregory F. Lawler, "A self avoiding walk", Duke Math. J. 47 (1980), 655-694. "The original definition and a proof of the Markov property".
*Gregory F. Lawler, "The logarithmic correction for loop-erased random walk in four dimensions", Proceedings of the Conference in Honor of Jean-Pierre kahane (Orsay, 1993). Special issue of J. Fourier Anal. Appl., 347-362.
*Gregory F. Lawler, Oded Schramm, Wendelin Werner, "Conformal invariance of planar loop-erased random walks and uniform spanning trees", Ann. Probab. 32:1B (2004), 939-995, [http://arxiv.org/abs/math.PR/0112234 online version] .
*Robin Pemantle, "Choosing a spanning tree for the integer lattice uniformly", Ann. Probab. 19:4 (1991), 1559-1574.
* Oded Schramm, "Scaling limits of loop-erased random walks and uniform spanning trees", Israel J. Math. 118 (2000), 221-288.
* David Bruce Wilson, "Generating random spanning trees more quickly than the cover time", Proceedings of the Twenty-eighth Annual ACM Symposium on the Theory of Computing (Philadelphia, PA, 1996), 296-303, ACM, New York, 1996.


Wikimedia Foundation. 2010.

Look at other dictionaries:

  • Random walk — A random walk, sometimes denoted RW, is a mathematical formalization of a trajectory that consists of taking successive random steps. The results of random walk analysis have been applied to computer science, physics, ecology, economics and a… …   Wikipedia

  • Schramm–Loewner evolution — In probability theory, Schramm–Loewner evolution, also known as stochastic Loewner evolution or SLE, is a conformally invariant stochastic process. It is a family of random planar curves that are generated by solving Charles Loewner s… …   Wikipedia

  • Oded Schramm — (hebräisch ‏עודד שרם‎; * 10. Dezember 1961 in Jerusalem; † 1. September 2008 am Guye Peak in Washington) war ein israelischer Mathematiker und mathematischer Physiker, der sich vor allem mit mathematischer statistischer P …   Deutsch Wikipedia

  • List of probability topics — This is a list of probability topics, by Wikipedia page. It overlaps with the (alphabetical) list of statistical topics. There are also the list of probabilists and list of statisticians.General aspects*Probability *Randomness, Pseudorandomness,… …   Wikipedia

  • List of mathematics articles (L) — NOTOC L L (complexity) L BFGS L² cohomology L function L game L notation L system L theory L Analyse des Infiniment Petits pour l Intelligence des Lignes Courbes L Hôpital s rule L(R) La Géométrie Labeled graph Labelled enumeration theorem Lack… …   Wikipedia

  • List of stochastic processes topics — In the mathematics of probability, a stochastic process can be thought of as a random function. In practical applications, the domain over which the function is defined is a time interval ( time series ) or a region of space ( random field… …   Wikipedia

  • Arbre (probabilité) — Pour les articles homonymes, voir Arbre (homonymie). En théorie des probabilités un arbre aléatoire est un arbre défini en utilisant une loi de probabilité sur un ensemble d arbres (au sens de graphe). Par exemple, un arbre aléatoire à n nœuds… …   Wikipédia en Français

  • Charles Loewner — in 63 Born 29 May 1893(1893 05 29) Lány …   Wikipedia

  • Charles Loewner — C. Loewner à droite (inconnu à gauche) en 1927 (coll. MFO) Naissance …   Wikipédia en Français

  • Oded Schramm — Born December 10, 1961(1961 12 10) Jerusalem, Israel …   Wikipedia