ConceptComplete

Special Functions - Examples and Constructions

Special functions appear naturally in solutions to physical problems with specific geometries and boundary conditions, providing exact expressions for wave propagation, diffusion, and quantum states.

Bessel Functions in Wave Problems

ExampleAcoustic Modes in a Circular Drum

A circular membrane clamped at radius aa vibrating at frequency ω\omega has displacement:

u(r,θ,t)=Jm(krr)cos⁡(mθ)cos⁡(ωt)u(r,\theta,t) = J_m(k_r r)\cos(m\theta)\cos(\omega t)

where kr=ω/ck_r = \omega/c and the boundary condition u(a,θ,t)=0u(a,\theta,t) = 0 requires Jm(kra)=0J_m(k_r a) = 0. The allowed frequencies are:

ωmn=cαmna\omega_{mn} = \frac{c\alpha_{mn}}{a}

where αmn\alpha_{mn} is the nn-th zero of JmJ_m. The fundamental mode (m=0,n=1)(m=0, n=1) has frequency ω01=2.405c/a\omega_{01} = 2.405c/a. Higher modes create nodal circles and radial lines.

ExampleBessel Beam in Optics

A Bessel beam maintains its profile while propagating:

E(r,z)=J0(k⊥r)eikzzE(r,z) = J_0(k_\perp r)e^{ik_z z}

where k⊥2+kz2=k2k_\perp^2 + k_z^2 = k^2. Unlike Gaussian beams, Bessel beams are diffraction-free (in the paraxial approximation) but require infinite energy. In practice, they're approximated over finite regions, useful in optical manipulation and microscopy.

Legendre Functions and Electrostatics

ExampleConducting Sphere in Uniform Field

A grounded conducting sphere of radius aa in external field E0z^E_0\hat{\mathbf{z}} has potential:

Φ(r,θ)={−E0rcos⁡θ+E0a3r2cos⁡θr>a0r≤a\Phi(r,\theta) = \begin{cases}-E_0 r\cos\theta + \frac{E_0 a^3}{r^2}\cos\theta & r > a \\ 0 & r \leq a\end{cases}

Using cos⁡θ=P1(cos⁡θ)\cos\theta = P_1(\cos\theta), this is:

Φ(r,θ)=−E0(r−a3r2)P1(cos⁡θ)\Phi(r,\theta) = -E_0\left(r - \frac{a^3}{r^2}\right)P_1(\cos\theta)

The induced surface charge density is:

σ(θ)=−ϵ0∂Φ∂r∣r=a=3ϵ0E0cos⁡θ\sigma(\theta) = -\epsilon_0\frac{\partial\Phi}{\partial r}\bigg|_{r=a} = 3\epsilon_0 E_0\cos\theta

ExampleGravitational Field of Oblate Spheroid

Earth's gravitational potential (oblate spheroid) can be expanded:

Φ(r,θ)=−GMr[1−J2(Rr)2P2(cos⁡θ)+⋯ ]\Phi(r,\theta) = -\frac{GM}{r}\left[1 - J_2\left(\frac{R}{r}\right)^2 P_2(\cos\theta) + \cdots\right]

where J2≈1.08×10−3J_2 \approx 1.08 \times 10^{-3} is the quadrupole moment coefficient and P2(x)=(3x2−1)/2P_2(x) = (3x^2-1)/2. This causes satellite orbit precession, crucial for GPS accuracy.

Quantum Mechanical Applications

ExampleAngular Momentum Eigenst ates

The angular part of the Schrödinger equation in 3D separates:

[1sin⁡θ∂∂θ(sin⁡θ∂∂θ)+1sin⁡2θ∂2∂ϕ2]Y(θ,ϕ)=−ℓ(ℓ+1)Y(θ,ϕ)\left[\frac{1}{\sin\theta}\frac{\partial}{\partial\theta}\left(\sin\theta\frac{\partial}{\partial\theta}\right) + \frac{1}{\sin^2\theta}\frac{\partial^2}{\partial\phi^2}\right]Y(\theta,\phi) = -\ell(\ell+1)Y(\theta,\phi)

Solutions are spherical harmonics Yℓm(θ,ϕ)Y_{\ell m}(\theta,\phi). For hydrogen:

  • n=1,ℓ=0n=1, \ell=0: ψ100∝e−r/a0Y00\psi_{100} \propto e^{-r/a_0}Y_{00} (1s orbital, spherical)
  • n=2,ℓ=1n=2, \ell=1: ψ21m∝re−r/(2a0)Y1m\psi_{21m} \propto re^{-r/(2a_0)}Y_{1m} (2p orbitals, lobed)
  • n=2,ℓ=0n=2, \ell=0: ψ200∝(2−r/a0)e−r/(2a0)Y00\psi_{200} \propto (2-r/a_0)e^{-r/(2a_0)}Y_{00} (2s orbital, radial node)

The angular momentum quantum number ℓ\ell determines the orbital's shape through the spherical harmonic.

ExampleScattering from a Spherical Potential

The radial Schrödinger equation for energy EE and angular momentum ℓ\ell:

d2udr2+[k2−ℓ(ℓ+1)r2−2mℏ2V(r)]u=0\frac{d^2u}{dr^2} + \left[k^2 - \frac{\ell(\ell+1)}{r^2} - \frac{2m}{\hbar^2}V(r)\right]u = 0

where u(r)=rR(r)u(r) = rR(r) and k2=2mE/ℏ2k^2 = 2mE/\hbar^2. For V(r)=0V(r) = 0, solutions are spherical Bessel functions jℓ(kr)j_\ell(kr) and nℓ(kr)n_\ell(kr):

jℓ(x)=π2xJℓ+1/2(x),nℓ(x)=π2xYℓ+1/2(x)j_\ell(x) = \sqrt{\frac{\pi}{2x}}J_{\ell+1/2}(x), \quad n_\ell(x) = \sqrt{\frac{\pi}{2x}}Y_{\ell+1/2}(x)

The scattering cross section involves phase shifts δℓ\delta_\ell determined by matching boundary conditions.

RemarkConfluent Hypergeometric Functions

Many special functions (Laguerre, Hermite polynomials) are special cases of confluent hypergeometric functions 1F1(a;b;z){}_1F_1(a;b;z). For example, hydrogen atom radial functions involve:

Rnℓ(r)∝rℓe−ρ/2Ln−ℓ−12ℓ+1(ρ),ρ=2r/(na0)R_{n\ell}(r) \propto r^\ell e^{-\rho/2}L_{n-\ell-1}^{2\ell+1}(\rho), \quad \rho = 2r/(na_0)

where LnαL_n^{\alpha} are associated Laguerre polynomials. This unified framework reveals deep connections between different physical systems.

These examples demonstrate how special functions provide exact, analytically tractable solutions to fundamental problems in wave physics, electrostatics, and quantum mechanics.