---
title: Feller–Wentzell Interface Diffusions
url: https://www.emergentmind.com/topics/feller-wentzell-interface-diffusions
type: topic
---

# Feller–Wentzell Interface Diffusions

Feller–Wentzell interface diffusions generalize classical diffusion processes by incorporating nonlocal interface (or boundary) conditions that couple local diffusion dynamics with probabilistic transitions—such as partial reflection, jump-like exits, absorption, and transmission—at specified interfaces. These processes arise as semigroups associated with generators whose domains are specified not only by local regularity but by Feller–Wentzell-type conditions at moving or fixed interfaces. The theory unifies classical and nonlocal (e.g., fractional) models and offers a rigorous framework for modeling transport and random motion across heterogeneous media, including phenomena observed in interfaces with complex reflection, absorption, or transmission behaviors.

## 1. Model Setup and Interface Domains

The foundational setting for Feller–Wentzell interface diffusions is a Markov process on a state space $\mathbb{R}$ partitioned at a (potentially moving) interface, or "membrane," $x = h(s)$ for $0 \leq s \leq T$. The space is decomposed as follows, with time-dependent domains:
- $D_1(s) = (-\infty, h(s))$ and $D_2(s) = (h(s), +\infty)$,
- $D(s) = D_1(s) \cup D_2(s)$.

On each half-line $D_i(s)$, the process follows a uniformly parabolic diffusion with time- and space-dependent generator:
\[
L_i u(s, x) = b_i(s, x) \partial_x^2 u(s, x) + a_i(s, x) \partial_x u(s, x), \qquad i=1,2,
\]
where $b_i(s, x) > 0$, $a_i$, and $b_i$ are bounded and Hölder continuous of order $\alpha \in (0, 1)$ in the domain. The interface curve $h(s)$ is continuous, Hölder of order $>\frac{1}{2}$, ensuring the nondegeneracy of subdomains in the space–time plane [1912.12433].

In superdiffusive kinetic scaling, the bulk equation is a fractional heat equation $(−Δ)^{\alpha/2}$ with $\alpha = 1 + 1/\beta$ (for scattering rates $R(k) \sim |k|^\beta$) on each side of a static interface, e.g., $x = 0$ [1905.10586].

## 2. Infinitesimal Generators, Domain, and Interface (Feller–Wentzell) Conditions

The infinitesimal generator $A$ of these interface diffusions acts as $L_1$ or $L_2$ away from the interface. The critical point is the domain $D(A)$: functions $f$ in $C(\mathbb{R}) \cap C^2(D_1(s)) \cap C^2(D_2(s))$ satisfying matched value and a generalized Feller–Wentzell conjugation condition at $x = h(s)$:

\[
\begin{cases}
f(h(s)-0) = f(h(s)+0) & \text{(continuity)} \\
q_1(s)\, \partial_x f(h(s)-0) - q_2(s)\, \partial_x f(h(s)+0) + \int_{D(s)} [f(h(s)) - f(y)]\, p(s,dy) = 0 & \text{(nonlocal interface condition)}
\end{cases}
\]

Here, $q_1(s), q_2(s) \geq 0$, $q_1(s) + q_2(s) > 0$, and $p(s, \cdot)$ is a Borel measure satisfying integrability and regularity constraints. The nonlocal term models jump-like exits from the interface into the domain with a distribution $p(s, dy)$. The alias notation,
\[
\alpha(s) \,\partial_x u(h(s)-0) + \beta(s)\, \partial_x u(h(s)+0) + \int K_s(dy)\, u(y)  = 0,
\]
absorbs the $u(h(s), \cdot)$-term into the measure $K_s$, where $\alpha(s) = q_1(s)$, $\beta(s) = -q_2(s)$ [1912.12433].

For nonlocal (e.g., fractional) equations, the generator at $x \neq 0$ is
\[
L f(x) = c \int_{y \in \mathbb{R}\setminus\{0\}} [f(x + y) - f(x)] q_\beta(y) dy,
\]
with domain determined by Feller–Wentzell interface conditions enforcing “jump-flux” balance at $x = 0$ [1905.10586].

## 3. Construction of the Semigroup and Transition Kernel

Given smooth initial data, the nonlocal parabolic problem
\[
\partial_s u + L_i u = 0,\quad u(t,x,t) = \varphi(x),
\]
coupled with Feller–Wentzell interface conditions at $x = h(s)$, admits a unique classical solution $u(s, x, t)$ for $0 \leq s < t \leq T$. The associated two-parameter semigroup $\{T_{s,t}\}$ is defined by $(T_{s,t} \varphi)(x) = u(s, x, t)$ and forms a strongly continuous, positive contraction family [1912.12433].

The transition probability kernel $P(s, x, t, dy)$ admits representation via Poisson (bulk) potentials and parabolic single-layer potentials centered at the moving membrane. The latter are determined by a $2 \times 2$ Volterra integral equation system encoding the interface's nonlocal behavior [1912.12433].

