Distributions and Sobolev Spaces

⚡ Key Takeaways

Distributions (generalized functions) are continuous linear functionals on a space of test functions, allowing differentiation of objects far rougher than classical functions. Sobolev spaces $W^{k,p}$ are Banach spaces of functions whose weak (distributional) derivatives up to order $k$ lie in $L^p$ . Together they form the foundational language of modern PDE theory, the calculus of variations, and numerical analysis. Introduced by Sergei Sobolev (1930s) and systematized by Laurent Schwartz (1940s–50s, Fields Medal 1950).

💡 The Core Insight

Stop asking "what is the value of $f$ at a point?" and instead ask "how does $f$ act when integrated against a smooth bump?" Integration by parts then becomes a definition rather than a theorem: the derivative of a distribution $T$ is the distribution that pairs with $\varphi$ the way $T$ pairs with $-\varphi'$ . Because we threw the derivative onto the infinitely-smooth test function, everything is infinitely differentiable in the distributional sense. Sobolev spaces are then the natural "in-between" world: rough enough to contain solutions, smooth enough ( $L^p$ derivatives) to carry a usable Banach/Hilbert structure.

DefinitionTest Functions and the Space

\mathcal{D}(\Omega)

Let $\Omega \subseteq \mathbb{R}^n$ be open. The space of test functions is

1

\quad \mathcal{D}(\Omega) = C_c^\infty(\Omega) = \{\varphi \in C^\infty(\Omega) : \operatorname{supp}\varphi \subset\subset \Omega\},

Derivation hint: $\operatorname{supp}\varphi = \overline{\{x : \varphi(x)\neq 0\}}$ , and $\subset\subset$ means the support is a compact subset of $\Omega$ . equipped with the following notion of convergence: $\varphi_j \to \varphi$ in $\mathcal{D}(\Omega)$ iff (i) there is a single compact $K \subset \Omega$ containing all $\operatorname{supp}\varphi_j$ , and (ii) $\partial^\alpha \varphi_j \to \partial^\alpha \varphi$ uniformly for every multi-index $\alpha$ . This is the locally convex inductive-limit (LF) topology; $\mathcal{D}(\Omega)$ is not metrizable.

DefinitionDistributions

A distribution is a linear functional $T:\mathcal{D}(\Omega)\to\mathbb{C}$ that is continuous: $\varphi_j\to\varphi$ in $\mathcal{D}$ implies $\langle T,\varphi_j\rangle\to\langle T,\varphi\rangle$ . The space of distributions is $\mathcal{D}'(\Omega)$ . Equivalently (sequential continuity = continuity here), $T\in\mathcal{D}'(\Omega)$ iff for each compact $K\subset\Omega$ there exist $C_K>0$ and $N_K\in\mathbb{N}$ with

2

\quad |\langle T,\varphi\rangle| \le C_K \sum_{|\alpha|\le N_K} \sup_{x\in K} |\partial^\alpha\varphi(x)|, \qquad \forall\,\varphi\in C_c^\infty(K).

Derivation hint: this is the seminorm characterization of continuity on the LF space; if a single $N_K = N$ works for all $K$ , the distribution has finite order $N$ . Every locally integrable function $f\in L^1_{loc}(\Omega)$ embeds as the regular distribution $T_f$ via

3

\quad \langle T_f,\varphi\rangle = \int_\Omega f(x)\,\varphi(x)\,dx.

Derivation hint: $T_f$ is well-defined because $\varphi$ has compact support, so the integral is over a compact set; the map $f\mapsto T_f$ is injective by the fundamental lemma of the calculus of variations.

DefinitionDistributional (Weak) Derivative

For $T\in\mathcal{D}'(\Omega)$ and multi-index $\alpha$ , the distributional derivative $\partial^\alpha T\in\mathcal{D}'(\Omega)$ is defined by

4

\quad \langle \partial^\alpha T,\varphi\rangle = (-1)^{|\alpha|}\,\langle T,\partial^\alpha\varphi\rangle, \qquad \forall\,\varphi\in\mathcal{D}(\Omega).

