Distributions and Sobolev Spaces Introduction

⚡ Key Takeaways

Distributions (generalized functions) extend the notion of a function so that every distribution is infinitely differentiable, while Sobolev spaces $W^{k,p}$ measure the integrability of a function together with its weak (distributional) derivatives up to order $k$ . Together they form the analytic backbone of modern PDE theory: weak formulations, existence/uniqueness via Hilbert-space methods, and regularity estimates. They apply wherever derivatives must be interpreted past the limits of classical calculus — shocks, point sources, finite elements, quantum mechanics. Foundational work by Sergei Sobolev (1930s) and Laurent Schwartz (1940s, Fields Medal 1950).

💡 The Core Insight

Differentiation is a chore for functions but trivial for distributions. Classical calculus forbids differentiating non-smooth functions; distribution theory dodges this by defining the derivative through the identity $\langle u', \varphi\rangle = -\langle u, \varphi'\rangle$ , which always makes sense because the test function $\varphi$ is smooth. Sobolev spaces then quantify how many such weak derivatives an object has and how integrable they are — turning "regularity" into a measurable, topological property and converting PDE problems into geometry in Banach/Hilbert spaces.

DefinitionTest Functions and Distributions

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 \text{ is compact}, \ \operatorname{supp}\varphi \subset \Omega\}.

Derivation hint: A function is in $C_c^\infty$ iff it is infinitely differentiable and vanishes outside a compact set; the standard bump $\varphi(x)=e^{-1/(1-|x|^2)}\mathbf 1_{|x|<1}$ exhibits non-emptiness. A sequence $\varphi_j \to \varphi$ in $\mathcal{D}(\Omega)$ iff there is a single compact $K$ containing all $\operatorname{supp}\varphi_j$ and $\partial^\alpha \varphi_j \to \partial^\alpha \varphi$ uniformly for every multi-index $\alpha$ . A distribution $u$ is a continuous linear functional $u : \mathcal{D}(\Omega) \to \mathbb{R}$ (or $\mathbb{C}$ ); the space of all distributions is $\mathcal{D}'(\Omega)$ . Continuity is equivalent to: for each compact $K \subset \Omega$ there exist $C>0$ and $N \in \mathbb{N}$ with

2

\quad |\langle u, \varphi\rangle| \le C \sum_{|\alpha|\le N} \sup_{x\in K} |\partial^\alpha \varphi(x)| \qquad \forall\, \varphi \in \mathcal{D}(\Omega),\ \operatorname{supp}\varphi\subset K.

Derivation hint: This seminorm bound is the explicit form of sequential continuity on the LF-space $\mathcal{D}(\Omega)$ ; $N$ is the local order of $u$ .

DefinitionDistributional (Weak) Derivative

For $u \in \mathcal{D}'(\Omega)$ and multi-index $\alpha = (\alpha_1,\dots,\alpha_n)$ , the derivative $\partial^\alpha u \in \mathcal{D}'(\Omega)$ is defined by

3

\quad \langle \partial^\alpha u, \varphi\rangle = (-1)^{|\alpha|}\langle u, \partial^\alpha \varphi\rangle, \qquad |\alpha| = \alpha_1+\cdots+\alpha_n.

Derivation hint: Repeated integration by parts on a smooth $u$ , with boundary terms killed by compact support of $\varphi$ , forces the sign $(-1)^{|\alpha|}$ ; one then promotes this identity to a definition for arbitrary $u$ . A locally integrable $u \in L^1_{\mathrm{loc}}(\Omega)$ is identified with the distribution $\varphi \mapsto \int_\Omega u\varphi\,dx$ . When $\partial^\alpha u$ (distributional) is again represented by an $L^1_{\mathrm{loc}}$ function $g$ , we call $g$ the weak derivative:

4

\quad \int_\Omega g\,\varphi\,dx = (-1)^{|\alpha|}\int_\Omega u\,\partial^\alpha\varphi\,dx \qquad \forall \varphi\in C_c^\infty(\Omega).

