Appendix a Dirac Delta Function
Total Page:16
File Type:pdf, Size:1020Kb
Appendix A Dirac Delta Function In 1880 the self-taught electrical scientist Oliver Heaviside introduced the following function 1forx > 0 (x) = (A.1) 0forx < 0 which is now called Heaviside step function. This is a discontinous function, with a discontinuity of first kind (jump) at x = 0, which is often used in the context of the analysis of electric signals. Moreover, it is important to stress that the Haviside step function appears also in the context of quantum statistical physics. In fact, the Fermi-Dirac function (or Fermi-Dirac distribution) 1 F (x) = , (A.2) β eβ x + 1 proposed in 1926 by Enrico Fermi and Paul Dirac to describe the quantum statistical distribution of electrons in metals, where β = 1/(kB T ) is the inverse of the absolute temperature T (with kB the Boltzmann constant) and x = − μ is the energy of the electron with respect to the chemical potential μ, becomes the function (−x) in the limit of very small temperature T , namely 0forx > 0 lim Fβ(x) = (−x) = . (A.3) β→+∞ 1forx < 0 Inspired by the work of Heaviside, with the purpose of describing an extremely localized charge density, in 1930 Paul Dirac investigated the following “function” +∞ for x = 0 δ(x) = (A.4) 0forx = 0 imposing that +∞ δ(x) dx = 1 . (A.5) −∞ © Springer International Publishing AG 2017 229 L. Salasnich, Quantum Physics of Light and Matter, UNITEXT for Physics, DOI 10.1007/978-3-319-52998-1 230 Appendix A: Dirac Delta Function Unfortunately, this property of δ(x) is not compatible with the definition (A.4). In fact, from Eq. (A.4) it follows that the integral must be equal to zero. In other words, it does not exist a function δ(x) which satisfies both Eqs. (A.4) and (A.5). Dirac suggested that a way to circumvent this problem is to interpret the integral of Eq. (A.5)as +∞ +∞ δ(x) dx = lim δ(x) dx , (A.6) −∞ →0+ −∞ where δ(x) is a generic function of both x and such that +∞ for x = 0 lim δ(x) = , (A.7) →0+ 0forx = 0 +∞ δ(x) dx = 1 . (A.8) −∞ Thus, the Dirac delta function δ(x) is a “generalized function” (but, strictly-speaking, not a function) which satisfy Eqs. (A.4) and (A.5) with the caveat that the integral in Eq. (A.5) must be interpreted according to Eq. (A.6) where the functions δ(x) satisfy Eqs. (A.7) and (A.8). There are infinite functions δ(x) which satisfy Eqs. (A.7) and (A.8). Among them there is, for instance, the following Gaussian 1 2 2 δ (x) = √ e−x / , (A.9) π which clearly satisfies Eq. (A.7) and whose integral is equal to 1 for any value of . Another example is the function 1 for |x|≤/2 δ (x) = , (A.10) 0for|x| > /2 which again satisfies Eq. (A.7) and whose integral is equal to 1 for any value of . In the following we shall use Eq. (A.10) to study the properties of the Dirac delta function. According to the approach of Dirac, the integral involving δ(x) must be interpreted as the limit of the corresponding integral involving δ(x), namely +∞ +∞ δ(x) f (x) dx = lim δ(x) f (x) dx , (A.11) −∞ →0+ −∞ for any function f (x). It is then easy to prove that +∞ δ(x) f (x) dx = f (0). (A.12) −∞ Appendix A: Dirac Delta Function 231 by using Eq. (A.10) and the mean value theorem. Similarly one finds +∞ δ(x − c) f (x) dx = f (c). (A.13) −∞ Several other properties of the Dirac delta function δ(x) follow from its definition. In particular δ(−x) = δ(x), (A.14) 1 δ(ax) = δ(x) with a = 0 , (A.15) |a| 1 δ( f (x)) = δ(x − x ) with f (x ) = 0 . (A.16) | f (x )| i i i i Up to now we have considered the Dirac delta function δ(x) with only one variable x. It is not difficult to define a Dirac delta function δ(D)(r) in the case of a D- D D dimensional domain R , where r = (x1, x2, ..., xD) ∈ R is a D-dimensional vector: ( ) +∞ for r = 0 δ D (r) = (A.17) 0forr = 0 and δ(D)(r) d Dr = 1 . (A.18) RD Notice that sometimes δ(D)(r) is written using the simpler notation δ(r). Clearly, also in this case one must interpret the integral of Eq. (A.18)as δ(D)(r) d Dr = lim δ(D)(r) d Dr , (A.19) + RD →0 RD (D)( ) where δ r is a generic function of both r and such that ( ) +∞ for r = 0 lim δ D (r) = , (A.20) → + 0forr = 0 0 (D)( ) D = . lim δ r d r 1 (A.21) →0+ Several properties of δ(x) remain valid also for δ(D)(r). Nevertheless, some prop- erties of δ(D)(r) depend on the space dimension D. For instance, one can prove the remarkable formula 1 ∇2 (ln |r|) for D = 2 (D)( ) = 2π , δ r − 1 ∇2 1 ≥ (A.22) D−2 for D 3 D(D−2)VD |r| 232 Appendix A: Dirac Delta Function ∇2 = ∂2 + ∂2 +···+ ∂2 = D/2/( + / ) where 2 2 2 and VD π 1 D 2 is the volume of ∂x1 ∂x2 ∂xD a D-dimensional ipersphere of unitary radius, with (x) the Euler Gamma function. In the case D = 3 the previous formula becomes 1 1 δ(3)(r) =− ∇2 , (A.23) 4π |r| which can be used to transform the Gauss law of electromagnetism from its integral form to its differential form. Appendix B Fourier Transform It was known from the times of Archimedes that, in some cases, the infinite sum of decreasing numbers can produce a finite result. But it was only in 1593 that the mathematician Francois Viete gave the first example of a function, f (x) = 1/(1 − x), written as the infinite sum of power functions. This function is nothing else than the geometric series, given by ∞ 1 = xn , for |x| < 1 . (B.1) 1 − x n=0 In 1714 Brook Taylor suggested that any real function f (x) which is infinitely differentiable in x0 and sufficiently regular can be written as a series of powers, i.e. ∞ n f (x) = cn (x − x0) , (B.2) n=0 where the coefficients cn are given by 1 c = f (n)(x ), (B.3) n n! 0 with f (n)(x) the n-th derivative of the function f (x). The series (B.2) is now called Taylor series and becomes the so-called Maclaurin series if x0 = 0. Clearly, the geometric series (B.1) is nothing else than the Maclaurin series, where cn = 1. We observe that it is quite easy to prove the Taylor series: it is sufficient to suppose that Eq. (B.2) is valid and then to derive the coefficients cn by calculating the derivatives of f (x) at x = x0; in this way one gets Eq. (B.3). In 1807 Jean Baptiste Joseph Fourier, who was interested on wave propagation and periodic phenomena, found that any sufficiently regular real function function f (x) which is periodic, i.e. such that © Springer International Publishing AG 2017 233 L. Salasnich, Quantum Physics of Light and Matter, UNITEXT for Physics, DOI 10.1007/978-3-319-52998-1 234 Appendix B: Fourier Transform f (x + L) = f (x), (B.4) where L is the periodicity, can be written as the infinite sum of sinusoidal functions, namely ∞ a 2π 2π f (x) = 0 + a cos n x + b sin n x , (B.5) 2 n L n L n=1 where 2 L/2 2π an = f (y) cos n y dy , (B.6) L − / L L 2 2 L/2 2π bn = f (y) sin n y dy . (B.7) L −L/2 L It is quite easy to prove also the series (B.5), which is now called Fourier series. In fact, it is sufficient to suppose that Eq. (B.5) is valid and then to derive the coeffi- 2π 2π cients an and bn by multiplying both side of Eq. (B.5) by cos n L x and cos n L x respectively and integrating over one period L; in this way one gets Eqs. (B.6) and (B.7). It is important to stress that, in general, the real variable x of the function f (x) can represent a space coordinate but also a time coordinate. In the former case L gives the spatial periodicity and 2π/L is the wavenumber, while in the latter case L is the time periodicity and 2π/L the angular frequency. Taking into account the Euler formula in 2π x 2π 2π e L = cos n x + i sin n x (B.8) L L √ with i = −1 the imaginary unit, Fourier observed that his series (B.5) can be re-written in the very elegant form +∞ in 2π x f (x) = fn e L , (B.9) n=−∞ where L/2 1 −in 2π y fn = f (y) e L dy (B.10) L −L/2 are complex coefficients, with f0 = a0/2, fn = (an − ibn)/2ifn > 0 and fn = ( + )/ < ∗ = a−n ib−n 2ifn 0, thus fn f−n. The complex representation (B.9) suggests that the function f (x) can be periodic but complex, i.e. such that f : R → C. Moreover, one can consider the limit L → +∞ of infinite periodicity, i.e. a function which is not periodic.