Derivation hint: motivated by repeated integration by parts on $T_f$ with $f\in C^{|\alpha|}$ — the boundary terms vanish since $\varphi$ has compact support, leaving the sign $(-1)^{|\alpha|}$ . A function $g\in L^1_{loc}(\Omega)$ is the weak derivative $\partial^\alpha f$ of $f\in L^1_{loc}(\Omega)$ if $T_g = \partial^\alpha T_f$ , i.e.

5

\quad \int_\Omega g\,\varphi\,dx = (-1)^{|\alpha|}\int_\Omega f\,\partial^\alpha\varphi\,dx \quad\forall\varphi\in\mathcal{D}(\Omega).

Derivation hint: combine [3] and [4]; weak derivatives are unique up to sets of measure zero by injectivity of $f\mapsto T_f$ .

DefinitionSobolev Spaces

W^{k,p}

and

H^k

Let $1\le p\le\infty$ and $k\in\mathbb{N}_0$ . The Sobolev space is

6

\quad W^{k,p}(\Omega) = \{ f\in L^p(\Omega) : \partial^\alpha f\in L^p(\Omega)\ \text{(weakly)}\ \forall\,|\alpha|\le k\},

Derivation hint: $\partial^\alpha f$ here means the weak derivative of [5]; membership requires every such derivative to exist as an $L^p$ function. with norm (for $p<\infty$ )

7

\quad \|f\|_{W^{k,p}(\Omega)} = \Bigg(\sum_{|\alpha|\le k}\|\partial^\alpha f\|_{L^p(\Omega)}^p\Bigg)^{1/p},

Derivation hint: equivalent to $\sum_{|\alpha|\le k}\|\partial^\alpha f\|_{L^p}$ ; the $\ell^p$ -sum of finitely many norms is itself a norm by Minkowski. and for $p=2$ we write $H^k(\Omega) = W^{k,2}(\Omega)$ , a Hilbert space with inner product

8

\quad \langle f,g\rangle_{H^k} = \sum_{|\alpha|\le k}\int_\Omega \partial^\alpha f\,\overline{\partial^\alpha g}\,dx.

Derivation hint: completeness follows because $L^2$ is complete and the weak-derivative operators are closed, so a Cauchy sequence in $H^k$ converges componentwise with limits matching by closedness. The closure of $C_c^\infty(\Omega)$ in the $W^{k,p}$ norm is denoted $W_0^{k,p}(\Omega)$ ( $H_0^k$ for $p=2$ ); it encodes vanishing boundary conditions.

TheoremCompleteness

For every $1\le p\le\infty$ , $k\in\mathbb{N}_0$ , and open $\Omega$ , $\big(W^{k,p}(\Omega),\|\cdot\|_{W^{k,p}}\big)$ is a Banach space; it is reflexive for $1<p<\infty$ , separable for $1\le p<\infty$ , and a Hilbert space for $p=2$ .

TheoremSobolev Embedding Theorem

Let $\Omega\subset\mathbb{R}^n$ be bounded with Lipschitz boundary, $1\le p<\infty$ . Set the Sobolev conjugate $p^* = \tfrac{np}{n-kp}$ when $kp<n$ . Then: 1. Subcritical ( $kp<n$ ): $W^{k,p}(\Omega)\hookrightarrow L^{q}(\Omega)$ continuously for $q\le p^*$ , and compactly for $q<p^*$ (Rellich–Kondrachov). 2. Critical/borderline ( $kp=n$ ): $W^{k,p}(\Omega)\hookrightarrow L^q(\Omega)$ for all $q<\infty$ (but not $L^\infty$ in general). 3. Supercritical ( $kp>n$ ): $W^{k,p}(\Omega)\hookrightarrow C^{m,\gamma}(\overline\Omega)$ where $m = k - \lfloor n/p\rfloor - 1$ and $\gamma$ is the Hölder exponent, with $\gamma = \lfloor n/p\rfloor + 1 - n/p$ if $n/p\notin\mathbb{Z}$ . The Gagliardo–Nirenberg–Sobolev inequality underlies case (1):

9

