9 Bessel Functions
9.1 Bessel functions of the first kind Friedrich Bessel (1784–1846) invented functions for problems with circular symmetry. The most useful ones are defined for any integer n by the series zn z2 z4 J (z)= 1 + ... n 2nn! 2(2n + 2) 2 4(2n + 2)(2n + 4) · (9.1) z n 1 ( 1)m z 2m = . 2 m!(m + n)! 2 m=0 ⇣ ⌘ X ⇣ ⌘ The first term of this series tells us that for small z 1 | |⌧ zn J (z) . (9.2) n ⇡ 2nn! The alternating signs in (9.1) make the waves plotted in Fig. 9.1, and we have for big z 1 the approximation (Courant and Hilbert, 1955, chap. VII) | | 2 n⇡ ⇡ 3/2 J (z) cos z + O( z ). (9.3) n ⇡ ⇡z 2 4 | | r ⇣ ⌘ The Jn(z) are entire transcendental functions. They obey Bessel’s equation d2J 1 dJ n2 n + n + 1 J = 0 (9.4) dz2 z dz z2 n ✓ ◆ (6.393) as one may show (exercise 9.1) by substituting the series (9.1) into the di↵erential equation (9.4). Their generating function is
z 1 exp (u 1/u) = un J (z) (9.5) 2 n n= h i X 1 9.1 Bessel functions of the first kind 383
The Bessel Functions J0(ρ), J1(ρ), and J2(ρ) 1
0.5
0
−0.5 0 2 4 6 8 10 12 ρ
The Bessel Functions J3(ρ), J4(ρ), and J5(ρ)
0.4
0.2
0
−0.2
−0.4 0 2 4 6 8 10 12 ρ
Figure 9.1 Top: Plots of J0(⇢) (solid curve), J1(⇢) (dot-dash), and J2(⇢) (dashed) for real ⇢. Bottom: Plots of J3(⇢) (solid curve), J4(⇢) (dot-dash), and J5(⇢)(dashed). The points at which Bessel functions cross the ⇢-axis are called zeros or roots; we use them to satisfy boundary conditions. from which one may derive (exercise 9.5) the series expansion (9.1) and (exercise 9.6) the integral representation (5.51)
⇡ 1 n Jn(z)= cos(z sin ✓ n✓) d✓ = J n( z)=( 1) J n(z) (9.6) ⇡ Z0 for all complex z. For n = 0, this integral is (exercise 9.7) more simply
1 2⇡ 1 2⇡ J (z)= eiz cos ✓ d✓ = eiz sin ✓ d✓. (9.7) 0 2⇡ 2⇡ Z0 Z0 These integrals (exercise 9.8) give J (0) = 0 for n = 0, and J (0) = 1. n 6 0 By di↵erentiating the generating function (9.5) with respect to u and 384 Bessel Functions identifying the coe cients of powers of u, one finds the recursion relation 2n Jn 1(z)+Jn+1(z)= Jn(z). (9.8) z Similar reasoning after taking the z derivative gives (exercise 9.10)
Jn 1(z) Jn+1(z)=2Jn0 (z). (9.9) By using the gamma function (section 5.15), one may extend Bessel’s equation (9.4) and its solutions Jn(z) to non-integral values of n
z ⌫ 1 ( 1)m z 2m J (z)= . (9.10) ⌫ 2 m! (m + ⌫ + 1) 2 m=0 ⇣ ⌘ X ⇣ ⌘ The self-adjoint form (6.392) of Bessel’s equation (9.4) with z = kx is (exercise 9.11) d d n2 x J (kx) + J (kx)=k2xJ (kx). (9.11) dx dx n x n n ✓ ◆ In the notation of equation (6.365), p(x)=x, k2 is an eigenvalue, and ⇢(x)=x is a weight function. To have a self-adjoint system (section 6.30) on an interval [0,b], we need the boundary condition (6.325) b b 0= p(J v0 J 0 v) = x(J (kx)v0(kx) J 0 (kx)v(kx)) (9.12) n n 0 n n 0 for all functions⇥ v(kx) in⇤ the⇥ domain D of the system. The definition⇤ (9.1) of Jn(z) implies that J0(0) = 1, and that Jn(0) = 0 for integers n>0. Thus since p(x)=x, the terms in this boundary condition vanish at x = 0 as long as the domain consists of functions v(kx) that are twice di↵erentiable on the interval [0,b]. To make these terms vanish at x = b, we require that Jn(kb) = 0 and that v(kb) = 0. Thus kb must be a zero zn,m of Jn(z), and so Jn(kb)=Jn(zn,m) = 0. With k = zn,m/b, Bessel’s equation (9.11) is d d n2 z2 x J (z x/b) + J (z x/b)= n,m xJ (z x/b) . (9.13) dx dx n n,m x n n,m b2 n n,m ✓ ◆ 2 2 2 For fixed n, the eigenvalue k = zn,m/b is di↵erent for each positive integer m. Moreover as m , the zeros z of J (x) rise as m⇡ as one might !1 n,m n expect since the leading term of the asymptotic form (9.3) of Jn(x) is pro- portional to cos(x n⇡/2 ⇡/4) which has zeros at m⇡+(n+1)⇡/2+⇡/4. It follows that the eigenvalues k2 (m⇡)2/b2 increase without limit as m ⇡ !1 in accordance with the general result of section 6.36. It follows then from the argument of section 6.37 and from the orthogonality relation (6.404) that for every fixed n, the eigenfunctions Jn(zn,mx/b), one for each zero, are 9.1 Bessel functions of the first kind 385 complete in the mean, orthogonal, and normalizable on the interval [0,b] with weight function ⇢(x)=x
b 2 2 zn,mx zn,m0 x b 2 b 2 xJn Jn dx = m,m0 Jn0 (zn,m)= m,m0 Jn+1(zn,m) 0 b b 2 2 Z ⇣ ⌘ ⇣ ⌘ (9.14) and a normalization constant (exercise 9.12) that depends upon the first derivative of the Bessel function or the square of the next Bessel function at the zero. Because they are complete, sums of Bessel functions Jn(zn,kx/b) can rep- resent Dirac’s delta function on the interval [0,b] as in the sum (6.451)
↵ 1 ↵ 2 1 Jn(zn,kx/b)Jn(zn,ky/b) (x y)= 2 x y /b . (9.15) J 2 (z ) k=1 n+1 n,k X For n = 0 and b = 1, this sum is
↵ 1 ↵ 1 J0(z0,kx)J0(z0,ky) (x y)=2x y 2 . (9.16) J1 (z0,k) Xk=1 Figure 9.2 plots the sum of the first 10,000 terms of this series for ↵ =0 and y =0.47, for ↵ =1/2 and y =1/2, and for ↵ = 1 and y =0.53. The integrals of these 10,000-term sums from 0 to 1 respectively are 0.9966, 0.9999, and 0.9999. The expansion (9.15) of the delta function expresses the completeness of the eigenfunctions Jn(zn,kx/b) for fixed n and lets us write any twice di↵erentiable function f(x) that vanishes at x = b as the Fourier series
↵ 1 f(x)=x ak Jn(zn,kx/b) (9.17) Xk=1 where b 2 1 ↵ a = y J (z y/b) f(y) dy. (9.18) k b2 J 2 (z ) n n,k n+1 n,k Z0 The orthogonality relation on an infinite interval is