ConceptComplete

Special Functions - Core Definitions

Special functions arise naturally as solutions to differential equations in physics, particularly when solving PDEs in curvilinear coordinates. They encode the geometric and symmetry properties of physical systems.

Bessel Functions

DefinitionBessel Functions of the First Kind

The Bessel function Jn(x)J_n(x) of order nn is a solution to Bessel's equation:

x2d2ydx2+xdydx+(x2−n2)y=0x^2\frac{d^2y}{dx^2} + x\frac{dy}{dx} + (x^2 - n^2)y = 0

The series representation is:

Jn(x)=∑k=0∞(−1)kk!(n+k)!(x2)n+2kJ_n(x) = \sum_{k=0}^{\infty}\frac{(-1)^k}{k!(n+k)!}\left(\frac{x}{2}\right)^{n+2k}

For integer nn, Jn(x)=(−1)nJ−n(x)J_n(x) = (-1)^nJ_{-n}(x). These functions appear in problems with cylindrical symmetry.

ExampleVibrating Circular Membrane

The displacement u(r,θ,t)u(r,\theta,t) of a circular drumhead satisfies:

∇2u=1c2∂2u∂t2\nabla^2 u = \frac{1}{c^2}\frac{\partial^2 u}{\partial t^2}

Separating variables u(r,θ,t)=R(r)Θ(θ)T(t)u(r,\theta,t) = R(r)\Theta(\theta)T(t) leads to:

r2R′′+rR′+(k2r2−m2)R=0r^2R'' + rR' + (k^2r^2 - m^2)R = 0

The solution is R(r)=Jm(kr)R(r) = J_m(kr). Boundary condition u(a,θ,t)=0u(a,\theta,t) = 0 requires Jm(ka)=0J_m(ka) = 0, giving discrete allowed frequencies Ήmn=cÎąmn/a\omega_{mn} = c\alpha_{mn}/a where Îąmn\alpha_{mn} are zeros of JmJ_m.

DefinitionBessel Functions of the Second Kind

The Neumann function (or Weber function) Yn(x)Y_n(x) is a second linearly independent solution to Bessel's equation, singular at x=0x = 0:

Yn(x)=Jn(x)cos⁥(nĪ€)−J−n(x)sin⁥(nĪ€)Y_n(x) = \frac{J_n(x)\cos(n\pi) - J_{-n}(x)}{\sin(n\pi)}

For integer nn, this is defined by taking the limit. The general solution is y(x)=AJn(x)+BYn(x)y(x) = AJ_n(x) + BY_n(x).

Legendre Polynomials

DefinitionLegendre Polynomials

The Legendre polynomials Pn(x)P_n(x) are solutions to Legendre's equation:

ddx[(1−x2)dydx]+n(n+1)y=0\frac{d}{dx}\left[(1-x^2)\frac{dy}{dx}\right] + n(n+1)y = 0

They can be generated by Rodrigues' formula:

Pn(x)=12nn!dndxn(x2−1)nP_n(x) = \frac{1}{2^n n!}\frac{d^n}{dx^n}(x^2 - 1)^n

The first few are P0(x)=1P_0(x) = 1, P1(x)=xP_1(x) = x, P2(x)=(3x2−1)/2P_2(x) = (3x^2 - 1)/2, P3(x)=(5x3−3x)/2P_3(x) = (5x^3 - 3x)/2.

Legendre polynomials form a complete orthogonal set on [−1,1][-1,1]:

âˆĢ−11Pm(x)Pn(x)dx=22n+1δmn\int_{-1}^1 P_m(x)P_n(x)dx = \frac{2}{2n+1}\delta_{mn}

ExampleGravitational and Electrostatic Multipole Expansion

The potential due to a localized charge distribution at position r′\mathbf{r}' observed at r\mathbf{r} with r>r′r > r':

1âˆŖr−râ€˛âˆŖ=∑ℓ=0∞r′ℓrℓ+1Pℓ(cosâĄÎŗ)\frac{1}{|\mathbf{r} - \mathbf{r}'|} = \sum_{\ell=0}^{\infty}\frac{r'^{\ell}}{r^{\ell+1}}P_\ell(\cos\gamma)

where Îŗ\gamma is the angle between r\mathbf{r} and r′\mathbf{r}'. This gives:

  • ℓ=0\ell = 0: monopole ∝1/r\propto 1/r
  • ℓ=1\ell = 1: dipole ∝1/r2\propto 1/r^2
  • ℓ=2\ell = 2: quadrupole ∝1/r3\propto 1/r^3

Essential for analyzing gravity and electromagnetic fields.

Spherical Harmonics

DefinitionSpherical Harmonics

The spherical harmonics Yℓm(θ,Ī•)Y_{\ell m}(\theta,\phi) are eigenfunctions of the angular momentum operator in quantum mechanics:

Yℓm(θ,Ī•)=(2ℓ+1)(ℓ−m)!4Ī€(ℓ+m)!Pℓm(cos⁥θ)eimĪ•Y_{\ell m}(\theta,\phi) = \sqrt{\frac{(2\ell+1)(\ell-m)!}{4\pi(\ell+m)!}}P_\ell^m(\cos\theta)e^{im\phi}

where PℓmP_\ell^m are associated Legendre functions. They satisfy:

∇angular2Yℓm=−ℓ(ℓ+1)Yℓm\nabla^2_{\text{angular}}Y_{\ell m} = -\ell(\ell+1)Y_{\ell m}

Spherical harmonics are orthonormal on the unit sphere:

âˆĢ02Ī€âˆĢ0Ī€Yℓm∗(θ,Ī•)Yℓ′m′(θ,Ī•)sin⁡θ dθ dĪ•=δℓℓ′δmm′\int_0^{2\pi}\int_0^{\pi}Y_{\ell m}^*(\theta,\phi)Y_{\ell' m'}(\theta,\phi)\sin\theta\,d\theta\,d\phi = \delta_{\ell\ell'}\delta_{mm'}

ExampleHydrogen Atom Wave Functions

The angular part of hydrogen atom wave functions are spherical harmonics. For example:

  • Y00=1/4Ī€Y_{00} = 1/\sqrt{4\pi} (s orbital, spherically symmetric)
  • Y10=3/(4Ī€)cos⁥θY_{10} = \sqrt{3/(4\pi)}\cos\theta (pz_z orbital)
  • Y1Âą1=∓3/(8Ī€)sin⁥θeÂąiĪ•Y_{1\pm 1} = \mp\sqrt{3/(8\pi)}\sin\theta e^{\pm i\phi} (px_x, py_y orbitals)

The quantum numbers ℓ\ell (orbital angular momentum) and mm (magnetic quantum number) label the symmetry of the wave function.

RemarkAddition Theorem

The addition theorem for spherical harmonics:

Pℓ(cosâĄÎŗ)=4Ī€2ℓ+1∑m=−ℓℓYℓm∗(θ′,Ī•â€˛)Yℓm(θ,Ī•)P_\ell(\cos\gamma) = \frac{4\pi}{2\ell+1}\sum_{m=-\ell}^{\ell}Y_{\ell m}^*(\theta',\phi')Y_{\ell m}(\theta,\phi)

where Îŗ\gamma is the angle between directions (θ,Ī•)(\theta,\phi) and (θ′,Ī•â€˛)(\theta',\phi'). This connects Legendre polynomials to spherical harmonics and is crucial in scattering theory.

These special functions provide the mathematical language for describing systems with spherical and cylindrical symmetry throughout classical and quantum physics.