In the fractional stable setting,
\[
W(t, x) = \mathbb{E}[W_0(\mathcal{X}(t; x)) \, 1_{t < \tau_{\text{kill}}}],
\]
where $\mathcal{X}(t; x)$ is a Lévy-type process with probabilistic reflection, transmission, or absorption at the interface, gives both a Duhamel (series) and a Feynman–Kac-type representation of solutions [1905.10586].

## 4. Interface and Boundary Condition Taxonomy

The general Feller–Wentzell interface condition encompasses:

- **Partial Reflection:** Parameters $q_1, q_2 > 0$, $p(\cdot) = 0$, leading to matched fluxes on either side—standard reflection.
- **Absorption:** One of $q_1$ or $q_2$ is zero, the corresponding flux vanishes, modeling absorption at the interface.
- **Jump-like Exit:** $q_1 = q_2 = 0$, $p(s, \cdot) > 0$, the process jumps on contact, distributed according to $p$.
- **Combination (General Case):** All terms present; the process reflects, absorbs, or jumps according to specified weights.
- **Nonlocal (Fractional) Transmission and Reflection:** Three exit channels—reflection, transmission, absorption—via weights $R$, $T$, $A$, in nonlocal fractional models [1905.10586], paralleling the “jump-flux” balance.

These interface conditions interpolate between classical boundary conditions (Neumann, Dirichlet, Robin, Wentzell) and nonlocal stochastic behaviors.

| Case                   | $q_1$, $q_2$     | $p(\cdot)$ | Interpretation              |
|------------------------|------------------|------------|----------------------------|
| Reflecting             | $q_1 = q_2 > 0$  | $= 0$      | Common flux (reflection)   |
| Absorbing              | $q_1 = 0$ or $q_2 = 0$ | $= 0$ | Dirichlet at interface     |
| Pure jump-exit         | $= 0$            | $> 0$      | Jump distributed by $p$    |
| Partial reflection + jump | $>0$          | $>0$       | General Feller–Wentzell    |

## 5. Green's Functions, Resolvents, and Spectral Structure

The resolvent operator,
\[
R_\lambda \varphi(x) = \int_s^T e^{-\lambda(t-s)} T_{s, t} \varphi(x) dt,
\]
satisfies $(\lambda - L_i) R_\lambda \varphi(x) = \varphi(x)$ in the bulk, together with interface conditions $B_1(R_\lambda \varphi) = B_2(R_\lambda \varphi) = 0$ at the membrane. The Green's function $G_\lambda(x, y)$ is constructed via the free-space resolvent and parabolic single- and double-layer potentials at the interface; see Section V.2 in Friedman (1964) for analogous constructions [1912.12433].

Under stationary or time-homogeneous conditions (e.g., fixed membrane), spectral expansions in eigenmodes subject to Wentzell boundary conditions yield explicit representations for the transition semigroup. In non-stationary (moving interface) cases, such spectral decompositions are only available via “freezing” $s$ [1912.12433].

## 6. Existence and Uniqueness of Solutions

Under the hypotheses of uniform parabolicity, regularity (Hölder continuity) of data, and positivity for coefficients ($q_i$, $p(s, \cdot)$), existence and uniqueness of bounded classical solutions to the nonlocal interface problem are established via:

- Construction of a Poisson potential for initial data and a simple-layer potential (unknown density $V_i$).
- Use of interface conditions to derive a coupled $2 \times 2$ Volterra system for $V_i$.
- Transformation and solution of the Volterra equations by successive approximations and kernel estimates to ensure convergence and boundedness.
- Reconstruction of the solution and verification of compatibility with the PDE and interface conditions.
- Uniqueness follows by reduction to a homogeneous Volterra system and application of energy estimates or Grönwall-type inequalities [1912.12433]; [1905.10586].

In probabilistic (fractional) frameworks, uniqueness in the finite-energy class is ensured by decay of the $L^2$ energy and absence of a boundary layer due to the interface condition [1905.10586].

## 7. Connections to Microscopic Models and Macroscopic Limits

In kinetic models with reflection, transmission, and absorption at a microscopic interface, the macroscopic limit under heavy-tailed jump kernels leads to fractional diffusion equations with Feller–Wentzell interface conditions. The reflection, transmission, and absorption weights arise precisely as the low-frequency limits of the microscopic probabilities, and the fractional index $\alpha$ encodes the degree of superdiffusivity determined by microscopic scattering [1905.10586].

This demonstrates that Feller–Wentzell interface diffusions serve as universal scaling limits for a broad class of physical and stochastic systems with interfacial mechanisms—unifying standard and anomalous diffusion, local and nonlocal transitions, and connecting mesoscopic (kinetic) rules directly to macroscopic interface laws. The analytical machinery of boundary integrals, Volterra equations, and careful semigroup construction underpins this theory, situating Feller–Wentzell interface diffusions as a central model class for rigorous analysis of processes with nontrivial interface or boundary phenomena [1912.12433]; [1905.10586].

Source: https://www.emergentmind.com/topics/feller-wentzell-interface-diffusions