---
title: Caffarelli–Silvestre Extension Method
url: https://www.emergentmind.com/topics/caffarelli-silvestre-extension
type: topic
---

# Caffarelli–Silvestre Extension Method

The Caffarelli–Silvestre extension is a degenerate-elliptic technique that realizes nonlocal fractional powers of second-order differential operators, such as the fractional Laplacian $(-\Delta)^s$ for $0 < s < 1$, as Dirichlet-to-Neumann maps for local PDEs in one higher dimension. This extension method underpins the analytic, variational, probabilistic, and numerical frameworks for fractional elliptic operators and is foundational in the analysis of nonlocal PDEs, inverse problems, and stochastic processes.

## 1. Classical Formulation and Weighted Extension Problem

Given $u \in H^s(\mathbb{R}^n)$, the Caffarelli–Silvestre extension $U(x,y)$ is defined as the unique solution to the weighted degenerate elliptic equation
\[
\begin{cases}
\operatorname{div}(y^{1-2s}\nabla U) = 0, & (x,y) \in \mathbb{R}^n \times (0,\infty), \\
U(x,0) = u(x), & x \in \mathbb{R}^n,
\end{cases}
\]
where the variable $y > 0$ is called the extension variable or "bulk" direction. The Dirichlet datum is imposed on $y = 0$, which acts as the "boundary" corresponding to the original domain.

The key property is the recovery of the fractional Laplacian via a weighted Neumann trace:
\[
-(\lim_{y \to 0^+} y^{1-2s} \partial_y U(x,y)) = C_s\, (-\Delta)^s u(x),
\]
in the sense of $H^{-s}$–$H^s$ duality, with $C_s = 2^{1-2s}\Gamma(1-s)/\Gamma(s)$ [2305.04227]. This provides a realization of the nonlocal operator $(-\Delta)^s$ as the Dirichlet-to-Neumann map of a local PDE in $\mathbb{R}^{n+1}_+$.

## 2. Analytic Structure and Weighted Sobolev Spaces

The extension problem is formulated in the weighted Sobolev space $\dot H^1(\mathbb{R}^{n+1}_+, y^{1-2s})$, which consists of functions $U$ with weak derivatives and finite energy
\[
\|U\|_{\dot H^1(\mathbb{R}^{n+1}_+, y^{1-2s})}^2 = \int_{\mathbb{R}^{n+1}_+} y^{1-2s} |\nabla U(x,y)|^2\, dx\, dy < \infty,
\]
and trace $U(\cdot,0) = u$ belonging to $H^s(\mathbb{R}^n)$. The mapping
\[
U \mapsto U|_{y=0}
\]
is continuous and surjective from $\dot H^1(\mathbb{R}^{n+1}_+, y^{1-2s})$ to $H^s(\mathbb{R}^n)$ [2206.04353].

The energy of the extension exactly captures the Gagliardo semi-norm of the fractional Sobolev space:
\[
E[U] = \int_{\mathbb{R}^{n+1}_+} y^{1-2s} |\nabla U|^2\, dx\, dy,
\]
is equivalent, up to a constant, to
\[
\iint_{\mathbb{R}^n \times \mathbb{R}^n} \frac{|u(x) - u(\xi)|^2}{|x - \xi|^{n + 2s}} \, dx\, d\xi.
\]
This yields a variational characterization and enables minimization principles for nonlocal equations.

## 3. Dirichlet-to-Neumann Operators and Generalizations

For variable-coefficient elliptic operators, the extension builds on a degenerate-elliptic operator in higher dimensions. Given a symmetric, positive-definite matrix $a(x)$, the extension problem is
\[
\operatorname{div}(y^{1-2s}\tilde a(x,y)\nabla \tilde U) = 0, \quad \tilde a(x,y) = \begin{pmatrix} a(x) & 0 \\ 0 & 1 \end{pmatrix},
\]
with appropriate boundary conditions. The resulting Dirichlet-to-Neumann map retrieves fractional elliptic operators of the form $(-\operatorname{div}(a \nabla))^s$ [2305.04227, 2309.00858].

The extension can equivalently be formulated for more general self-adjoint operators $L$, where the extension PDE becomes
\[
\partial_{yy} U + \frac{1-2s}{y} \partial_y U = -L U, \qquad U|_{y=0} = f,
\]
and the Dirichlet-to-Neumann map recovers $L^s f$ up to normalization. The abstract framework extends to generators of strongly continuous semigroups and Dirichlet forms, as in [1207.7203, 2403.18984, 2101.11305], and to more general Banach-space settings via functional calculi based on complete Bernstein functions.

## 4. Applications: Inverse Problems, Regularity, and Numerical Methods

A notable application is in inverse problems for nonlocal elliptic operators. In the fractional Calderón problem, the Caffarelli–Silvestre extension provides a constructive pathway to reduce uniqueness proofs for nonlocal inverse problems to their local analogues [2305.04227, 2309.00858]. This connection is explicit in the mapping of fractional Dirichlet-to-Neumann data to the local counterpart via vertical averaging of the extension solution:
\[
v(x') = \int_0^\infty y^{1-2s} \tilde U(x', y) \, dy,
\]
which solves a local elliptic PDE and determines the local Dirichlet-to-Neumann map.