\quad \|f\|_{L^{p^*}(\mathbb{R}^n)} \le C(n,p)\,\|\nabla f\|_{L^p(\mathbb{R}^n)}, \qquad p^* = \tfrac{np}{n-p},\ \ 1\le p<n.

Derivation hint: for $p=1$ start from $|f(x)|\le \int_{-\infty}^{x_i}|\partial_i f|\,dt$ in each coordinate, multiply the $n$ resulting bounds, integrate, and apply the generalized Hölder inequality; for $p>1$ apply the $p=1$ case to $|f|^\gamma$ with $\gamma=\tfrac{p(n-1)}{n-p}$ .

ProofSobolev's Inequality for

p=1

,

n=2

(model case)

We prove $\|f\|_{L^{2}(\mathbb{R}^2)} \le \tfrac12\|\nabla f\|_{L^1(\mathbb{R}^2)}$ for $f\in C_c^\infty(\mathbb{R}^2)$ (here $p^*=2$ ). 1. Line integration. For each point $(x_1,x_2)$ ,

10

\quad |f(x_1,x_2)| \le \int_{-\infty}^{\infty}|\partial_1 f(t,x_2)|\,dt =: A(x_2), \qquad |f(x_1,x_2)| \le \int_{-\infty}^{\infty}|\partial_2 f(x_1,s)|\,ds =: B(x_1).

Derivation hint: $f$ vanishes at $\pm\infty$ , so $f(x_1,x_2)=\int_{-\infty}^{x_1}\partial_1 f\,dt$ and $|f|\le \int|\partial_1 f|\,dt$ ; same in the other variable. 2. Multiply the two bounds. $|f(x_1,x_2)|^2 \le A(x_2)\,B(x_1)$ . 3. Integrate in $x_1$ then $x_2$ .

11

\quad \int_{\mathbb{R}^2}|f|^2\,dx \le \int_{\mathbb{R}}\!\!\int_{\mathbb{R}} A(x_2)B(x_1)\,dx_1\,dx_2 = \Big(\int_{\mathbb{R}} B(x_1)\,dx_1\Big)\Big(\int_{\mathbb{R}} A(x_2)\,dx_2\Big).

Derivation hint: $A$ depends only on $x_2$ , $B$ only on $x_1$ , so the double integral factors (Fubini–Tonelli). 4. Recognize the gradient. Each factor equals $\int_{\mathbb{R}^2}|\partial_i f|\,dx \le \|\nabla f\|_{L^1}$ . Thus

12

\quad \|f\|_{L^2}^2 \le \|\partial_1 f\|_{L^1}\,\|\partial_2 f\|_{L^1} \le \tfrac14\big(\|\partial_1 f\|_{L^1}+\|\partial_2 f\|_{L^1}\big)^2 \le \tfrac14\|\nabla f\|_{L^1}^2.

Derivation hint: AM–GM gives $ab\le\frac14(a+b)^2$ ; taking square roots yields the claim. Density of $C_c^\infty$ in $W^{1,1}$ extends it to all of $W^{1,1}(\mathbb{R}^2)$ . $\blacksquare$

∎

ProofWeak Derivative of the Absolute Value

f(x)=|x|

on

\mathbb{R}

1. $|x|\in L^1_{loc}$ , hence defines $T_{|x|}$ . For $\varphi\in\mathcal{D}(\mathbb{R})$ ,

13

\quad \langle T_{|x|}', \varphi\rangle = -\int_{\mathbb{R}}|x|\,\varphi'(x)\,dx = -\int_{-\infty}^0(-x)\varphi'\,dx -\int_0^\infty x\,\varphi'\,dx.

Derivation hint: apply definition [4] with $|\alpha|=1$ and split at the corner $x=0$ . 2. Integrate each piece by parts (boundary terms at $\pm\infty$ vanish, and the $x=0$ terms cancel since both equal $0\cdot\varphi(0)$ ):

14

\quad \langle T_{|x|}',\varphi\rangle = -\int_{-\infty}^0(-1)\varphi\,dx - \int_0^\infty(1)\varphi\,dx = \int_{\mathbb{R}}\operatorname{sgn}(x)\,\varphi(x)\,dx.

