Bessel Functions of the First and Second Kind

⚡ Key Takeaways

Bessel functions $J_\nu(x)$ (first kind) and $Y_\nu(x)$ (second kind) are the two linearly independent solutions of Bessel's differential equation, the canonical equation arising from the Laplacian in cylindrical coordinates. The key result: $J_\nu$ is finite at the origin and oscillates with slowly decaying amplitude, while $Y_\nu$ diverges to $-\infty$ at the origin. They appear wherever waves, heat, or potentials live on circular or cylindrical domains. Introduced by Daniel Bernoulli and Euler, systematized by Friedrich Bessel in 1824.

💡 The Core Insight

Bessel's equation is "the radial wave equation." When you demand that a vibrating circular membrane or cylindrical field be a product $R(r)\Theta(\theta)Z(z)e^{i\omega t}$ , the angular and axial parts give sinusoids, but the radial part is forced to satisfy Bessel's equation. The two kinds, $J_\nu$ and $Y_\nu$ , are simply the "regular at center" and "singular at center" solutions — analogous to $\cos$ and $\sin$ , but adapted to the geometric fact that area grows like $r\,dr$ in the plane.

DefinitionBessel's Differential Equation and the Two Kinds

Bessel's equation of order $\nu$ (with $\nu \ge 0$ real) is

1

\quad x^2\frac{d^2 y}{dx^2} + x\frac{dy}{dx} + (x^2 - \nu^2)\,y = 0.

Derivation hint: Separate the Helmholtz equation $\nabla^2 u + k^2 u = 0$ in cylindrical coordinates; the radial ODE, after substituting $x = kr$ , takes exactly this form. Here $x$ is the (dimensionless) independent variable, $y(x)$ the unknown, and $\nu$ the order governing the angular periodicity. The general solution is

2

\quad y(x) = A\,J_\nu(x) + B\,Y_\nu(x),

Derivation hint: Equation [1] is a second-order linear ODE, so its solution space is two-dimensional; $J_\nu$ and $Y_\nu$ are a fundamental set since their Wronskian is nonzero. where $A,B$ are constants fixed by boundary conditions. The Bessel function of the first kind is defined by its power series

3

\quad J_\nu(x) = \sum_{m=0}^{\infty} \frac{(-1)^m}{m!\,\Gamma(m+\nu+1)}\left(\frac{x}{2}\right)^{2m+\nu},

Derivation hint: Substitute the Frobenius ansatz $y = \sum_m a_m x^{m+\nu}$ into [1], obtain the indicial equation $\nu^2 = \nu^2$ and the two-step recurrence $a_m = -a_{m-2}/[m(m+2\nu)]$ , then normalize with the Gamma function. The Bessel function of the second kind (Neumann function) is defined, for non-integer $\nu$ , as

4

\quad Y_\nu(x) = \frac{J_\nu(x)\cos(\nu\pi) - J_{-\nu}(x)}{\sin(\nu\pi)},

Derivation hint: For non-integer $\nu$ , $J_\nu$ and $J_{-\nu}$ are independent; this particular combination is constructed so the limit $\nu \to n$ (integer) exists and yields a solution independent of $J_n$ . For integer order $n$ , equation [4] is taken as a limit $\nu \to n$ , giving a logarithmic term.

TheoremIndependence and the Wronskian

For all real $\nu$ and $x>0$ , the functions $J_\nu(x)$ and $Y_\nu(x)$ are linearly independent, with Wronskian

5

\quad W\{J_\nu, Y_\nu\}(x) = J_\nu(x)\,Y_\nu'(x) - J_\nu'(x)\,Y_\nu(x) = \frac{2}{\pi x}.

