Question regarding the Dirac delta function

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
7 replies · 4K views
Kevin McHugh
Messages
319
Reaction score
165
Given the definition:

δ(x) = 0 for all x ≠ 0
∞ for x = 0

-∞δ(x)dx = 1

I don't understand how the integral can equal unity. The integral from -∞ to zero is zero, and the integral from 0 to ∞ is zero. How can one integrate the discontinuity at zero to equal one? What am I missing? TIA for your insight.
 
Physics news on Phys.org
The Dirac delta isn't really a function in the strict mathematical sense. One way to think of it is as the limit of a sequence of functions that becomes more and more sharply peaked. For example, consider the function f(x) = 0 for x < -a; 1/2a for -a ≤ x ≤ +a; 0 for x > +a. Its graph is a rectangle with area (integral) = 1. Now imagine making a smaller and smaller. The rectangle becomes a narrower and narrower and taller and taller "spike." Now take the limit as a --> 0.

It turns out you can make this work in a mathematically well-defined way, so long as you keep this "function" inside an integral, and you can derive properties like $$\int^{+\infty}_{-\infty} f(x)\delta(x) dx = f(0)$$ and $$\int^{+\infty}_{-\infty} f(x)\delta(x-a) dx = f(a)$$ Maybe a mathematician will come along and make this description more precise, but I think this is the way most physicists think of it.
 
  • Like
Likes   Reactions: vanhees71
Kevin McHugh said:
Given the definition:

δ(x) = 0 for all x ≠ 0
∞ for x = 0

-∞δ(x)dx = 1

I don't understand how the integral can equal unity. The integral from -∞ to zero is zero, and the integral from 0 to ∞ is zero. How can one integrate the discontinuity at zero to equal one? What am I missing? TIA for your insight.

There are two different ways to think of the delta function: (1) As a limit of functions, and (2) as a distribution.

As a limit of functions:
One approach to dealing with the delta function that almost always works is to think of it as an ordinary function that has a parameter, [itex]W[/itex] that in some sense measures how sharp the function is. The delta function is sort of (in a sense to be explained) a limit as [itex]W \rightarrow 0[/itex]

Some examples:
  1. [itex]\delta(x) \approx \frac{1}{W \sqrt{\pi}} e^{-\frac{x^2}{W^2}}[/itex]
  2. [itex]\delta(x) \approx \frac{1}{\pi} \frac{sin(\frac{x}{W})}{x}[/itex]
  3. [itex]\delta(x) \approx \frac{1}{2W}[/itex], if [itex]-W \leq x \leq +W[/itex]), and zero otherwise.
In all three cases, the function tends to zero when [itex]|x| \gg W[/itex], and the integral over all space is equal to 1.

These functions don't actually have a limit as [itex]W \rightarrow 0[/itex], but in most cases of interest, integrals involving those functions do have limits as [itex]W \rightarrow 0[/itex], so you can solve the problem under the assumption that [itex]W[/itex] is small, but nonzero, and then let [itex]W \rightarrow 0[/itex] at the end.

As a distribution:
This is the mathematically rigorous way to reason about delta functions.

Let [itex]\mathcal{H}[/itex] be a set of square-integrable complex-valued functions of the real numbers. Then we can say that [itex]F[/itex] is a "distribution" for [itex]\mathcal{H}[/itex] if
  1. [itex]F[/itex] takes any element of [itex]\mathcal{H}[/itex] and returns a complex number. That is, if [itex]f[/itex] is a function in [itex]\mathcal{H}[/itex], then [itex]F(f)[/itex] is a complex number.
  2. [itex]F[/itex] is linear: For any two functions [itex]f(x)[/itex] and [itex]g(x)[/itex], and for any complex number [itex]\alpha[/itex], if [itex]h(x) = f(x) + \alpha g(x)[/itex], then [itex]F(h) = F(f) + \alpha F(g)[/itex].
So any function satisfying 1 and 2 is a distribution. We can then use integrals to define some particular distributions: For any function [itex]\phi[/itex] in [itex]H[/itex], there is a corresponding distribution [itex]F_\phi[/itex] defined via:

[itex]F_\phi(f) = \int_{-\infty}^{+\infty} \phi^*(x) f(x) dx[/itex] (where [itex]^*[/itex] means complex-conjugate)

But there are distributions that don't correspond to any actual function [itex]\phi[/itex]. For example:

[itex]F(f) \equiv f(0)[/itex]

That's a distribution, since it takes any function and returns a complex number, and it's linear. But it doesn't correspond to [itex]F_\phi[/itex] for any actual function [itex]\phi[/itex]. But we can pretend that it is of that form, by writing: [itex]F(f) = \int_{-\infty}^{+\infty} \delta(x) f(x) dx[/itex]. The right-hand side is just suggestive notation; it's not really an integral, although for many purposes, it works as if it were an integral.
 
  • Like
Likes   Reactions: vanhees71, ShayanJ and PeroK
Thanks Steven, I did read the wiki definitions before I asked the question. I think where i am having difficulty is at x = 0. By definition, f(x) is infinity at zero. I see that as a discontinuity. There are no limits of integration. Where is the Δx? The limit of a function view implies there is a Δx→ 0. (Or the W parameter in the example). I might be being pig headed about this.
 
Yeah, you are being a tad stubborn Kevin. The delta function is not defined as infinity at x=0 and zero elsewhere. Like you say, a function defined that way is a discontinuity, and it's integral is not defined.
Instead, one definition is that the area under the delta function, δ(x) dx = 1. (Imagine the integral is from -∞ to ∞. I can't get the limits to print right.)
Another definition is that:

-∞ f(x) δ(x) dx = f(0) for any suitable f(x).

The way I see it, according to the 1st definition, δ(x) dx = 1⋅δ(x) dx = 1. The latter puts the integral in the form required by the 2nd definition as long as f(x) = f(0) = 1 for all x. In other words, f(x)⋅δ(x) dx = f(0) = 1. The 2 definitions are equivalent.

 
Mark Harder said:
Yeah, you are being a tad stubborn Kevin. The delta function is not defined as infinity at x=0 and zero elsewhere. Like you say, a function defined that way is a discontinuity, and it's integral is not defined.
Instead, one definition is that the area under the delta function, δ(x) dx = 1. (Imagine the integral is from -∞ to ∞. I can't get the limits to print right.)
Another definition is that:

-∞ f(x) δ(x) dx = f(0) for any suitable f(x).

The way I see it, according to the 1st definition, δ(x) dx = 1⋅δ(x) dx = 1. The latter puts the integral in the form required by the 2nd definition as long as f(x) = f(0) = 1 for all x. In other words, f(x)⋅δ(x) dx = f(0) = 1. The 2 definitions are equivalent.
The first definition above is not correct since the constant function is not an appropriate test function. The second one is correct. You always need to use an appropriate test function with the ##\delta## distribution under the integral.