Derivation hint: $\partial_x(-x)=-1$ and $\partial_x x=1$ produce the sign function as the weak derivative. 3. Hence $\,|x|' = \operatorname{sgn}(x)\,$ weakly — an $L^\infty$ function, so $|x|\in W^{1,p}_{loc}$ for all $p$ . Differentiating once more gives $\operatorname{sgn}' = 2\delta_0$ , not a function: $|x|\notin W^{2,p}$ . $\blacksquare$

∎

PD

Paul Dirac & Laurent Schwartz

1902–1984 / 1915–2002

Dirac introduced the "delta function" $\delta(x)$ in his 1930 Principles of Quantum Mechanics as a symbol with $\int\delta(x)\varphi(x)\,dx=\varphi(0)$ — mathematically scandalous, since no ordinary function does this. For two decades it was a brilliant but unjustified formal trick. Laurent Schwartz (Fields Medal, 1950) resolved the scandal: $\delta$ is not a function but a distribution, the evaluation functional $\langle\delta,\varphi\rangle=\varphi(0)$ , perfectly rigorous within $\mathcal{D}'$ . Schwartz's Théorie des distributions (1950–51) unified Dirac's delta, Heaviside's operational calculus, Sobolev's weak derivatives, and Hadamard's finite parts into a single coherent framework. Sergei Sobolev, independently and earlier (1936–38), had introduced both weak derivatives and the function spaces now bearing his name while studying hyperbolic equations.

🧠 Intuition

Think of a distribution as a measuring device: you cannot read its pointwise value, but you can probe it with any smooth bump $\varphi$ and read off a number $\langle T,\varphi\rangle$ . The delta is a perfect "spike detector" centered at the origin. Differentiation is "passing the derivative to the probe": you can't differentiate a sharp corner, but you can differentiate the smooth $\varphi$ that touches it, and the minus sign bookkeeps the integration-by-parts. For Sobolev spaces, picture a tug-of-war between roughness budget and dimension. The exponent inequality $kp$ vs $n$ decides who wins:

If you have plenty of derivatives relative to dimension ( $kp>n$ ), your function is forced to be continuous — even Hölder continuous.
If derivatives are scarce ( $kp<n$ ), the best you get is improved integrability $L^{p^*}$ , not continuity. In $\mathbb{R}^3$ , $H^1$ functions can blow up at points (e.g. like $|x|^{-1/2+\epsilon}$ ), which is exactly why hydrogen wavefunctions, lying in $H^1$ , can have cusps yet finite energy.

The "energy" picture is the deepest: $\|\nabla f\|_{L^2}^2$ is literally the Dirichlet energy of a membrane or electrostatic field, and minimizing it over $H_0^1$ is solving Laplace's equation.

Figure

Sobolev embedding regimes: roughness budget

kp

versus dimension

n

ExampleHeaviside step and its derivative

Given: the Heaviside function $H(x)=1$ for $x>0$ , $H(x)=0$ for $x<0$ , on $\mathbb{R}$ .

Find: its distributional derivative.

Steps:

15

\quad \langle H',\varphi\rangle = -\langle H,\varphi'\rangle = -\int_0^\infty \varphi'(x)\,dx = -[\varphi(x)]_0^\infty = \varphi(0).

Derivation hint: $\varphi$ has compact support so $\varphi(\infty)=0$ ; the definition [4] supplies the leading minus. Since $\langle H',\varphi\rangle=\varphi(0)=\langle\delta_0,\varphi\rangle$ for all test functions,

