---
title: Regional Fractional Laplacian
url: https://www.emergentmind.com/topics/regional-fractional-laplacian
type: topic
---

# Regional Fractional Laplacian

The regional fractional Laplacian is a bounded-domain nonlocal operator defined by restricting the singular integral of the fractional Laplacian to interior interactions. For a bounded open or Lipschitz domain $\Omega\subset\mathbb{R}^n$ and $0<s<1$, a standard form is
\[
(-\Delta)^s_{\Omega}u(x)=C_{n,s}\,\mathrm{P.V.}\!\int_{\Omega}\frac{u(x)-u(y)}{|x-y|^{n+2s}}\,dy,\qquad x\in\Omega,
\]
with the normalization chosen so that on $\mathbb{R}^n$ one recovers the Fourier symbol $|\xi|^{2s}$ [2403.00357]. In probabilistic terms it is the generator of the censored $2s$-stable process, and analytically it differs from the whole-space, restricted, and spectral fractional Laplacians by the way the boundary enters both the operator and the admissible data [1711.06916].

## 1. Definition and functional framework

The operator is naturally paired with the Gagliardo seminorm
\[
[u]_{H^s(\Omega)}^2=\iint_{\Omega\times\Omega}\frac{|u(x)-u(y)|^2}{|x-y|^{n+2s}}\,dx\,dy,
\]
and with the Sobolev-type spaces
\[
H^s(\Omega)=\{u\in L^2(\Omega):[u]_{H^s(\Omega)}<\infty\},\qquad
H_0^s(\Omega)=\overline{C_c^\infty(\Omega)}^{\,H^s},
\]
which support the bilinear form associated with $(-\Delta)^s_\Omega$ [1509.05838]. In variational form, the operator is characterized by
\[
Q_{n,s,\Omega}(u,v)=\frac{C_{n,s}}{2}\int_\Omega\int_\Omega
\frac{(u(x)-u(y))(v(x)-v(y))}{|x-y|^{n+2s}}\,dx\,dy,
\]
and weak solutions are defined by duality against $H_0^s(\Omega)$ or $H^s(\Omega)$ test functions, depending on the boundary regime [2112.06272].

Because the integral is restricted to $\Omega$, boundary geometry enters through the domain itself rather than through an exterior tail over $\mathbb{R}^n\setminus\Omega$. One paper formulates this by saying that the operator only “sees” the interior $\Omega$, so boundary geometry appears via an effective nonlocal “Hardy” term near $\partial\Omega$ [2403.00357]. In the linear elliptic setting, coercivity follows from a Hardy–Sobolev inequality and a Poincaré inequality on $H_0^\alpha(\Omega)$ for $\alpha\in(\tfrac12,1)$ [1509.05838].

A recurrent point in the literature is that the boundary formulation is regime-dependent. Some treatments regard the appropriate Dirichlet condition as the local one $u=0$ on $\partial\Omega$ because the singular integral involves only $y\in\Omega$ [1711.06916]. Variational treatments instead encode Dirichlet data by requiring $u\in H_0^s(\Omega)$; for $s\le \tfrac12$ in bounded Lipschitz domains, one has $H_0^s(\Omega)=H^s(\Omega)$, so the Dirichlet and Neumann variational formulations coincide at the level of the energy space [2007.04808].

## 2. Position among bounded-domain fractional operators

On bounded domains, the regional fractional Laplacian is one of several inequivalent fractional Laplace operators. Their distinction is structural rather than notational.

| Operator | Integral domain / prescription | Distinctive feature |
|---|---|---|
| Regional fractional Laplacian | Integral over $\Omega$ only | Generator of the censored $2s$-stable process |
| Restricted fractional Laplacian | Integral over $\mathbb{R}^n$ with $u\equiv0$ on $\mathbb{R}^n\setminus\Omega$ | Exterior volume constraint |
| Spectral fractional Laplacian | Functional calculus of the Dirichlet Laplacian | Uses Laplace eigenpairs on $\Omega$ |

For the spectral operator, if $(\mu_k,\varphi_k)$ is the $k$th Dirichlet eigenpair of $-\Delta$, then $\lambda_k^{(\mathrm{spec})}=\mu_k^{\alpha/2}$ and $\phi_k^{(\mathrm{spec})}=\varphi_k$ [1711.06916]. By contrast, the regional operator has its own spectrum and its own boundary regularity class. In one dimension on $\Omega=(-1,1)$, numerical experiments showed
\[
\lambda_k^{(\mathrm{regional})}
<
\lambda_k^{(\mathrm{restricted})}
<
\lambda_k^{(\mathrm{spectral})},
\qquad k=1,2,\dots,\ \alpha\in(0,2),
\]
and all three families converge to the classical Dirichlet Laplacian eigenvalues as $\alpha\to2$ [1711.06916].

