Pólya Web: Coalescing Random Walks on ℕ²

• Pólya Web is a system of coalescing random walks on ℕ² where each walk follows state-dependent probabilities based on its current coordinates.
• It establishes key connections between local bias, negative association via the BKR inequality, and the strong law showing cluster count scaling as approximately √(πn).
• The model’s edge-scaling limit transitions to the Yule Web, where discrete Pólya walks converge to continuum processes with Gamma-distributed martingale limits.

The Pólya Web is a system of coalescing random walks founded on the classical Pólya urn process, situated on the integer lattice $\mathbb{N}^2$. It serves as an analogue to the up-right oriented web of coalescing random walks introduced by Tóth and Werner (1998), but replaces simple symmetric random walks with Pólya walks as the fundamental constituents. The Pólya Web reveals deep connections between the geometry of random walks with state-dependent bias, negative dependence structures, strong laws for clustering phenomena, determinantal formulas for joint limit distributions, and continuum scaling limits leading to the Yule Web of coalescing Poissonian processes (Urbán, 17 Jan 2026).

1. Construction of the Pólya Web

A single Pólya walk is defined via an inhomogeneous Markov process on $\mathbb{N}^2$. Given initial state $(X_0, Y_0) = (x_0, y_0)$, the transitions are

$$P\bigl((X_{n+1}, Y_{n+1}) = (x+1, y)\mid (X_n, Y_n) = (x, y)\bigr) = \frac{x}{x+y},$$

$$P\bigl((X_{n+1}, Y_{n+1}) = (x, y+1)\mid (X_n, Y_n) = (x, y)\bigr) = \frac{y}{x+y},$$

where the walk increment is either in the $x$-direction or $y$-direction with probability proportional to the current coordinate values.

A lattice-arrow representation is available: assign to each site $\lambda = (k, \ell)$ an independent arrow $\omega_\lambda \in \{(1,0), (0,1)\}$, with

$$P(\omega_\lambda = (1,0)) = \frac{k}{k+\ell}, \qquad P(\omega_\lambda = (0,1)) = \frac{\ell}{k+\ell}.$$

Given $\mathbb{N}^2$0, the walk $\mathbb{N}^2$1 born at time $\mathbb{N}^2$2 starts at $\mathbb{N}^2$3 and follows the arrows: $\mathbb{N}^2$4.

In the Pólya Web, a Pólya walk is initiated at every lattice site. Trajectories that intersect at any site coalesce and continue together thereafter, forming an up-right web of coalescing Pólya walks.

2. Negative Association and the van den Berg–Kesten–Reimer Inequality

Consider a finite region $\mathbb{N}^2$5 and product space $\mathbb{N}^2$6 under the natural product measure. For two neighboring start-points $\mathbb{N}^2$7 and $\mathbb{N}^2$8 at level $\mathbb{N}^2$9, define the indicator

$(X_0, Y_0) = (x_0, y_0)$0

As $(X_0, Y_0) = (x_0, y_0)$1, $(X_0, Y_0) = (x_0, y_0)$2never coalesce$(X_0, Y_0) = (x_0, y_0)$3.

The geometry of the dual web ensures that “no-meet” events for disjoint sets $(X_0, Y_0) = (x_0, y_0)$4 of indices are certified on disjoint subsets of arrow variables, enabling the use of the van den Berg–Kesten–Reimer (BKR) inequality: $(X_0, Y_0) = (x_0, y_0)$5 for any two events $(X_0, Y_0) = (x_0, y_0)$6.

In the current context, for disjoint sets $(X_0, Y_0) = (x_0, y_0)$7,

$(X_0, Y_0) = (x_0, y_0)$8

Letting $(X_0, Y_0) = (x_0, y_0)$9 yields negative association for the family $$P\bigl((X_{n+1}, Y_{n+1}) = (x+1, y)\mid (X_n, Y_n) = (x, y)\bigr) = \frac{x}{x+y},$$0: for any two increasing functions $$P\bigl((X_{n+1}, Y_{n+1}) = (x+1, y)\mid (X_n, Y_n) = (x, y)\bigr) = \frac{x}{x+y},$$1 of disjoint subsets,

