advanced25 min studyLesson 390 in Mathematics

Fractional Sobolev spaces and interpolation

Measure nonlocal differences at intermediate smoothness and connect the Gagliardo seminorm with Fourier and interpolation viewpoints.

On this page

What you’ll understand

  • Compute a fractional difference seminorm with explicit s and p.
  • Distinguish whole-space Fourier Hˢ from general-p fractional conventions.
  • State real interpolation and zero-extension claims with their domain and endpoint restrictions.

Intermediate smoothness cannot be expressed by taking half of a classical derivative pointwise. Fractional Sobolev spaces instead measure how function values differ across pairs of points, with a distance weight that emphasizes short scales. For square integrability, Fourier weights give a complementary frequency interpretation. Interpolation provides a third viewpoint linking function norms and derivative norms. These descriptions agree under precise hypotheses, while boundary choices and general integrability exponents can lead to genuinely different spaces.

A nonlocal difference energy

For an open set Ω⊂Rn\Omega\subset\mathbb R^n, 0<s<10<s<1, and 1≤p<∞1\le p<\infty, define the Gagliardo seminorm by integrating ∣u(x)−u(y)∣p/∣x−y∣n+sp|u(x)-u(y)|^p/|x-y|^{n+sp} over all ordered pairs in the domain. The fractional space Wˢ,ᵖ consists of Lᵖ functions with finite seminorm, with full norm formed by adding the pth power of the Lᵖ norm. Constants have zero seminorm, so the order-zero term is again necessary. The singular kernel does not make every smooth function fail: a Lipschitz difference contributes a compensating power of the distance when s is below one.

The energy is nonlocal because points interact across arbitrary separations. On a bounded domain, the intrinsic seminorm only counts pairs inside the domain. Extending by zero and computing a whole-space seminorm also counts pairs crossing the boundary, producing an additional boundary-layer cost. These are not interchangeable definitions. Even a constant function has zero intrinsic seminorm but its zero extension can have a singular interaction near the boundary. State whether a PDE uses intrinsic differences, an exterior condition, or an interpolation space before comparing its energy to an integer-order norm.

[u]Ws,p(Ω)p=∫Ω∫Ω∣u(x)−u(y)∣p∣x−y∣n+sp dx dy,∥u∥Ws,pp=∥u∥pp+[u]Ws,pp[u]_{W^{s,p}(\Omega)}^p=\int_\Omega\int_\Omega\frac{|u(x)-u(y)|^p}{|x-y|^{n+sp}}\,dx\,dy,\qquad \|u\|_{W^{s,p}}^p=\|u\|_p^p+[u]_{W^{s,p}}^p
The ordered-pair convention counts both orientations; this matters in exact seminorm calculations.↔ Scroll wide equations horizontally.

Fourier weights and interpolation

For p equal to two on the whole space and a unitary Fourier transform, the squared Gagliardo seminorm is a positive dimensional constant times ∫∣ξ∣2s∣u^(ξ)∣2 dξ\int|\xi|^{2s}|\widehat u(\xi)|^2\,d\xi. The full norm is equivalent to the Bessel-potential Hˢ norm with weight (1+∣ξ∣2)s(1+|\xi|^2)^s. It is equivalent rather than numerically identical with the convention above unless the normalization constant is included. For general p, Slobodeckij spaces and Bessel-potential spaces are not automatically the same; using Hˢ language for the Hilbert case avoids silently identifying different scales.

On a bounded Lipschitz domain and for 1<p<∞1<p<\infty, real interpolation gives Ws,p(Ω)=(Lp(Ω),W1,p(Ω))s,pW^{s,p}(\Omega)=(L^p(\Omega),W^{1,p}(\Omega))_{s,p} with equivalent norms. The K-functional measures the cheapest decomposition of a function into a part small in Lᵖ and a part controlled in W¹,ᵖ, with a parameter weighting the latter. Integrating these costs across scales creates fractional regularity. Replacing W¹,ᵖ by its zero-trace subspace changes the boundary behavior of the interpolation space; at critical boundary exponents extra conditions can appear. It is unsafe to identify all spaces with a subscript zero just because their names look similar.

K(t,u;X0,X1)=inf⁡u=v+w(∥v∥X0+t∥w∥X1),∥u∥(X0,X1)s,p=(∫0∞[t−sK(t,u)]pdtt)1/pK(t,u;X_0,X_1)=\inf_{u=v+w}\bigl(\|v\|_{X_0}+t\|w\|_{X_1}\bigr),\qquad \|u\|_{(X_0,X_1)_{s,p}}=\left(\int_0^\infty[t^{-s}K(t,u)]^p\frac{dt}{t}\right)^{1/p}
Real interpolation weighs decomposition costs at every scale; the stated Sobolev identification assumes a bounded Lipschitz domain and finite p strictly above one.↔ Scroll wide equations horizontally.

Jumps, embeddings, and boundary thresholds

