---
title: Fractional Discrete Laplacians
url: https://www.emergentmind.com/topics/fractional-discrete-laplacians
type: topic
---

# Fractional Discrete Laplacians

A fractional discrete Laplacian is a nonlocal operator extending the concept of fractional Laplacians from continuous domains to discrete settings such as lattices, periodic chains, or general graphs. These operators form the cornerstone for diverse research fields including analysis on graphs, nonlocal PDEs, harmonic analysis, fractional quantum mechanics, and discrete modeling of anomalous diffusion. Fractional discrete Laplacians interpolate between the identity and local (integer-order) Laplacians and capture the signature nonlocal, scale-free behavior characteristic of their continuous counterparts.

## 1. Algebraic and Spectral Definitions

On the uniform lattice $\mathbb{Z}_h = h\mathbb{Z}$, the standard discrete Laplacian, $\Delta_h u_j = (u_{j+1} - 2u_j + u_{j-1})/h^2$, admits fractional powers defined spectrally via the semigroup:
\[
(-\Delta_h)^s u_j = \frac{1}{\Gamma(-s)} \int_0^\infty (e^{t \Delta_h} u_j - u_j)\, t^{-1-s}\, dt, \quad 0 < s < 1.
\]
On finite or infinite graphs, the operator is analogously defined by spectral calculus in terms of the Laplacian's (generalized) eigenpairs; for a finite graph $(V,E)$ with Laplacian $-\Delta$, the spectral definition is
\[
(-\Delta)^s f = \sum_{k=0}^{n-1} \lambda_k^s \langle f, \phi_k\rangle \phi_k,
\]
where $\{\phi_k\}$ and $\{\lambda_k\}$ are the eigenfunctions and eigenvalues of $-\Delta$ [2403.19987]. The heat semigroup formula generalizes to weighted graphs and locally finite graphs:
\[
(-\Delta)^s u(x) = \frac{1}{\Gamma(-s)} \int_0^\infty (e^{t\Delta}u(x) - u(x)) t^{-1-s} dt,
\]
where $e^{t\Delta}$ is the graph heat semigroup [2408.02902].

For translation-invariant cases, the Fourier representation yields the symbol $[2\sin(h\xi/2)/h]^{2s}$ [1507.04986]. In $\mathbb{Z}^d$ with anisotropic orders $\vec s$, the operator is most naturally viewed as
\[
\Delta_{\mathbb{Z}^d}^{\vec s} f(n) = \sum_{j=1}^d \Delta_j^{s_j} f(n),
\]
where each $\Delta_j^{s_j}$ is defined by the spectral calculus for the shift-invariant 1D Laplacian acting along $j$ [2510.10680, 2509.06117].

## 2. Explicit Kernel Formulas and Nonlocal Series Representations

Fractional discrete Laplacians admit pointwise, global (nonlocal) series representations. In one dimension,
\[
(-\Delta_h)^s u_j = \sum_{m\in\mathbb{Z}} (u_j-u_m) K_s(j-m),
\]
where the kernel $K_s(m)$ has explicit expressions:
\[
K_s(m) = \frac{4^s\Gamma(1/2+s)\Gamma(|m|-s)}{\sqrt{\pi} |\Gamma(-s)| \Gamma(|m|+1+s)} \, h^{-2s}.
\]
This kernel decays as $|m|^{-1-2s}$, so the operator is genuinely nonlocal but with algebraically decaying tails [1608.08913, 1507.04986, 2101.03629].

On general graphs or weighted graphs, the kernel $K_s(x,y)$ is built from the heat kernel $p_t(x,y)$ via
\[
K_s(x,y) = \frac{1}{|\Gamma(-s)|} \int_0^\infty p_t(x,y) t^{-1-s} dt,
\]
and the operator writes as a discrete “Riesz-type” integral:
\[
(-\Delta)^s u(x) = \sum_{y\in V} (u(x)-u(y)) K_s(x,y).
\]
Such symmetric, positive, rapidly decaying kernels underpin the nonlocal character of discrete fractional Laplacians on arbitrary (possibly infinite) graphs [2408.02902].

For the periodic chain of length $L=Nh$, the operator is a matrix function of the Born–von Karman Laplacian, with explicit entries using generalized binomial coefficients:
\[
f_N^{(\alpha)}(m) = \Omega_\alpha^2\sum_{n=-\infty}^\infty f^{(\alpha)}_\infty(|m+nN|), \quad
f^{(\alpha)}_\infty(k) = \Omega_\alpha^2 (-1)^k \binom{\alpha}{\frac{\alpha}{2}+k} [1412.5904].
\]

## 3. Analytical Properties and Limiting Cases

Fractional discrete Laplacians interpolate continuously between the identity and the standard Laplacian:
- $\lim_{s\to 1^-} (-\Delta)^s = -\Delta$
- $\lim_{s\to 0^+} (-\Delta)^s = \mathrm{Id}$ on zero-mean functions, for finite graphs [2403.19987]
- For even-integer $\alpha=2m$ in chains, the kernel recovers the $2m$-point central difference formula [1412.5904, 1501.01942].
- As $h\to 0$ (continuum limit), $(-\Delta_h)^s$ converges to the continuous space Riesz fractional Laplacian, both in the pointwise and spectral sense [1608.08913, 1412.5904].

