ConceptComplete

Sturm-Liouville Theory

Sturm-Liouville problems arise throughout mathematical physics whenever separation of variables is applied to PDEs. The theory provides a complete orthogonal basis of eigenfunctions, generalizing Fourier analysis to arbitrary domains and potentials.


The Sturm-Liouville Problem

Definition5.4Regular Sturm-Liouville Problem

A regular Sturm-Liouville problem on [a,b][a,b] is: βˆ’ddx[p(x)dydx]+q(x)y=Ξ»w(x)y-\frac{d}{dx}\left[p(x)\frac{dy}{dx}\right] + q(x)y = \lambda w(x) y with boundary conditions Ξ±1y(a)+Ξ±2yβ€²(a)=0\alpha_1 y(a) + \alpha_2 y'(a) = 0, Ξ²1y(b)+Ξ²2yβ€²(b)=0\beta_1 y(b) + \beta_2 y'(b) = 0, where p(x)>0p(x) > 0, w(x)>0w(x) > 0 (weight function), and q(x)q(x) is continuous on [a,b][a,b]. The operator Ly=βˆ’1w[(pyβ€²)β€²βˆ’qy]\mathcal{L}y = -\frac{1}{w}[(py')' - qy] is self-adjoint with respect to the inner product ⟨f,g⟩=∫abf(x)g(x)w(x) dx\langle f, g \rangle = \int_a^b f(x)g(x)w(x)\,dx.

ExampleSturm-Liouville in Quantum Mechanics

The time-independent Schrodinger equation βˆ’β„22mΟˆβ€²β€²+V(x)ψ=Eψ-\frac{\hbar^2}{2m}\psi'' + V(x)\psi = E\psi on [0,L][0, L] with ψ(0)=ψ(L)=0\psi(0) = \psi(L) = 0 is a Sturm-Liouville problem with p=ℏ2/(2m)p = \hbar^2/(2m), q=V(x)q = V(x), w=1w = 1, Ξ»=E\lambda = E. For the infinite square well (V=0V = 0): eigenvalues En=n2Ο€2ℏ2/(2mL2)E_n = n^2\pi^2\hbar^2/(2mL^2), eigenfunctions ψn=2/Lsin⁑(nΟ€x/L)\psi_n = \sqrt{2/L}\sin(n\pi x/L).


Eigenvalue Theory

Definition5.5Sturm-Liouville Eigenvalue Properties

For a regular Sturm-Liouville problem: (1) eigenvalues form an increasing sequence Ξ»1<Ξ»2<β‹―β†’βˆž\lambda_1 < \lambda_2 < \cdots \to \infty; (2) eigenfunctions Ο•n\phi_n corresponding to distinct eigenvalues are orthogonal: βŸ¨Ο•m,Ο•n⟩w=0\langle\phi_m, \phi_n\rangle_w = 0 for mβ‰ nm \neq n; (3) Ο•n\phi_n has exactly nβˆ’1n - 1 zeros in (a,b)(a,b) (Sturm oscillation theorem); (4) the eigenfunctions form a complete orthonormal basis in Lw2([a,b])L^2_w([a,b]).

RemarkGreen's Functions for Sturm-Liouville

The Green's function G(x,ΞΎ;Ξ»)G(x, \xi; \lambda) solves LGβˆ’Ξ»wG=Ξ΄(xβˆ’ΞΎ)\mathcal{L}G - \lambda wG = \delta(x - \xi) and provides the resolvent: (Lβˆ’Ξ»)βˆ’1f(x)=∫abG(x,ΞΎ;Ξ»)f(ΞΎ)w(ΞΎ) dΞΎ(\mathcal{L} - \lambda)^{-1}f(x) = \int_a^b G(x,\xi;\lambda)f(\xi)w(\xi)\,d\xi. The eigenfunction expansion of GG: G(x,ΞΎ;Ξ»)=βˆ‘nΟ•n(x)Ο•n(ΞΎ)Ξ»nβˆ’Ξ»G(x,\xi;\lambda) = \sum_n \frac{\phi_n(x)\phi_n(\xi)}{\lambda_n - \lambda}, with poles at the eigenvalues, connecting the spectral theory to integral operators.