---
title: Higher-Order Fractional Laplacians
url: https://www.emergentmind.com/topics/higher-order-fractional-laplacians
type: topic
---

# Higher-Order Fractional Laplacians

Higher-order fractional Laplacians generalize both the classical Laplacian and its integer iterates to non-integer powers, providing a unified framework for nonlocal operators of arbitrary real order. They are fundamental in the analysis of partial differential equations (PDEs) involving nonlocality and play crucial roles in conformal geometry, potential theory, elliptic regularity, and various applications in physics, such as anomalous diffusion and elasticity. Their study encompasses operator theory, variational methods, explicit kernel representations, maximum principle failures, boundary regularity, and computational discretization schemes.

## 1. Definitions and Operator Formulations

Several equivalent definitions of the higher-order fractional Laplacian $(-\Delta)^s$ for $s>0$ are standard:

- **Spectral (Fourier) Definition:** For $u\in\mathcal{S}(\mathbb{R}^n)$,
  $$
  (-\Delta)^s u(x) = \mathcal{F}^{-1}\left(|\xi|^{2s} \widehat{u}(\xi)\right)(x),
  $$
  where $\mathcal{F}$ denotes the Fourier transform. This definition immediately extends to tempered distributions and to Sobolev/Bessel potential spaces by duality [1702.02078].

- **Integral (Hypersingular) Formulation (Riesz, pointwise):** For non-integer $s>0$, write $s=m+\sigma$ with $m\in\mathbb{N}_0$, $\sigma\in(0,1)$. For smooth $u$
  $$
  (-\Delta)^s u(x) = C_{n,\sigma} \,\mathrm{P.V.} \int_{\mathbb{R}^n} \frac{\Delta^m u(x) - \Delta^m u(y)}{|x-y|^{n+2\sigma}} \,dy,
  $$
  with normalization so that the Fourier symbol is $|\xi|^{2s}$ [1607.00929, 1810.08435, 1702.02078].

- **Finite-Difference ("Grünwald-Letnikov"-type) Formula:** Alternative representation using a finite number of function evaluations, typically as
  $$
  (-\Delta)^s u(x) = C_{n,s} \int_{\mathbb{R}^n} \sum_{k=0}^{m} (-1)^k \binom{m}{k} u\big(x + (\tfrac{m}{2}-k)y\big) |y|^{-n-2s} dy,
  $$
  which recovers classical polyharmonic operators for integer $s$ [1810.08435, 2205.12610].

- **Spectral Fractional Laplacian in Domains:** For a bounded domain $\Omega$ with Dirichlet Laplacian eigenpairs $\{(\lambda_k, \phi_k)\}$,
  $$
  (-\Delta)^s u = \sum_{k=1}^\infty \lambda_k^s U_k \phi_k, \qquad U_k = (u, \phi_k)_{L^2(\Omega)},
  $$
  mapping the Sobolev domain $H^s_\Delta(\Omega)$ to its dual [2603.06945].

## 2. Fundamental Solutions, Integral Representations, and Extension Problems

- **Riesz Potential and Green Kernels:** For $n>2s$, the fundamental solution is $F_{n,s}(x) = K_{n,s}|x|^{2s-n}$, satisfying $(-\Delta)^s F_{n,s} = \delta_0$ distributionally [1607.00929]. For the Dirichlet problem in the unit ball $B$, explicit representations via Boggio's formula are available, giving $G_s(x,y)$ with precise boundary decay and regularity [1810.08435, 1707.03603].

- **Integral Equation Equivalence:** For PDEs of the form $(-\Delta)^s u = f(u)$ on $\mathbb{R}^n$, solutions often coincide with solutions of the nonlinear Riesz potential equation
  $$
  u(x) = C_{n,s} \int_{\mathbb{R}^n} \frac{f(u(y))}{|x-y|^{n-2s}} dy,
  $$
  under suitable decay or integrability [1609.04105, 2109.09441, 1905.04300].

- **Extension Problem (Polyharmonic and Weighted):** The Caffarelli-Silvestre extension for $s\in(0,1)$ and its generalizations for $s>1$ realize $(-\Delta)^s$ as a Dirichlet-to-Neumann map for a $k$-th order degenerate elliptic equation in $\mathbb{R}^{n+1}_+$, with $y^b$-weighted Laplacian: $\mathcal{A}_b^k U = 0$ and boundary data prescribed at $y=0$, where $b=1-2\{s\}$, $k=1+[s]$ [2106.11669, 1302.4413]. This framework provides variational characterizations and explicit extension kernels.

