Distribution

**Distributions** (or **generalized functions**) are objects that generalize the classical notion of functions in mathematical analysis. Distributions make it possible to differentiate functions whose derivatives do not exist in the classical sense. In particular, any locally integrable function has a distributional derivative. Distributions are widely used in the theory of partial differential equations, where it may be easier to establish the existence of distributional solutions than classical solutions, or appropriate classical solutions may not exist. Distributions are also important in physics and engineering where many problems naturally lead to differential equations whose solutions or initial conditions are distributions, such as the Dirac delta function.

The practical use of distributions can be traced back to the use of Green functions in the 1830s to solve
ordinary differential equations, but was not formalized until much later. According to Kolmogorov & Fomin (1957), generalized functions originated in the work of Sergei Sobolev (1936) on second-order hyperbolic partial differential equations, and the ideas were developed in somewhat extended form by Laurent Schwartz in the late 1940s. According to his autobiography, Schwartz introduced the term "distribution" by analogy with a distribution of electrical charge, possibly including not only point charges but also dipoles and so on. Gårding (1997) comments that although the ideas in the transformative book by Schwartz (1951) were not entirely new, it was Schwartz's broad attack and conviction that distributions would be useful almost everywhere in analysis that made the difference.

Distribution theory reinterprets functions as linear functionals acting on a space of **test functions**. Standard functions act by integration against a test function, but many other linear functionals do not arise in this way, and these are the "generalized functions". There are different possible choices for the space of test functions, leading to different spaces of distributions. The basic space of test function consists of smooth functions with compact support, leading to standard distributions. Use of the space of smooth, rapidly (faster than any polynomial increases) decreasing test functions (these functions are called Schwartz functions) gives instead the tempered distributions, which are important because they have a well-defined distributional Fourier transform. Every tempered distribution is a distribution in the normal sense, but the converse is not true: in general the larger the space of test functions, the more restrictive the notion of distribution. On the other hand, the use of spaces of analytic test functions leads to Sato's theory of hyperfunctions; this theory has a different character from the previous ones because there are no analytic functions with non-empty compact support.

## Basic idea

Distributions are a class of linear functionals that map a set of *test functions* (conventional and well-behaved functions) into the set of real numbers. In the simplest case, the set of test functions considered is D(**R**), which is the set of functions *φ

{\displaystyle \varphi}* : **R** → **R** having two properties:

* *φ

{\displaystyle \varphi}* is smooth (infinitely differentiable);
* *φ

{\displaystyle \varphi}* has compact support (is identically zero outside some bounded interval).

A distribution *T* is a linear mapping *T* : D(**R**) → **R**. Instead of writing *T*(*φ

{\displaystyle \varphi}*), it is conventional to write

⟨
T
,
φ
⟩

{\displaystyle \langle T,\varphi \rangle}

 for the value of *T* acting on a test function *φ

{\displaystyle \varphi}*. A simple example of a distribution is the Dirac delta *δ*, defined by

⟨

δ
,
φ

⟩

=
φ
(
0
)
,

{\displaystyle \left\langle \delta, \varphi \right\rangle = \varphi(0),}

meaning that *δ* evaluates a test function at 0. Its physical interpretation is as the density of a point source.

As described next, there are straightforward mappings from both locally integrable functions and Radon measures to corresponding distributions, but not all distributions can be formed in this manner.

### Functions and measures as distributions

Suppose that *f* : **R** → **R** is a locally integrable function. Then a corresponding distribution *Tf* may be defined by

⟨

T

f

,
φ

⟩

=

∫

R

f
(
x
)
φ
(
x
)

d
x

for

φ
∈
D
(

R

)
.

{\displaystyle \left\langle T\_{f}, \varphi \right\rangle = \int\_\mathbf{R} f(x) \varphi(x) \,dx \quad \text{for} \quad \varphi\in D(\mathbf{R}).}

This integral is a real number which depends linearly and continuously on

φ

{\displaystyle \varphi}

. Conversely, the values of the distribution *Tf* on test functions in D(**R**) determine the pointwise almost everywhere values of the function *f* on **R**. In a conventional abuse of notation, *f* is often used to represent both the original function *f* and the corresponding distribution *Tf*. This example suggests the definition of a distribution as a linear and, in an appropriate sense, continuous functional on the space of test functions D(**R**).

Similarly, if μ is a Radon measure on **R**, then a corresponding distribution *R*μ may be defined by

⟨

R
μ

,
φ

⟩

=

∫

R

φ

d
μ

for

φ
∈
D
(

R

)
.

{\displaystyle \left\langle R\_\mu, \varphi \right\rangle = \int\_{\mathbf{R}} \varphi\, d\mu \quad \text{for} \quad \varphi\in D(\mathbf{R}). }

This integral also depends linearly and continuously on

φ

{\displaystyle \varphi}

, so that *R*μ is a distribution. If μ is absolutely continuous with respect to Lebesgue measure with density *f* and *d*μ = *f* *dx*, then this definition for *R*μ is the same as the previous one for *Tf*, but if μ is not absolutely continuous, then *R*μ is a distribution that is not associated with a function. For example, if *P* is the point-mass measure on **R** that assigns measure one to the singleton set {0} and measure zero to sets that do not contain zero, then

∫

R

φ

d
# P
φ
(
0
)
,

{\displaystyle \int\_{\mathbf{R}} \varphi\, dP = \varphi(0), }

so that *R**P* = *δ* is the Dirac delta.

### Adding and multiplying distributions

Distributions may be multiplied by real numbers and added together, so they form a real vector space.
Distributions may also be multiplied by infinitely differentiable functions, but
it is not possible to define a product of general distributions that extends the usual pointwise product of functions and has the same algebraic properties. This result was shown by Schwartz (1954), and is usually referred to as the
*Schwartz' Impossibility Theorem*.

### Derivatives of distributions

It is desirable to choose a definition for the derivative of a distribution which, at least for distributions derived from smooth functions, has the property that

T
f
′

=

T