Derivation hint: Specialize [3] to a regular distribution and demand that the result be represented by an $L^1_{\mathrm{loc}}$ kernel.

DefinitionSobolev Spaces

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

5

\quad W^{k,p}(\Omega) = \{ u \in L^p(\Omega) : \partial^\alpha u \in L^p(\Omega)\ \text{for all } |\alpha| \le k\},

with norm

6

\quad \|u\|_{W^{k,p}(\Omega)} = \Big(\sum_{|\alpha|\le k} \|\partial^\alpha u\|_{L^p(\Omega)}^p\Big)^{1/p} \quad (1\le p<\infty), \qquad \|u\|_{W^{k,\infty}} = \max_{|\alpha|\le k}\|\partial^\alpha u\|_{L^\infty}.

Derivation hint: Each $\partial^\alpha u$ is the weak derivative from [4]; the norm is the $\ell^p$ -aggregate of the component $L^p$ norms, making $W^{k,p}$ a closed subspace-style construction. When $p=2$ we write $H^k(\Omega) := W^{k,2}(\Omega)$ , a Hilbert space with inner product

7

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

Derivation hint: Polarize the $p=2$ norm; completeness follows from completeness of $L^2$ on each derivative component. Define $W^{k,p}_0(\Omega)$ as the closure of $C_c^\infty(\Omega)$ in the $W^{k,p}$ norm (zero boundary trace), and $H^{-k}(\Omega) = (H^k_0(\Omega))'$ for negative-order spaces.

TheoremCompleteness, Embedding, Density

1. (Banach/Hilbert) $W^{k,p}(\Omega)$ is a Banach space for $1 \le p \le \infty$ ; $H^k(\Omega)$ is a separable Hilbert space for $1\le p<\infty$ . 2. (Meyers–Serrin, " $H=W$ ") For $1\le p<\infty$ , $C^\infty(\Omega)\cap W^{k,p}(\Omega)$ is dense in $W^{k,p}(\Omega)$ . 3. (Sobolev Embedding, $\Omega=\mathbb{R}^n$ , $kp<n$ ) With $p^* = \frac{np}{n-kp}$ (the Sobolev conjugate),

8

\quad W^{k,p}(\mathbb{R}^n) \hookrightarrow L^{p^*}(\mathbb{R}^n), \qquad \|u\|_{L^{p^*}} \le C(n,k,p)\,\|u\|_{W^{k,p}}.

Derivation hint: For $k=1$ the Gagliardo–Nirenberg–Sobolev inequality follows from the fundamental theorem of calculus along each axis plus the generalized Hölder inequality; iterate for higher $k$ . 4. (Morrey, $kp>n$ ) $W^{k,p}(\mathbb{R}^n) \hookrightarrow C^{0,\gamma}(\mathbb{R}^n)$ with Hölder exponent $\gamma = k - n/p$ (mod 1), so functions become continuous (even Hölder). 5. (Rellich–Kondrachov) If $\Omega$ is bounded with Lipschitz boundary and $kp<n$ , then $W^{k,p}(\Omega) \hookrightarrow\hookrightarrow L^{q}(\Omega)$ compactly for all $1\le q<p^*$ .

ProofSobolev Inequality for

W^{1,1}(\mathbb R^n)

— the prototype

We prove the sharp scaling case $p=1$ , $k=1$ : $\|u\|_{L^{n/(n-1)}} \le C\|\nabla u\|_{L^1}$ for $u\in C_c^\infty(\mathbb R^n)$ .

Step 1 — One-dimensional control. For each coordinate $i$ , the fundamental theorem of calculus gives

9

\quad |u(x)| \le \int_{-\infty}^{\infty} |\partial_i u(x_1,\dots,t,\dots,x_n)|\,dt =: g_i(\hat x_i),

