[Paper Review] From the long jump random walk to the fractional Laplacian

Enrico Valdinoci|ArXiv.org|Jan 21, 2009

Advanced Mathematical Modeling in Engineering37 references202 citations

TL;DR

This paper establishes a rigorous connection between long jump Lévy flights—discrete random walks with heavy-tailed jumps—and the fractional Laplacian via a continuous limit. By analyzing the scaling limit of such walks with power-law jump kernels $|y|^{-(n+\alpha)}$, it derives the fractional Laplacian as a singular integral operator through Fourier multipliers, showing that the generator of the process converges to $-(-\Delta)^{\alpha/2}u$ in the limit, with $\alpha \in (0,2)$.

ABSTRACT

This note illustrates how a simple random walk with possibly long jumps is related to fractional powers of the Laplace operator. The exposition is elementary and self-contained.

Motivation & Objective

• To establish a probabilistic and analytical link between long jump random walks and nonlocal operators, particularly the fractional Laplacian.
• To demonstrate how singular integral operators naturally emerge as the continuous limit of discrete long jump processes.
• To provide a self-contained, elementary derivation of the fractional Laplacian as the generator of a Lévy process with infinite variance.
• To clarify the relationship between the Fourier multiplier of the operator and the jump kernel via the formula $S(\xi) = \int (\cos(\xi \cdot y) - 1) K(y) \, dy$.
• To show that the fractional Laplacian $(-\Delta)^{\alpha/2}$ can be equivalently defined via both Fourier multipliers and singular integral representations.

Proposed method

• Define a discrete long jump random walk on the lattice $h\mathbb{Z}^n$ with jump probabilities governed by a symmetric, homogeneous kernel $K(k) = |k|^{-(n+\alpha)}$.
• Derive the discrete evolution equation $u(x,t+\tau) - u(x,t) = \sum_{k \in \mathbb{Z}^n} K(k) \left[ u(x + hk, t) - u(x,t) \right]$ with $\tau = h^\alpha$.
• Take the continuous limit as $h \to 0^+$, transforming the sum into a Riemann sum that converges to the singular integral $\partial_t u(x,t) = \int_{\mathbb{R}^n} \frac{u(x+y,t) - u(x,t)}{|y|^{n+\alpha}} \, dy$.
• Use Fourier analysis to relate the jump kernel $K(y)$ to the Fourier multiplier $S(\xi)$ via $S(\xi) = \int_{\mathbb{R}^n} (\cos(\xi \cdot y) - 1) K(y) \, dy$, establishing the symbol of the operator.
• Prove the equivalence between the Fourier definition $(-\Delta)^{\alpha/2}u = \mathcal{F}^{-1}( |\xi|^\alpha \mathcal{F}u )$ and the singular integral representation $(-\Delta)^{\alpha/2}u = -\int_{\mathbb{R}^n} \frac{u(x+y) + u(x-y) - 2u(x)}{|y|^{n+\alpha}} \, dy$.
• Show that the normalization-free identity $|\xi|^\alpha = \int_{\mathbb{R}^n} \frac{1 - \cos(\xi \cdot y)}{|y|^{n+\alpha}} \, dy$ holds up to constants, using rotational invariance and scaling.

Experimental results

Research questions

• RQ1 How does a long jump random walk with heavy-tailed jumps converge to a nonlocal diffusion process governed by the fractional Laplacian?
• RQ2 What is the precise relationship between the jump kernel $K(y) = |y|^{-(n+\alpha)}$ and the Fourier multiplier of the resulting operator?
• RQ3 Why is the singular integral $\int \frac{u(x+y) - u(x)}{|y|^{n+\alpha}} \, dy$ well-defined as a principal value for $\alpha \in (0,2)$ and smooth $u$?
• RQ4 How can the fractional Laplacian be rigorously derived as the generator of a Lévy process with infinite variance via discrete approximations?
• RQ5 What is the mathematical justification for equating the Fourier multiplier $|\xi|^\alpha$ with the integral $\int \frac{1 - \cos(\xi \cdot y)}{|y|^{n+\alpha}} \, dy$?

Key findings

• The continuous limit of a long jump random walk with jump kernel $K(y) = |y|^{-(n+\alpha)}$ yields the evolution equation $\partial_t u = \int_{\mathbb{R}^n} \frac{u(x+y,t) - u(x,t)}{|y|^{n+\alpha}} \, dy$, which is a nonlocal diffusion equation.
• The singular integral in the limit is well-defined as a principal value when $\alpha \in (0,2)$, due to cancellation of the linear term in the Taylor expansion.
• The fractional Laplacian $(-\Delta)^{\alpha/2}$ is equivalent to the singular integral $-\int_{\mathbb{R}^n} \frac{u(x+y) + u(x-y) - 2u(x)}{|y|^{n+\alpha}} \, dy$ up to normalization constants.
• The Fourier multiplier of the operator is $S(\xi) = |\xi|^\alpha$, which matches the symbol of $(-\Delta)^{\alpha/2}$, confirming the equivalence of the Fourier and integral definitions.
• The identity $|\xi|^\alpha = \int_{\mathbb{R}^n} \frac{1 - \cos(\xi \cdot y)}{|y|^{n+\alpha}} \, dy$ holds up to normalization, proven via rotational invariance and scaling.
• The process has infinite variance since the $\beta$-moment $\sum_k |k|^\beta K(k)$ diverges for $\beta \geq \alpha$, confirming it is in the domain of attraction of an $\alpha$-stable law.