Bessel Functions

Bessel functions are the canonical solutions

J_{ν} (x)

and

Y_{ν} (x)

of Bessel's differential equation

x^{2} y^{''} + x y^{'} + (x^{2} - ν^{2}) y = 0

. They are the "cylindrical analogues" of sine and cosine, governing oscillations and fields with circular or cylindrical symmetry. Applies to: vibrating drumheads, heat flow in cylinders, electromagnetic waves in fibres, planetary motion, frequency modulation. Origin: Daniel Bernoulli (1730s), formalized by Friedrich Wilhelm Bessel (1817–1824).

A Bessel function is a solution to Bessel's equation, a second-order linear ODE that appears whenever the Laplacian or wave operator is separated in cylindrical or polar coordinates. The functions $J_{ν}$ (first kind) are bounded at the origin and oscillate with slowly decaying amplitude; $Y_{ν}$ (second kind) blow up at the origin. They matter because an enormous class of physical problems — anything with a circular or cylindrical boundary — reduces to them, just as Fourier series reduce rectangular problems to sines and cosines. A surprising fact: although $J_{ν} (x)$ looks like a damped sinusoid, its zeros are not evenly spaced, and the spacing only approaches $π$ asymptotically — the irregularity encodes the geometry of the disk.

The Core Insight

Separate the wave or Laplace equation in cylindrical coordinates and the radial part is forced to obey Bessel's equation. Bessel functions are therefore the radial eigenfunctions of the Laplacian on a disk — the natural "sines and cosines of the round world." Their oscillation comes from the $+ x^{2}$ term (a restoring force); their amplitude decay like $x^{- 1/2}$ comes from energy spreading over an ever-larger circumference.

Formal Statement

Bessel's Differential Equation

For a parameter $ν \geq 0$ (the order), Bessel's equation is

x^{2} \frac{d ^{2} y}{d x ^{2}} + x \frac{d y}{d x} + (x^{2} - ν^{2}) y = 0.

⟳separate

\nabla^{2} u + k^{2} u = 0

in cylindrical coordinates; the radial factor

R (r)

with

x = k r

obeys this equation, with

ν

the angular separation constant.

Here $x$ is the (dimensionless) radial variable and $ν$ is the order, typically an integer $n$ or half-integer for physical problems.

Bessel Function of the First Kind

The solution bounded at $x = 0$ is the Bessel function of the first kind, defined by its series

J_{ν} (x) = m = 0 \sum \infty \frac{( - 1 ) ^{m}}{m ! Γ ( m + ν + 1 )} (\frac{x}{2})^{2 m + ν} .

⟳substitute the Frobenius ansatz

y = \sum_{m} a_{m} x^{m + r}

into [1]; the indicial equation gives

r = \pm ν

, and the recurrence

a_{m} = - a_{m - 2} / [m (m + 2 ν)]

yields the coefficients with

Γ

extending the factorial.

Here $Γ$ is the Gamma function, $Γ (n + 1) = n!$ for integers. For integer order, $J_{- n} (x) = (- 1)^{n} J_{n} (x)$ .

Bessel Function of the Second Kind

The second, linearly independent solution (singular at the origin) is the Bessel function of the second kind (Weber's function):

Y_{ν} (x) = \frac{J _{ν} ( x ) cos ( ν π ) - J _{- ν} ( x )}{sin ( ν π )},

⟳for non-integer

ν

J_{ν}

and

J_{- ν}

are independent; this combination is constructed so the limit

ν \to n

exists (via L'Hôpital) and produces a logarithmic second solution.

For integer $n$ , $Y_{n} (x)$ contains a $ln x$ term and diverges like $- \frac{2}{π} ln (2/ x)$ as $x \to 0^{+}$ .

General Solution

Every solution of Bessel's equation [1] is a linear combination

y (x) = A J_{ν} (x) + B Y_{ν} (x), A, B \in C .

⟳[1] is a second-order linear ODE, so its solution space is two-dimensional;

J_{ν}

and

Y_{ν}

have nonzero Wronskian

W = 2/ (π x)

, hence are independent.

Why It Holds / Derivation

Frobenius Construction of $J_\nu$

Step 1. The point

x = 0

is a regular singular point of [1], so seek

y = \sum_{m \geq 0} a_{m} x^{m + r}

with

a_{0} \neq = 0

. Step 2. Substitute into [1]. The lowest-power term

x^{r}

gives the indicial equation

r^{2} - ν^{2} = 0

, so

r = ν

r = - ν

. Step 3. Equating the coefficient of

x^{m + r}

to zero gives the recurrence

a_{m} [(m + r)^{2} - ν^{2}] + a_{m - 2} = 0 \Rightarrow a_{m} = \frac{- a _{m - 2}}{m ( m + 2 r )} (taking r = ν) .

⟳collect like powers after differentiating the series twice and multiplying by

x^{2}

x

as in [1].

Step 4. All odd-index

a_{m}

vanish (since

a_{- 1} = 0

). For even

m = 2 k

, iterate:

a_{2 k} = \frac{( - 1 ) ^{k} a _{0}}{2 ^{2 k} k ! ( ν + 1 ) ( ν + 2 ) \dots ( ν + k )} .

⟳repeatedly apply [5] with

r = ν

; the denominator telescopes into a rising factorial.

Step 5. Choose the conventional normalization

a_{0} = 1/ [2^{ν} Γ (ν + 1)]

. Using

Γ (ν + k + 1) = (ν + k) \dots (ν + 1) Γ (ν + 1)

, the coefficients collapse exactly into the series [2].

■

Friedrich Wilhelm Bessel | 1784–1846

A self-taught German astronomer and mathematician, Bessel was director of the Königsberg Observatory. He encountered these functions in 1817 while studying Kepler's equation for planetary perturbations — the Fourier coefficients of the eccentric anomaly turned out to be exactly $J_{n}$ . His 1824 memoir gave the first systematic treatment, including recurrence relations and the integral representation. He is also famous for the first reliable measurement of stellar parallax (61 Cygni, 1838) — the same precision instincts that drove his function theory.

Building Intuition

Think of a circular drumhead. A radial slice must vibrate, but unlike a guitar string, the "length" of string available at radius $r$ shrinks toward the center and the energy is spread around a circle of circumference $2 π r$ . So as the wave travels inward it gets compressed (no clean nodes near the center) and outward it gets stretched out (amplitude decays). That is exactly $J_{0}$ : an oscillation whose wavelength is roughly constant but whose amplitude falls like $1/ r$ .

The asymptotic form makes this precise:

J_{ν} (x) \sim \frac{2}{π x} cos (x - \frac{ν π}{2} - \frac{π}{4}), x \to \infty.

⟳apply the method of stationary phase or the WKB approximation to the integral representation; the

x^{- 1/2}

prefactor reflects conservation of energy over an expanding wavefront.

So far from the origin, $J_{ν}$ is just a phase-shifted, amplitude-damped cosine — the cylindrical version of a plane wave.

$J_0$, $J_1$, $J_2$: oscillation with decaying amplitude and unequal zeros

Worked Examples

Fundamental Frequency of a Circular Drum

Given: A circular drumhead of radius

a = 0.30 m

, wave speed

c = 100 m/s

(set by tension and density). The lowest mode is axisymmetric,

u \propto J_{0} (k r) cos (ω t)

. Step 1. Boundary condition: the rim is fixed,

u (a) = 0

, so

J_{0} (k a) = 0

. The smallest positive root is

j_{0, 1} = 2.40483

. Step 2. Thus

k a = j_{0, 1} \Rightarrow k = 2.40483/0.30 = 8.016 m^{- 1}

. Step 3. Frequency:

ω = c k = 100 \times 8.016 = 801.6 rad/s

, so

f = \frac{ω}{2 π} = \frac{801.6}{6.2832} .

f \approx 127.6 Hz

Note: unlike a string, the overtones (next roots $j_{0, 2} = 5.520$ , $j_{0, 3} = 8.654$ ) are not integer multiples — which is why drums sound inharmonic rather than pitched.

Series Evaluation of $J_0$ at a Point

Given: Compute

J_{0} (1)

to four decimals using the series [2] with

ν = 0

. Step 1. With

ν = 0

Γ (m + 1) = m!

, so

J_{0} (x) = m = 0 \sum \infty \frac{( - 1 ) ^{m}}{( m ! ) ^{2}} (\frac{x}{2})^{2 m} .

Step 2. At

x = 1

(x /2)^{2 m} = (1/4)^{m}

. Tabulate terms:

$m = 0 : 1$
$m = 1 : - \frac{1}{( 1 ) ^{2}} (0.25) = - 0.25$
$m = 2 : + \frac{1}{( 2 ) ^{2}} (0.0625) = + 0.015625$
$m = 3 : - \frac{1}{( 6 ) ^{2}} (0.015625) = - 0.000434$
$m = 4 : + \frac{1}{( 24 ) ^{2}} (0.003906) = + 0.0000068$

Step 3. Sum:

1 - 0.25 + 0.015625 - 0.000434 + 0.0000068 = 0.7651978

J_{0} (1) \approx 0.7652

The rapid factorial-squared decay means just five terms give full four-digit accuracy.

Common Misconceptions

$J_{ν}$ is NOT periodic, and its zeros are NOT equally spaced. Students often treat

J_{0}

like a cosine. The zeros

j_{ν, k}

approach equal spacing

π

only asymptotically (from [7]); near the origin they are noticeably farther apart. Likewise the amplitude is not constant — it decays as

x^{- 1/2}

. Treating Bessel functions as exact sinusoids will give wrong eigenvalues and wrong resonant frequencies. "Both solutions are usable on a full disk." Consider heat flow on a solid disk

0 \leq r \leq a

. One might keep

Y_{0}

in the general solution [4]. But

Y_{0} (r) \to - \infty

r \to 0

, which is physically impossible for a finite temperature at the center. Hence

B = 0

is forced: only

J_{0}

survives.

Y_{ν}

is essential only on annular domains (a washer shape) where

r = 0

is excluded — there both kinds are needed to satisfy boundary conditions at the inner and outer radii.

Connections

Bessel functions sit at a crossroads of analysis and physics. They are a special case of the Hypergeometric Function (confluent type), and the modified Bessel functions $I_{ν}, K_{ν}$ solve the same equation with $x^{2} \to - x^{2}$ , appearing in diffusion and decaying fields. They are eigenfunctions in Sturm–Liouville Theory, giving rise to Fourier–Bessel series that generalize Fourier Series to circular domains. The half-integer orders reduce to elementary spherical Bessel functions $j_{n} (x)$ (e.g. $j_{0} = sin x / x$ ), central to Spherical Harmonics and quantum scattering. Their generating function $exp [\frac{x}{2} (t - 1/ t)] = \sum_{n} J_{n} (x) t^{n}$ links them directly to the Laplace Equation and to Generating Functions.

Real-World Applications

Vibrating membranes & musical acoustics — drumheads, cymbals, and the resonance of circular plates use zeros $j_{n, k}$ to set mode frequencies.
Optical fibres & waveguides — the transverse field of guided modes in step-index fibres is exactly $J_{n}$ inside the core and $K_{n}$ in the cladding.
Frequency modulation (FM) radio — the sidebands of an FM-modulated carrier have amplitudes $J_{n} (β)$ where $β$ is the modulation index; Bessel zeros give "carrier null" calibration.
Heat conduction in cylinders — cooling of wires, fuel rods, and pipes solves the cylindrical heat equation with Fourier–Bessel expansions.
Diffraction (Airy disk) — the intensity pattern of light through a circular aperture is $[2 J_{1} (x) / x]^{2}$ , setting the resolution limit of telescopes and microscopes.
Electromagnetic scattering — Mie scattering and antenna theory expand fields in cylindrical Bessel modes.
Quantum mechanics — the radial wavefunction of a free particle and partial-wave scattering amplitudes use spherical Bessel functions.

Historical Development

1732 — Daniel Bernoulli encounters $J_{0}$ studying the oscillation of a hanging chain (the first appearance).
1764 — Leonhard Euler analyzes the vibrating circular membrane, deriving the functions for integer order.
1770 — Lagrange uses series of this type in his theory of planetary orbits.
1817–1824 — Friedrich Bessel, working on Kepler's equation, systematizes the functions, deriving recurrences and integral forms; the 1824 Berlin memoir is the definitive founding paper.
1867 — Carl Neumann and Hermann Hankel develop the second kind and the Hankel functions $H_{ν}^{(1, 2)} = J_{ν} \pm i Y_{ν}$ for travelling cylindrical waves.
1878 — Lord Rayleigh applies Bessel functions extensively in The Theory of Sound.
1879 — Heinrich Weber gives the modern definition of $Y_{ν}$ used in [3].
1922 — G. N. Watson publishes A Treatise on the Theory of Bessel Functions, still the encyclopedic reference.
1964 — Abramowitz & Stegun's Handbook tabulates zeros and values, standardizing computation.

Subtleties & Limits

Several deep questions remain active. The distribution of Bessel zeros $j_{ν, k}$ — their precise spacing and their dependence on $ν$ — connects to random matrix theory and the still-open conjectures about eigenvalue statistics. The transcendence and irrationality of individual zeros and of values like $J_{0} (1)$ is largely unknown; only partial results exist (e.g. positive zeros of $J_{ν}$ are transcendental for rational $ν$ , by Siegel). Efficient uniform asymptotics valid across the transition region near $x \approx ν$ (where the function changes from monotone to oscillatory) require Airy-function corrections and remain a subtle numerical challenge for large order. Finally, the analytic continuation of $ν$ into the complex plane underlies open problems in the theory of integral transforms.

Bessel Functions

The Core Insight

Formal Statement

Bessel's Differential Equation

Bessel Function of the First Kind

Bessel Function of the Second Kind

General Solution

Why It Holds / Derivation

Frobenius Construction of $J_\nu$

Friedrich Wilhelm Bessel | 1784–1846

Building Intuition

$J_0$, $J_1$, $J_2$: oscillation with decaying amplitude and unequal zeros

Worked Examples

Fundamental Frequency of a Circular Drum

Series Evaluation of $J_0$ at a Point

Common Misconceptions

Connections

Real-World Applications

Historical Development

Subtleties & Limits

See Also

Before this

Where this leads