## 3. Boundary Value Problems and Regularity

- **Dirichlet and Nonlocal Boundary Conditions:** For fractional orders, the correct notion of Dirichlet problem is to prescribe $u$ on $\mathbb{R}^n\setminus\Omega$. For $s\in(m,m+1)$, corresponding traces ("Edenhofer kernels") correspond to consecutive boundary operators recovering standard normal derivatives as $s\to\mathbb{N}$ [1707.03603, 1810.08435].

- **Boundary Regularity:** Solutions $u$ to
  $$
  (-\Delta)^s u = f \text{ in }\Omega, \quad u=0\text{ in }\mathbb{R}^n\setminus\Omega
  $$
  satisfy $u/d^s \in C^{s-\epsilon}(\bar{\Omega})$ for all $\epsilon>0$ ($d(x)$ distance to $\partial\Omega$), with $u/d^s\in C^\infty(\bar{\Omega})$ if $f\in C^\infty$ [1406.1107, 1607.00929].

- **Positivity and Maximum Principle:** The maximum principle holds for $s\in(0,1)$, but fails for $s>1$ in general, as sign-changing Green functions and eigenfunctions emerge in nontrivial domains (especially disconnected or nonconvex), with explicit characterizations via kernel sign patterns [1607.00929, 2205.12610, 2005.09514].

## 4. Liouville Theorems, Radial Symmetry, and Classification

- **Liouville-type Theorems:** For subcritical powers $p$ in semilinear equations with $(-\Delta)^s u = u^p$ or Hardy-Hénon weights, all nonnegative entire solutions must vanish, with the critical exponent given by $p^*_{s,\alpha} = \frac{n+2s-2\alpha}{n-2s}$ for $(-\Delta)^s u = |x|^{-\alpha} u^p$ [2109.09441, 1609.04105, 1905.04300].

- **Symmetry and Moving Spheres/Planes:** In subcritical or critical cases, all nonnegative solutions are radially symmetric and monotone (up to conformal transformation). These are established both via the method of moving spheres (integral forms) and via splitting systems for the operator factorization [1609.04105, 2109.09441, 1905.04300].

- **Super Polyharmonic Properties:** Nonnegative solutions $u$ of higher-order equations satisfy $(−\Delta)^{i+\alpha/2}u \geq 0$ for all $i=0,\ldots,m-1$ for $s=m+\alpha/2$ [1905.04300].

## 5. Oscillatory Phenomena, Maximum Principle Failure, and Rearrangement Inequalities

- **Oscillatory Behavior and Eigenfunction Sign Changes:** For $s\in(k,k+1),$ the kernel and first Dirichlet eigenfunction may change sign, especially for disconnected or large aspect-ratio domains. Explicit domains (unions of balls, ellipsoids) realize this anti-maximum behavior [2205.12610, 2005.09514, 1607.00929].

- **Breakdown of Rearrangement Inequalities:** For $s>1$, polarization and Pólya–Szegő type inequalities (reduction of Dirichlet energy by symmetric decreasing rearrangement) can fail, and may even reverse direction in certain ranges of $s$ [2205.12610].

- **Exceptions in One Dimension:** Despite these oscillations, the Faber–Krahn inequality holds for $s>1$ in dimension one: the interval always minimizes the fundamental gap among open sets of fixed length [2205.12610].

## 6. Variational Properties, Integration by Parts, and Pohozaev Identities

- **Fractional Energy Forms:** The natural energy is $u\mapsto \mathcal{E}_s(u,u) = \int_{\mathbb{R}^n} |\xi|^{2s} |\widehat u(\xi)|^2 d\xi$. Minimizing this over $u$ vanishing outside $\Omega$ yields the unique weak solution to $(-\Delta)^s u = f$ [1607.00929, 1810.08435].

- **Integration by Parts and Boundary Terms:** For $u,v\in H^s(\mathbb{R}^n),u=v=0$ on $\mathbb{R}^n\setminus\Omega$, and smooth $(-\Delta)^s u, (-\Delta)^s v$,
  $$
  \int_\Omega (-\Delta)^s u\, \partial_i v\, dx = -\int_\Omega \partial_i u\,(-\Delta)^s v\, dx - \Gamma(1+s)^2 \int_{\partial \Omega} \frac{u}{d^s}\frac{v}{d^s}\nu_i\, d\sigma,
  $$
  with corresponding Pohozaev-type identities for fractional orders [1406.1107].

