Community Hub0 Subscribers
Community hub
Legendre function
View on Wikipediafrom Wikipedia
 | This article includes a list of general references, but it lacks sufficient corresponding inline citations. (January 2013) |
In physical science and mathematics, the Legendre functions Pλ, Qλ and associated Legendre functions Pμ
λ, Qμ
λ, and Legendre functions of the second kind, Qn, are all solutions of Legendre's differential equation. The Legendre polynomials and the associated Legendre polynomials are also solutions of the differential equation in special cases, which, by virtue of being polynomials, have a large number of additional properties, mathematical structure, and applications. For these polynomial solutions, see the separate Wikipedia articles.
Legendre's differential equation
[edit]The general Legendre equation reads
where the numbers λ and μ may be complex, and are called the degree and order of the relevant function, respectively. The polynomial solutions when λ is an integer (denoted n), and μ = 0 are the Legendre polynomials Pn; and when
λ is an integer (denoted n), and μ = m is also an integer with |m| < n are the associated Legendre polynomials. All other cases of λ and μ can be discussed as one, and the solutions are written Pμ
λ, Qμ
λ. If μ = 0, the superscript is omitted, and one writes just Pλ, Qλ. However, the solution Qλ when λ is an integer is often discussed separately as Legendre's function of the second kind, and denoted Qn.
![{\displaystyle \left(1-x^{2}\right)y''-2xy'+\left[\lambda (\lambda +1)-{\frac {\mu ^{2}}{1-x^{2}}}\right]y=0,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bbb35fc0e28f0018a0188e4b157b25a13cc7d5f3)
This is a second order linear equation with three regular singular points (at 1, −1, and ∞). Like all such equations, it can be converted into a hypergeometric differential equation by a change of variable, and its solutions can be expressed using hypergeometric functions.
Solutions of the differential equation
[edit]Since the differential equation is linear, homogeneous (the right hand side =zero) and of second order, it has two linearly independent solutions, which can both be expressed in terms of the hypergeometric function, . With being the gamma function, the first solution is and the second is
![{\displaystyle P_{\lambda }^{\mu }(z)={\frac {1}{\Gamma (1-\mu )}}\left[{\frac {z+1}{z-1}}\right]^{\mu /2}\,_{2}F_{1}\left(-\lambda ,\lambda +1;1-\mu ;{\frac {1-z}{2}}\right),\qquad {\text{for }}\ |1-z|<2,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/311143e4f03184fdec4c1f8ae7f9e619a890012e)
_with_n%3D0.5_in_the_complex_plane_from_-2-2i_to_2%2B2i_with_colors_created_with_Mathematica_13.1_function_ComplexPlot3D.svg)
These are generally known as Legendre functions of the first and second kind of noninteger degree, with the additional qualifier 'associated' if μ is non-zero. A useful relation between the P and Q solutions is Whipple's formula.
Positive integer order
[edit]For positive integer the evaluation of above involves cancellation of singular terms. We can find the limit valid for as[1]
![{\displaystyle P_{\lambda }^{m}(z)=\lim _{\mu \to m}P_{\lambda }^{\mu }(z)={\frac {(-\lambda )_{m}(\lambda +1)_{m}}{m!}}\left[{\frac {1-z}{1+z}}\right]^{m/2}\,_{2}F_{1}\left(-\lambda ,\lambda +1;1+m;{\frac {1-z}{2}}\right),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/270c3b94e8cddf4350ca5ac598ce668b198c2ce9)
with the (rising) Pochhammer symbol.
Legendre functions of the second kind (Qn)
[edit]
The nonpolynomial solution for the special case of integer degree , and , is often discussed separately. It is given by
This solution is necessarily singular when .
The Legendre functions of the second kind can also be defined recursively via Bonnet's recursion formula
Associated Legendre functions of the second kind
[edit]The nonpolynomial solution for the special case of integer degree , and is given by
Integral representations
[edit]The Legendre functions can be written as contour integrals. For example, where the contour winds around the points 1 and z in the positive direction and does not wind around −1. For real x, we have
Legendre function as characters
[edit]The real integral representation of are very useful in the study of harmonic analysis on where is the double coset space of (see Zonal spherical function). Actually the Fourier transform on is given by where
Singularities of Legendre functions of the first kind (Pλ) as a consequence of symmetry
[edit]Legendre functions Pλ of non-integer degree are unbounded at the interval [-1, 1] . In applications in physics, this often provides a selection criterion. Indeed, because Legendre functions Qλ of the second kind are always unbounded, in order to have a bounded solution of Legendre's equation at all, the degree must be integer valued: only for integer degree, Legendre functions of the first kind reduce to Legendre polynomials, which are bounded on [-1, 1]. It can be shown[2] that the singularity of the Legendre functions Pλ for non-integer degree is a consequence of the mirror symmetry of Legendre's equation. Thus there is a symmetry under the selection rule just mentioned.
See also
[edit]References
[edit]- ^ Creasey, Peter E.; Lang, Annika (2018). "Fast generation of isotropic Gaussian random fields on the sphere". Monte Carlo Methods and Applications. 24 (1): 1–11. arXiv:1709.10314. Bibcode:2018MCMA...24....1C. doi:10.1515/mcma-2018-0001. S2CID 4657044.
- ^ van der Toorn, Ramses (4 April 2022). "The Singularity of Legendre Functions of the First Kind as a Consequence of the Symmetry of Legendre's Equation". Symmetry. 14 (4): 741. Bibcode:2022Symm...14..741V. doi:10.3390/sym14040741. ISSN 2073-8994.
- Abramowitz, Milton; Stegun, Irene Ann, eds. (1983) [June 1964]. "Chapter 8". Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables. Applied Mathematics Series. Vol. 55 (Ninth reprint with additional corrections of tenth original printing with corrections (December 1972); first ed.). Washington D.C.; New York: United States Department of Commerce, National Bureau of Standards; Dover Publications. p. 332. ISBN 978-0-486-61272-0. LCCN 64-60036. MR 0167642. LCCN 65-12253.
- Courant, Richard; Hilbert, David (1953), Methods of Mathematical Physics, Volume 1, New York: Interscience Publisher, Inc.
- Dunster, T. M. (2010), "Legendre and Related Functions", in Olver, Frank W. J.; Lozier, Daniel M.; Boisvert, Ronald F.; Clark, Charles W. (eds.), NIST Handbook of Mathematical Functions, Cambridge University Press, ISBN 978-0-521-19225-5, MR 2723248.
- Ivanov, A.B. (2001) [1994], "Legendre function", Encyclopedia of Mathematics, EMS Press
- Snow, Chester (1952) [1942], Hypergeometric and Legendre functions with applications to integral equations of potential theory, National Bureau of Standards Applied Mathematics Series, No. 19, Washington, D.C.: U. S. Government Printing Office, hdl:2027/mdp.39015011416826, MR 0048145
- Whittaker, E. T.; Watson, G. N. (1963), A Course in Modern Analysis, Cambridge University Press, ISBN 978-0-521-58807-2
{{citation}}: ISBN / Date incompatibility (help)
External links
[edit]- Legendre function P on the Wolfram functions site.
- Legendre function Q on the Wolfram functions site.
- Associated Legendre function P on the Wolfram functions site.
- Associated Legendre function Q on the Wolfram functions site.
| International | |
|---|---|
| National | |
| Other | |
Legendre function
View on Grokipediafrom Grokipedia
Mathematical Foundations
Legendre's Differential Equation
The Legendre differential equation is a second-order linear ordinary differential equation that defines the Legendre functions. In its general form, known as the associated Legendre equation, it is given by
where denotes the degree and the order, both typically complex parameters, though often taken as real or integers in applications.[1] When , this reduces to the standard Legendre equation
which serves as the foundational case for the unassociated Legendre functions.[1]
This equation was introduced by Adrien-Marie Legendre in 1782 as part of his work on the gravitational attraction of ellipsoids, where the functions arose in the series expansion of the Newtonian potential for spheroidal mass distributions.[4] The equation exhibits regular singular points at and at infinity, with indicial exponents at the finite singularities and at infinity, classifying it as a Fuchsian equation.[1]
The Legendre equation can be transformed into the Gauss hypergeometric differential equation via the substitution , which maps the interval to and facilitates expressing solutions in terms of hypergeometric functions. The two linearly independent solutions to the associated equation are conventionally denoted and .[2]
General Solutions via Hypergeometric Functions
The Legendre differential equation, along with its associated form, possesses regular singular points at and at infinity. To derive general analytic solutions, the method of Frobenius is applied around these points, particularly at infinity via the substitution to expand near . This yields an indicial equation whose roots are and , determining the leading behaviors of the solutions as . These exponents reflect the Fuchsian nature of the equation and facilitate the construction of series expansions that converge in appropriate domains.[6] The resulting series solutions can be expressed in terms of the Gauss hypergeometric function . For the associated Legendre function of the first kind, one solution is given by
up to normalization constants that ensure standard conventions. This representation arises from the transformation , which maps the interval to and converts the differential equation into the standard hypergeometric form. The second linearly independent solution is the associated Legendre function of the second kind , which involves a more complex combination of hypergeometric functions to account for the logarithmic singularity at the branch points. The general solution to the associated Legendre equation is thus , where and are arbitrary constants.[2]
The hypergeometric series converges absolutely for , corresponding to in the transformed variable, with radius of convergence 1 in the -plane. For non-integer and , the functions exhibit branch points at , necessitating analytic continuation beyond the principal domain using connection formulas or integral representations to define single-valued branches on the complex plane, often with cuts along . This continuation preserves linear independence between and across the Riemann surface.
In the specific case , the associated form reduces to the standard Legendre equation, yielding the Legendre function of the first kind , with the second solution ensuring completeness. For integer , this specializes to polynomial solutions, but the hypergeometric form remains valid for general .
Legendre Functions of the First Kind
Definition and Explicit Forms for Integer Order
The Legendre functions of the first kind for integer order, denoted with an integer, are polynomial solutions of degree to Legendre's differential equation with azimuthal order , uniquely determined by the normalization condition .[7] These polynomials form an orthogonal basis on the interval with respect to the weight function , and their leading coefficient is , where denotes the Pochhammer symbol.[7] One explicit representation is given by Rodrigues' formula:
[8] This formula generates the polynomials directly from repeated differentiation of the binomial , ensuring the required degree and normalization. An equivalent form uses with an adjusted sign factor, but the version above is conventional for deriving further properties.[9]
The generating function for the sequence is
valid for and .[10] This closed-form expression facilitates the expansion of functions in terms of Legendre polynomials and highlights their role in potential theory and spherical harmonics.
For low orders, the polynomials take simple forms: , , and .[9] These examples illustrate the progression from constants to higher-degree terms with increasing powers of .
The Legendre polynomials exhibit definite parity: , making them even functions for even and odd for odd .[9] Additionally, they are bounded on the interval of orthogonality, satisfying for all , with equality at the endpoints due to the normalization.[8]
General Order and Branch Cuts
The Legendre function of the first kind for general real or complex order is defined via analytic continuation of the hypergeometric representation, extending beyond integer degrees where it reduces to polynomials. For , it is given by
where denotes the Gauss hypergeometric function.[2] This expression provides the principal branch, which is real-valued when and are real. As approaches a non-negative integer , limits to the Legendre polynomial .
The domain of the principal branch is the interval , with analytic continuation to the complex plane excluding branch cuts along the real axis from to 1. This cut structure arises from the singularities at the endpoints and at infinity, ensuring single-valuedness within the principal sheet. For , values are obtained via connection formulas that account for the monodromy around the branch points.[11]
For associated Legendre functions, the Ferrers function of the first kind serves as a variant suited to the interval , defined as
where is the regularized hypergeometric function and . This form, often expressed in terms of in some conventions, incorporates a prefactor that ensures real values for real arguments in the principal domain and facilitates applications in toroidal coordinates.[2] The Ferrers function relates to the standard associated Legendre function by a transformation involving , highlighting normalization differences: Legendre functions emphasize polynomial behavior for integer orders, while Ferrers variants prioritize boundedness and reality on .[2]
For large , asymptotic approximations of are obtained using the Mehler-Dirichlet integral representation:
valid for . This integral allows evaluation via methods such as stationary phase, yielding leading-order behavior proportional to near and oscillatory decay elsewhere, which is crucial for high-frequency wave problems.[12] Uniqueness of these functions for general holds up to the choice of normalization, with Legendre and Ferrers conventions differing primarily in phase factors and prefactors for associated cases.[2]
Legendre Functions of the Second Kind
Definition and Explicit Forms
The Legendre functions of the second kind, denoted $ Q_\nu(x) $, provide the second linearly independent solution to Legendre's differential equation, complementing the functions of the first kind $ P_\nu(x) $. For $ x > 1 $, an explicit representation is given by
where $ W_{\nu-1}(x) $ is a series expressible in terms of hypergeometric functions or other special functions, ensuring the form captures the singular behavior.[13]
For integer orders $ n = 0, 1, 2, \dots $, the expression simplifies, with $ W_{n-1}(x) $ becoming a polynomial of degree $ n-1 $:
This sum arises from the connection formulas and ensures $ Q_n(x) $ is well-defined for computation.[14]
These functions are analytic in the domain $ |x| > 1 $, where the principal branch of the logarithm is taken, but exhibit logarithmic singularities at the branch points $ x = \pm 1 $, reflecting their role in solutions with singular behavior.[15]
A key normalization occurs for $ n = 0 $:
which follows directly from the general formula with $ P_0(x) = 1 $ and $ W_{-1}(x) = 0 $.[16]
The functions $ Q_\nu(x) $ and $ P_\nu(x) $ are linearly independent over the complex plane, as their Wronskian is nonzero, forming a fundamental pair for solving Legendre's equation.
Behavior at Singular Points
The Legendre functions of the second kind, $ Q_\nu(x) $, exhibit singular behavior at the regular singular points $ x = \pm 1 $ of Legendre's differential equation, characterized by logarithmic divergences that distinguish them from the functions of the first kind. Near $ x = 1^- $, for fixed $ \nu \not\in {-1, -2, \dots} $, the leading asymptotic behavior is $ Q_\nu(x) \sim -\frac{1}{2} \ln(1 - x) $, with higher-order terms involving constants such as Euler's constant $ \gamma $ and the digamma function $ \psi(\nu + 1) $. This logarithmic singularity arises from the integral representation or hypergeometric series expansion of $ Q_\nu(x) $, ensuring the function is unbounded as the argument approaches the endpoint of the interval $ (-1, 1) $. A similar logarithmic divergence occurs near $ x = -1^+ $, where the coefficient depends on $ \nu $ through connection formulas relating values across the branch cut, typically yielding $ Q_\nu(x) \sim \frac{1}{2} P_\nu(-1) \ln(1 + x) + $ subleading terms, with $ P_\nu(-1) = \cos(\pi \nu) $ providing the parity adjustment for integer orders. This behavior reflects the symmetric placement of singular points in the differential equation, with the precise coefficient incorporating $ \nu $ to maintain consistency with recurrence relations. At infinity, for large $ |x| $ with $ x > 1 $, $ Q_\nu(x) $ decays algebraically as $ Q_\nu(x) \sim \frac{\sqrt{\pi}}{\Gamma(\nu + 3/2) (2x)^{\nu+1}} $, providing the dominant term in the expansion that ensures convergence in applications requiring solutions at large distances. This asymptotic form highlights the polynomial-like decay modulated by $ \nu $. The singular behaviors at $ x = \pm 1 $ are linked to the monodromy of solutions around these branch points, as described by the Riemann P-symbol for Legendre's equation, which encodes the local exponents $ 0, 0 $ at each finite singular point and the branching structure for non-integer $ \nu $.[17] This monodromy analysis reveals how encircling the points $ \pm 1 $ induces logarithmic phase shifts in $ Q_\nu(x) $, contrasting with the single-valued nature of $ P_\nu(x) $ near these points. In comparison, the Legendre functions of the first kind $ P_\nu(x) $ remain finite at $ x = \pm 1 $ for all $ \nu $, with $ P_\nu(1) = 1 $ and $ P_\nu(-1) = \cos(\pi \nu) $, thus completing an independent basis of solutions without singularities in the physical domain $ [-1, 1] $ for integer orders.Associated Legendre Functions
Definitions for First and Second Kind
The associated Legendre functions generalize the Legendre functions of the first and second kinds by introducing a non-zero order parameter , which arises naturally in problems with azimuthal dependence, such as those in spherical coordinates. These functions satisfy the associated Legendre differential equation,
where is the degree and is the order, typically taken as integers in many applications.[1]
For integer , the associated Legendre function of the first kind is defined as
where is the Legendre function of the first kind. This definition ensures the function is regular at for appropriate and , and it reduces to the standard Legendre function when . For integer degree and order with , becomes a polynomial of degree , often referred to as an associated Legendre polynomial. A representative example is .[18][19]
The associated Legendre function of the second kind, for integer , is similarly defined as
where is the Legendre function of the second kind. This function exhibits logarithmic singularities at and is used to form the complete set of solutions to the associated Legendre equation. An illustrative example is .[13]
In spherical coordinates, the associated Legendre functions of the first kind are integral to the definition of spherical harmonics, where , providing the -dependent part of the angular solutions to Laplace's equation in three dimensions.[20]
Normalization and Symmetry Properties
The normalization of associated Legendre functions of the first kind, , often incorporates the Condon-Shortley phase in physics contexts to ensure orthonormality when combined with azimuthal factors in spherical harmonics. Specifically, for integer orders and , the convention is , where is the Legendre polynomial of degree . This phase factor simplifies matrix elements in angular momentum calculations and aligns with the unit normalization of spherical harmonics over the unit sphere, . Symmetry properties of these functions include parity relations that facilitate their use in even-odd decompositions. For nonnegative integers and with , . The relation to negative order is given by , ensuring consistency across positive and negative . Alternative normalizations appear in specialized fields, such as geophysics, where Ferrers functions are defined without the Condon-Shortley phase as for , often paired with Schmidt quasi-normalization for surface spherical harmonics to handle real-valued geomagnetic data efficiently.[21] For imaginary orders ( real), Ferrers functions extend to toroidal coordinates, aiding in geophysical modeling of Earth's magnetic field. A key reflection formula connects functions of degree and : , which is essential for analytic continuation and reciprocity in potential theory.[22] In quantum applications, half-integer orders arise in contexts like the Dirac equation for relativistic particles, where functions such as exhibit behaviors tied to hypergeometric representations, maintaining orthogonality under specific boundary conditions while preserving the overall symmetry structure.Representations and Expansions
Series Expansions
The Legendre functions of the first kind possess a series expansion expressible through the Gauss hypergeometric function, facilitating numerical computation within appropriate domains. Specifically,
P_\nu(x) = \, _2F_1\left(-\nu, \nu+1; 1; \frac{1-x}{2}\right),
where the hypergeometric function expands as the power series
with Pochhammer symbols and . Substituting the parameters yields
This series converges for , equivalent to , and extends analytically to the complex plane excluding the branch cut . For integer orders , the series terminates at due to the vanishing Pochhammer symbol for , reducing to a polynomial of degree .
For the Legendre functions of the second kind , the expansion incorporates a logarithmic singularity reflecting the branch point at . For ,
Q_\nu(x) = \frac{\sqrt{\pi} \, \Gamma(\nu+1)}{2^{\nu+1} \Gamma\left(\nu + \frac{3}{2}\right)} x^{-\nu-1} \, _2F_1\left( \frac{\nu+1}{2}, \frac{\nu+2}{2}; \nu + \frac{3}{2}; \frac{1}{x^2} \right),
which provides a power series in converging for . For integer orders and ,
where the finite sum involves polynomials of lower degree, enabling efficient evaluation. This form arises from the limiting case of the general hypergeometric representation and highlights the non-polynomial nature of , diverging logarithmically as . For , the principal value uses the Ferrers function of the second kind, which includes similar logarithmic and hypergeometric components but requires careful branch handling.
Associated Legendre functions of the first kind also admit hypergeometric series expansions. For integer orders and with ,
P_n^m(x) = (-1)^m \left(1 - x^2\right)^{m/2} \frac{(n+m)!}{2^m m! (n-m)!} \, _2F_1\left(-n+m, n+m+1; m+1; \frac{1-x}{2}\right),
converging under the same condition . For general , the form generalizes accordingly, with termination occurring when parameters cause Pochhammer symbols to vanish. Similar expansions hold for associated functions of the second kind , often combining hypergeometric terms with logarithmic factors.
Integral Representations
Integral representations of Legendre functions trace their origins to the late 18th-century developments in potential theory by Adrien-Marie Legendre and Pierre-Simon Laplace, who employed such forms to expand the gravitational potential of spheroids and ellipsoids in series and integrals for celestial mechanics problems. These early Laplace integrals facilitated the solution of Poisson's equation in spherical coordinates, laying the groundwork for modern analytic expressions that enable evaluation and continuation beyond integer orders. A prominent real integral representation for the Legendre function of the first kind, known as the Mehler-Dirichlet integral, is given by
valid for and complex with appropriate convergence conditions. This formula, derived by Dirichlet in 1836 and generalized by Mehler in 1847, proves especially effective for analytic continuation to non-integer and for deriving asymptotic behaviors in the complex plane.
For general complex arguments, the Schläfli contour integral offers a powerful representation for :
where the closed contour encircles the branch cut in the positive sense, assuming outside the cut and suitable branch choices for the multi-valued functions. Introduced by Ludwig Schläfli in 1850, this integral circumvents singularities and supports analytic continuation across branch cuts, making it ideal for computational purposes in regions away from the real interval .
The Legendre function of the second kind possesses a real-line integral representation linking it directly to :
for , interpreted as a Cauchy principal value when approaches the cut. This form, arising from the theory of singular integrals and orthogonal expansions, facilitates the computation of via known values of and is instrumental in asymptotic expansions for large .
For associated Legendre functions, the Schläfli integral generalizes to
with the contour encircling and branches chosen consistently; a similar form holds for negative order . This representation extends the utility of contour integrals to associated cases, aiding in the study of toroidal and conical functions in potential problems.
These integral forms are particularly valuable for deriving asymptotic approximations in high-order or large-argument regimes.
Key Properties
Recurrence Relations and Differentiation Formulas
Legendre functions satisfy a variety of recurrence relations that enable the computation of higher-order functions from lower-order ones, facilitating numerical evaluation and analytical manipulations. For the Legendre polynomials of integer degree , the fundamental three-term recurrence relation is
with initial conditions and . This relation, derived from the hypergeometric representation or the differential equation, allows recursive generation of the polynomials and holds for .
Differentiation formulas for express the derivative in terms of the polynomials themselves. One such relation is
valid for , which follows from combining the recurrence with the Legendre differential equation. An equivalent form is .
For Legendre functions of general (non-integer) degree , a key relation involving derivatives is
This identity arises from logarithmic differentiation of the hypergeometric series expansion and is useful in asymptotic analysis and solutions to boundary value problems.[23]
The associated Legendre functions obey similar recurrences adjusted for the order . A prominent differentiation formula is
which generalizes the integer-degree case and supports derivations in spherical coordinate systems. This can be paired with the recurrence for computational stability.
Sums of Legendre functions are connected via the Christoffel-Darboux identity, a summation formula for orthogonal polynomials. For Legendre polynomials, it states
for , with the confluent form at involving derivatives: . This identity, fundamental for kernel representations and quadrature, extends to associated functions with appropriate weights.
Orthogonality and Completeness for Integer Orders
The Legendre polynomials for nonnegative integers form an orthogonal set on the interval with respect to the constant weight function 1. Specifically, their orthogonality relation is given by
where is the Kronecker delta, equal to 1 if and 0 otherwise.[7]
This property extends to the associated Legendre functions of integer order, where are integers. For fixed , the functions are orthogonal on with
The Legendre polynomials constitute a complete orthogonal basis for the Hilbert space equipped with the inner product . Consequently, any function admits a Fourier-Legendre series expansion
where the coefficients are , and the series converges to in the norm.[24]
Associated with this completeness is Parseval's identity, which equates the norm of to a sum over the squared coefficients:
$$
\int_{-1}^{1} |f(x)|^2 , dx = \sum_{n=0}^{\infty} \frac{2n+1}{2} |a_n|^2.[25]
In contrast, the Legendre functions of the second kind exhibit logarithmic singularities at the endpoints , precluding standard orthogonality relations on in the sense.[15]
Applications
Potential Theory and Multipole Expansions
In potential theory, Legendre functions play a fundamental role in describing axisymmetric solutions to Laplace's equation, particularly for problems involving gravitational or electrostatic potentials. Adrien-Marie Legendre introduced these functions in his 1782 memoir on the gravitational attraction of ellipsoids, where he developed methods to compute the potential due to homogeneous ellipsoidal mass distributions at exterior points, laying the groundwork for their use in celestial mechanics and geophysics. A key application arises in the expansion of the reciprocal distance between two points, which represents the fundamental solution to Laplace's equation in three dimensions. For points and with magnitudes and , and angle between them, the expansion is given by [ \frac{1}{|\mathbf{r} - \mathbf{r}'|} = \sum_{n=0}^{\infty} \frac{r_<^n}{r_>^{n+1}} P_n(\cos \gamma), $$ where and , and are the Legendre polynomials.[26] This series converges for and enables the decomposition of potentials from distributed sources into separable angular and radial components./04:_Series_Solutions/4.05:_Legendre_Polynomials) For axisymmetric boundary value problems, solutions to Laplace's equation in spherical coordinates take the form , assuming azimuthal independence.[27] The terms with positive powers of describe interior solutions (regular at the origin), while the inverse powers handle exterior solutions (vanishing at infinity).[28] Coefficients and are determined by boundary conditions, such as specified potentials on spheres or other surfaces of revolution. In electrostatics, this framework underpins the multipole expansion of the potential due to a localized charge distribution. The leading term for is the monopole, proportional to the total charge and scaling as , representing the net Coulomb potential.[29] The dipole term involves the dipole moment and falls as , capturing the potential asymmetry for neutral distributions with separated charges.[30] Higher-order terms, such as the quadrupole, describe further deviations and decay as , essential for precise modeling of molecular or planetary fields. Legendre series also facilitate numerical solutions to boundary value problems in potential theory, where the potential on a boundary is expanded as , and coefficients are found via orthogonality integrals to approximate the full solution inside or outside the domain. This approach is particularly effective for axisymmetric geometries, reducing the problem to solving for radial coefficients after angular decomposition. For non-axisymmetric cases, associated Legendre functions extend this method to include azimuthal dependence.Quantum Mechanics and Spherical Harmonics
In quantum mechanics, associated Legendre functions play a central role in describing the angular dependence of wave functions for systems with spherical symmetry, particularly through their incorporation into spherical harmonics. These functions arise as solutions to the angular part of the Schrödinger equation in spherical coordinates, where the separation of variables leads to the associated Legendre equation for the polar angle θ. The spherical harmonics $ Y_l^m(\theta, \phi) $, which form a complete orthonormal basis for functions on the sphere, are defined as
where $ P_l^m $ are the associated Legendre functions, $ l $ is the orbital angular momentum quantum number ($ l = 0, 1, 2, \dots $), and $ m $ is the magnetic quantum number ($ -l \leq m \leq l $). This normalization ensures the orthogonality relation $ \int Y_l^{m*}(\theta, \phi) Y_{l'}^{m'}(\theta, \phi) , d\Omega = \delta_{ll'} \delta_{mm'} $, with the integral over the solid angle $ d\Omega = \sin\theta , d\theta , d\phi $.[31]
A prime example is the hydrogen atom, where the total wave function $ \psi_{nlm}(r, \theta, \phi) = R_{nl}(r) Y_l^m(\theta, \phi) $ separates into radial and angular parts. The radial function $ R_{nl}(r) $ involves associated Laguerre polynomials, while the angular part relies on the associated Legendre functions via the spherical harmonics to capture the quantum states' dependence on direction. This structure reflects the conservation of angular momentum, with $ L^2 Y_l^m = \hbar^2 l(l+1) Y_l^m $ and $ L_z Y_l^m = \hbar m Y_l^m $, enabling the classification of atomic orbitals (e.g., s, p, d) by $ l $ and $ m \theta = 0, \pi $) for $ m \neq 0 $, maintaining physical boundary conditions.[32][33]
The addition theorem for spherical harmonics connects Legendre polynomials to multipole expansions in quantum contexts, such as scattering or perturbation theory:
where $ \cos \gamma $ is the angle between directions $ (\theta, \phi) $ and $ (\theta', \phi') $. This relation facilitates the decomposition of rotationally invariant operators and is crucial for calculating matrix elements in angular momentum coupling. The orthogonality of spherical harmonics directly implies selection rules for quantum transitions, such as electric dipole radiation, where $ \Delta l = \pm 1 $ and $ \Delta m = 0, \pm 1 $ arise from the integral $ \int Y_{l'}^{m'*} , \mathbf{r} , Y_l^m , d\Omega $ vanishing unless these conditions hold, governing allowed spectral lines in atomic spectra.[5][34][35]
Advanced Topics
Legendre Functions as Characters in Harmonic Analysis
In the representation theory of the Lie group SL(2,ℝ), Legendre functions arise as zonal spherical functions on the symmetric space SL(2,ℝ)/SO(2), which is isometric to the hyperbolic plane ℍ². Specifically, for the principal series representations parameterized by ν ∈ ℝ, the spherical function φ_ν(g) is given by P_{-1/2 + iν}(cosh t), where t is the hyperbolic distance from the identity coset e·SO(2) to g·SO(2), and P_μ denotes the Legendre function of the first kind. These functions satisfy the spherical function equation derived from the Casimir operator of the group algebra, confirming their role as the Legendre equation solutions in this context.[36] As characters, P_ν(cosh t) traces the irreducible unitary representations induced from the Borel subgroup, providing the radial part of the matrix coefficients in the K-biinvariant functions on the group.[37] In Fourier analysis on SL(2,ℝ)/SO(2), the Legendre transform serves as the Harish-Chandra transform for radial functions, decomposing them into spherical harmonics via integration against P_ν(cosh t). The Plancherel formula for this decomposition states that for a K-biinvariant function f on SL(2,ℝ),
where \hat{f}(ν) is the Legendre transform \hat{f}(ν) = \int_0^{\infty} f(t) P_{-1/2 + iν}(cosh t) \sinh t , dt, and the measure involves |P_ν|^2 implicitly through the inversion.[38] This formula arises from the decomposition of L²(SL(2,ℝ)/SO(2)) into principal series representations, with the continuous spectrum parameterized by ν.[38]
The connection to spherical functions on general symmetric spaces follows Harish-Chandra's integral formula, which for rank-one spaces like SL(2,ℝ)/SO(2) specializes to the explicit form involving Legendre functions. Harish-Chandra's c-function c(λ) = Γ(1/2 + iλ) Γ(1/2 - iλ) / [2π Γ(iλ) Γ(-iλ)] appears in the inversion, yielding the density |c(λ)|^{-2} in the Plancherel measure.[39] This framework extends Legendre functions to higher-rank semisimple groups, where they generalize to Jacobi functions, but retain the Legendre form for the SL(2,ℝ) case.[39]