Spectral decompositions and kernel representations reveal that discrete fractional Laplacians are self-adjoint on the $\ell^2$ space, have absolutely continuous spectrum in the infinite lattice setting, and inherit the nonlocality and scaling properties of their continuous counterparts.

## 4. Boundary Conditions, Compression, and Operator Extension

For bounded domains or compressed settings, forms of the fractional discrete Laplacian are defined by restriction/projection:
- The half-lattice discrete fractional Laplacian is given by
\[
\Delta_{\mathbb{Z}_+^d}^{\vec s}= P_+ \Delta_{\mathbb{Z}^d}^{\vec s} P_+ + K_{\vec s}
\]
where $P_+$ is the projection to the half-lattice and $K_{\vec s}$ is a compact "boundary correction," ensuring that the essential spectrum matches that of the infinite lattice operator and that interior spectral/dynamical results transfer between full- and half-lattice settings [2510.10680].

Generalizations to finite graphs employ the spectral functional calculus; for graphs with boundary conditions, Dirichlet, Neumann, and Robin extensions have been systematically discretized via finite element and quadrature-based methods [1708.03602]. For negative fractional powers (fractional integrals), similar constructions yield discrete versions of Riesz potentials, preserving mapping properties and Hardy–Littlewood–Sobolev inequalities [1608.08913, 2603.01039].

## 5. Analytic, Probabilistic, and Variational Structures

Discrete fractional Laplacians on graphs generate regular Dirichlet forms:
\[
Q^s(f) = \frac{1}{2} \sum_{x,y} K_s(x,y) (f(x) - f(y))^2
\]
which serve as discrete energy functionals and underlie variational approaches to nonlocal equations. The fractional discrete Laplacian is the Dirichlet-to-Neumann map for a semidiscrete extension problem, mirroring the Caffarelli–Silvestre extension in the continuum [1608.08913].

s-Harmonic functions satisfy a nonlocal mean-value property, which has a (fractional) random walk interpretation: the probability to "jump" is proportional to $|j-m|^{-1-2s}$ for large $|j-m|$ [1608.08913]. On graphs, fractional Sobolev spaces $W^{s,2}(V)$ are defined, yielding embeddings and variational principles for related nonlinear equations, such as the discrete fractional Schrödinger equation [2408.02902].

For general weighted graphs, the structure of positive-critical Hardy weights for fractional Laplacians has been characterized, allowing for optimal Hardy inequalities and sharp control of ground-state behaviors, with explicit asymptotics on Cayley, curvature, and fractal graphs [2605.19423].

## 6. Spectral Theory, Scattering, and Propagation Estimates

On hypercubic lattices, the anisotropic fractional discrete Laplacian’s spectral theory, Mourre estimates, and threshold analysis are developed in detail:
- The spectrum is determined by the sum of powers of the discrete Laplacians along each axis, with thresholds at critical values of the group velocity.
- A strict Mourre estimate is established away from thresholds, yielding a Limiting Absorption Principle (LAP), absence of singular continuous spectrum, and finiteness of eigenvalues on spectral windows [2509.06117, 2510.10680].
- Scattering theory for long-range perturbations involves stationary representations and Birman–Krein formulae for the scattering matrix, ensuring completeness and spectral-shift function representation [2509.06117].

Propagation and time-decay estimates (weighted and minimal velocity bounds) have been derived for fractional discrete Laplacians, matching sharp results from the continuum theory.

## 7. Numerical Discretization, Algorithms, and Applications

Multiple high-accuracy and structure-preserving discretizations for fractional Laplacians and their nonlinear variants have been constructed:
- Finite-difference schemes (spectral, regularized, tent, quadratic) yield discretizations with explicit positive weights, controlled accuracy, and convergent to the correct continuum limit [1611.00164, 1311.7691].
- Factorization methods and Toeplitz structures enable efficient multidimensional matrix-vector operations via FFT [2103.03658, 1802.03770]. Monte Carlo and fast preconditioned solvers are effective for large-scale elliptic and parabolic problems.
- Discretizations on general bounded domains are resolved using heat-semigroup quadrature and finite element temporal-spatial discretization, with proven convergence for spectral and integral fractional operators under Dirichlet, Neumann, or Robin boundary conditions [1708.03602].

Operator and kernel regularity, mean-value formulas, and tail corrections are central to precise numerical treatment in practice. The stability, maximal principle, and convergence are well-understood for the schemes constructed using semigroup or quadrature-based frameworks [2401.09926].

---

**References**:  
- [1412.5904], [1501.01942], [1507.04986], [1608.08913], [1611.00164], [1708.03602], [1802.03770], [2101.03629], [2103.03658], [2401.09926], [2403.19987], [2408.02902], [2509.06117], [2510.10680], [2603.01039], [2605.19423].

Source: https://www.emergentmind.com/topics/fractional-discrete-laplacians