f
′

{\displaystyle T'\_f = T\_{f'}}

. If

φ

{\displaystyle \varphi}

 is a test function, we can use integration by parts to see that

⟨

f
′

,
φ

⟩

=

∫

R

f
′

φ

d
# x

[

f
(
x
)
φ
(
x
)

]

−
∞

∞

−

∫

R

f

φ
′

d
# x
−

⟨

f
,

φ
′

⟩

{\displaystyle \left\langle f', \varphi\right\rangle = \int\_{\mathbf{R}} f'\varphi \,dx = \Big[ f(x) \varphi(x) \Big]\_{-\infty}^\infty - \int\_{\mathbf{R}} f\varphi' \,dx = -\left\langle f, \varphi' \right\rangle}

where the last equality follows from the fact that

φ

{\displaystyle \varphi}

 has compact support, so is zero outside of a bounded set. This suggests that if

T

{\displaystyle T}

 is a *distribution*, we should define its derivative

T
′

{\displaystyle T'}

 by

⟨

T
′

,
φ

⟩

=
−

⟨

T
,

φ
′

⟩

.

{\displaystyle \left\langle T', \varphi \right\rangle = - \left\langle T, \varphi' \right\rangle.}

It turns out that this is the proper definition; it extends the ordinary definition of derivative, every distribution becomes infinitely differentiable and the usual properties of derivatives hold.

**Example:** Recall that the Dirac delta (so-called Dirac delta function) is the distribution defined by the equation

⟨

δ
,
φ

⟩

=
φ
(
0
)
.

{\displaystyle \left\langle \delta, \varphi \right\rangle = \varphi(0).}

It is the derivative of the distribution corresponding to the Heaviside step function *H*: For any test function *φ

{\displaystyle \varphi}

,*

⟨

H
′

,
φ

⟩

=
−

∫

−
∞

∞

H
(
x
)

φ
′

(
x
)

d
# x
−
φ
(
∞
)
φ
(
0
# )

⟨

δ
,
φ

⟩

,

{\displaystyle \left\langle H', \varphi \right\rangle = - \int\_{-\infty}^\infty H(x) \varphi'(x) \, dx = - \varphi(\infty) +\varphi(0) = \left\langle \delta, \varphi \right\rangle,}

so *H*′ = *δ*. Note,

φ

{\displaystyle \varphi}

(∞) = 0 because

φ

{\displaystyle \varphi}

 has compact support by our definition of a test function. Similarly, the derivative of the Dirac delta is the distribution defined by the equation

⟨

δ
′

,
φ
# ⟩
−

φ
′

(
0
)
.

{\displaystyle \langle\delta',\varphi\rangle= -\varphi'(0).}

This latter distribution is an example of a distribution that is not derived from a function or a measure. Its physical interpretation is the density of a dipole source. Just as the Dirac impulse can be realized in the weak limit as a sequence of various kinds of constant norm bump functions of ever increasing amplitude and narrowing support, its derivative can by definition be realized as the weak limit of the negative derivatives of said functions, which are now antisymmetric about the eventual distribution's point of singular support.

## Test functions and distributions

In the following, real-valued distributions on an open subset *U* of **R***n* will be formally defined. With minor modifications, one can also define complex-valued distributions, and one can replace **R***n* by any (paracompact) smooth manifold.

The first object to define is the space D(*U*) of test functions on *U*. Once this is defined, it is then necessary to equip it with a topology by defining the limit of a sequence of elements of D(*U*). The space of distributions will then be given as the space of continuous linear functionals on D(*U*).

### Test function space

The space D(*U*) of **test functions** on *U* is defined as follows. A function *φ

{\displaystyle \varphi}* : *U* → **R** is said to have compact support if there exists a compact subset *K* of *U* such that *φ

{\displaystyle \varphi}*(*x*) = 0 for all *x* in *U* \ *K*. The elements of D(*U*) are the infinitely differentiable functions *φ

{\displaystyle \varphi}* : *U* → **R** with compact support – also known as bump functions. This is a real vector space. It can be given a topology by defining the limit of a sequence of elements of D(*U*). A sequence (*φ

{\displaystyle \varphi}**k*) in D(*U*) is said to converge to *φ

{\displaystyle \varphi}* ∈ D(*U*) if the following two conditions hold:

* There is a compact set *K* ⊂ *U* containing the supports of all *φ

{\displaystyle \varphi}**k*:

⋃
k

supp
⁡
(

φ
k

)
⊂
K
.

{\displaystyle \bigcup\nolimits\_k \operatorname{supp}(\varphi\_k)\subset K.}

* For each multi-index α, the sequence of partial derivatives

∂
α

φ
k

{\displaystyle \partial^\alpha \varphi\_k}

 tends uniformly to

∂
α

φ

{\displaystyle \partial^\alpha\varphi}

.

With this definition, D(*U*) becomes a complete locally convex topological vector space satisfying the Heine–Borel property.

This topology can be placed in the context of the following general construction: let

# X

⋃
i

X
i

{\displaystyle X = \bigcup\nolimits\_i X\_i}

be a countable increasing union of locally convex topological vector spaces and ι*i* : *Xi* → *X* be the inclusion maps. In this context, the inductive limit topology, or final topology, τ on *X* is the finest locally convex vector space topology making all the inclusion maps

ι
i

{\displaystyle \iota\_i}

 continuous. The topology τ can be explicitly described as follows: let *β* be the collection of convex balanced subsets *W* of *X* such that *W* ∩ *Xi* is open for all *i*. A base for the inductive limit topology τ then consists of the sets of the form *x* + *W*, where *x* in *X* and *W* in *β*.

The proof that τ is a vector space topology makes use of the assumption that each *Xi* is locally convex. By construction, *β* is a local base for *τ*. That any locally convex vector space topology on *X* must necessarily contain *τ* means it is the weakest one. One can also show that, for each *i*, the subspace topology *Xi* inherits from τ coincides with its original topology. When each *Xi* is a Fréchet space, (*X*, τ) is called an LF space.

Now let *U* be the union of *Ui* where {*Ui*} is a countable nested family of open subsets of *U* with compact closures *Ki* = *U**i*. Then we have the countable increasing union

D

(
U
# )

⋃
i

D

K
i

{\displaystyle \mathrm{D}(U) = \bigcup\nolimits\_i \mathrm{D}\_{K\_i} }

where D*Ki* is the set of all smooth functions on *U* with support lying in *Ki*. On each D*Ki*, consider the topology given by the seminorms

‖
φ

‖

α

=

max

x
∈

K
i

∂

α

φ

,

i.e. the topology of uniform convergence of derivatives of arbitrary order. This makes each D*Ki* a Fréchet space. The resulting LF space structure on D(*U*) is the topology described in the beginning of the section.

On D(*U*), one can also consider the topology given by the seminorms

‖
φ

‖

α
,

K
i

=

max

x
∈

K
i

∂

α

φ

.

However, this topology has the disadvantage of not being complete. On the other hand, because of the particular features of D*Ki*'s, a set this bounded with respect to τ if and only if it lies in some D*Ki*'s. The completeness of (*D*(*U*), τ) then follow from that of D*Ki*'s.

The topology τ is not metrizable by the Baire category theorem, since D(*U*) is the union of subspaces of the first category in D(*U*).

### Distributions

A **distribution** on *U* is a continuous linear functional *T* : D(*U*) → **R** (or *T* : D(*U*) → **C**). That is, a distribution *T* assigns to each test function *φ

{\displaystyle \varphi}* a real (or complex) scalar *T*(*φ

{\displaystyle \varphi}*) such that

T
(

c
1

φ
1

c
2

φ
2

# )

c
1

T
(

φ
1

)

c
2

T
(

φ
2

)

{\displaystyle T(c\_1\varphi\_1 + c\_2\varphi\_2) = c\_1 T(\varphi\_1) + c\_2 T(\varphi\_2)}

for all test functions *φ

{\displaystyle \varphi}*1, *φ

{\displaystyle \varphi}*2 and scalars c1, c2.
Moreover, *T* is continuous if and only if

lim

k
→
∞

T
(

φ
k

# )
T

(

lim

k
→
∞

φ
k

)

{\displaystyle \lim\_{k\to\infty}T(\varphi\_k)= T\left(\lim\_{k\to\infty}\varphi\_k\right)}

for every convergent sequence *φ

{\displaystyle \varphi}**k* in D(*U*). (Even though the topology of D(*U*) is not metrizable, a linear functional on D(*U*) is continuous if and only if it is sequentially continuous.) Equivalently, *T* is continuous if and only if for every compact subset *K* of *U* there exists a positive constant *CK* and a non-negative integer *NK* such that

T
(
φ
)

≤

C
K

sup
{

∂
α

φ
(
x
)

∣
x
∈
K
,

α

≤

N
K

}

for all test functions *φ

{\displaystyle \varphi}* with support contained in *K*.

The space of distributions on *U* is denoted by D′(*U*) and it is the continuous dual space of D(*U*). No matter what dual topology is placed on D′(*U*), a *sequence* of distributions converges in this topology if and only if it converges pointwise (although this need not be true of a net), which is why the topology is sometimes defined to be the weak-\* topology. But often the topology of bounded convergence, which in this case is the same as the topology of uniform convergence on compact sets, is placed on D′(*U*) since it is with this topology that D′(*U*) becomes a nuclear Montel space and it is with this topology that the kernels theorem of Schwartz holds. No matter which topology is chosen, D′(*U*) will be a non-metrizable, locally convex topological vector space.

The duality pairing between a distribution *T* in D′(*U*) and a test function *φ

{\displaystyle \varphi}* in D(*U*) is denoted using angle brackets by

{

D

′

(
U
)
×

D

(
U
)
→

R

(
T
,
φ
)
↦
⟨
T
,
φ
⟩
,

{\displaystyle \begin{cases}
\mathrm{D}'(U) \times \mathrm{D}(U) \to \mathbf{R} \\
(T, \varphi) \mapsto \langle T, \varphi \rangle,
\end{cases}}

so that ⟨*T*,*φ

{\displaystyle \varphi}*⟩ = *T*(*φ

{\displaystyle \varphi}*). One interprets this notation as the distribution *T* acting on the test function *φ

{\displaystyle \varphi}* to give a scalar, or symmetrically as the test function *φ

{\displaystyle \varphi}* acting on the distribution *T*.

A sequence of distributions (*Tk*) converges with respect to the weak-\* topology on D′(*U*) to a distribution *T* if and only if

⟨

T
k

,
φ
⟩
→
⟨
T
,
φ
⟩

{\displaystyle \langle T\_k, \varphi\rangle \to \langle T, \varphi\rangle}

for every test function *φ

{\displaystyle \varphi}* in D(*U*). For example, if *fk* : **R** → **R** is the function

f
k

(
x
# )

{

k

if

0
≤
x
≤
1

/

k

0

otherwise

{\displaystyle f\_k(x) = \begin{cases} k & \text{if}\ 0\le x \le 1/k \\ 0 & \text{otherwise} \end{cases}}

and *Tk* is the distribution corresponding to *fk*, then

⟨

T
k

,
φ
# ⟩
k

∫
0

1

/

k

φ
(
x
)

d
x
→
φ
(
0
# )
⟨
δ
,
φ
⟩

{\displaystyle \langle T\_k, \varphi\rangle = k\int\_0^{1/k} \varphi(x)\, dx \to \varphi(0) = \langle \delta, \varphi\rangle}

as *k* → ∞, so *T**k* → *δ* in D′(**R**). Thus, for large *k*, the function *f**k* can be regarded as an approximation of the Dirac delta distribution.

### Functions as distributions

The function *f* : *U* → **R** is called **locally integrable** if it is Lebesgue integrable over every compact subset *K* of *U*. This is a large class of functions which includes all continuous functions and all *Lp* functions. The topology on D(*U*) is defined in such a fashion that any locally integrable function *f* yields a continuous linear functional on D(*U*) – that is, an element of D′(*U*) – denoted here by *Tf*, whose value on the test function *φ

{\displaystyle \varphi}* is given by the Lebesgue integral:

⟨

T
f

,
φ
# ⟩

∫
U

f
φ

d
x
.

{\displaystyle \langle T\_f,\varphi \rangle = \int\_U f\varphi\,dx.}

Conventionally, one abuses notation by identifying *Tf* with *f*, provided no confusion can arise, and thus the pairing between *Tf* and *φ

{\displaystyle \varphi}* is often written

⟨
f
,
φ
# ⟩
⟨

T
f

,
φ
⟩
.

{\displaystyle \langle f, \varphi\rangle = \langle T\_f,\varphi\rangle.}

If *f* and *g* are two locally integrable functions, then the associated distributions *Tf* and *Tg* are equal to the same element of D′(*U*) if and only if *f* and *g* are equal almost everywhere (see, for instance, Hörmander (1983, Theorem 1.2.5)). In a similar manner, every Radon measure μ on *U* defines an element of D′(*U*) whose value on the test function *φ

{\displaystyle \varphi}* is ∫*φ

{\displaystyle \varphi}* *dμ*. As above, it is conventional to abuse notation and write the pairing between a Radon measure *μ* and a test function *φ

{\displaystyle \varphi}* as

⟨
μ
,
φ
⟩

{\displaystyle \langle \mu, \varphi \rangle}

. Conversely, as shown in a theorem by Schwartz (similar to the Riesz representation theorem), every distribution which is non-negative on non-negative functions is of this form for some (positive) Radon measure.

The test functions are themselves locally integrable, and so define distributions. As such they are dense in D′(*U*) with respect to the topology on D′(*U*) in the sense that for any distribution *T* ∈ D′(*U*), there is a sequence *φ

{\displaystyle \varphi}**n* ∈ D(*U*) such that

⟨

φ
n

,
ψ
⟩
→
⟨
T
,
ψ
⟩

{\displaystyle \langle\varphi\_n,\psi\rangle\to \langle T,\psi\rangle}

for all *Ψ* ∈ D(*U*). This fact follows from the Hahn–Banach theorem, since the dual of D′(*U*) with its weak-\* topology is the space D(*U*), and it can also be proven more constructively by a convolution argument.

## Operations on distributions

Many operations which are defined on smooth functions with compact support can also be defined for distributions. In general, if *A* : D(*U*) → D(*U*) is a linear mapping of vector spaces which is continuous with respect to the weak-\* topology, then it is possible to extend *A* to a mapping *A* : D′(*U*) → D′(*U*) by passing to the limit. (This approach works for non-linear mappings as well, provided they are assumed to be uniformly continuous.)

In practice, however, it is more convenient to define operations on distributions by means of the transpose. If *A* : D(*U*) → D(*U*) is a continuous linear operator, then the transpose is an operator *At* : D(*U*) → D(*U*) such that

∫
U

A
φ
(
x
)
⋅
ψ
(
x
)

d
# x

∫
U

φ
(
x
)
⋅

A
t

ψ
(
x
)

d
x

for all

φ
,
ψ
∈
D
(
U
)
.

{\displaystyle \int\_U A\varphi(x)\cdot \psi(x) \,dx = \int\_U \varphi(x) \cdot A^t\psi(x)\, dx\qquad \text{for all}\ \varphi,\psi\in D(U).}

(For operators acting on spaces of complex-valued test functions, the transpose *At* differs from the adjoint *A\** in that it does not include a complex conjugate.)

If such an operator *At* exists and is continuous on D(*U*), then the original operator *A* may be extended to D′(*U*) by defining *AT* for a distribution *T* as

⟨
A
T
,
φ
# ⟩
⟨
T
,

A
t

φ
⟩

for all

φ
∈
D
(
U
)
.

{\displaystyle \langle AT, \varphi\rangle = \langle T, A^t\varphi\rangle\qquad \text{for all}\ \varphi\in D(U).}

### Differentiation

Suppose *A* : D(*U*) → D(*U*) is the partial derivative operator

A
# φ

∂
φ

∂

x
k

.

{\displaystyle A\varphi = \frac{\partial\varphi}{\partial x\_k}.}

If *φ

{\displaystyle \varphi}* and *ψ* are in D(*U*), then an integration by parts gives

∫
U

∂
φ

∂

x
k

ψ

d
# x
−

∫
U

φ

∂
ψ

∂

x
k

d
x
,

{\displaystyle \int\_U \frac{\partial\varphi}{\partial x\_k} \psi \, dx = -\int\_U\varphi \frac{\partial\psi}{\partial x\_k}\, dx,}

so that *At* = −*A*. This operator is a continuous linear transformation on D(*U*). So, if *T* ∈ D′(*U*) is a distribution, then the partial derivative of *T* with respect to the coordinate *xk* is defined by the formula

⟨

∂
T

∂

x

k

,
φ

⟩

=
−

⟨

T
,

∂
φ

∂

x

k

⟩

for all

φ
∈
D
(
U
)
.

{\displaystyle \left\langle \frac{\partial T}{\partial x\_{k}}, \varphi \right\rangle = - \left\langle T, \frac{\partial \varphi}{\partial x\_{k}} \right\rangle
\qquad \text{for all}\ \varphi\in D(U).}

With this definition, every distribution is infinitely differentiable, and the derivative in the direction *xk* is a linear operator on D′(*U*).

More generally, if α = (α1, ..., α*n*) is an arbitrary multi-index and ∂α is the associated partial derivative operator, then the partial derivative ∂α*T* of the distribution *T* ∈ D′(*U*) is defined by

⟨

∂

α

T
,
φ

⟩

=
(
−
1

)

α

⟨

T
,

∂

α

φ

⟩

 for all

φ
∈

D

(
U
)
.

{\displaystyle \left\langle \partial^{\alpha} T, \varphi \right\rangle = (-1)^{| \alpha |} \left\langle T, \partial^{\alpha} \varphi \right\rangle \mbox{ for all } \varphi \in \mathrm{D}(U).}

Differentiation of distributions is a continuous operator on D′(*U*); this is an important and desirable property that is not shared by most other notions of differentiation.

### Multiplication by a smooth function

If *m* : *U* → **R** is an infinitely differentiable function and *T* is a distribution on *U*, then the product m*T* is defined by

⟨
m
T
,
φ
# ⟩
⟨
T
,
m
φ
⟩

for all

φ
∈
D
(
U
)
.

{\displaystyle \langle mT, \varphi\rangle = \langle T, m\varphi \rangle\qquad \text{for all}\ \varphi\in D(U).}

This definition coincides with the transpose definition since if *M* : D(*U*) → D(*U*) is the operator of multiplication by the function *m* (i.e., *M

φ

{\displaystyle \varphi}* = *m* *φ

{\displaystyle \varphi}*), then

∫
U

M
φ
(
x
)
⋅
ψ
(
x
)

d
# x

∫
U

m
(
x
)
φ
(
x
)
⋅
ψ
(
x
)

d
# x

∫
U

φ
(
x
)
⋅
m
(
x
)
ψ
(
x
)

d
# x

∫
U

φ
(
x
)
⋅
M
ψ
(
x
)

d
x
,

{\displaystyle \int\_U M\varphi(x)\cdot \psi(x)\,dx = \int\_U m(x)\varphi(x)\cdot \psi(x)\,dx = \int\_U \varphi(x)\cdot m(x)\psi(x)\,dx = \int\_U \varphi(x)\cdot M\psi(x)\,dx,}

so that *Mt* = *M*.

Under multiplication by smooth functions, D′(*U*) is a module over the ring C∞(*U*). With this definition of multiplication by a smooth function, the ordinary product rule of calculus remains valid. However, a number of unusual identities also arise. For example, if *δ* is the Dirac delta distribution on **R**, then *mδ* = *m*(0)*δ*, and if *δ*′ is the derivative of the delta distribution, then

m

δ
′

=
m
(
0
)

δ
′

−

m
′

# δ
m
(
0
)

δ
′

−

m
′

(
0
)
δ
.

{\displaystyle m\delta' = m(0)\delta' - m'\delta = m(0)\delta' - m'(0)\delta.\,}

These definitions of differentiation and multiplication also make it possible to define the operation of a linear differential operator with smooth coefficients on a distribution. A linear differential operator *P* takes a distribution *T* ∈ D′(*U*) to another distribution *PT* given by a sum of the form

P
# T

∑

α

≤
k

p
α

∂
α

T
,

{\displaystyle PT = \sum\nolimits\_{|\alpha|\le k} p\_\alpha \partial^\alpha T,}

where the coefficients *p*α are smooth functions on *U*. The action of the distribution *PT* on a test function *φ

{\displaystyle \varphi}* is given by

⟨

∑

α

≤
k

p
α

∂
α

T
,
φ

⟩

=

⟨

T
,

∑

α

≤
k

(
−
1

)

α

∂
α

(

p
α

φ
)

⟩

.

The minimum integer *k* for which such an expansion holds for every distribution *T* is called the **order** of *P*.
The space D′(*U*) is a D-module with respect to the action of the ring of linear differential operators.

### Composition with a smooth function

Let *T* be a distribution on an open set *U* ⊂ **R***n*. Let *V* be an open set in **R***n*, and *F* : *V* → *U*. Then provided *F* is a submersion, it is possible to define

T
∘
F
∈

D

′

(
V
)
.

{\displaystyle T\circ F \in \mathrm{D}'(V).}

This is the **composition** of the distribution *T* with *F*, and is also called the **pullback** of *T* along *F*, sometimes written

F
♯

:
T
↦

F
♯

# T
T
∘
F
.

{\displaystyle F^\sharp : T\mapsto F^\sharp T = T\circ F.}

The pullback is often denoted *F\**, although this notation should not be confused with the use of '\*' to denote the adjoint of a linear mapping.

The condition that *F* be a submersion is equivalent to the requirement that the Jacobian derivative *dF*(*x*) of *F* is a surjective linear map for every *x* ∈ *V*. A necessary (but not sufficient) condition for extending *F*# to distributions is that *F* be an open mapping. The inverse function theorem ensures that a submersion satisfies this condition.

If *F* is a submersion, then *F*# is defined on distributions by finding the transpose map. Uniqueness of this extension is guaranteed since *F*# is a continuous linear operator on D(*U*). Existence, however, requires using the change of variables formula, the inverse function theorem (locally) and a partition of unity argument.

In the special case when *F* is a diffeomorphism from an open subset *V* of **R***n* onto an open subset *U* of **R***n* change of variables under the integral gives

∫
V

φ
∘
F
(
x
)
ψ
(
x
)

d
# x

∫
U

φ
(
x
)
ψ

(

F

−
1

(
x
)

)

det
d

F

−
1

(
x
)

d
x
.

{\displaystyle \int\_V\varphi\circ F(x) \psi(x)\,dx = \int\_U\varphi(x) \psi \left (F^{-1}(x) \right ) \left |\det dF^{-1}(x) \right |\,dx.}

In this particular case, then, *F*# is defined by the transpose formula:

⟨

F
♯

T
,
φ

⟩

=

⟨

T
,

det
d
(

F

−
1

)

φ
∘

F

−
1

⟩

.

{\displaystyle \left \langle F^\sharp T,\varphi \right \rangle = \left \langle T, \left |\det d(F^{-1}) \right | \varphi\circ F^{-1} \right \rangle.}

## Localization of distributions

There is no way to define the value of a distribution in D′(*U*) at a particular point of *U*. However, as is the case with functions, distributions on *U* restrict to give distributions on open subsets of *U*. Furthermore, distributions are *locally determined* in the sense that a distribution on all of *U* can be assembled from a distribution on an open cover of *U* satisfying some compatibility conditions on the overlap. Such a structure is known as a sheaf.

### Restriction

Let *U* and *V* be open subsets of **R***n* with *V* ⊂ *U*. Let *EVU* : D(*V*) → D(*U*) be the operator which *extends by zero* a given smooth function compactly supported in *V* to a smooth function compactly supported in the larger set *U*. Then the restriction mapping ρ*VU* is defined to be the transpose of *EVU*. Thus for any distribution *T* ∈ D′(*U*), the restriction ρ*VU**T* is a distribution in the dual space D′(*V*) defined by

⟨

ρ

V
U

T
,
φ
# ⟩
⟨
T
,

E

V
U

φ
⟩

{\displaystyle \langle \rho\_{VU}T,\varphi\rangle = \langle T, E\_{VU}\varphi\rangle}

for all test functions *φ

{\displaystyle \varphi}* ∈ D(*V*).

Unless *U* = *V*, the restriction to *V* is neither injective nor surjective. Lack of surjectivity follows since distributions can blow up towards the boundary of *V*. For instance, if *U* = **R** and *V* = (0, 2), then the distribution

T
(
x

# n
1

∞

n

δ

(

x
−

1
n

)

{\displaystyle T(x) = \sum\_{n=1}^\infty n\,\delta\left(x-\frac{1}{n}\right)}

is in D′(*V*) but admits no extension to D′(*U*).

### Support of a distribution

Let *T* ∈ D′(*U*) be a distribution on an open set *U*. Then *T* is said to vanish on an open set *V* of *U* if *T* lies in the kernel of the restriction map ρ*VU*. Explicitly *T* vanishes on *V* if

⟨
T
,
φ
# ⟩
0

{\displaystyle \langle T,\varphi\rangle = 0}

for all test functions *φ

{\displaystyle \varphi}* ∈ C∞(*U*) with support in *V*. Let *V* be a maximal open set on which the distribution *T* vanishes; i.e., *V* is the union of every open set on which *T* vanishes. The **support** of *T* is the complement of *V* in *U*. Thus

supp

# T
U
∖
⋃

{

V
∣

ρ

V
U

# T
0

}

.

{\displaystyle \operatorname{supp}\,T = U \setminus \bigcup\left\{V \mid \rho\_{VU}T = 0\right\}.}

The distribution *T* has **compact support** if its support is a compact set. Explicitly, *T* has compact support if there is a compact subset *K* of *U* such that for every test function *φ

{\displaystyle \varphi}* whose support is completely outside of *K*, we have *T*(*φ

{\displaystyle \varphi}*) = 0. Compactly supported distributions define continuous linear functionals on the space C∞(*U*); the topology on C∞(*U*) is defined such that a sequence of test functions *φ

{\displaystyle \varphi}**k* converges to 0 if and only if all derivatives of *φ

{\displaystyle \varphi}**k* converge uniformly to 0 on every compact subset of *U*. Conversely, it can be shown that every continuous linear functional on this space defines a distribution of compact support. The embedding of Cc∞(*U*) into C∞(*U*), where the spaces are given their respective topologies, is continuous and has dense image. Thus compactly supported distributions can be identified with those distributions that can be extended from Cc∞(*U*) to C∞(*U*).

## Tempered distributions and Fourier transform

By using a larger space of test functions *S*(**R***n*),
one can define the space of **tempered distributions**
S′(**R***n*), a subspace of D′(**R***n*). These distributions are useful if one studies the Fourier transform: all tempered distributions have a Fourier transform, but not all distributions in D′(**R***n*) have one.

The space of test functions employed here, the so-called Schwartz space *S*(**R***n*), is the function space of all infinitely differentiable functions that are rapidly decreasing at infinity along with all partial derivatives. Thus *φ* : **R***n* → **R** is in the Schwartz space provided that any derivative of *φ

p

α
,
β

(
φ
# )

sup

x
∈

R

n

x
α

D
β

φ
(
x
)

{\displaystyle p\_{\alpha , \beta} (\varphi) = \sup\_{x \in \mathbf{R}^n} | x^\alpha D^\beta \varphi(x)| }

for *α*, *β* multi-indices of size *n*. Then *φ

{\displaystyle \varphi}* is a Schwartz function if all the values satisfy

p

α
,
β

(
φ
)
<
∞
.

{\displaystyle p\_{\alpha, \beta} (\varphi) < \infty.}

The family of seminorms *p**α*, *β* defines a locally convex topology on the Schwartz space. The seminorms are, in fact, norms on the Schwartz space. The Schwartz space is metrizable and complete. Because the Fourier transform changes differentiation by *x**α* into multiplication by *x**α* and vice versa, this symmetry implies that the Fourier transform of a Schwartz function is also a Schwartz function.

The space of **tempered distributions** is defined as the (continuous) dual of the Schwartz space. In other words, a distribution *T* is a tempered distribution if and only if

lim

m
→
∞

T
(

φ
m

# )

{\displaystyle \lim\_{m\to\infty} T(\varphi\_m)=0. }

is true whenever

lim

m
→
∞

p

α
,
β

(

φ
m

# )
0

{\displaystyle \lim\_{m\to\infty} p\_{\alpha , \beta} (\varphi\_m) = 0 }

holds for all multi-indices *α*, *β*.

The derivative of a tempered distribution is again a tempered distribution. Tempered distributions generalize the bounded (or slow-growing) locally integrable functions; all distributions with compact support and all square-integrable functions are tempered distributions. More generally, all functions that are products of polynomials with elements of *Lp*(**R***n*) for *p* ≥ 1 are tempered distributions.

The *tempered distributions* can also be characterized as *slowly growing*, meaning that each derivative of *T* grows at most as fast as some polynomial. This characterization is dual to the *rapidly falling* behaviour of the derivatives of a function in the Schwartz space, where each derivative of *φ

{\displaystyle \varphi}* decays faster than every inverse power of *|x|*. An example of a rapidly falling function is

x

n

⋅
exp
⁡
(
−
λ

x

β

)

 for any positive *n*, *λ*, *β*.

To study the Fourier transform, it is best to consider *complex*-valued test functions and complex-linear distributions. The ordinary continuous Fourier transform *F* yields then an automorphism of Schwartz function space, and we can define the **Fourier transform** of the tempered distribution *T* by (*FT*)(*Ψ*) = *T*(*Fψ*) for every Schwartz function *Ψ*. *FT* is thus again a tempered distribution. The Fourier transform is a continuous, linear, bijective operator from the space of tempered distributions to itself. This operation is compatible with differentiation in the sense that

F

d
T

d
x

=
i
x
F
T

{\displaystyle F\dfrac{dT}{dx}=ixFT}

and also with convolution: if *T* is a tempered distribution and *Ψ* is a *slowly increasing* infinitely differentiable function on **R***n*, then *ψT* is again a tempered distribution and

F
(
ψ
T
# )
F
ψ
∗
F
T

{\displaystyle F(\psi T)=F\psi\*FT\,}

is the convolution of *FT* and *Fψ*. In particular, the Fourier transform of the constant function equal to 1 is the *δ* distribution.

## Convolution

Under some circumstances, it is possible to define the convolution of a function with a distribution, or even the convolution of two distributions.

Convolution of a test function with a distribution
If *f* ∈ D(**R***n*) is a compactly supported smooth test function, then convolution with *f*,

{

C
f

:

D

(

R

n

)
→

D

(

R

n

)

C
f

:
g
↦
f
∗
g

{\displaystyle \begin{cases}
C\_f : \mathrm{D}(\mathbf{R}^n)\to \mathrm{D}(\mathbf{R}^n) \\
C\_f : g \mapsto f \* g
\end{cases}}

defines a linear operator which is continuous with respect to the LF space topology on D(**R***n*).

Convolution of *f* with a distribution *T* ∈ D′(**R***n*) can be defined by taking the transpose of *Cf* relative to the duality pairing of D(**R***n*) with the space D′(**R***n*) of distributions . If *f*, *g*, *φ

{\displaystyle \varphi}* ∈ D(**R***n*), then by Fubini's theorem

⟨

C
f

g
,
φ

⟩

=

∫

R

n

φ
(
x
)

∫

R

n

f
(
x
−
y
)
g
(
y
)

d
y

d
# x

⟨

g
,

C

f
~

φ

⟩

{\displaystyle \left \langle C\_fg, \varphi \right \rangle = \int\_{\mathbf{R}^n}\varphi(x)\int\_{\mathbf{R}^n}f(x-y) g(y) \, dy \, dx = \left \langle g, C\_{\widetilde{f}}\varphi \right \rangle}

where

f
~

(
x
# )
f
(
−
x
)

{\displaystyle \scriptstyle{\widetilde{f}(x) = f(-x)}}

. Extending by continuity, the convolution of *f* with a distribution *T* is defined by

⟨
f
∗
T
,
φ
# ⟩

⟨

T
,

f
~

∗
φ

⟩

{\displaystyle \langle f\*T, \varphi\rangle = \left \langle T, \widetilde{f}\*\varphi \right \rangle}

for all test functions *φ

{\displaystyle \varphi}* ∈ D(**R***n*).

An alternative way to define the convolution of a function *f* and a distribution *T* is to use the translation operator τ*x* defined on test functions by

τ
x

φ
(
y
# )
φ
(
y
−
x
)

{\displaystyle \tau\_x \varphi(y) = \varphi(y-x)}

and extended by the transpose to distributions in the obvious way. The convolution of the compactly supported function *f* and the distribution *T* is then the function defined for each *x* ∈ **R***n* by

(
f
∗
T
)
(
x
# )

⟨

T
,

τ
x

f
~

⟩

.

{\displaystyle (f\*T)(x) = \left \langle T, \tau\_x\widetilde{f} \right \rangle.}

It can be shown that the convolution of a smooth, compactly supported function and a distribution is a smooth function. If the distribution *T* has compact support as well, then *f*∗*T* is a compactly supported function, and the Titchmarsh convolution theorem implies that

ch
⁡
(
supp
⁡
(
f
∗
T
)
# )
ch
⁡
supp
⁡
f
ch
⁡
supp
⁡
T

{\displaystyle \operatorname{ch}(\operatorname{supp}(f\*T)) = \operatorname{ch}\operatorname{supp}f + \operatorname{ch} \operatorname{supp} T}

where *ch* denotes the convex hull and supp denotes the support.

Distribution of compact support
It is also possible to define the convolution of two distributions *S* and *T* on **R***n*, provided one of them has compact support. Informally, in order to define *S*∗*T* where *T* has compact support, the idea is to extend the definition of the convolution ∗ to a linear operation on distributions so that the associativity formula

S
∗
(
T
∗
φ
# )
(
S
∗
T
)
∗
φ

{\displaystyle S\*(T\*\varphi) = (S\*T)\*\varphi}

continues to hold for all test functions *φ

{\displaystyle \varphi}*.

It is also possible to provide a more explicit characterization of the convolution of distributions . Suppose that it is *T* that has compact support. For any test function *φ

{\displaystyle \varphi}* in D(**R***n*), consider the function

ψ
(
x
# )
⟨
T
,

τ

−
x

φ
⟩
.

{\displaystyle \psi(x) = \langle T, \tau\_{-x} \varphi\rangle.}

It can be readily shown that this defines a smooth function of *x*, which moreover has compact support. The convolution of *S* and *T* is defined by

⟨
S
∗
T
,
φ
# ⟩
⟨
S
,
ψ
⟩
.

{\displaystyle \langle S \* T,\varphi\rangle = \langle S, \psi\rangle.}

This generalizes the classical notion of convolution of functions and is compatible with differentiation in the following sense:

∂
α

(
S
∗
T
# )
(

∂
α

S
)
∗
# T
S
∗
(

∂
α

T
)
.

{\displaystyle \partial^\alpha(S\*T)=(\partial^\alpha S)\*T=S\*(\partial^\alpha T).}

This definition of convolution remains valid under less restrictive assumptions about *S* and *T*.

## Distributions as derivatives of continuous functions

The formal definition of distributions exhibits them as a subspace of a very large space, namely the topological dual of D(*U*) (or S(**R***d*) for tempered distributions). It is not immediately clear from the definition how exotic a distribution might be. To answer this question, it is instructive to see distributions built up from a smaller space, namely the space of continuous functions. Roughly, any distribution is locally a (multiple) derivative of a continuous function. A precise version of this result, given below, holds for distributions of compact support, tempered distributions, and general distributions. Generally speaking, no proper subset of the space of distributions contains all continuous functions and is closed under differentiation. This says that distributions are not particularly exotic objects; they are only as complicated as necessary.

### Tempered distributions

If *f* ∈ *S*′(**R***n*) is a tempered distribution, then there exists a constant *C* > 0, and positive integers *M* and *N* such that for all Schwartz functions *φ

{\displaystyle \varphi}* ∈ *S*(**R***n*)

⟨
f
,
φ
⟩
≤
C

∑

α

≤
N
,

β

≤
M

sup

x
∈

R

n

x
α

D
β

φ
(
x
)

=
C

∑

α

≤
N
,

β

≤
M

p

α
,
β

(
φ
)
.

This estimate along with some techniques from functional analysis can be used to show that there is a continuous slowly increasing function *F* and a multi-index *α* such that

# f

D
α

F
.

{\displaystyle f=D^\alpha F.\,}

### Restriction of distributions to compact sets

If *f* ∈ D′(**R***n*), then for any compact set *K* ⊂ **R***n*, there exists a continuous function *F*  compactly supported
in **R***n* (possibly on a larger set than *K* itself) and a multi-index *α* such that *f* = *D**α**F* on Cc∞(*K*).
This follows from the previously quoted result on tempered distributions by means of a localization argument.

### Distributions with point support

If *f* has support at a single point {*P*}, then *f* is in fact a finite linear combination of distributional derivatives of the *δ* function at *P*. That is, there exists an integer *m* and complex constants *a**α* for multi-indices |*α*| ≤ *m* such that

# f

∑

α

≤
m

a
α

D
α

(

τ
P

δ
)

{\displaystyle f = \sum\nolimits\_{|\alpha|\le m} a\_\alpha D^\alpha(\tau\_P\delta)}

where τ*P* is the translation operator.

### General distributions

A version of the above theorem holds locally in the following sense . Let *T* be a distribution on *U*, then one can find for every multi-index *α* a continuous function *g**α* such that

# T

∑
α

D
α

g
α

{\displaystyle \displaystyle T = \sum\nolimits\_\alpha D^\alpha g\_\alpha}

and that any compact subset *K* of *U* intersects the supports of only finitely many *g**α*; therefore, to evaluate the value of *T* for a given smooth function *f* compactly supported in *U*, we only need finitely many *g**α*; hence the infinite sum above is well-defined as a distribution. If the distribution *T* is of finite order, then one can choose *g**α* in such a way that only finitely many of them are nonzero.

## Using holomorphic functions as test functions

The success of the theory led to investigation of the idea of hyperfunction, in which spaces of holomorphic functions are used as test functions. A refined theory has been developed, in particular Mikio Sato's algebraic analysis, using sheaf theory and several complex variables. This extends the range of symbolic methods that can be made into rigorous mathematics, for example Feynman integrals.

## Problem of multiplication

It is easy to define the product of a distribution with a smooth function, or more generally the product of two distributions whose singular supports are disjoint. With more effort it is possible to define a well-behaved product of several distributions provided their wave front sets at each point are compatible.
A limitation of the theory of distributions (and hyperfunctions) is that there is no associative product of two distributions extending the product of a distribution by a smooth function, as has been proved by Laurent Schwartz in the 1950s. For example, if p.v. 1/*x* is the distribution obtained by the Cauchy principal value

(

p
.
v
.

⁡

1
x

)

[
φ
# ]

lim

ε
→

0

∫

x

≥
ε

φ
(
x
)

x

d
x

{\displaystyle \left(\operatorname{p.v.}\frac{1}{x}\right)[\varphi] = \lim\_{\varepsilon\to 0^+} \int\_{|x|\ge\varepsilon} \frac{\varphi(x)}{x}\, dx}

for all *φ

{\displaystyle \varphi}* ∈ *S*(**R**), and *δ* is the Dirac delta distribution then

(

δ
×
x

)

×

p
.
v
.

⁡

1
x

=
0

{\displaystyle \left(\delta \times x \right) \times \operatorname{p.v.} \frac{1}{x} = 0}

but

δ
×

(

x
×

p
.
v
.

⁡

1
x

)

=
δ

{\displaystyle \delta \times \left( x \times \operatorname{p.v.} \frac{1}{x} \right) = \delta}

so the product of a distribution by a smooth function (which is always well defined) cannot be extended to an associative product on the space of distributions.

Thus, nonlinear problems cannot be posed in general and thus not solved within distribution theory alone. In the context of quantum field theory, however, solutions can be found. In more than two spacetime dimensions the problem is related to the regularization of divergences. Here Henri Epstein and Vladimir Glaser developed the mathematically rigorous (but extremely technical) *causal perturbation theory*. This does not solve the problem in other situations. Many other interesting theories are non linear, like for example the Navier–Stokes equations of fluid dynamics.

Several not entirely satisfactory theories of algebras of generalized functions have been developed, among which Colombeau's (simplified) algebra is maybe the most popular in use today.

Inspired by Lyons' rough path theory, Martin Hairer proposed a consistent way of multiplying distributions with certain structure (regularity structures ), available in many examples from stochastic analysis, notably stochastic partial differential equations. See also Gubinelli–Imkeller–Perkowski (2015) for a related development based on Bony's paraproduct from Fourier analysis.