Derivation hint: $u(x)=\int_{-\infty}^{x_i}\partial_i u\,dt$ and $u$ has compact support, so bound the integral over all of $\mathbb R$ ; $g_i$ does not depend on $x_i$ .

Step 2 — Multiply the $n$ estimates. Raise to the power $1/(n-1)$ and multiply over $i$ :

10

\quad |u(x)|^{\frac{n}{n-1}} \le \prod_{i=1}^n \Big(\int_{\mathbb R}|\partial_i u|\,dt\Big)^{\frac{1}{n-1}}.

Derivation hint: Each factor $|u|^{1/(n-1)}\le g_i^{1/(n-1)}$ ; multiplying all $n$ of them produces $|u|^{n/(n-1)}$ on the left.

Step 3 — Integrate iteratively with generalized Hölder. Integrate in $x_1$ , pulling out the factor independent of $x_1$ , and apply Hölder with $n-1$ exponents each equal to $n-1$ . Repeating in $x_2,\dots,x_n$ yields

11

\quad \int_{\mathbb R^n}|u|^{\frac{n}{n-1}}\,dx \le \prod_{i=1}^n \Big(\int_{\mathbb R^n}|\partial_i u|\,dx\Big)^{\frac{1}{n-1}}.

Derivation hint: At each integration the term lacking the integration variable comes out; the remaining $n-1$ terms are paired against each other via Hölder’s inequality with conjugate exponents $(n-1,\dots,n-1)$ .

Step 4 — AM–GM finish. Take the $(n-1)/n$ power and apply the arithmetic–geometric mean inequality to the product:

12

\quad \|u\|_{L^{n/(n-1)}} \le \prod_{i=1}^n \|\partial_i u\|_{L^1}^{1/n} \le \frac1n \sum_{i=1}^n \|\partial_i u\|_{L^1} \le \|\nabla u\|_{L^1}.

Derivation hint: Geometric mean $\le$ arithmetic mean of the $n$ quantities $\|\partial_i u\|_{L^1}$ ; absorb constants. Density of $C_c^\infty$ extends the inequality to all of $W^{1,1}$ . The general $W^{1,p}$ case applies this to $|u|^{\gamma}$ with $\gamma=\frac{(n-1)p}{n-p}$ and uses Hölder once more. $\qquad\blacksquare$

∎

LS

Laurent Schwartz

1915–2002

Schwartz unified decades of "improper" objects — Heaviside’s operational calculus, Dirac’s delta, Hadamard’s finite parts — into the rigorous theory of distributions in Théorie des distributions (1950–51). His decisive idea was duality: define generalized functions as continuous linear functionals on a carefully topologized space of test functions, so that operations on the singular object are dual to operations on smooth functions. For this he received the Fields Medal (1950). He emphasized that the "trick" was less in the objects than in choosing the right test-function topology, the inductive-limit space $\mathcal D$ , whose dual is large enough to host delta and its derivatives.

🧠 Intuition

Think of a distribution as a measuring device, not a pointwise object. You never ask "what is $u(x)$ ?"; you ask "what is the weighted average $\langle u,\varphi\rangle$ against this smooth probe $\varphi$ ?" The Dirac delta is the device that returns $\varphi(0)$ — it reads a value, not an average, which no ordinary function can do. Differentiation works by passing the derivative to the probe. Since you can’t differentiate $u$ directly, you let the smooth $\varphi$ absorb the derivative, paying a sign $(-1)^{|\alpha|}$ as the price of integration by parts. Because $\varphi$ is infinitely smooth, this never fails — hence every distribution is infinitely differentiable. Sobolev spaces add a bookkeeping ledger: not just "does the derivative exist as a distribution," but "is it integrable to power $p$ ?" The number $k$ counts derivatives; $p$ counts integrability. The Sobolev embedding theorem is the exchange rate: spend enough derivatives (large $kp$ relative to dimension $n$ ) and you buy back classical continuity. The competition is always against dimension $n$ — high dimension makes regularity expensive.