The discrepancy is especially pronounced in the strongly nonlocal regime. For $\alpha\ll2$, the regional operator neglects jumps from $x\in\Omega$ to $y\notin\Omega$, while the restricted fractional Laplacian still “sees” all of $\mathbb{R}^n$; in the one-dimensional formula recorded for $\Omega=(-l,l)$, the difference between the two operators contains an explicit boundary-singular factor $\bigl[(l+x)^{-\alpha}+(l-x)^{-\alpha}\bigr]u(x)$ [1711.06916]. This suggests that replacing the restricted operator by the regional one is benign only near the local limit $\alpha\to2$, not uniformly over the full nonlocal range.

A further variant appears when homogeneous nonlocal Neumann conditions are imposed on the whole-space fractional Laplacian. Under the condition $\widetilde{\mathcal N}_s u=0$ on $\mathbb{R}^n\setminus\Omega$, the whole-space operator can be rewritten on $\Omega$ as a purely regional operator with a modified kernel $K_s(x,y)$ that encodes the exterior jumps; the correction term has logarithmic blow-up near the boundary, so the resulting kernel is not stable-like in the sense of kernels bounded above and below on $\Omega\times\Omega$ [1712.00320].

## 3. Boundary regimes and regularity theory

The boundary behavior of the regional operator depends sharply on whether $s$ is above or below $\tfrac12$. For $\alpha\in(0,\tfrac12]$, the censored symmetric $2\alpha$-stable process is conservative and never approaches the boundary, and the corresponding PDE theory exhibits strong nonexistence phenomena [2007.05775]. In particular, for bounded $C^2$ domains, the equation
\[
(-\Delta)^\alpha_\Omega u=1\quad\text{in }\Omega
\]
admits no viscosity solution bounded from above or bounded from below, and the Lane–Emden problem
\[
(-\Delta)^\alpha_\Omega u=u^p\quad\text{in }\Omega,\qquad u=0\quad\text{on }\partial\Omega
\]
admits no nonnegative nontrivial solution for any $p>0$ [2007.05775].

For $s\in(\tfrac12,1)$, the boundary rate $\delta^{\,2s-1}$ becomes fundamental. A Hopf boundary lemma shows that if $u$ is a pointwise or weak super-solution of
\[
(-\Delta)^s_\Omega u=c(x)u\quad\text{in }\Omega,
\]
then either $u\equiv0$ or
\[
\liminf_{x\to z,\ x\in\Omega}\frac{u(x)}{\mathrm{dist}(x,\partial\Omega)^{2s-1}}>0
\]
at every boundary point $z\in\partial\Omega$; the same paper identifies the torsion function as having two-sided bounds of order $\mathrm{dist}(x,\partial\Omega)^{2s-1}$, so the exponent $2s-1$ is sharp [2112.09522].

Boundary regularity has been developed in both Hölder and Sobolev scales. For bounded $C^{1,\sigma}$ domains and $f\in L^p(\Omega)$, zero-Neumann solutions satisfy $w\in C^{2s-N/p}(\overline\Omega)$, and if $2s-N/p>1$ then $w\in C^{1,2s-N/p-1}(\overline\Omega)$; for the Dirichlet problem with $s\in(\tfrac12,1)$ and $p>N$, one obtains
\[
\frac{w}{\delta^{2s-1}}\in C^{1-N/p}(\overline\Omega)
\]
up to the boundary [2007.04808]. The proofs combine half-space Liouville classifications, blow-up arguments, compactness, and kernel-freezing.

A more refined Sobolev-scale boundary theory is available in flat geometry. On the upper half-space $\mathbb{R}^n_+$, if $\tfrac12<s<1$, $p\in(2,\infty)$, and $u$ solves
\[
L_s^{\mathbb{R}^n_+}u=g,
\]
then for any $\tilde s$ with $s<\tilde s<\min(1,2s-t)$ one has a global tangential Calderón–Zygmund estimate
\[
[u]_{W_T^{\tilde s,p}(\mathbb{R}^n_+)}\le C\|G\|_{L^p(\mathbb{R}^n)},
\]
and the $n$-dimensional boundary regularity problem can be reduced to a one-dimensional regional fractional Laplace equation on $(0,\infty)$ [2209.07939].

## 4. Spectral, asymptotic, and variational structure

