| Legendre polynomial | |
|---|---|
| Name | Legendre polynomials |
| Notation | P_n(x) |
| Developer | Adrien-Marie Legendre |
| Field | Mathematics; applications in Quantum mechanics |
| Introduced | 1782 |
Legendre polynomial
Legendre polynomials are a sequence of real polynomials P_n(x) that solve Legendre's differential equation and form an orthogonal basis on the interval [-1,1]. They are fundamental in expanding angular parts of solutions to the Schrödinger equation in spherical problems and appear throughout quantum physics in contexts such as atomic orbitals and scattering theory. Their mathematical structure underpins analytical and numerical techniques used across physics and engineering.
Legendre polynomials P_n(x) are defined as the polynomial solutions of degree n to Legendre's differential equation, (1 - x^2) y'' - 2x y' + n(n + 1) y = 0, with regularity on [-1,1]. The sequence was studied by Adrien-Marie Legendre and later by Carl Friedrich Gauss in relation to potential theory and spherical harmonics. Key algebraic properties include parity P_n(-x) = (-1)^n P_n(x), explicit expressions via Rodrigues' formula, P_n(x) = (1 / (2^n n!)) d^n/dx^n [(x^2 - 1)^n], and closed-form representations using the hypergeometric function {}_2F_1. Connections exist to classical orthogonal polynomials in the Askey scheme and to special functions studied by Nathaniel Bernstein and others.
A standard generating function for Legendre polynomials is (1 - 2xt + t^2)^{-1/2} = sum_{n=0}^\infty P_n(x) t^n, which is central for deriving identities used in theoretical and computational work. Recurrence relations such as (n+1) P_{n+1}(x) = (2n+1) x P_n(x) - n P_{n-1}(x) enable efficient construction of higher-order polynomials and are exploited in algorithms developed at institutions like Los Alamos National Laboratory and computational projects like GNU Scientific Library. These recurrences relate to three-term recurrences studied in numerical linear algebra (e.g., in work by Gene H. Golub).
Legendre polynomials satisfy orthogonality ∫_{-1}^1 P_n(x) P_m(x) dx = 2/(2n+1) δ_{nm}, which underpins their use in expanding square-integrable functions on the sphere. In quantum mechanics, orthogonality ensures angular eigenfunctions are mutually independent and supports spectral decompositions in Hilbert spaces central to Paul Dirac's formalism. Normalized forms feed directly into construction of spherical harmonics Y_{l}^{m} used in describing eigenstates of the angular momentum operator and in solutions of the hydrogen atom problem as formulated by Erwin Schrödinger and refined by P. A. M. Dirac.
When the time-independent Schrödinger equation is separated in spherical coordinates for central potentials (e.g., the Coulomb potential in the hydrogen atom), the polar-angle equation reduces to Legendre's equation for the m=0 case; more generally it yields associated Legendre equations. Solutions in terms of P_l(cos θ) multiply radial functions R_{nl}(r) to produce atomic orbitals described in textbooks by L. D. Landau and E. M. Lifshitz. Legendre polynomials also enter partial-wave expansions in scattering theory developed in the context of nuclear and particle physics at institutions like CERN and Brookhaven National Laboratory.
Associated Legendre functions P_l^m(x) generalize P_l(x) and appear in the construction of spherical harmonics Y_l^m(θ,φ), which are simultaneous eigenfunctions of the squared angular momentum operator L^2 and L_z. These functions are central to the representation theory of the rotation group SO(3) and its cover SU(2), with historical contributions from Eugène Wigner and modern applications in quantum chemistry (e.g., at IUPAC-referenced computational codes) and condensed matter physics. The ladder operator formalism used by Julian Schwinger and others provides algebraic routes to these functions, linking to selection rules and quantum numbers observed in spectroscopy experiments at facilities like National Institute of Standards and Technology (NIST).
Computing Legendre and associated Legendre functions reliably is crucial for numerical solutions of the Schrödinger equation, multipole expansions in electrodynamics, and simulations in quantum chemistry packages such as Gaussian (software) and Quantum ESPRESSO. Stable algorithms include upward/downward recurrence, Clenshaw summation, and use of orthogonal polynomial libraries developed in projects at Argonne National Laboratory and by researchers like Walter Gautschi. Careful handling of floating-point errors, underflow/overflow, and normalization is required for high angular momenta in scattering calculations and for eigenstate expansions used in density functional theory (DFT) and many-body methods.
Access to robust implementations of Legendre polynomials affects equity in education and research: open-source libraries (e.g., SciPy, GNU Scientific Library) democratize tools, while proprietary codes and paywalled textbooks can create barriers for students in under-resourced institutions. Initiatives by organizations like the Open Source Initiative and programs at universities such as Massachusetts Institute of Technology (MIT OpenCourseWare) and University of Cape Town aim to broaden access to computational curricula. Ensuring inclusive pedagogy around mathematical topics—highlighting contributions from diverse scientists and providing multilingual resources—aligns with broader scientific justice goals promoted by groups like Association for Women in Science and community-driven repositories that implement Legendre-related routines for global researchers.
Category:Special functions Category:Mathematical physics Category:Quantum mechanics