Figure

Duality of differentiation: the weak derivative passes

d/dx

onto the smooth test function

Worked Examples

ExampleIs

u(x)=|x|^{-s}

in

H^1

of the unit ball in

\mathbb R^3

?

Given: $\Omega = B(0,1)\subset\mathbb R^3$ ( $n=3$ ), $u(x)=|x|^{-s}$ , $s>0$ . Determine the threshold $s$ for which $u\in H^1(\Omega)=W^{1,2}$ .

Step 1 — $L^2$ condition. Using spherical coordinates $dx = 4\pi r^2\,dr$ :

\int_\Omega |u|^2\,dx = 4\pi\int_0^1 r^{-2s}\,r^2\,dr = 4\pi\int_0^1 r^{2-2s}\,dr,

finite iff $2-2s>-1$ , i.e. $s < 3/2$ .

Step 2 — Gradient. $|\nabla u| = s\,|x|^{-s-1}$ , so

\int_\Omega |\nabla u|^2\,dx = 4\pi s^2\int_0^1 r^{-2s-2}\,r^2\,dr = 4\pi s^2\int_0^1 r^{-2s}\,dr,

finite iff $-2s>-1$ , i.e. $s < 1/2$ .

Step 3 — Intersect. The binding constraint is the gradient:

\boxed{u=|x|^{-s}\in H^1(B_1\subset\mathbb R^3) \iff 0 < s < \tfrac12.}

Cross-check with embedding. Here $kp = 1\cdot 2 = 2 < n=3$ , so $H^1$ does not embed into $L^\infty$ in 3D — unbounded functions like $|x|^{-s}$ , $s<1/2$ , are legitimately in $H^1$ . The Sobolev conjugate is $p^*=\frac{3\cdot2}{3-2}=6$ ; indeed $|x|^{-s}\in L^6$ iff $6s<3$ iff $s<1/2$ , consistent.

⚠️ Common mistake

"Weak derivative = classical derivative wherever the function is differentiable." False in general. The weak derivative is determined by behavior on sets of positive measure plus singular structure. The Cantor–Lebesgue (devil's staircase) function $c(x)$ has classical derivative $0$ almost everywhere, yet its distributional derivative is the singular Cantor measure, not the zero function. So $c \notin W^{1,1}$ even though $c'=0$ a.e. — the pointwise-a.e. derivative ignores the measure concentrated on a null set.

⚠️ Common mistake

"Sobolev functions have pointwise values." Generally false. Elements of $W^{k,p}$ are equivalence classes modulo a.e. equality. A "value at a point" makes sense only after the trace theorem (on lower-dimensional sets) or Morrey embedding (when $kp>n$ gives a continuous representative). Asking $u(0)$ for $u\in H^1(\mathbb R^3)$ is meaningless without further structure.

CounterexampleMultiplication of distributions is impossible in general

There is no associative product on $\mathcal D'(\mathbb R)$ extending pointwise multiplication and the usual calculus (Schwartz's impossibility). Naively, $\delta \cdot H$ would need to satisfy contradictory rules. Consider, with $H^2 = H$ and $H' = \delta$ :

\delta = (H)' = (H^2)' = 2HH' = 2H\delta \implies \delta = 2H\delta.

Yet also $H^3=H$ gives $\delta = (H^3)' = 3H^2 H' = 3H\delta$ , so $2H\delta = 3H\delta$ , forcing $H\delta=0$ and hence $\delta=0$ — absurd. Moral: products like $\delta^2$ or $H\delta$ are undefined; one needs Colombeau algebras or microlocal wavefront-set conditions to multiply selectively.

connection

Distributions and Sobolev spaces sit at the crossroads of analysis:

