Special Functions - Key Properties

Special functions satisfy remarkable recurrence relations and integral properties that facilitate their use in solving physical problems and evaluating complex integrals.

Recurrence Relations for Bessel Functions

TheoremBessel Function Recurrence Relations

Bessel functions satisfy several key recurrence relations:

$J_{n-1}(x) + J_{n+1}(x) = \frac{2n}{x}J_n(x)$

$J_{n-1}(x) - J_{n+1}(x) = 2J_n'(x)$

$\frac{d}{dx}\left[x^n J_n(x)\right] = x^n J_{n-1}(x)$

$\frac{d}{dx}\left[x^{-n} J_n(x)\right] = -x^{-n} J_{n+1}(x)$

These allow efficient computation of Bessel functions of different orders and their derivatives.

ExampleElectromagnetic Waveguides

In a cylindrical waveguide, the electric field component $E_z(r,\phi,z) = J_m(k_c r)e^{im\phi}e^{ikz}$ must satisfy boundary conditions. The recurrence relations help evaluate field derivatives needed for Maxwell's equations:

$\frac{\partial E_z}{\partial r} = k_c J_m'(k_c r)e^{im\phi}e^{ikz} = -k_c J_{m+1}(k_c r)e^{im\phi}e^{ikz}$

using the recurrence relation $J_n'(x) = -J_{n+1}(x) + (n/x)J_n(x)$ specialized to the boundary.

Generating Functions

DefinitionGenerating Function for Legendre Polynomials

The Legendre polynomials have generating function:

$\frac{1}{\sqrt{1 - 2xt + t^2}} = \sum_{n=0}^{\infty}P_n(x)t^n, \quad |t| < 1$

Expanding around $t = 0$ generates each $P_n(x)$ successively.

This generating function has direct physical meaning: it's the potential at distance $r$ from a unit point charge at distance $r'$ when $x = \cos\gamma$ is the angle between them.

ExampleNewtonian Potential Expansion

For $r > r'$ :

$\frac{1}{|\mathbf{r} - \mathbf{r}'|} = \frac{1}{\sqrt{r^2 - 2rr'\cos\gamma + r'^2}} = \frac{1}{r}\sum_{n=0}^{\infty}\left(\frac{r'}{r}\right)^n P_n(\cos\gamma)$

This is the foundation of multipole expansions in electrostatics and gravitation.

Orthogonality and Completeness

TheoremOrthogonality of Bessel Functions

For zeros $\alpha_m$ and $\alpha_n$ of $J_\nu$ :

$\int_0^a rJ_\nu(\alpha_m r/a)J_\nu(\alpha_n r/a)dr = \frac{a^2}{2}[J_{\nu+1}(\alpha_m)]^2\delta_{mn}$

This orthogonality on $[0,a]$ with weight $r$ allows expansion of functions in Fourier-Bessel series.

ExampleHeat Conduction in a Cylinder

Temperature $T(r,t)$ in a cylinder of radius $a$ satisfies:

$\frac{\partial T}{\partial t} = \alpha\nabla^2 T = \alpha\left(\frac{\partial^2 T}{\partial r^2} + \frac{1}{r}\frac{\partial T}{\partial r}\right)$

The solution is:

$T(r,t) = \sum_{n=1}^{\infty} A_n J_0(\alpha_n r/a)e^{-\alpha(\alpha_n/a)^2 t}$

where $J_0(\alpha_n) = 0$ and:

$A_n = \frac{2}{a^2[J_1(\alpha_n)]^2}\int_0^a r T(r,0)J_0(\alpha_n r/a)dr$

Asymptotics and Approximations

TheoremAsymptotic Behavior of Bessel Functions

For large argument $x \gg n$ :

$J_n(x) \sim \sqrt{\frac{2}{\pi x}}\cos\left(x - \frac{n\pi}{2} - \frac{\pi}{4}\right)$

$Y_n(x) \sim \sqrt{\frac{2}{\pi x}}\sin\left(x - \frac{n\pi}{2} - \frac{\pi}{4}\right)$

These show that Bessel functions behave like damped oscillations for large arguments.

RemarkModified Bessel Functions

For imaginary arguments, we define modified Bessel functions:

$I_n(x) = i^{-n}J_n(ix), \quad K_n(x) = \frac{\pi}{2}i^{n+1}[J_n(ix) + iY_n(ix)]$

These are exponentially growing/decaying rather than oscillatory, appearing in diffusion and heat conduction problems with cylindrical symmetry.

These properties make special functions powerful computational tools, enabling both exact solutions and reliable approximations in mathematical physics.