Fundamentals of Bessel Functions

Bessel functions are canonical solutions to the second-order ordinary differential equation known as Bessel’s equation. Named after the German astronomer and mathematician Friedrich Bessel (1784–1846), these functions arise naturally in problems possessing cylindrical or spherical symmetry. Their role in electromagnetic wave propagation is fundamental: from the analysis of circular waveguides and optical fibers to the computation of electromagnetic scattering by cylinders and spheres, Bessel functions provide the mathematical backbone for describing wave behavior in curvilinear coordinates.

The Bessel Differential Equation

The standard form of Bessel’s differential equation is:

x² y'' + x y' + (x² – ν²) y = 0

where ν is a real or complex parameter called the order of the equation. The solutions to this equation are the Bessel functions of the first kind, Jν(x), and the Bessel functions of the second kind, Yν(x) (also known as Neumann functions or Weber functions). For integer orders, Jn(x) and Yn(x) are the most commonly used. A comprehensive reference for these functions and their properties can be found at Wolfram MathWorld.

Modified Bessel Functions

When the argument of the Bessel equation becomes imaginary (i.e., the equation is of the form x² y'' + x y' – (x² + ν²) y = 0), the solutions are the modified Bessel functions of the first kind, Iν(x), and of the second kind, Kν(x). These are essential in analyzing evanescent fields and cutoff behavior in waveguides, as they describe fields that decay exponentially rather than oscillate.

Wave Equations in Cylindrical Coordinates

Maxwell’s equations in a source-free, homogeneous, isotropic medium can be combined to yield the vector wave equation for the electric and magnetic fields. In cylindrical coordinates (ρ, φ, z), the wave equation separates into three ordinary differential equations via the method of separation of variables. Assuming a time-harmonic dependence ejωt and propagation in the z-direction, the fields can be expressed as:

E(ρ, φ, z) = R(ρ) Φ(φ) e–jβz

Substituting into the wave equation yields a radial equation that is exactly Bessel’s differential equation. Specifically, the radial function R(ρ) satisfies:

ρ² R'' + ρ R' + (kρ² ρ² – n²) R = 0

where kρ² = ω²με – β² and n is the azimuthal mode number (an integer). The solutions are Jn(kρρ) and Yn(kρρ) for traveling waves, or In(kρρ) and Kn(kρρ) when the transverse wavenumber becomes imaginary.

Separation of Variables Procedure

To illustrate, consider the scalar wave equation for the z-component of the electric field in cylindrical coordinates:

∇²Ez + k²Ez = 0 where k = ω√(με).

Using separation of variables, we assume Ez(ρ, φ, z) = R(ρ) Φ(φ) Z(z). This substitution leads to ordinary differential equations for each coordinate. The azimuthal equation yields sinusoidal functions, and the longitudinal equation gives exponential propagation factors. The radial equation reduces to Bessel’s equation with parameter ν = n (the separation constant from the azimuthal equation). The general solution for the radial part is thus a linear combination of Bessel functions of the first and second kinds.

Application in Circular Waveguides

Circular waveguides are a classic application of Bessel functions in electromagnetics. In a metallic cylindrical waveguide with radius a, the boundary condition at the perfectly conducting wall requires the tangential electric field to vanish. For TE (transverse electric) modes, this translates into the condition that the derivative of the Bessel function of the first kind, Jn' (kρ a) = 0. For TM (transverse magnetic) modes, the condition is Jn(kρ a) = 0. The roots of these equations determine the cutoff wavenumbers and thus the propagation constants of the waveguide modes.

Field Components and Mode Patterns

For TEnm and TMnm modes, the field components are expressed directly in terms of Jn(kρρ) and its derivative. The integer n denotes the number of full-period variations in the azimuthal direction, while m indicates the m-th root of the Bessel function or its derivative, corresponding to the number of radial variations. Detailed tables of these roots are widely available; an authoritative source is the NIST Digital Library of Mathematical Functions.

For example, the dominant TE11 mode in a circular waveguide has a cutoff frequency determined by the first non-zero root of J1'(x) = 0. The radial electric field is proportional to J1(kρρ), and the pattern exhibits a single variation in the radial direction. These modes are critical in the design of microwave components, antennas, and feedhorns.

Cutoff Frequencies and Dispersion

Applying the boundary condition gives the eigenvalues kρ mn = pnm/a for TM modes and kρ mn = p'nm/a for TE modes, where pnm are the zeros of Jn(x) and p'nm are the zeros of Jn'(x). The cutoff frequency is fc = (kρ mn) / (2π√(με)). Below cutoff, the propagation constant becomes imaginary, leading to evanescent fields described by modified Bessel functions In and Kn. This behavior is extensively used in waveguide filters and directional couplers.

Hollow Waveguides vs. Dielectric Waveguides

In dielectric waveguides such as optical fibers, the fields are not confined by a perfect conductor but by total internal reflection. The radial equation in the core (where the refractive index is higher) yields oscillatory Bessel functions of the first kind, while in the cladding the solutions are modified Bessel functions that decay exponentially. The matching of fields at the core-cladding interface leads to the characteristic equation that determines the modal propagation constants. Bessel functions of the first kind Jn(uρ/a) appear in the core, and modified Bessel functions of the second kind Kn(wρ/a) appear in the cladding, where u and w are normalized transverse wavenumbers.

Spherical Bessel Functions and Spherical Harmonics

When the problem exhibits spherical symmetry, the Helmholtz equation is separated in spherical coordinates (r, θ, φ). The radial equation becomes a spherical Bessel differential equation:

r² R'' + 2r R' + [k²r² – l(l+1)] R = 0

This equation can be transformed into the ordinary Bessel equation by the substitution R(r) = (1/√r) S(r), leading to solutions involving half-integer-order Bessel functions. These are the spherical Bessel functions jl(kr), yl(kr), and their combinations hl(1) and hl(2) (spherical Hankel functions). Spherical Bessel functions describe the radial dependence of electromagnetic fields in resonant cavities, scattering by spheres, and radiation from spherical antennas.

Electromagnetic Waves in Spherical Cavities

In spherical cavities, the boundary conditions at the conducting surface r = a lead to transcendental equations involving spherical Bessel functions. For TM modes, the condition is jl(ka) = 0; for TE modes, it is (d/dr)[r jl(kr)] evaluated at r = a equals zero. The lowest order modes are used in applications ranging from satellite communication filters to particle accelerators. The mathematical connection between spherical Bessel functions and ordinary Bessel functions of half-integer order is given by jl(x) = √(π/(2x)) Jl+1/2(x).

Advanced Applications in Electromagnetic Wave Propagation

Scattering by Cylinders

Bessel functions are indispensable in the Mie theory of scattering by cylinders and spheres. When a plane wave impinges on a dielectric cylinder, the scattered field is expanded in a series of Bessel and Hankel functions. The internal fields involve Jn(k1ρ) while the external scattered field uses Hn(2)(k0ρ) (the Hankel function of the second kind, which represents outgoing waves). The coefficients of the expansion are determined by the continuity of tangential fields at the cylinder surface. This analysis is crucial in understanding radar cross sections and optical scattering from fibers.

Bessel Beams and Non-Diffracting Fields

An intriguing recent development is the generation of Bessel beams — electromagnetic fields whose transverse profile follows a Bessel function of the first kind. These beams are solutions to the Helmholtz equation that exhibit non-diffracting propagation over a certain distance. They are used in optical tweezers, laser machining, and imaging. The fundamental Bessel beam of order zero has a field proportional to J0(kρρ) e–jβz, and its intensity profile shows a central spot surrounded by concentric rings.

Conical Waveguides and Horn Antennas

Horn antennas with conical geometries also employ spherical Bessel functions. The fields inside a conical horn are expressed using spherical wave functions, with the radial dependence given by spherical Hankel functions. The radiation pattern is then computed from the aperture field. This approach is standard in the design of feedhorns for satellite dishes and radio telescopes.

Computational Techniques and Software

Modern computational electromagnetics software relies on accurate calculations of Bessel functions. Many programming environments (Matlab, Python with SciPy, Mathematica) provide built-in functions for computing Jν, Yν, Iν, Kν, and their spherical counterparts. For large arguments or orders, asymptotic expansions are used to maintain numerical stability. Understanding the underlying mathematics enables engineers and physicists to choose appropriate modes and correctly interpret simulation results. The SciPy documentation provides comprehensive routines for Bessel functions and their zeros.

Conclusion

Bessel functions remain an indispensable mathematical tool in the analysis of electromagnetic wave propagation in geometries with cylindrical or spherical symmetry. From fundamental waveguide theory to advanced topics like non-diffracting beams and Mie scattering, these special functions enable exact analytical solutions that form the foundation of engineering design. Their properties — orthogonality, recurrence relations, asymptotic behavior, and connection to other special functions — allow for systematic modal expansions and closed-form field representations. For any engineer or physicist working with wave phenomena in curvilinear coordinates, a solid grasp of Bessel functions is not merely useful but essential. A thorough exploration of these functions and their role in electromagnetic theory can be found in classical texts such as those by Stratton and Balanis. Mastery of Bessel functions unlocks the ability to analyze and design a wide range of radio-frequency, microwave, and optical systems.