$$P\bigl((X_{n+1}, Y_{n+1}) = (x+1, y)\mid (X_n, Y_n) = (x, y)\bigr) = \frac{x}{x+y},$$2

This structure is pivotal in establishing rigorous large deviations, variance bounds, and limit laws.

3. Strong Law for the Number of Clusters

For walks started at level $$P\bigl((X_{n+1}, Y_{n+1}) = (x+1, y)\mid (X_n, Y_n) = (x, y)\bigr) = \frac{x}{x+y},$$3, coalescence partitions the $$P\bigl((X_{n+1}, Y_{n+1}) = (x+1, y)\mid (X_n, Y_n) = (x, y)\bigr) = \frac{x}{x+y},$$4 start-points $$P\bigl((X_{n+1}, Y_{n+1}) = (x+1, y)\mid (X_n, Y_n) = (x, y)\bigr) = \frac{x}{x+y},$$5 into clusters, with $$P\bigl((X_{n+1}, Y_{n+1}) = (x+1, y)\mid (X_n, Y_n) = (x, y)\bigr) = \frac{x}{x+y},$$6 denoting the total number. One has

$$P\bigl((X_{n+1}, Y_{n+1}) = (x+1, y)\mid (X_n, Y_n) = (x, y)\bigr) = \frac{x}{x+y},$$7

where $$P\bigl((X_{n+1}, Y_{n+1}) = (x+1, y)\mid (X_n, Y_n) = (x, y)\bigr) = \frac{x}{x+y},$$8 indicates that $$P\bigl((X_{n+1}, Y_{n+1}) = (x+1, y)\mid (X_n, Y_n) = (x, y)\bigr) = \frac{x}{x+y},$$9 and $$P\bigl((X_{n+1}, Y_{n+1}) = (x, y+1)\mid (X_n, Y_n) = (x, y)\bigr) = \frac{y}{x+y},$$0 never meet.

Using the joint law of limiting Beta variables (see Section 4), it is shown that

$$P\bigl((X_{n+1}, Y_{n+1}) = (x, y+1)\mid (X_n, Y_n) = (x, y)\bigr) = \frac{y}{x+y},$$1

Large deviations for negatively associated $$P\bigl((X_{n+1}, Y_{n+1}) = (x, y+1)\mid (X_n, Y_n) = (x, y)\bigr) = \frac{y}{x+y},$$2-$$P\bigl((X_{n+1}, Y_{n+1}) = (x, y+1)\mid (X_n, Y_n) = (x, y)\bigr) = \frac{y}{x+y},$$3 variables (by Dubhashi–Ranjan 1998) imply the exponentially small probability of significant deviation from the mean.

Consequently, by the strong law of large numbers,

$$P\bigl((X_{n+1}, Y_{n+1}) = (x, y+1)\mid (X_n, Y_n) = (x, y)\bigr) = \frac{y}{x+y},$$4

and thus almost surely $$P\bigl((X_{n+1}, Y_{n+1}) = (x, y+1)\mid (X_n, Y_n) = (x, y)\bigr) = \frac{y}{x+y},$$5.

4. Joint Limiting Law for Normalized Coordinates

For a Pólya walk started at $$P\bigl((X_{n+1}, Y_{n+1}) = (x, y+1)\mid (X_n, Y_n) = (x, y)\bigr) = \frac{y}{x+y},$$6, the classical martingale gives

$$P\bigl((X_{n+1}, Y_{n+1}) = (x, y+1)\mid (X_n, Y_n) = (x, y)\bigr) = \frac{y}{x+y},$$7

For points related by $$P\bigl((X_{n+1}, Y_{n+1}) = (x, y+1)\mid (X_n, Y_n) = (x, y)\bigr) = \frac{y}{x+y},$$8 (i.e., $$P\bigl((X_{n+1}, Y_{n+1}) = (x, y+1)\mid (X_n, Y_n) = (x, y)\bigr) = \frac{y}{x+y},$$9): $x$0

