← All research articles

Research article · Mathematics

Why randomness returns in 2D—but may escape in 3D

At each step, a walker randomly chooses one neighbouring direction. On a line and on a plane it returns to its starting point almost surely. In space there is a positive probability that it never returns. One extra dimension changes the answer.

One-sentence conclusion

For simple symmetric random walk on Zd\mathbb{Z}^d, the probability of being at the origin after 2n steps decays as n−d/2n^{-d/2}; its sum diverges in dimensions 1 and 2 but converges from d=3d = 3 onward, exactly the criterion separating almost-sure return from possible escape.

1. The problem and precise definitions

The walk starts at S0=0S_0 = 0 on the integer lattice Zd\mathbb{Z}^d. At each time it adds, with equal probability, one of the 2d vectors ±e1,…,±ed\pm e_1,\ldots,\pm e_d: exactly one coordinate changes by one unit. Steps are independent. The model is “simple” because only nearest neighbours are allowed and “symmetric” because no direction is favoured.

The first return time is τ0=inf⁡{n≥1:Sn=0}\tau_0 = \inf\{n \ge 1 : S_n = 0\}. The walk is recurrent if P(τ0<∞)=1\mathbb{P}(\tau_0 < \infty) = 1 and transient if this probability is strictly below 1. “Almost surely” permits exceptional paths of probability zero and gives no maximum return time.

The question is therefore sharp: in which dimensions do the opportunities for return remain abundant enough to make a return certain? Pólya answered it in 1921: recurrence for d=1,2d = 1,2; transience for d≥3d \ge 3.[1][9]

2. Search and source selection

The search was conducted on 11 September 2026 across EuDML, DOI and publisher pages, the Dartmouth, Chicago and Duke author archives, and the Cambridge, Wiley and AMS catalogues. Queries combined “Pólya random walk recurrence”, “lattice Green function”, “local central limit theorem”, “Chung–Fuchs criterion”, “multiple returns” and “electric networks recurrence”.

Ten central references were retained: the original paper, the general Chung–Fuchs criterion, two historical works on lattice integrals and multiple returns, two classic monographs, three modern or pedagogical treatments, and one extension to infinite graphs. Popular accounts without proofs, irreproducible simulations and variants remote from Zd\mathbb{Z}^d were excluded.

Bibliographic pages were used to verify editions, DOIs and scope; fully accessible texts were used to check definitions, the Green-function equivalence, the asymptotic proof and the electrical analogy. No numerical result is borrowed from a third-party simulation.

3. From returns to expected visits

Set pn=P(Sn=0)p_n = \mathbb{P}(S_n = 0) and G=∑n=0∞pnG = \sum_{n=0}^{\infty} p_n. The quantity G, the Green function at the origin, is the expected number of visits to the origin including time zero. It adds the probabilities of being there at each exact time n; it does not assume that those events are independent.

Let F=P(τ0<∞)F = \mathbb{P}(\tau_0 < \infty), the probability of returning at least once. After every return, the Markov property restarts the same problem. The number of visits therefore has the geometric structure 1+F+F2+⋯1 + F + F^2 + \cdots: if F<1F < 1, then G=11−FG = \frac{1}{1-F}; if F=1F = 1, G diverges. Proving recurrence is thus exactly the same as testing whether ∑pn\sum p_n diverges.[8][9]

This avoids a trap: divergence of ∑pn\sum p_n does not generally justify a naïve Borel–Cantelli argument because return events are dependent. The Markov restart is what makes the equivalence exact here.