Derivation hint: Write [1] in self-adjoint Sturm–Liouville form $(xy')' + (x - \nu^2/x)y = 0$ ; Abel's identity gives $W = C/x$ , and the constant $C = 2/\pi$ is fixed using the small- $x$ behaviour of $J_\nu$ and $Y_\nu$ . Since $W \neq 0$ for $x>0$ , equation [2] is the complete general solution.

ProofFrobenius Construction of

J_\nu

Step 1 — Indicial equation. Insert $y = \sum_{m\ge 0} a_m x^{m+s}$ into [1]. The lowest-power term ( $x^s$ ) yields $(s^2 - \nu^2)a_0 = 0$ , so the indicial roots are $s = +\nu$ and $s = -\nu$ .

Step 2 — Recurrence. Matching the coefficient of $x^{m+s}$ gives

6

\quad \big[(m+s)^2 - \nu^2\big]a_m + a_{m-2} = 0 \quad\Longrightarrow\quad a_m = \frac{-a_{m-2}}{m(m+2s)}\ \ (s=\nu).

Derivation hint: $(m+\nu)^2 - \nu^2 = m(m+2\nu)$ ; odd-indexed coefficients vanish since $a_1=0$ forced by the $m=1$ relation.

Step 3 — Solve the recurrence. Only even indices survive. Writing $m = 2k$ ,

7

\quad a_{2k} = \frac{(-1)^k\,a_0}{2^{2k}\,k!\,(\nu+1)(\nu+2)\cdots(\nu+k)}.

Derivation hint: Iterate [6] $k$ times; the product in the denominator telescopes into a ratio of Gamma functions.

Step 4 — Normalize. Choose $a_0 = 1/[2^\nu \Gamma(\nu+1)]$ so the Gamma product becomes $\Gamma(k+\nu+1)$ , producing series [3].

Step 5 — Second solution. Repeating with $s = -\nu$ gives $J_{-\nu}$ . For non-integer $\nu$ these are independent. For integer $n$ , however, $J_{-n}(x) = (-1)^n J_n(x)$ , so a genuinely independent second solution must be built — this is $Y_n$ , obtained from [4] by L'Hôpital's rule, which introduces a $\ln x$ singularity. $\blacksquare$

∎

FW

Friedrich Wilhelm Bessel

1784–1846

A German astronomer and mathematician, Bessel was director of the Königsberg Observatory. While studying the perturbations of planetary orbits (specifically the problem of computing elliptical motion via Kepler's equation), he encountered and systematically analyzed the functions $J_n$ in his 1824 memoir Untersuchung des Theils der planetarischen Störungen. Though Daniel Bernoulli (vibrating chains, 1730s) and Euler (membranes, 1764) had met special cases earlier, Bessel gave the first comprehensive treatment of integer-order functions and their recurrence relations. He is equally famous for measuring the first stellar parallax (61 Cygni, 1838).

🧠 Intuition

Think of $\cos x$ and $\sin x$ as solutions of $y'' + y = 0$ — undamped oscillations. Bessel's equation [1] differs by the friction-like first-derivative term $xy'$ . Dividing [1] by $x^2$ :

y'' + \frac{1}{x}y' + \left(1 - \frac{\nu^2}{x^2}\right)y = 0.

The $\tfrac{1}{x}y'$ term acts like a position-dependent damping that weakens as $x$ grows, and the $\nu^2/x^2$ term is a centrifugal barrier near the origin. For large $x$ , the damping fades and the equation looks like $y'' + y = 0$ but with an amplitude that bleeds off as $1/\sqrt{x}$ . Hence the asymptotic form

8

\quad J_\nu(x) \sim \sqrt{\frac{2}{\pi x}}\cos\!\left(x - \frac{\nu\pi}{2} - \frac{\pi}{4}\right),\qquad x\to\infty.

Derivation hint: Substitute $y = u(x)/\sqrt{x}$ into the divided equation; the $1/x$ term cancels and you get $u'' + (1 - (\nu^2-\tfrac14)/x^2)u = 0$ , which $\to u''+u=0$ as $x\to\infty$ . So Bessel functions are "sinusoids on a stretched circular drumhead": they oscillate, but each successive ring of the drum carries less energy density, decaying like $1/\sqrt{x}$ .

Figure

J_0

,

J_1

(first kind) and

Y_0

,

Y_1

(second kind).

J

's are finite at

x=0

;

Y

's plunge to

-\infty

.

```svg

```

ExampleVibrating Circular Drumhead — Fundamental Frequency

Given: A circular drum of radius $a = 0.30$ m, with the membrane fixed at the edge. Wave speed on the membrane $c = 100$ m/s. Find the fundamental (lowest) vibration frequency.

Setup. For the axisymmetric (radial, $\nu=0$ ) mode, displacement $\propto J_0(kr)$ with $k = \omega/c$ . The fixed-edge condition $u(a)=0$ requires

J_0(ka) = 0.

Step 1 — First zero. The smallest positive root of $J_0$ is $j_{0,1} = 2.4048$ . So $ka = 2.4048$ .

Step 2 — Wavenumber. $k = 2.4048 / 0.30 = 8.016\ \text{m}^{-1}$ .

Step 3 — Angular frequency. $\omega = ck = 100 \times 8.016 = 801.6\ \text{rad/s}$ .

Step 4 — Frequency. $f = \omega/(2\pi) = 801.6 / 6.2832 = 127.6\ \text{Hz}$ .

\boxed{f_{0,1} \approx 1.28\times 10^{2}\ \text{Hz}}

(Note: unlike a string, the overtones $j_{0,2}=5.520,\ j_{0,3}=8.654$ are not integer multiples — drums sound inharmonic.)

ExampleSeries Evaluation of

J_0(1)

Given: Compute $J_0(1)$ from series [3] to 4 significant figures.

Step 1 — Write terms ( $\nu=0$ , so $\Gamma(m+1)=m!$ ):

J_0(1) = \sum_{m=0}^{\infty}\frac{(-1)^m}{(m!)^2}\left(\frac12\right)^{2m}.

Step 2 — Tabulate with $(1/2)^{2m} = 4^{-m}$ :

$m=0$ : $1/1 = 1.000000$
$m=1$ : $-\dfrac{1}{(1)^2\cdot 4} = -0.250000$
$m=2$ : $+\dfrac{1}{(2)^2\cdot 16} = +0.015625$
$m=3$ : $-\dfrac{1}{(6)^2\cdot 64} = -0.000434$
$m=4$ : $+\dfrac{1}{(24)^2\cdot 256} = +0.0000068$

Step 3 — Sum: $1 - 0.25 + 0.015625 - 0.000434 + 0.0000068 = 0.765198$ .

\boxed{J_0(1) \approx 0.7652}

(The tabulated value is $0.76519769$ — agreement confirms the rapid convergence of [3].)

⚠️ Common mistake

$J_\nu$ and $Y_\nu$ are NOT periodic. Their zeros only approach uniform spacing $\pi$ as $x\to\infty$ . Treating $J_0$ as if its zeros were equally spaced (like $\cos$ ) gives wrong overtone ratios for drums and wrong cutoff frequencies for waveguides. Also note: $Y_\nu$ must be discarded whenever the domain includes the origin (e.g., a solid disk), because $Y_\nu(x)\to -\infty$ as $x\to 0^+$ . It is retained only for annular regions ( $r_{\min}>0$ ) where the origin is excluded.

Counterexample

One might guess that for integer order, $J_{-n}$ supplies a second independent solution (as $J_{-\nu}$ does for non-integer $\nu$ ). It does not. From series [3], one can show

J_{-n}(x) = (-1)^n\, J_n(x),

so $J_{-n}$ is just a scalar multiple of $J_n$ — not independent. This degeneracy is precisely why $Y_n$ (with its logarithmic singularity) had to be invented: it is the only way to complete the solution basis at integer order.

connection

Bessel functions sit at a crossroads of analysis and physics. They are the cylindrical members of the family of Special Functions alongside Legendre and Hermite functions. They arise from Separation of Variables applied to the Helmholtz Equation and Laplace's Equation in cylindrical coordinates. Their orthogonality on $[0,a]$ with weight $r$ makes them a complete set for Fourier–Bessel Series, the cylindrical analogue of Fourier Series. Modified Bessel functions $I_\nu, K_\nu$ (solutions when the sign of the $x^2$ term flips) connect to Sturm–Liouville Theory and decay/growth problems. The integer-order $J_n$ are also the Fourier coefficients of $e^{ix\sin\theta}$ , linking them to the Generating Functions formalism.

🔬 Where this appears

⏳ History

- 1732–1734 — Daniel Bernoulli analyzes the oscillations of a heavy hanging chain, implicitly encountering

J_0

.

- 1764 — Leonhard Euler studies the vibrations of a stretched circular membrane and writes down series for

J_n

of arbitrary integer order.

- 1770 — Joseph-Louis Lagrange uses Bessel-type series in his work on planetary orbital perturbations.

- 1816 — Bessel begins systematic analysis of the functions in connection with Kepler's equation.

- 1824 — Friedrich Bessel publishes his foundational memoir, giving recurrence relations and integral representations for integer-order

J_n

; the functions take his name.

- 1867 — Carl Neumann introduces the function of the second kind

Y_n

(hence "Neumann function," symbol

N_n

), supplying the missing second solution.

- 1868 — Hermann Hankel develops the functions

H_\nu^{(1,2)} = J_\nu \pm iY_\nu

for radiating/incoming cylindrical waves.

- 1871–1874 — Eugen Lommel writes the first comprehensive treatise, Studien über die Bessel'schen Functionen, unifying the theory.

- 1922 — G. N. Watson publishes A Treatise on the Theory of Bessel Functions, still the definitive reference.

open-problem

The distribution of Bessel zeros remains a source of deep results and open refinements. While McMahon's asymptotic expansion gives the $m$ -th zero of $J_\nu$ as $j_{\nu,m} \approx \beta - \frac{4\nu^2-1}{8\beta}$ with $\beta = (m + \tfrac{\nu}{2} - \tfrac14)\pi$ , the precise gap statistics and monotonicity properties (e.g., whether $j_{\nu,m}$ is convex in $\nu$ , conjectures of Elbert) are still actively studied. Connections between Bessel-zero spacings and random-matrix eigenvalue statistics (the "Bessel kernel" of random matrix theory near a hard edge) link this classical topic to the modern theory of universality. Whether closed-form expressions exist for any nontrivial $j_{\nu,m}$ is essentially answered in the negative, but optimal bounds remain a research frontier.

Bessel Functions of the First and Second Kind

Formal Statement

Why It Holds / Derivation

Building Intuition

Worked Examples

Common Misconceptions

Connections

Real-World Applications

Historical Development

Subtleties & Limits

See Also