The regional fractional Laplacian has a discrete Dirichlet spectrum on bounded Lipschitz domains, with eigenvalues characterized by min–max principles in $H_0^s(\Omega)$ [2112.08856]. As $s\to0^+$, the operator admits an expansion
\[
(-\Delta)^s_\Omega u=s\,L_0u+o(s),
\]
where
\[
L_0u(x)=|S^{n-1}|\,\mathrm{P.V.}\!\int_\Omega\frac{u(x)-u(y)}{|x-y|^n}\,dy
\]
is the regional logarithmic Laplacian [2112.08856]. In the same limit, for each fixed $k$, $\lambda_k(s)\to\lambda_k(0)$ and $u_k^s\to u_k^0$ in $L^2(\Omega)$ and even in $C^0(\Omega)$ [2112.08856].

The dependence on the order parameter is differentiable at the level of weak solutions. For the free Poisson problem with zero-average datum,
\[
(-\Delta)^s_\Omega u_s=f,\qquad \int_\Omega f=0,
\]
the map $s\mapsto u_s$ belongs to $C^1((0,1);L^2(\Omega))$ [2112.09547]. For the first nontrivial eigenvalue $\lambda_{1,s}$ on the zero-mean space, right differentiability is established, and the derivative is expressed through a logarithmic correction to the energy form [2112.09547].

Critical variational problems reveal a strong domain dependence absent in the full-space theory. For
\[
S_{n,s}(\Omega)=\inf_{u\in H_0^s(\Omega)\setminus\{0\}}
\frac{[u]_{H^s(\Omega)}^2}{\|u\|_{L^{2_s^*}(\Omega)}^2},
\qquad 2_s^*=\frac{2n}{n-2s},
\]
the best constant genuinely depends on $\Omega$ [2403.00357]. If $\mathbb{R}^n\setminus\Omega$ has nonempty interior, then $S_{n,s}(\Omega)<S_{n,s}(\mathbb{R}^n)$; for $s>\tfrac12$ and suitable smooth domains with a flat boundary patch, the infimum is attained, and a curvature expansion shows how $S_{n,s}(\Omega)$ departs from its full-space counterpart [2403.00357].

Existence of minimizers has also been obtained through concentration-compactness and geometric test-function arguments. For bounded $C^1$ domains there exists $s_0\in(\tfrac12,1)$ such that $S_{n,s}(\Omega)$ is achieved for every $s\in(\tfrac12,s_0)$; in the unit ball, positive radial minimizers exist for all $s\in(\tfrac12,1)$ when $n\ge4s$, while the low-dimensional case $n<4s$ is handled under a positive mass condition [2112.06272]. At the shape-optimization level, a Rayleigh–Faber–Krahn-type problem for the regional energy admits a compactly supported nonnegative minimizer, but symmetrization techniques may fail, and the regularity and geometry of the positivity set remain open [2108.07967].

## 5. Boundary-value problems, nonlinear equations, and identities

For the linear Dirichlet problem on bounded $C^2$ domains with $\alpha\in(\tfrac12,1)$,
\[
(-\Delta)^\alpha_\Omega u=f\quad\text{in }\Omega,\qquad u=g\quad\text{on }\partial\Omega,
\]
existence and uniqueness are available at three levels of source data when $g\equiv0$: $f\in L^2(\Omega)$, $f\in L^1(\Omega,\rho^\beta dx)$, and $\nu\in\mathcal M(\Omega,\rho^\beta)$, where $\rho(x)=\mathrm{dist}(x,\partial\Omega)$ and $\beta=2\alpha-1$ [1509.05838]. The same work provides an integration by parts formula involving the nonlocal normal derivative
\[
D_\beta g(x)=\lim_{t\to0^+}\frac{g(x)-g(x-t\,\mathbf n(x))}{t^\beta}
\]
for classical solutions with general boundary data [1509.05838].

Nonlinear equations display boundary phenomena specific to the regional operator. For boundary blow-up problems
\[
(-\Delta)^a_\Omega u+f(u)=0\quad\text{in }\Omega,\qquad \lim_{p(x)\to0}u(x)=+\infty,
\]
with $p(x)=\mathrm{dist}(x,\partial\Omega)$, the pure power nonlinearity $f(s)=s^p$ has the sharp threshold $p_c=1+2a$: no blow-up solution exists if $p\le p_c$, while a minimal blow-up solution exists if $p>p_c$ [1602.02869]. In the existence regime, two-sided boundary estimates of power type are obtained, and the proof uses Perron’s method, Green-operator estimates, and barriers of the form $V_\tau(x)=p(x)^\tau$ [1602.02869].

The low-order regime $s\in(0,\tfrac12]$ is usually hostile to boundary-vanishing positive solutions, but special existence results are available on the unit ball for equations of the form
\[
(-\Delta)^s_{B_1}u+u=h(x,u),\qquad u\in C_0(B_1),
\]
when the nonlinearity is either $h(x,t)=f(x)$ or $h(x,t)=h_1(x)t^p+\epsilon h_2(x)$ with radially symmetric decreasing data under suitable conditions [2211.10554]. The resulting solutions are positive, radial, and strictly decreasing [2211.10554].

