Bessel Functions
Bessel Functions
Bessel Functions
Jennifer Niedziela
+ xy
+ (x
2
2
)y = 0 (1)
By re-writing this equation as:
x(xy
+ (x
2
2
)y = 0 (2)
and employing the use of a generalized power series, we
re-write the terms of (2) in terms of the series:
y =
n=0
a
n
x
n+s
y
n=0
a
n
(n + s)x
n+s1
xy
n=0
a
n
(n + s)x
n+s
(xy
n=0
a
n
(n + s)
2
x
n+s1
(xy
n=0
a
n
(n + s)
2
x
n+s
When the coecients of the powers of x are organized, we
nd that the coecient on x
s
gives the indicial equation
s
2
2
= 0, = s = , and we develop the general
formula for the coecient on the x
s+n
term:
a
n
=
a
n2
(n + s)
2
2
(3)
In the case s = :
a
n
=
a
n2
n(n + 2)
(4)
and since a
1
= 0, a
n
= 0 for all n = odd integers. Coe-
cients for even powers of n are found:
a
2n
=
a
2n2
2
2
n(n + )
(5)
Recalling that for the gamma function:
( + 2) = ( + 1)( + 1),
( + 3) = ( + 2)( + 2) = ( + 2)( + 1)( + 1),
we can write the coecients:
a
2
=
a
0
2
2
(1 + )
=
(1 + )
2
2
(2 + )
a
2n
=
a
0
(1 + )
n!2
2n
(n + 1 + )
Which allows us to write the terms of the series :
y = J
(x) =
n=0
(1)
n
(n + 1)(n + + 1)
_
x
2
_
2n+
(6)
Where J
(x) =
cosJ
(x) J
(x)
sin
(7)
For integral values of , the expression of N
(x) has an
indeterminate form, and N
(x)|
x=0
= . Nevertheless
the limit of this function for x = 0, the expression for N
(x) + BN
(x) (8)
with A and B as arbitrary constants determined from
boundary conditions.
Bessel functions of the rst and second kind are the
most commonly found forms of the Bessel function in ap-
plications. Many applications in hydrodynamics, elastic-
ity, and oscillatory systems have solutions that are based
on the Bessel functions. One such example is that of a
uniform density chain xed at one end undergoing small
oscillations. The dierential equation of this situation is:
d
2
u
dz
2
+
1
z
du
dz
+
k
2
u
z
= 0 (9)
where z references a point on the chain, k
2
=
p
2
g
, with p
as the frequency of small oscillations at that point, and
g the gravitational constant of acceleration. Eq. (9) is a
form of eq. (1), and solution is:
u = AJ
0
(2kz
1
2
) + BY
0
(2kz
1
2
), (10)
where the A and B are determined by the boundary con-
ditions.
Modied Bessel Functions
Modied Bessel functions are found as solutions to the
modied Bessel equation
x
2
y
+ xy
(x
2
2
)y = 0 (11)
which transforms into eq. (1) when x is replaced with
ix. However, this leaves the general solution of eq. (1)
a complex function of x. To avoid dealing with complex
solutions in practical applications [2], the solutions to
(11) are expressed in the form:
I
(x) = e
i
2
J
(xe
i
2
) (12)
The I
(x):
y = AI
(x) BI
(x) (13)
where again A and B are determined from the boundary
conditions.
A solution for non-integer orders of is found:
K
(x) =
2
I
(x) I
(x)
sin
(14)
The functions K
(x))is
shown in gure (2). A plot of the Modied Second Kind
functions (K
n
(x)) is shown in g. (3).
Modied Bessel functions appear less frequently in ap-
plications, but can be found in transmission line studies,
non-uniform beams, and the statistical treatment of a
relativistic gas in statistical mechanics.
Zeroes of Bessel Functions
The zeroes of Bessel functions are of great importance
in applications [5]. The zeroes, or roots, of the Bessel
functions are the values of x where value of the Bessel
function goes to zero (J
(x))
functions do not have zeroes.
Bessel function zeros are exploited in frequency modu-
lated (FM) radio transmission. FM transmission is math-
ematically represented by a harmonic distribution of a
sine wave carrier modulated by a sine wave signal which
can be represented with Bessel Functions. The carrier
or sideband frequencies disappear when the modulation
index (the peak frequency deviation divided by the mod-
ulation frequency) is equal to the zero crossing of the
function for the n
th
sideband.
APPLICATION - SOLUTION TO
SCHROEDINGERS EQUATION IN A
CYLINDRICAL WELL
Consider a particle of mass m placed into a two-
dimensional potential well, where the potential is zero
inside of the radius of the disk, innite outside of the
radius of the disk. In polar coordinates using r, as
representatives of the system, the Laplacian is written:
2
=
1
r
r
_
r
r
_
+
1
r
2
2
. (15)
Which in the Schroedinger equation presents:
h
2
2m
_
1
r
r
_
r
r
_
+
1
r
2
2
_
= E. (16)
Using the method of separation of variables with a pro-
posed solution = R(r)T() in (16), produces
h
2
2m
_
T()
1
R(r)
r
_
r
R(r)
r
_
+
1
r
2
R(r)
2
T
2
_
= ER(r)T()
(17)
and then dividing by :
_
1
R
1
r
r
_
r
R
r
_
+
1
r
2
1
T
2
T
2
_
=
2mE
h
2
(18)
setting
2mE
h
2
= k
2
and multiplying through by r
2
pro-
duces
r
R
d
dr
_
r
dR
dr
_
+ k
2
r
2
+
1
T
d
2
T
d
2
= 0 (19)
which is fully separated in r and . To solve, the
dependent portion is set to m
2
, yielding the harmonic
oscillator equation in T(), which presents the solution:
T() = Ae
im
(20)
Where A is a constant determined via proper normaliza-
tion in :
_
2
0
A
2
T()T() d = 1 =A =
_
1
2
(21)
Leaving the dependent portion T() =
_
1
2
e
im
Working now with the r dependent portion of the sep-
arated equation, multiplying the r dependent portion of
(19) by r
2
, and setting equal to m
2
one obtains:
r
2
R
d
2
R
dr
2
+
r
R
dR
dr
+ k
2
r
2
= m
2
(22)
4
which when rearranged:
r
2
d
2
R
dr
2
+ r
dR
dr
+ (k
2
r
2
m
2
) = 0 (23)
which is of the same form as eq. (1), Bessels dierential
equation. The general solution to eq. (23) is of the form
of eq. (8), and we write that general solution:
R(r) = AJ
m
(kr) + BN
m
(kr), (24)
where J
m
(kr) and N
m
(kr) are respectively the Bessel
and Neumann functions of order m, and A and B are
constants to be determined via application of the bound-
ary conditions. As the solution must be nite at x = 0,
and as N
m
(kr) as x 0, this means that the coef-
cient of N
m
(kr) = B = 0, leaving R(r) to be expressed:
R(r) = AJ
m
(kr) (25)
Using the boundary condition that = 0 at the radius of
the disk, we have the condition that J
m
(kr
b
) = 0, which
implicitly requires the argument of J
m
to be a zero of the
Bessel function. As noted earlier, these zeroes must be
calculated individually in numerical fashion. Requiring
that kr
b
=
m,n
, which is the nth zero of the m
th
or-
der Bessel function[6], the energy of the system is solved
by expressing k in terms of
m,n
in (eq. k =
_
2mE
h
2
),
arriving at:
E
m,n
=
2
m,n
h
2
2mr
2
b
h (26)
The full solution for is thus:
m
(r, ) = AJ
m
(
m,n
r
r
b
)e
im
(27)
Because the Bessel function zeroes cannot be determined
apriori, it is dicult to nd a closed solution to express
the normalization constant A. We select an order for m
to continue with the determination of the normalization
constant, and arbitrarily choose m = 2, which has a zero
at r = 5.13562, which we will set to be the radius of the
circle. Given the preceding, the normalization can be for
the m = 2, n = 1 case can be found:
_
r
boundary
0
A
2
J
2
(
2,1
r
r
b
)J
2
(
2,1
r
r
b
) dr = 1 (28)
which for r
boundary
= 5.13562 = A =
_
1
0.510377
, (nu-
merical values obtained using numerical integration in
Mathematica). Thus, we can express the full solution for
the m = 2 scenario:
(r, ) =
_
1
0.510377
_
1
2
J
2
(
2,1
r
r
b
)e
im
(29)
And since weve eectively set
m,n
= r
b
,
(r, ) =
_
1
0.510377
_
1
2
J
2
(r)e
im
(30)
Admittedly, this solution is somewhat contrived, but it
shows the importance of working with the zeroes of the
Bessel function to generate the particular solution using
the boundary conditions.
CONCLUSION
The Bessel functions appear in many diverse scenar-
ios, particularly situations involving cylindrical symme-
try. The most dicult aspect of working with the Bessel
function is rst determining that they can be applied
through reduction of the system equation to Bessels dif-
ferential or modied equation, and then manipulating
boundary conditions with appropriate application of ze-
roes, and the coecient values on the argument of the
Bessel function. This topic can be greatly expanded
upon, and the reader is highly encouraged to review the
applications and development presented in [2].
REFERENCES
jniedziela@ornl.gov
[1] J. J. O Connor and R. E. F., Friedrich Wilhelm Bessel
(School of Mathematics and Statistics University of St An-
drews Scotland, 1997).
[2] F. E. Relton, Applied Bessel Functions (Blackie and Son
Limited, 1946).
[3] H. J. Arfken, G. B., Weber, Mathematical Methods for
Physicists (Elsevier Academic Press, 2005).
[4] E. W. Weissten, Modied Bessel Function of the Second
Kind (Eric Weissteins World of Physics, 2008).
[5] M. Boas, Mathematical Methods for the Physical Sciences
(Wiley, 1983).
[6] P. F. Newhouse and K. C. McGill, Journal of Chemical
Education 81, 424 (2004).