\boxed{\,H' = \delta_0\,}

Consequence (numbers): Is $H\in W^{1,1}(0,1)$ ? No — its weak derivative $\delta_0$ is not an $L^1$ function. But $H\in L^p(0,1)$ for all $p$ , and $H\in W^{1,p}$ only on intervals not containing $0$ , where $H'=0$ . The corner-vs-jump distinction is sharp: a corner ( $|x|$ ) keeps you in $W^{1,p}$ ; a jump ( $H$ ) does not.

ExampleIs

f(x)=|x|^{-\alpha}

in

H^1

of the unit ball in

\mathbb{R}^3

?

Given: $f(x)=|x|^{-\alpha}$ on $B = \{|x|<1\}\subset\mathbb{R}^3$ , $\alpha>0$ .

Find: the range of $\alpha$ for which $f\in H^1(B)=W^{1,2}(B)$ .

Step 1 — $L^2$ membership. In spherical coordinates $dx = r^2\,dr\,d\sigma$ :

16

\quad \int_B |f|^2\,dx = 4\pi\int_0^1 r^{-2\alpha}\,r^2\,dr = 4\pi\int_0^1 r^{2-2\alpha}\,dr <\infty \iff 2-2\alpha>-1 \iff \alpha<\tfrac32.

Derivation hint: $\int_0^1 r^\beta\,dr$ converges iff $\beta>-1$ ; here $\beta = 2-2\alpha$ .

Step 2 — gradient in $L^2$ . $|\nabla f| = \alpha\,|x|^{-\alpha-1}$ , so

17

\quad \int_B |\nabla f|^2\,dx = 4\pi\alpha^2\int_0^1 r^{-2\alpha-2}\,r^2\,dr = 4\pi\alpha^2\int_0^1 r^{-2\alpha}\,dr <\infty \iff -2\alpha>-1 \iff \alpha<\tfrac12.

Derivation hint: same convergence test with $\beta=-2\alpha$ ; the gradient condition is strictly more restrictive. The gradient condition dominates. (One also checks no singular distributional part appears at the origin for $\alpha<\tfrac12$ , which holds.)

\boxed{\,f=|x|^{-\alpha}\in H^1(B)\subset\mathbb{R}^3 \iff 0<\alpha<\tfrac12\,}

Interpretation: For $\alpha$ just under $\tfrac12$ , $f$ is unbounded at $0$ yet has finite Dirichlet energy — concretely illustrating that $H^1(\mathbb{R}^3)\not\hookrightarrow C^0$ (here $kp=2<3=n$ , the subcritical regime).

⚠️ Common mistake

"Weak derivatives are just classical derivatives almost everywhere." False in general. The Cantor function (Devil's staircase) is continuous, increasing, with classical derivative $0$ a.e. — yet its distributional derivative is the singular Cantor measure, not the zero function. The two notions agree only when the classical derivative is itself locally integrable and no concentrated mass (jumps, measures) is hidden. Always verify there is no singular part.

⚠️ Common mistake

" $H^1$ functions are continuous." Only when $n=1$ (where $W^{1,p}\hookrightarrow C^{0,1-1/p}$ for $p>1$ ). In $n\ge 2$ , $H^1$ functions can be unbounded, as Example 2 shows. The critical case $kp=n$ is treacherous: $W^{1,n}(\mathbb{R}^n)$ embeds into every $L^q$ , $q<\infty$ , and even into exponential-integrability (Trudinger–Moser) classes — but not into $L^\infty$ .

CounterexampleDistributions cannot always be multiplied

There is no associative product on $\mathcal{D}'$ extending pointwise multiplication and the action on smooth functions. Consider

18

\quad (\delta_0 \cdot x)\cdot \operatorname{pv}\tfrac1x = 0\cdot \operatorname{pv}\tfrac1x = 0, \qquad \delta_0\cdot\big(x\cdot\operatorname{pv}\tfrac1x\big) = \delta_0\cdot 1 = \delta_0.

Derivation hint: $x\,\delta_0 = 0$ (since $\langle x\delta_0,\varphi\rangle=\varphi(0)\cdot 0$ ), while $x\cdot\operatorname{pv}\frac1x = 1$ as distributions; associativity would force $0=\delta_0$ . This Schwartz impossibility theorem is why nonlinear PDEs require extra structure (Sobolev embeddings to give products meaning, Colombeau algebras, or paraproducts/rough-path methods) — you cannot naïvely square a distribution.

connection

Distributions and Sobolev spaces sit at the crossroads of analysis:

<a class="wiki-link" href="/wiki/fourier-transform">Fourier Transform</a>: On the Schwartz space $\mathcal{S}$ , the Fourier transform is a homeomorphism, and it extends to tempered distributions $\mathcal{S}'$ . This gives the clean characterization $H^s(\mathbb{R}^n) = \{f : (1+|\xi|^2)^{s/2}\hat f\in L^2\}$ , valid for real $s$ , defining fractional Sobolev spaces.
Lebesgue Integration and $L^p$ Spaces: Sobolev spaces are built atop $L^p$ ; completeness and density arguments rest on Lebesgue theory.
Weak Solutions of PDEs: The Lax–Milgram theorem produces solutions in $H_0^1$ to elliptic problems by reformulating $-\Delta u = f$ as a bounded coercive bilinear form.
<a class="wiki-link" href="/wiki/spectral-theory-of-self-adjoint-operators">Spectral Theory of Self-Adjoint Operators</a>: $H^k$ are the natural domains; $-\Delta$ is self-adjoint on $H^2\cap H_0^1$ .
<a class="wiki-link" href="/wiki/calculus-of-variations">Calculus of Variations</a>: The direct method minimizes energy functionals over $H^1$ using weak lower semicontinuity and Rellich compactness.

🔬 Where this appears

⏳ History

- 1893–1903 — Heaviside & Hadamard. Oliver Heaviside's operational calculus and Jacques Hadamard's "finite parts" (parties finies) of divergent integrals foreshadow generalized functions without rigor.

- 1907 — Beppo Levi & 1918 — Friedrichs precursors. Early notions of generalized derivatives appear in variational problems.

- 1930 — Dirac. Introduces

\delta(x)

in quantum mechanics as a formal device.

- 1936–1938 — Sergei Sobolev. Defines weak (distributional) derivatives and the spaces

W^{k,p}

while studying the Cauchy problem for hyperbolic equations; proves the first embedding theorems. Independently introduces "functionals" close to distributions.

- 1944–1948 — Friedrichs & mollifiers. Kurt Friedrichs develops mollification (smoothing by convolution), the workhorse for density theorems (

C^\infty

dense in

W^{k,p}

, Meyers–Serrin "

H=W

", 1964).

- 1945 — Laurent Schwartz. Conceives the general theory of distributions as duals of test-function spaces; full Théorie des distributions published 1950–1951.

- 1950 — Fields Medal awarded to Schwartz, largely for distribution theory.

- 1955–1959 — Aronszajn, Gagliardo, Slobodeckij. Define fractional-order Sobolev spaces

W^{s,p}

via Gagliardo seminorms; Nirenberg and Gagliardo prove the interpolation/embedding inequalities.

- 1961 — Lars Hörmander. Linear Partial Differential Operators and the analysis of wavefront sets / microlocal analysis push distribution theory deep into PDE regularity.

open-problem

Nonlinear distribution theory remains incomplete. The Schwartz impossibility result (no associative multiplication) means deeply nonlinear PDEs — Navier–Stokes, nonlinear wave/Schrödinger equations, stochastic PDEs driven by white noise — lack a canonical product of rough objects. Partial remedies have emerged:

Colombeau algebras embed $\mathcal{D}'$ into a differential algebra, but at the cost of a non-canonical product and "association" relations.
Paraproducts (Bony, 1981) and paradifferential calculus decompose products by frequency to control low–high interactions.
Rough paths (Lyons, 1998) and Regularity Structures (Hairer, 2014, Fields Medal) provide a renormalization framework giving meaning to equations like the KPZ and $\Phi^4_3$ models, where products of distributions appear and must be "renormalized" by subtracting infinities.

Open frontiers: a fully canonical, dimension-independent theory of nonlinear operations on distributions; sharp regularity for 3D Navier–Stokes (a Clay Millennium Problem, posed in $H^s$ / $L^p$ -based spaces); and optimal function-space settings (critical Besov/Triebel–Lizorkin spaces) for global well-posedness of dispersive equations.

Formal Statement

Why It Holds / Derivation

Building Intuition

Worked Examples

Common Misconceptions

Connections

Real-World Applications

Historical Development

Subtleties & Limits