For nonlinear nonlocal operators modeled on the regional fractional $p$-Laplacian in the subquadratic regime $1<p<2$, weak solutions enjoy local self-improvement. If
\[
L_{s,\Omega}u=f
\]
with measurable symmetric kernel comparable to $|x-y|^{-(n+sp)}$ and $f\in (W_0^{s-(p-1),p}(\Omega))^*$, then there exists $\varepsilon_0>0$ such that $u\in W^{s+\varepsilon,p}(\Omega_1)$ for every $\Omega_1\Subset\Omega$, together with a fractional Sobolev gain in integrability [2606.26610]. The proof uses nonlinear commutator estimates, localization, and a Gehring-type iteration [2606.26610].

A recent development is a Pohozaev-like identity for the regional operator on bounded $C^2$ domains. For $s\in(\tfrac12,1)$, weak solutions of
\[
(-\Delta)^s_\Omega u=f(u)\quad\text{in }\Omega,\qquad u=0\quad\text{on }\partial\Omega
\]
satisfy an identity with an explicit remainder term and a boundary contribution involving $(u/\delta^{2s-1})^2(x\cdot\nu)$ [2507.21962]. The same work establishes a new integration-by-parts formula for $(-\Delta)^s_\Omega$ and discusses applications to eigenvalue problems in the unit ball and to boundary-type unique continuation [2507.21962].

## 6. Stochastic, numerical, and data-driven realizations

The stochastic interpretation is central. The regional fractional Laplacian is the generator of the censored $\alpha$-stable process, and several papers also describe it as arising from jumps that are allowed only inside the domain [1711.06916]. Under the nonlocal Neumann condition of Dipierro–Ros-Oton–Valdinoci, the induced regional kernel corresponds instead to a “reflected-and-redistributed” process: when a particle lands outside $\Omega$, it immediately jumps back into $\Omega$ with density proportional to $|x-y|^{-n-2s}$; the logarithmic boundary singularity of the kernel prevents the process from lying in the stable-like class [1712.00320].

The operator also appears as a scaling limit in interacting particle systems. On $(0,1)$, a boundary-driven exclusion process with long jumps and asymmetric rates converges, depending on the asymmetry strength, either to stationary solutions of the Ornstein–Uhlenbeck equation or to energy solutions of the stochastic Burgers equation with characteristic operator given by the regional fractional Laplacian $L^{\gamma/2}$ [2412.10068]. In that framework, the energy form
\[
\ll f,g\gg_{\gamma/2}
\]
makes $L^{\gamma/2}$ symmetric, and the martingale problem is formulated on test-function spaces whose boundary conditions depend on $\gamma$ and on the reservoir regime [2412.10068].

Numerically, the regional operator is not a universally accurate surrogate for the restricted fractional Laplacian. Extensive experiments comparing the fractional, spectral, regional, and peridynamic operators show that all collapse to the classical Laplacian as $\alpha\to2$, but for $\alpha\ll2$ the regional fractional Laplacian generally gives results inconsistent with the fractional Laplacian, whereas the peridynamic model approximates the fractional Laplacian at rate $O(\delta^{-\alpha})$ when the horizon $\delta$ is sufficiently large [1711.06916].

In data analysis, nonlocal diffusion maps with polynomial-tail kernels provide a discrete approximation of the regional fractional Laplacian on manifolds with boundary. The construction uses graph distances to approximate geodesic distances, a two-stage normalization to remove sampling bias, and a graph generator
\[
L_\epsilon^\beta u(x_i)=\frac{1}{t}\Bigl[u(x_i)-\sum_jP_{ij}u(x_j)\Bigr],
\]
which converges in the large-data, small-bandwidth limit to $(-\Delta)^s_\Omega$ with $s=\beta/2$ [1810.03952]. Because the point cloud is confined to $\Omega$, the limiting integral is automatically regional rather than whole-space [1810.03952].

Taken together, these results place the regional fractional Laplacian at the intersection of censored jump processes, boundary-sensitive nonlocal PDE, domain-dependent critical phenomena, and computational schemes that exploit the operator’s strictly in-domain character. The present body of work also isolates several persistent open problems: precise Green’s function estimates, sharper boundary regularity in low-order regimes, Weyl-type asymptotics beyond currently known cases, the geometry of shape minimizers, and the analytic status of the logarithmically corrected kernels arising from nonlocal Neumann reformulations [1711.06916].

Source: https://www.emergentmind.com/topics/regional-fractional-laplacian