G=∑n=0∞P(Sn=0)F=P(τ0<∞)G=∑k=0∞Fk={11−F,F<1,+∞,F=1.\begin{aligned}G &= \sum_{n=0}^{\infty}\mathbb{P}(S_n=0) \\[0.6em] F &= \mathbb{P}(\tau_0<\infty) \\[0.6em] G &= \sum_{k=0}^{\infty}F^k = \begin{cases}\dfrac{1}{1-F}, & F<1, \\[0.5em]+\infty, & F=1.\end{cases}\end{aligned}

4. Why the threshold is two

The characteristic function of one step is ϕ(θ)=1d∑j=1dcos⁡θj\phi(\theta) = \frac{1}{d}\sum_{j=1}^{d}\cos\theta_j. Fourier inversion gives pn=(2π)−d∫[−π,π]dϕ(θ)n dθp_n = (2\pi)^{-d}\int_{[-\pi,\pi]^d}\phi(\theta)^n\,\mathrm{d}\theta. To sum safely, introduce 0<s<10 < s < 1 and then let s tend to one: G(s)=∑n=0∞snpn=(2π)−d∫[−π,π]ddθ1−sϕ(θ)G(s) = \sum_{n=0}^{\infty}s^np_n = (2\pi)^{-d}\int_{[-\pi,\pi]^d}\frac{\mathrm{d}\theta}{1-s\phi(\theta)}.[2][8]

The entire question concentrates near θ=0\theta = 0, where cos⁡x=1−x22+o(x2)\cos x = 1 - \frac{x^2}{2} + o(x^2). Consequently, 1−ϕ(θ)1 - \phi(\theta) behaves like ∥θ∥22d\frac{\lVert\theta\rVert^2}{2d}. The Green kernel therefore behaves like a constant divided by r2r^2, where r=∥θ∥r = \lVert\theta\rVert.

In radial coordinates, a shell of radius r contributes a volume factor proportional to rd−1 drr^{d-1}\,\mathrm{d}r. The decisive integral becomes ∫0εrd−3 dr\int_0^\varepsilon r^{d-3}\,\mathrm{d}r. It diverges for d=1d = 1 and d=2d = 2 but converges for d≥3d \ge 3. By the previous equivalence, this is exactly the recurrence–transience threshold.

1−ϕ(θ)∼∥θ∥22d(θ→0)C∫0εrd−1r2 dr=C∫0εrd−3 dr∫0εrd−3 dr=+∞  ⟺  d≤2\begin{gathered}1-\phi(\theta)\sim\frac{\lVert\theta\rVert^2}{2d}\quad(\theta\to0) \\[0.9em] C\int_0^\varepsilon\frac{r^{d-1}}{r^2}\,\mathrm{d}r = C\int_0^\varepsilon r^{d-3}\,\mathrm{d}r \\[0.9em] \int_0^\varepsilon r^{d-3}\,\mathrm{d}r=+\infty \iff d\le2\end{gathered}

5. Reproducible calculation in dimensions 1–3

Returns are possible only at even times. In one dimension, p2n=(2nn)4np_{2n} = \frac{\binom{2n}{n}}{4^n}. In two dimensions, the combinatorial identity p2n=[(2nn)4n]2p_{2n} = \left[\frac{\binom{2n}{n}}{4^n}\right]^2 yields the n−1/2n^{-1/2} and n−1n^{-1} decay rates. In three dimensions, an exact recurrence for multinomial coefficients computes p2np_{2n} without simulation.[4][5][9]

The script sums terms through n=100,000n = 100{,}000, or 200,000 steps. For d=3d = 3 the partial sum is 1.514910606. A first-order local-central-limit tail correction gives G3≈1.516386065G_3 \approx 1.516386065 and hence F3=1−1G3≈0.340537332F_3 = 1 - \frac{1}{G_3} \approx 0.340537332. The difference from the published constant 0.3405373296 is 2.36×10−92.36 \times 10^{-9}; this agreement is a numerical check, not the proof of transience.[3][6][9]

Logarithmic curves of return probabilities after 2n steps in dimensions one, two and three.
The slopes −12-\frac{1}{2}, −1-1 and −32-\frac{3}{2} reveal the change in the series: only the last is summable.
Probability of being exactly back at the origin at the stated time.
Total stepsd=1d = 1d=2d = 2d=3d = 3
200.1761970.03104540.00710835
2,0000.01783900.0003182300.00000737453
200,0000.001784120.000003183090.00000000737727
Reproduce or auditCalculation script (.mjs)Figure data (.csv)Method and results (.json)

6. Three intuitions to correct

  1. Almost-sure return does not mean quick return. In dimension two the sum diverges only logarithmically: returns become rare, but not fast enough for their total weight to remain finite.
  2. Recurrent does not mean confined. In dimensions one and two the walk reaches arbitrarily remote positions and still returns; the theorem concerns return, not bounded distance from the origin.
  3. Three recurrent coordinates do not make their coincidence recurrent. A three-dimensional return requires all three coordinates to be zero at the same time. That extra synchronisation changes the decay to n−3/2n^{-3/2}, which is summable.

7. What the theorem does not say

  • The proof targets a simple, symmetric, drift-free walk with independent steps on the integer lattice. Non-zero drift already makes a one-dimensional walk transient; heavy-tailed jumps can also alter the threshold.
  • Euclidean dimension alone does not classify every graph. Volume growth, resistance, isoperimetry and global geometry can change recurrence; a graph drawn in the plane is not automatically recurrent.[2][10]
  • The 200,000-step calculation uses exact recurrences in floating-point arithmetic and then an asymptotic tail extrapolation in dimension three. Its nine-decimal result is checked against the literature but is not presented as an interval-arithmetic certificate.
  • The stated problem is nevertheless closed: for simple symmetric random walk on Zd\mathbb{Z}^d, the theorem and proof completely determine the answer. Open questions begin with other graphs, dependencies or step laws.

8. Central references actually consulted

Each note identifies the source’s role. Links point to the text, DOI or official publisher page.

  1. G. Pólya (1921). Über eine Aufgabe der Wahrscheinlichkeitsrechnung betreffend die Irrfahrt im Straßennetz. Mathematische Annalen 84, 149–160.

    The original paper establishing recurrence on one- and two-dimensional lattices and transience from dimension three onward.

  2. K. L. Chung & W. H. J. Fuchs (1951). On the Distribution of Values of Sums of Random Variables. Memoirs of the American Mathematical Society, no. 6.

    Generalises recurrence beyond the simple walk and underpins the characteristic-function criterion.

  3. G. N. Watson (1939). Three Triple Integrals. Quarterly Journal of Mathematics os-10(1), 266–276.

    Evaluates lattice integrals that enter the cubic lattice Green constant.

  4. W. Feller (1968). An Introduction to Probability Theory and Its Applications, Volume I. 3rd ed., Wiley.

    Reference for random walks, recurrent events, generating functions, and Stirling’s approximation.

  5. C. Domb (1954). On Multiple Returns in the Random-Walk Problem. Proceedings of the Cambridge Philosophical Society 50(4), 586–591.

    Connects fixed-time return probabilities with asymptotic methods and multiple returns.

  6. E. W. Montroll (1956). Random Walks in Multidimensional Spaces, Especially on Periodic Lattices. Journal of the Society for Industrial and Applied Mathematics 4(4), 241–260.

    Develops Green-function and return-probability calculations for several periodic lattices.

  7. P. G. Doyle & J. L. Snell (1984/2006). Random Walks and Electric Networks. Carus Mathematical Monographs 22; open 2006 edition.

    Gives a second reading of the dimensional threshold through effective resistance from the origin to infinity.

  8. G. F. Lawler & V. Limic (2010). Random Walk: A Modern Introduction. Cambridge Studies in Advanced Mathematics 123.

    Covers the local central limit theorem, Green functions, and the equivalence between recurrence and divergent expected visit counts.

  9. R. Durrett (2019). Probability: Theory and Examples. 5th ed., Section 5.4.

    Provides a modern combinatorial proof, the dimension-one-to-three asymptotics, and the numerical three-dimensional return constant.

  10. W. Woess (2000). Random Walks on Infinite Graphs and Groups. Cambridge Tracts in Mathematics 138.

    Shows why graph geometry, rather than its Euclidean drawing alone, governs recurrence beyond regular lattices.

Educational article about an ideal mathematical model. Illustrated paths are not observed data.