Regularity results, Almgren-type monotonicity formulas, and optimal Hölder estimates for obstacle and free boundary problems are naturally formulated in the extension framework, facilitating the treatment of degenerate elliptic equations and providing access to weighted harmonic measure theory [1403.5015].

In computational mathematics, the extension transforms nonlocal PDEs into local, albeit higher-dimensional, variational problems, enabling the design of efficient numerical methods, including deep Ritz approaches in high dimensions [2108.11592]. The degeneracy and singularity near $y=0$ are addressed by specialized approximation spaces and integration schemes based on the structure of the extension.

## 5. Extensions: Higher-Order, Discrete, and Non-Euclidean Generalizations

Higher-order fractional Laplacians, $(-\Delta)^s$ with $s > 1$, admit an analogous extension, but the related PDE becomes higher-order in $y$ and requires additional boundary conditions, involving both Dirichlet and conormal data [1302.4413, 1809.09496]. The Dirichlet-to-Neumann characterization persists, though with increased complexity in the formulation.

On discrete structures, such as $\mathbb{Z}$ or graphs, a "discrete Caffarelli–Silvestre extension" realizes fractional discrete Laplacians as boundary traces of local operators defined on higher-dimensional graphs (e.g., $\mathbb{Z} \times \mathbb{N}$), enabling the analysis of discrete nonlocal models via local Gaussian free fields on layered structures [2312.04536].

For subelliptic and CR geometries (e.g., the Heisenberg group), the extension operates in the context of non-Euclidean metrics and sub-Laplacians, often linked to scattering theory and the geometry of the Siegel upper half-space. The Dirichlet-to-Neumann map for the extension again recovers conformally covariant fractional powers of subelliptic operators, with identical normalization to the Euclidean case [1312.3381].

## 6. Limitations, Obstructions, and Inverse Correspondence

While the extension provides a surjective reduction from fractional to local operators at the level of Dirichlet-to-Neumann data, the inverse is generally obstructed. In the Calderón-type framework, key obstructions arise from two independent one-dimensional averaging effects: (i) vertical averaging (losing $y$-dependence of the extension on the boundary), and (ii) tangential averaging related to density arguments for traces from "partial" to "full" boundary data. Consequently, the mapping is surjective only up to closure and cannot be continuously inverted in general [2305.04227].

## 7. Further Generalizations and Connections

Recent developments generalize the extension method to higher-order and logarithmic operators (via Hadamard regularization), discrete settings, general Dirichlet forms, and nonlocal problems involving drift terms or singular potentials [2502.10747, 1403.5015, 2403.18984]. The functional calculus based on extension operators characterizes a broad class of spectral multipliers for which the extension paradigm applies, including Phillips–Bochner subordination and complete Bernstein functions [2101.11305].

The extension also underpins a rich theory of trace inequalities, embeddings, and capacity theory in function spaces (e.g., Sobolev, Besov, and Lebesgue spaces) and interfaces with stochastic analysis, via its probabilistic realization as the boundary value of degenerate diffusions with suitable Bessel-type processes in the extension variable [2310.01070, 2201.00765, 2007.00713].

---

**References**:
- "A Reduction of the Fractional Calderón Problem to the Local Calderón Problem by Means of the Caffarelli-Silvestre Extension" [2305.04227]
- "Optimal regularity of solutions to the obstacle problem for the fractional Laplacian with drift" [1403.5015]
- "An extension problem of higher order operators and operators of logarithmic type via renormalization" [2502.10747]
- "Revisiting the Anisotropic Fractional Calderón Problem Using the Caffarelli-Silvestre Extension" [2309.00858]
- "On a family of integral operators on the ball" [2110.11913]
- "Extension problem and fractional operators: semigroups and wave equations" [1207.7203]
- "Functional Calculus via the extension technique: a first hitting time approach" [2101.11305]
- "A note on a stochastic approach to Caffarelli-Silvestre Theorem" [2310.01070]
- "Quantization of nonlocal fractional field theories via the extension problem" [1907.00733]
- "Invisibility of the integers for the discrete Gaussian chain via a Caffarelli-Silvestre extension of the discrete fractional Laplacian" [2312.04536]
- "Singularities of fractional Emden's equations via Caffarelli-Silvestre extension" [2206.04353]
- "Deep Ritz method for the spectral fractional Laplacian equation using the Caffarelli-Silvestre extension" [2108.11592]
- "Embeddings of Function Spaces via the Caffarelli-Silvestre Extension, Capacities and Wolff potentials" [2007.00713]
- "Unique continuation principles for a higher order fractional Laplace equation" [1809.09496]
- "Extension method in Dirichlet spaces with sub-Gaussian estimates and applications to regularity of jump processes on fractals" [2403.18984]
- "Fractional Besov Trace/Extension Type Inequalities via the Caffarelli-Silvestre extension" [2201.00765]
- "On higher order extensions for the fractional Laplacian" [1302.4413]
- "An extension problem for the CR fractional Laplacian" [1312.3381]
- "An extension problem related to inverse fractional operators" [1603.07988]

Source: https://www.emergentmind.com/topics/caffarelli-silvestre-extension