<a class="wiki-link" href="/wiki/fourier-transform">Fourier Transform</a> — On $\mathbb R^n$ , $H^s$ is characterized by $\int (1+|\xi|^2)^s |\hat u(\xi)|^2\,d\xi<\infty$ , defining fractional and negative order spaces $H^s$ for all real $s$ . The transform turns differentiation into multiplication by $i\xi$ , making smoothness a decay condition on $\hat u$ .
Lp Spaces — Sobolev spaces are built on top of $L^p$ ; the embedding theorems trade derivative count $k$ against integrability $p$ and dimension $n$ .
Weak Solutions of PDEs — The Lax–Milgram theorem on $H^1_0$ produces unique weak solutions of elliptic equations; this is the payoff of the whole machinery.
<a class="wiki-link" href="/wiki/hilbert-spaces">Hilbert Spaces</a> — $H^k$ inherits orthogonality, projections, and the Riesz representation theorem, enabling variational methods.
Tempered Distributions — The dual $\mathcal S'$ of the Schwartz space $\mathcal S$ is the natural home for the Fourier transform, intermediate between $\mathcal D'$ and finite-energy spaces.
Measure Theory — Radon measures are exactly the distributions of order $0$ , linking the abstract dual to classical integration.

🔬 Where this appears

⏳ History

- 1822 — Fourier introduces what becomes the delta as a limit of his integral kernel, used freely but without rigor.

- 1894–1920s — Heaviside develops operational calculus; physicists differentiate discontinuities pragmatically, drawing criticism from rigorists.

- 1908 — Hadamard introduces finite parts (parties finies) of divergent integrals — an early regularization that foreshadows distributions.

- 1926–1927 — Dirac formalizes the delta "function"

\delta

in quantum mechanics; mathematically illegitimate but operationally indispensable.

- 1935–1938 — Sergei Sobolev (Leningrad/Moscow) defines generalized functions and the spaces

W^{l,p}

(now bearing his name) to solve the wave/Cauchy problem; proves the embedding inequality. This is the genuine first rigorous theory.

- 1936 — Bochner & Friedrichs develop related notions of generalized derivatives and mollification ("Friedrichs mollifiers").

- 1944–1948 — Laurent Schwartz crystallizes the duality framework, choosing the topology on

\mathcal D

and identifying

\mathcal D'

.

- 1950–1951 — Schwartz publishes Théorie des distributions (2 vols); receives the Fields Medal (1950). Tempered distributions

\mathcal S'

and the distributional Fourier transform follow.

- 1955 — Lax / Sobolev-Kondrachov / Rellich mature compact embeddings; Meyers–Serrin (1964) prove "

H=W

", unifying two historical definitions of Sobolev spaces.

- 1960s–70s — Hörmander founds microlocal analysis and the wavefront set, refining where distributions can be multiplied and how singularities propagate.

open-problem

Multiplication and nonlinearity. Distribution theory is linear by design; nonlinear PDEs (Navier–Stokes, nonlinear wave/Schrödinger, stochastic PDEs) require multiplying rough objects, which is generically forbidden (see the impossibility counterexample). Partial resolutions — Colombeau algebras, paraproducts (Bony), rough paths (Lyons), and regularity structures (Hairer, Fields Medal 2014) — provide self-consistent calculi for specific equations, but no universal nonlinear theory of distributions exists.

Global regularity of Navier–Stokes (Clay Millennium) is the headline open question framed entirely in Sobolev language: do 3D weak (Leray–Hopf) solutions in $L^2_t H^1_x$ remain smooth for all time, or can the $H^1$ norm blow up?

Optimal constants and rigidity. The sharp constant in the Sobolev inequality [8] and the classification of extremizers (Aubin–Talenti functions) connect to the Yamabe problem in conformal geometry; sharp traces, fractional analogues ( $H^s$ , $W^{s,p}$ via Gagliardo seminorms), and stability of these inequalities remain active.

Distributions and Sobolev Spaces Introduction

Formal Statement

Why It Holds

Building Intuition

Worked Examples

Common Misconceptions

Connections

Real-World Applications

Historical Development

Subtleties & Limits

See Also