Dirac Delta Function

A. Dirac Delta Function

[Let us introduce a quantity 𝛿(x), which depends on x satisfying the following conditions:] with [This quantity is called the Dirac Delta function and is quite useful in physics. We can show that]

Proof: If f(x) varies relatively slowly in the range to we can replace it by an average of f(x) over that interval. Letting , we get

The delta function is useful since it may be operated on as though it were a real function. The only justification we give for this is that it gives the correct answer. The following useful relations may be easily proved:

The following figures are approximate representations of and .

A diagram shows a sharp peak for δ(x) at x=0
A diagram shows a sharp positive peak followed by a sharp negative peak for δ'(x) at x=0

Problem: (with integrals)

  1. Prove: for
  2. Find:

B. Step Function

Explanation

The function being described on the right-hand side is the Heaviside step function, often written as H(x) or θ(x)

  • Case 1 (x<0) The integral’s upper limit  x is negative. The integration range is from −∞ to x. This range does not include t=0 where the delta functon is non-zero. Therefore, you are integrating a function that is zero over the entire interval, and the result is 0.
  • Case 2 (x>0): The integral’s upper limit x is positive. The integration range from  −∞ to xdoes include t=0. Because the entire “spike” of the delta function is included in the range, the integral evaluates to the total area of the delta function, which is 1.

Notice that .

Consider the curve in the following figure. Its derivative may be expressed in terms of the delta function as follows:

A diagram shows a function f(x) with a step discontinuity of height c at x=a
Explanation

Let’s take any function f(x) that is smooth everywhere except for a single jump discontinuity of height c at x=a.

We can represent this function f(x) as the sum of two parts:

  1. A continuous smooth “base” function, let’s call it g(x).
  2. A step function that adds the jump at the right place.

The formula is:

How do we find g(x)? The function g(x) is simply the “pre-jump” part of the function extended across the entire domain.

  • For x < a, we have H(x-a)=0, so f(x) = g(x). This means g(x) is identical to f(x) before the jump.
  • For x > a, we have H(x-a)=1, so f(x) = g(x) + c. This correctly describes the function after the jump.

Differentiating the Representation

Now that we have expressed f(x) in a form that is a sum of well-behaved functions (a continuous function and a Heaviside function), we can differentiate it using the sum rule:

Using the relationship , we get:

Now, what is g’(x)? Since g(x) is the continuous version of f(x), its derivative g’(x) is exactly the same as the derivative of f(x) away from the jump, i.e., for .

So, we can replace g’(x) with the notation f’(x)_{x a} to get the final formula:

C. Delta function and the integral of cos(ßx)

Consider the integral: . We shall show that it is equal to .

Explanation

Notice that cos(ßx) dx does not converge in the usual sense because the cosine function oscillates forever. To give it a meaningful value, we treat it as a distribution by examining its behavior in a limit.

Explanation

It was previously shown that

For small this integral looks approximately as in the following figure.

A diagram shows a bell-shaped curve centered at ß = 0 with height 1/a

The area under the curve is: Note that the area is independent of the value of a.

Explanation

The function has an area of π for any value of a. As a→0, the function becomes an infinitely thin, infinitely high spike at β=0, while its total area remains π. This is precisely the behavior of πδ(β).

D. Use of the delta function in evaluating a definite integral

Example: Consider Substitute The general solution of such a differential equation is: where

Explanation

From differential equations, recall that if and are two solutions the equation , then a particular solution of is given by where is the Wronskian.

For the equation , we look for solutions of the form . The characteristic equation is , which has roots . This gives us two independent solutions of the corresponding homogeneous equation: The homogeneous solution is The Wronskian is According to the theorem mentioned above, a particular solution is given by and the total solution is

Solutions:

Explanation

Previously, we obtained If : and the total solution is If then and the total solution is

Example: Let us consider the integral [If we differentiate with respect to , we get the where we know is :] [Solving this differential equation, we get:] We note that (A graph shows a step function, going from -π for β<0 to +π for β>0)

Consider now the integral

Explanation

Notice that we previously studied this integral.

A graph of this integral looks as follows:

Graph of pi/2-arctan(alpha/beta)

In the lim , this integral reduces to the one considered previous to it (there is a factor 2 to be considered)


Reference: Principals of Quantum Mechanics, Paul Dirac, 4th ed., p. 58–61