# The Dirac Delta function

by arildno
Tags: delta, dirac, function
 PF Patron HW Helper Sci Advisor Thanks P: 11,935 1. INTRODUCTION Many students become frustrated when they first meet the Dirac Delta function, typically in a course involving electrostatics, or Laplace transforms. As it is commonly presented, the Dirac function seems totally meaningless: Either, it is "defined" as: $$\int_{-\infty}^{\infty}\delta(x)dx=1,\delta(x)=0,x\neq0$$ or, in conjunction with an arbitrary function f: $$\int_{I}\delta(x)f(x)dx=f(0),0\in{I},\int_{I}\delta(x)f(x)dx=0,0\notin{ I}$$ where I is any interval on the real axis. What sort of function can do these sort of things?? The student's confusion is not markedly lessened when: the lecturer praises the Dirac function's ability to sample a function value, then the confusion increases when a mathematician scoffs at the whole "definition" and says "it is just a formalism, you've got to think of it as a distribution", and when the lecturer is confronted by this, he backs off a bit and says "yeah, the mathematician is right, but we won't get into distributions in this course, and, anyways, the theory of distributions is unnecessary for PRACTICAL purposes, and that's what this course is all about.." What shall the poor student think of all of this? My aim in creating this thread is to make an ultra-short intro to the mathematical way of looking at it. It won't be very technical, but fairly rigorous, yet, I hope easily understood. In the next post, I'll look into functionals, and the last post will provide the link between the functional perspective and the common way of "defining" the Dirac function.
 PF Patron HW Helper Sci Advisor Thanks P: 11,935 A functional can be thought of as a function whose domain is a function space and some range, for example the real numbers. For example, when we compute "the definite integral" of some (integrable) function on an interval [a,b], we may think of that computation operation as an action done upon f by a functional J(f): $$J(f)=\int_{a}^{b}f(x)dx$$ That is, for any eligible f, the functional J associates f the real number given by the definite integral (which corresponds to the "function value/output" by comparing the action of J with the action of a normal function on a member of its domain) If we let the space of functions be those which are integrable on any interval I on the real axis, we may form a functional with two "arguments": $$J(f,I)=\int_{I}f(x)dx$$ Now, this particular functional is what is known as a LINEAR functional, because for any constants $$\alpha,\beta$$ and functions f,g,h satisfying $$h=\alpha{f}+\beta{g}$$ (i.e, h is a linear combination of f and g), we have: $$J(h,I)=J(\alpha{f}+\beta{g},I)=\int_{I}(\alpha{f}+\beta{g})dx=\alpha\in t_{I}fdx+\beta\int_{I}gdx=\alpha{J}(f,I)+\beta{J}(g,I)$$ by well-known properties of the definite integral. Must all functionals be defined by the aid of integrals, or must at least all LINEAR functionals be defined with integrals? By no means. Here's a counter- example: Let J(f) be defined by: $$J(f)=f(3)$$ That is, every function in our domain is mapped onto its own value at x=3 by the action of J. (Clearly then, functions in the domain of J must be defined at x=3) Now, by letting h be as before, we have: $$J(h)=h(3)=\alpha{f}(3)+\beta{g}(3)=\alpha{J}(f)+\beta{J}(g)$$ Thus, J is linear, but there's not the slightest mention of integrals here!! We are now ready to introduce the Dirac sampling functional: We want a functional D which to any f assigns it the value f(0), provided 0 is within a specified interval I (that's the sampling feature), but which assigns the value 0 to each f if 0 is not in the specified interval. Thus, we define in symbols: $$D(f,I)=f(0),0\in{I},D(f,I)=0,0\notin{I}$$ That's it! Period. This is what the Dirac functional "really is". There's nothing weird about it; it is just a functional among many others, it is trivial to compute its output, and it seems totally unrelated to any features of integrals. It does, however, have in common with the "integral functional " J the property of llinearity, which is easily proven. So, what's all this talk about integrals and how the Dirac functional naturally occurs in conjunction with them? That will be the topic in my next post.
 PF Patron HW Helper Sci Advisor Thanks P: 11,935 While we have seen that there's nothing really strange about the Dirac functional when we consider it simply as a functional, we have yet to see how we may link this object to integrals. In order to do this, we intoduce a set S of weight functions or comparisons functions with respect to some interval [-a,a] on the real axis. Let every member $$\phi$$of S be a non-negative, integrable function with the following basic specification: $$\phi(x)=0,|x|\geq{a},\int_{-a}^{a}\phi{dx}=1, a>0$$ Surely, there are many functions fulfilling this requirement. For later simplicity, we'll add the condition that S only includes those functions $$\phi$$ which are infinitely differentiable. That behaviour might seem a bit difficult to achieve, so here's an actual example: Let $$\phi(x)$$ be defined as: $$\phi(x)=K(a)e^{-\frac{x^{2}}{a^{2}-x^{2}}},|x| PF Patron HW Helper Sci Advisor Thanks P: 11,935 ## The Dirac Delta function If [tex][-a_{n},a_{n}]$$ is outside I for n greater than some N, then we evidently have:
$$D_{n\geq{N}}(f,I)=0\to\lim_{n\to\infty}D_{n}(f,I)=0$$

If the interval is strictly inside of I, we may write:
$$D_{n}(f,I)=\int_{I}f\phi_{n}dx=\int_{-a_{n}}^{a_{n}}f\phi_{n}dx$$

Now, since f was assumed to be continuous, there exist $$a_{m}$$, so that for any given $$\epsilon_{m}>0$$ we have:
$$|f(x)-f(0)|\leq\epsilon_{m},|x|<a_{m}$$
Thus, we have:
$$|D_{m}(f,I)-f(0)|=|\int_{-a_{m}}^{a_{m}}f\phi_{m}dx-f(0)\int_{-a_{m}}^{a_{m}}\phi_{m}dx|=|\int_{-a_{m}}^{a_{m}}(f(x)-f(0))\phi_{m}dx|\leq|=\int_{-a_{m}}^{a_{m}}|f(x)-f(0)|\phi_{m}dx\leq\epsilon_{m}$$
that is, we have $$\lim_{n\to\infty}D_{n}=f(0)$$, by letting the sequence $$\epsilon_{m}$$ go to zero as m goes to infinity.

But this means that we have managed to form a proper INTEGRAL representation of our sampling functional D(f,I)!!!

And, furthermore, since we have $$D_{n}\to{D}$$, we often say instead $$\phi_{n}\to\delta$$ where $$\delta$$ is the commonly "defined" Dirac Delta "function".
It must be emphasized, however, that it is the convergence of FUNCTIONALS which is mathematically meaningful.

In my next post, I'll broach the topic of how to define derivatives of the Dirac Delta function(al), and subtleties concerning choices of the set of comparison functions.
 PF Patron HW Helper Sci Advisor Thanks P: 11,935 Now, the example of a comparison function I posted has the additional nice property that all its derivatives vanish at $$x=\pm{a}$$ Let us now consider a sequence of such functions which has this property: $$\phi^{(m)}(\pm{a})=0, m=0,1,2,3..$$ We will find a useful definition of the first derivative of the functional D, by investigating the behaviour of a particular sequence of functionals, called $$D_{n}^{(1)}(f,I)$$ (We focus our attention upon the behaviour of the sequence on the non-trivial intervals $$[-a_{n},a_{n}]$$): $$D_{n}^{(1)}(f,I)=\int_{-a_{n}}^{a_{n}}f\phi^{(1)}dx=f\phi_{n}\mid_{x=a_{n}}-f\phi_{n}\mid_{x=-a_{n}}+\int_{-a_{n}}^{a_{n}}(-f^{(1)})\phi_{n}dx=\int_{-a_{n}}^{a_{n}}(-f^{(1)})\phi_{n}dx=D_{n}(-f^{(1)},I)$$ where we have assumed that f has a continuous first derivative. Thus, we DEFINE the first derivative of the D to be the functional: $$D^{(1)}(f,I)\equiv{D}(-f^{(1)},I)$$ By extension, we may then define the n'th derivative of D as the functional $$D^{(n)}(f,I))\equiv{D}((-1)^{n}f^{(n)},I)$$ Note that as long as the $$\phi's$$ are AT LEAST n times differentiable, we are therefore able to define in a consistent manner the n'th derivative, provided that the f is n times differentiable. But: Surely, there are lots of weight functions which would, for example give us the functional D, that is, could be regarded as a sequence "converging" to the delta function $$\delta$$, but which wasn't n times differentiable? For example: If we let $$\phi_{n}(x)=\frac{1}{2a_{n}},|x|  PF Patron HW Helper Sci Advisor Thanks P: 11,935 The importance of rate of decay: We have up till now used comparison functions which are strictly zero outside a finite interval. These therefore kill off any "dangerous" behaviour f exhibits outside that interval. However, there are lots of function sequences which might be used to approximate the delta function, which are non-zero on the whole axis, for example, the function sequence [tex]\phi_{n}(x)=\frac{e^{-(\frac{x}{a_{n}})^{2}}}{a_{n}\sqrt{\pi}}$$ will also do, as long as f does not diverge too fast towards infinity. For example, if $$f=Ae^{x^{3}}$$, then $$\int_{-\infty}^{\infty}f\phi_{n}dx$$ is divergent for each choice of n, that is, this particular choice of comparison sequence is wholly unsuitable to be used to give an integral representation of D(f,I) (the sampling functional). Again, it is important to realize that the FUNCTIONAL D(f,I) is perfectly well behaved for that choice of f; we have merely chosen a comparison sequence which won't yield an integral representation of that functional!!! In short then, my choice of S (infinitely differentiable vanishing outside a finite interval) can be regarded as the CORE SET for construction of valid integral representations of the Dirac sampling functional for continuous f's. Other choices might well work for some particular continuous f's, but not for other continuous f's, or those choices are unsuitable for defining derivatives of the functionals involved. Discontinuous f's: Let's look at the following ugly f: $$f(x)=0,x\neq{0},f(0)=1$$ Now, we have as before, equally unproblematically: $$D(f,I)=f(0)=1,0\in{I},D(f,I)=0,x\notin{I}$$ But how should we manage an INTEGRAL representation of D(f,I)??? If we try any of the comparison functions we've used so far, all we'd ever get out is 0. The simplest way of managing this, is to introduce a sequence of $$f_{n}$$-functions converging to f (in some sense of the word): $$f_{n}(x)=1, |x|  PF Patron HW Helper Sci Advisor Thanks P: 11,935 The last example with a discontinuous f shows us how we may form an integral representation of Dirac's sampling function for an extremely wide class of functions. Let's simply assume that f is defined on the whole real axis, and in particular, at the origin. Then the following sequence of functions converges pointwise to f: [tex]f_{n}(x)=f_{n}(0),|x| HW Helper Sci Advisor P: 1,593 Hello Arildno. I'm a few parts back and have the following question: The following relation (from Part 2) seems to be at the heart of the theory and I wish to understand it fully. You say:  If for any $\epsilon_m> 0$, we can find an $a_m$ such that: [tex]|f(x)-f(0)|\leq\epsilon_{m},|x| I wonder if there is a minor typo in that last relation: a stray absolute value sign and equal sign appearing right after the $\leq$ sign. Should it read: [tex]|D_{m}(f,I)-f(0)|=|\int_{-a_{m}}^{a_{m}}f\phi_{m}dx-f(0)\int_{-a_{m}}^{a_{m}}\phi_{m}dx|=|\int_{-a_{m}}^{a_{m}}(f(x)-f(0))\phi_{m}dx|\leq\int_{-a_{m}}^{a_{m}}|f(x)-f(0)|\phi_{m}dx\leq\epsilon_{m}$$

To me, that relation really shows why:

$$\lim_{n\to\infty}D_{n}[f,I]=f(0)$$

and I want to make sure I have it understood correctly. Can you confirm this suspicion?

Thanks,
Salty
 PF Patron HW Helper Sci Advisor Thanks P: 11,935 You're absolutely right; they shouldn't be there..
 PF Patron HW Helper Sci Advisor Thanks P: 11,935 While the technique I espoused in part 5 for the formation of an integral representation of the sampling functional is okay enough, I guess, coupling the limiting processes of $$\phi_{n},D_{n},f_{n}$$ isn't really necessary, and rather inelegant. Thus, I would prefer the following: Let [tex]f_{m}(x)=f(0),|x|
 P: 395 This thread was very helpful to me Arildno, thanks for taking the time to write it.
 P: 1 On the MIT OCW (Open CourseWare) the course on Differential Equations 18.03 includes Laplace transforms and a very clear* explanation of the Dirac delta function in Lecture #22. *It's clear to me and I am mathematically naive at the college level (just learning DiffEq and other stuff now.) Anyone not familiar with Laplace transforms should go back to Lecture #19 and proceed from 19-22. -- Herb Martin
 P: 56 When solving problems involving the dirac dela function, can you recommend any sort of method or steps to follow to evaluate intergrals involving the dirac delta function.
PF Patron
HW Helper