The joint density of ordered limits $x$1 admits a determinantal formula. Let $x$2 and $x$3 be the marginal Beta density and CDF for parameter $x$4: $x$5 The derivation applies the Karlin–McGregor determinant for noncrossing walks to the unnormalized coordinate process $x$6, divides by $x$7, and lets $x$8, incorporating coalescence via inclusion-exclusion on neighbors.

For two neighbors $x$9 at level $y$0,

$y$1

This exact formula recovers $y$2.

5. Edge-Scaling Limit and the Yule Web

At fixed level $y$3 and interval $y$4, consider

$y$5

the number of up-steps from $y$6. $y$7 is a sum of independent Bernoulli random variables with success probability $y$8. The sum converges in distribution to a Poisson: $y$9

The logarithmic change of variables $\lambda = (k, \ell)$0 renders the rate $\lambda = (k, \ell)$1 constant over the continuum, suggesting a scaling window to real-valued levels. In this scaling limit, replace each discrete level $\lambda = (k, \ell)$2 by $\lambda = (k, \ell)$3 and place a Poisson process of rate $\lambda = (k, \ell)$4 along $\lambda = (k, \ell)$5. The Yule walk $\lambda = (k, \ell)$6 at $\lambda = (k, \ell)$7 proceeds to level $\lambda = (k, \ell)$8 at the next Poisson point to the right; coalescence is inherited from the discrete Pólya Web.

In the Yule Web, the analogue of the Beta limit is given by the martingale

$\lambda = (k, \ell)$9

with

$\omega_\lambda \in \{(1,0), (0,1)\}$0

The Yule Web maintains negative association for “no-meet” indicators as in the discrete case. The number of clusters among $\omega_\lambda \in \{(1,0), (0,1)\}$1 edge-starts in a window of width $\omega_\lambda \in \{(1,0), (0,1)\}$2 satisfies

$\omega_\lambda \in \{(1,0), (0,1)\}$3

Joint distributions of the limiting Gamma variables are again given by a $\omega_\lambda \in \{(1,0), (0,1)\}$4 determinant, analogous to the Beta-case, but constructed from the Gamma densities $\omega_\lambda \in \{(1,0), (0,1)\}$5 and cumulative functions $\omega_\lambda \in \{(1,0), (0,1)\}$6.

6. Summary of Key Properties and Connections

1. The Pólya walk is a directed, inhomogeneous random walk on $\omega_\lambda \in \{(1,0), (0,1)\}$7, with local bias determined by site coordinates as $\omega_\lambda \in \{(1,0), (0,1)\}$8.
2. The coalescing system forms an up-right web where “no-meet” indicators for disjoint pairs of neighbors are negatively associated, owing to geometric separation in the dual web structure and the BKR–Reimer inequality.
3. The total number of clusters for walks born at a fixed level $\omega_\lambda \in \{(1,0), (0,1)\}$9 satisfies a strong law: $$P(\omega_\lambda = (1,0)) = \frac{k}{k+\ell}, \qquad P(\omega_\lambda = (0,1)) = \frac{\ell}{k+\ell}.$$0 almost surely, and $$P(\omega_\lambda = (1,0)) = \frac{k}{k+\ell}, \qquad P(\omega_\lambda = (0,1)) = \frac{\ell}{k+\ell}.$$1.
4. The joint law of normalized limiting coordinates $$P(\omega_\lambda = (1,0)) = \frac{k}{k+\ell}, \qquad P(\omega_\lambda = (0,1)) = \frac{\ell}{k+\ell}.$$2 for finite systems is given by a determinantal formula using the Karlin–McGregor method and inclusion-exclusion to account for coalescence.
5. Scaling to the edge by considering Poisson process approximations leads to the Yule Web. In this regime, normalized coordinate martingales converge almost surely to independent Gamma$$P(\omega_\lambda = (1,0)) = \frac{k}{k+\ell}, \qquad P(\omega_\lambda = (0,1)) = \frac{\ell}{k+\ell}.$$3 distributions, and joint laws are furnished by parallel determinantal constructions.

Further details, derivations, and all critical formulas are provided systematically in (Urbán, 17 Jan 2026).