- **Unique Continuation and Runge Approximation:** Solutions exhibit strong unique continuation from open sets and support Runge approximation results, providing the basis for inverse problem identifiability [2001.06210].

## 7. Numerical Methods, Discretization, and Implementation

- **Finite Difference & High-Order Schemes:** Accurate discretizations for $(-\Delta)^s$ employ either high-order central difference stencils (with composite quadrature or lattice-Boltzmann–Hermite expansions for isotropy), or Fourier-based weights, yielding $O(h^p)$ convergence for $p$-th order stencils [2601.21388, 2202.06333]. The discretized operators are Toeplitz and amenable to FFT acceleration.

- **PDE-Based Polyharmonic Extensions:** Extension-based approaches directly discretize the degenerate polyharmonic problem for the extension variable, coupling high-regularity finite elements in $x$ and spline approximations in $y$ to achieve algebraic convergence in both variables [2603.06945, 2106.11669].

- **Kernel Evaluation on General Domains:** Explicit formulas for $(-\Delta)^s$ acting on polynomials or radial profiles are available for balls and ellipsoids, but general domains require weighted $L^1$ spaces for integrability, and explicit kernel assembly becomes intractable outside balls [1707.03603, 2005.09514].

## 8. Sharp Inequalities and Extremal Sequences

- **Adams and Moser–Trudinger Inequalities:** For the critical Sobolev exponent $p=n/a$, the sharp exponential integrability bound holds for $u\in W^{a,p}(\mathbb{R}^n)$ with $\|u\|_{L^p} + \|(-\Delta)^{a/2}u\|_{L^p} \leq 1$,
  $$
  \int_{E} \exp\Bigl(B_1|u(x)|^{p'}\Bigr)dx \leq C(1+|E|),\qquad B_1 = n/\omega_{n-1},\quad p' = n/(n-a),
  $$
  sharpness demonstrated by explicit extremal ("bubble") sequences [1702.02078].

- **Extremal Solutions and Criticality:** While true maximizers do not exist, carefully constructed squeezing sequences achieve blow-up at the critical constant in Adams-type inequalities [1702.02078].

---

**References**

- [1702.02078] Fontana & Morpurgo, "Sharp exponential integrability for critical Riesz potentials and fractional Laplacians on ℝ^n"
- [2109.09441] Chen, Li & Wang, "Liouville-type theorems, radial symmetry and integral representation of solutions to Hardy-Hénon equations involving higher order fractional Laplacians"
- [1808.03241] Parisis & Aifantis, "Fractional Generalization of Higher-Order Diffusion"
- [1707.03603] Abatangelo & Jarohs, "Integral representation of solutions to higher-order fractional Dirichlet problems on balls"
- [2005.09514] Gal & Jarohs, "Fractional Laplacians on ellipsoids"
- [2001.06210] Ghosh, Salo, Uhlmann, "Unique continuation property and Poincaré inequality for higher order fractional Laplacians"
- [1609.04105] Zhuo & Li, "A Liouville Theorem for the Higher Order Fractional Laplacian"
- [1607.00929] Abatangelo, Jarohs & Saldaña, "On the maximum principle for higher-order fractional Laplacians"
- [1905.04300] Cao, Dai, Qin, "Super poly-harmonic properties, Liouville theorems and classification of nonnegative solutions to equations involving higher-order fractional Laplacians"
- [1810.08435] Saldaña, Abatangelo & Jarohs, "On fractional higher-order Dirichlet boundary value problems: between the Laplacian and the bilaplacian"
- [1302.4413] Stinga & Torrea, "On higher order extensions for the fractional Laplacian"
- [1406.1107] Ros-Oton & Serra, "Local integration by parts and Pohozaev identities for higher order fractional Laplacians"
- [2205.12610] Abatangelo & Jarohs, "Oscillatory Phenomena for Higher-Order Fractional Laplacians"
- [2106.11669] Cora & Musina, "The $s$-polyharmonic extension problem and higher-order fractional Laplacians"
- [2601.21388] Wang et al., "Higher-Order Finite Difference Methods for the Tempered Fractional Laplacian"
- [2202.06333] Ghaffari, "On Generalisation of Isotropic Central Difference for Higher Order Approximation of Fractional Laplacian"
- [2603.06945] Otárola & Salgado, "Approximation of higher-order powers of the spectral fractional Laplacian via polyharmonic extension"

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