A one-dimensional jump can belong to Wˢ,ᵖ when sp is below one. Near the interface, points on opposite sides give a kernel integral with the same convergence threshold as ∫0εr−sp dr\int_0^\varepsilon r^{-sp}\,dr. At sp equal to one it diverges logarithmically, and above one it also diverges. Thus fractional regularity can admit jumps excluded from first-order Sobolev spaces, but only below the appropriate exponent threshold. This is a local interface calculation, not a claim that every discontinuous function has finite fractional energy.

On a bounded Lipschitz domain, for 0<s<10<s<1 and 1<p<∞1<p<\infty, fractional embedding has conjugate exponent np/(n−sp)np/(n-sp) when sp<nsp<n, while sp>nsp>n gives a Hölder representative with exponent s−n/ps-n/p. Compactness is again subcritical under these bounded-domain hypotheses. Smooth-domain trace and zero-extension thresholds are linked to sp crossing one, but the equality case needs special care. For example, the constant function on an interval lies in its intrinsic H¹ᐟ² space while its whole-line zero extension has infinite H¹ᐟ² seminorm. Boundary conventions are part of the mathematical problem and cannot be inferred from an interior formula alone.

An exact half-order energy

Use the displayed ordered-pair convention on (0,1)(0,1) with s = one-half and p = two. Compute the fractional seminorm and full squared norm of u(x)=xu(x)=x.

  1. The numerator is ∣x−y∣2|x-y|^2, and the denominator has exponent n+sp=1+1=2n+sp=1+1=2. Away from the diagonal their quotient is exactly one. The diagonal has zero two-dimensional measure and does not change the integral.

  2. The double integral over the unit square is therefore one, so [u]W1/2,22=1[u]_{W^{1/2,2}}^2=1. This calculation includes both orientations of each pair as specified in the definition.

  3. The squared L² norm is ∫01x2=1/3\int_0^1x^2=1/3. Adding it gives ∥u∥W1/2,22=4/3\|u\|_{W^{1/2,2}}^2=4/3. This is the intrinsic interval norm, not the whole-space energy of an unspecified continuation.

The result

A simple affine function has a finite explicitly computed half-order energy, with the pair-counting and boundary conventions both visible.

Explore the model

Interactive lab

Interactive lab · spectral

Fractional energy reweights the same modes

Pause and Reset; set mode count N to 4 and coefficient decay p to 1. At time zero the periodic sine coefficients are 1/k1/k. For a fresh comparison Reset, set N to 4 and p to 2. On paper compare the homogeneous H-half squared energies π∑k=14k/k2p\pi\sum_{k=1}^4 k/k^{2p}; predict that faster decay reduces the weighted high-mode contribution.

Fractional energy reweights the same modesA finite sum of sin(nx − t)/nᵖ. Playback changes t; increasing the cutoff includes finer spatial scales. This is not a proof about a continuum limit or renormalization. Current readings: Sum at x = π/2: 0.867; Smallest wavelength: 2π / 5; Highest-mode coefficient: 0.2. All parameters are adjustable with the labeled controls below.01.573.144.716.28−2.56−1.2801.282.56x · radiansFinite sine sum
Finite sumHighest included mode
Sum at x = π/2
0.867
Smallest wavelength
2π / 5
Highest-mode coefficient
0.2
t = 0 s

A finite sum of sin(nx − t)/nᵖ. Playback changes t; increasing the cutoff includes finer spatial scales. This is not a proof about a continuum limit or renormalization. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forThe finite sine instrument lets you inspect the modes whose fractional Fourier weights are computed on paper. It does not evaluate a double-integral seminorm or establish infinite-series membership; its parameter p is coefficient decay, distinct from the integrability exponent.

Check your understanding

For the intrinsic Gagliardo seminorm of a one-dimensional step with one interior jump, which threshold is correct?

Try it yourself

Let v=1(1/2,1)v=\mathbf1_{(1/2,1)} on (0,1)(0,1), with p = two and s = one-quarter. Compute its squared Gagliardo seminorm and explain what happens if s increases to one-half.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Only opposite-side pairs contribute. Write 2∫01/2∫1/21(y−x)−3/2 dy dx2\int_0^{1/2}\int_{1/2}^1(y-x)^{-3/2}\,dy\,dx and integrate the inner variable first.

Compare with the solution

The inner integral is 2[(1/2−x)−1/2−(1−x)−1/2]2[(1/2-x)^{-1/2}-(1-x)^{-1/2}]. Integrating over x and multiplying by the outer factor two gives 8(2−1)8(\sqrt2-1). At s one-half, the kernel exponent becomes two; the corresponding corner integral has a logarithmic divergence. Thus v belongs to the intrinsic W¹ᐟ⁴,² interval space but not to W¹ᐟ²,², in agreement with the threshold sp < 1.

Take it with you

  • Fractional Sobolev energy measures weighted differences over pairs of points.
  • Fourier and interpolation characterizations require explicit exponent, normalization, and domain conventions.
  • Jumps and zero-extension boundary costs change at sp = 1, where endpoint care is essential.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.