In mathematics, the associated Legendre polynomials are the canonical solutions of the general Legendre equation

(1  x^2) \frac{d^2}{d x^2} P_\ell^m(x)  2 x \frac{d}{d x} P_\ell^m(x) + \left[ \ell (\ell + 1)  \frac{m^2}{1  x^2} \right] P_\ell^m(x) = 0,
or equivalently

\frac{d}{d x} \left[ (1  x^2) \frac{d}{d x} P_\ell^m(x) \right] + \left[ \ell (\ell + 1)  \frac{m^2}{1  x^2} \right] P_\ell^m(x) = 0,
where the indices ℓ and m (which are integers) are referred to as the degree and order of the associated Legendre polynomial respectively. This equation has nonzero solutions that are nonsingular on [−1, 1] only if ℓ and m are integers with 0 ≤ m ≤ ℓ, or with trivially equivalent negative values. When in addition m is even, the function is a polynomial. When m is zero and ℓ integer, these functions are identical to the Legendre polynomials. In general, when ℓ and m are integers, the regular solutions are sometimes called "associated Legendre polynomials", even though they are not polynomials when m is odd. The fully general class of functions with arbitrary real or complex values of ℓ and m are Legendre functions. In that case the parameters are usually labelled with Greek letters.
The Legendre ordinary differential equation is frequently encountered in physics and other technical fields. In particular, it occurs when solving Laplace's equation (and related partial differential equations) in spherical coordinates. Associated Legendre polynomials play a vital role in the definition of spherical harmonics.
Contents

Definition for nonnegative integer parameters ℓ and m 1

Alternative notations 1.1

Orthogonality 2

Negative m and/or negative ℓ 3

The first few associated Legendre functions 4

Recurrence formula 5

Gaunt's formula 6

Generalization via hypergeometric functions 7

Reparameterization in terms of angles 8

Applications in physics: spherical harmonics 9

Generalizations 10

See also 11

Notes and references 12

External links 13
Definition for nonnegative integer parameters ℓ and m
These functions are denoted P_\ell^{m}(x), where the superscript indicates the order, and not a power of P. Their most straightforward definition is in terms of derivatives of ordinary Legendre polynomials (m ≥ 0)

P_\ell^{m}(x) = (1)^m\ (1x^2)^{m/2}\ \frac{d^m}{dx^m}\left(P_\ell(x)\right)\,
The (−1)^{m} factor in this formula is known as the Condon–Shortley phase. Some authors omit it. The functions described by this equation satisfy the general Legendre differential equation with the indicated values of the parameters ℓ and m follows by differentiating m times the Legendre equation for P_{ℓ}:^{[1]}

(1x^2) \frac{d^2}{dx^2}P_\ell(x) 2x\frac{d}{dx}P_\ell(x)+ \ell(\ell+1)P_\ell(x) = 0.
Moreover, since by Rodrigues' formula,

P_\ell(x) = \frac{1}{2^\ell\,\ell!} \ \frac{d^\ell}{dx^\ell}\left[(x^21)^\ell\right],
the Pm
ℓ can be expressed in the form

P_\ell^{m}(x) = \frac{(1)^m}{2^\ell \ell!} (1x^2)^{m/2}\ \frac{d^{\ell+m}}{dx^{\ell+m}}(x^21)^\ell.
This equation allows extension of the range of m to: −ℓ ≤ m ≤ ℓ. The definitions of P_{ℓ}^{±m}, resulting from this expression by substitution of ±m, are proportional. Indeed, equate the coefficients of equal powers on the left and right hand side of

\frac{d^{\ellm}}{dx^{\ellm}} (x^21)^{\ell} = c_{lm} (1x^2)^m \frac{d^{\ell+m}}{dx^{\ell+m}}(x^21)^{\ell},
then it follows that the proportionality constant is

c_{lm} = (1)^m \frac{(\ellm)!}{(\ell+m)!} ,
so that

P^{m}_\ell(x) = (1)^m \frac{(\ellm)!}{(\ell+m)!} P^{m}_\ell(x).
Alternative notations
The following alternative notations are also used in literature:^{[2]}

P_{\ell m}(x) = (1)^m P_\ell^{m}(x)
Orthogonality
Assuming 0 ≤ m ≤ ℓ, they satisfy the orthogonality condition for fixed m:

\int_{1}^{1} P_k ^{m} P_\ell ^{m} dx = \frac{2 (\ell+m)!}{(2\ell+1)(\ellm)!}\ \delta _{k,\ell}
Where δ_{k, ℓ} is the Kronecker delta.
Also, they satisfy the orthogonality condition for fixed ℓ:

\int_{1}^{1} \frac{P_\ell ^{m} P_\ell ^{n}}{1x^2}dx = \begin{cases} 0 & \mbox{if } m\neq n \\ \frac{(\ell+m)!}{m(\ellm)!} & \mbox{if } m=n\neq0 \\ \infty & \mbox{if } m=n=0\end{cases}
Negative m and/or negative ℓ
The differential equation is clearly invariant under a change in sign of m.
The functions for negative m were shown above to be proportional to those of positive m:

P_\ell ^{m} = (1)^m \frac{(\ellm)!}{(\ell+m)!} P_\ell ^{m}
(This followed from the Rodrigues' formula definition. This definition also makes the various recurrence formulas work for positive or negative m.)
\textrm{If}\quad {\mid}m{\mid} > \ell\,\quad\mathrm{then}\quad P_\ell^{m} = 0.\,
The differential equation is also invariant under a change from ℓ to −ℓ − 1, and the functions for negative ℓ are defined by

P_{\ell} ^{m} = P_{\ell1} ^{m},\ (\ell=1,\,2,\, ...).
The first few associated Legendre functions
The first few associated Legendre functions, including those for negative values of m, are:

P_{0}^{0}(x)=1

P_{1}^{1}(x)=\begin{matrix}\frac{1}{2}\end{matrix}P_{1}^{1}(x)

P_{1}^{0}(x)=x

P_{1}^{1}(x)=(1x^2)^{1/2}

P_{2}^{2}(x)=\begin{matrix}\frac{1}{24}\end{matrix}P_{2}^{2}(x)

P_{2}^{1}(x)=\begin{matrix}\frac{1}{6}\end{matrix}P_{2}^{1}(x)

P_{2}^{0}(x)=\begin{matrix}\frac{1}{2}\end{matrix}(3x^{2}1)

P_{2}^{1}(x)=3x(1x^2)^{1/2}

P_{2}^{2}(x)=3(1x^2)

P_{3}^{3}(x)=\begin{matrix}\frac{1}{720}\end{matrix}P_{3}^{3}(x)

P_{3}^{2}(x)=\begin{matrix}\frac{1}{120}\end{matrix}P_{3}^{2}(x)

P_{3}^{1}(x)=\begin{matrix}\frac{1}{12}\end{matrix}P_{3}^{1}(x)

P_{3}^{0}(x)=\begin{matrix}\frac{1}{2}\end{matrix}(5x^33x)

P_{3}^{1}(x)=\begin{matrix}\frac{3}{2}\end{matrix}(5x^{2}1)(1x^2)^{1/2}

P_{3}^{2}(x)=15x(1x^2)

P_{3}^{3}(x)=15(1x^2)^{3/2}

P_{4}^{4}(x)=\begin{matrix}\frac{1}{40320}\end{matrix}P_{4}^{4}(x)

P_{4}^{3}(x)=\begin{matrix}\frac{1}{5040}\end{matrix}P_{4}^{3}(x)

P_{4}^{2}(x)=\begin{matrix}\frac{1}{360}\end{matrix}P_{4}^{2}(x)

P_{4}^{1}(x)=\begin{matrix}\frac{1}{20}\end{matrix}P_{4}^{1}(x)

P_{4}^{0}(x)=\begin{matrix}\frac{1}{8}\end{matrix}(35x^{4}30x^{2}+3)

P_{4}^{1}(x)=\begin{matrix}\frac{5}{2}\end{matrix}(7x^33x)(1x^2)^{1/2}

P_{4}^{2}(x)=\begin{matrix}\frac{15}{2}\end{matrix}(7x^21)(1x^2)

P_{4}^{3}(x)=  105x(1x^2)^{3/2}

P_{4}^{4}(x)=105(1x^2)^{2}
Recurrence formula
These functions have a number of recurrence properties:

(\ellm+1)P_{\ell+1}^{m}(x) = (2\ell+1)xP_{\ell}^{m}(x)  (\ell+m)P_{\ell1}^{m}(x)

2mxP_{\ell}^{m}(x)=\sqrt{1x^2}\left[P_{\ell}^{m+1}(x)+(\ell+m)(\ellm+1)P_{\ell}^{m1}(x)\right]

\frac{1}{\sqrt{1x^2}}P_\ell^m(x) = \frac{1}{2m} \left[ P_{\ell1}^{m+1}(x) + (\ell+m1)(\ell+m)P_{\ell1}^{m1}(x) \right]

\frac{1}{\sqrt{1x^2}}P_\ell^m(x) = \frac{1}{2m} \left[ P_{\ell+1}^{m+1}(x) + (\ellm+1)(\ellm+2)P_{\ell+1}^{m1}(x) \right]

\sqrt{1x^2}P_\ell^m(x) = \frac1{2\ell+1} \left[ (\ellm+1)(\ellm+2) P_{\ell+1}^{m1}(x)  (\ell+m1)(\ell+m) P_{\ell1}^{m1}(x) \right]

\sqrt{1x^2}P_\ell^m(x) = \frac1{2\ell+1} \left[  P_{\ell+1}^{m+1}(x) + P_{\ell1}^{m+1}(x) \right]

\sqrt{1x^2}P_\ell^{m+1}(x) = (\ellm)xP_{\ell}^{m}(x)  (\ell+m)P_{\ell1}^{m}(x)

\sqrt{1x^2}P_\ell^{m+1}(x) = (\ellm+1)P_{\ell+1}^m(x)  (\ell+m+1)xP_\ell^m(x)

\sqrt{1x^2}\frac{d}{dx}{P_\ell^m}(x) = \frac12 \left[ (\ell+m)(\ellm+1)P_\ell^{m1}(x)  P_\ell^{m+1}(x) \right]

(1x^2)\frac{d}{dx}{P_\ell^m}(x) = \frac1{2\ell+1} \left[ (\ell+1)(\ell+m)P_{\ell1}^m(x)  \ell(\ellm+1)P_{\ell+1}^m(x) \right]

(x^21)\frac{d}{dx}{P_{\ell}^{m}}(x) = {\ell}xP_{\ell}^{m}(x)  (\ell+m)P_{\ell1}^{m}(x)

(x^21)\frac{d}{dx}{P_{\ell}^{m}}(x) = \sqrt{1x^2}P_{\ell}^{m+1}(x) + mxP_{\ell}^{m}(x)

(x^21)\frac{d}{dx}{P_{\ell}^{m}}(x) = (\ell+m)(\ellm+1)\sqrt{1x^2}P_{\ell}^{m1}(x)  mxP_{\ell}^{m}(x)
Helpful identities (initial values for the first recursion):

P_{\ell +1}^{\ell +1}(x) =  (2\ell+1) \sqrt{1x^2} P_{\ell}^{\ell}(x)

P_{\ell}^{\ell}(x) = (1)^l (2\ell1)!! (1 x^2)^{(l/2)}

P_{\ell +1}^{\ell}(x) = x (2\ell+1) P_{\ell}^{\ell}(x)
with !! the double factorial.
Gaunt's formula
The integral over the product of three associated Legendre polynomials (with orders matching as shown below) is a necessary ingredient when developing products of Legendre polynomials into a series linear in the Legendre polynomials. For instance, this turns out to be necessary when doing atomic calculations of the Hartree–Fock variety where matrix elements of the Coulomb operator are needed. For this we have Gaunt's formula ^{[3]}
\frac{1}{2} \int_{1}^1 P_l^u(x) P_m^v(x) P_n^w(x) dx =

(1)^{smw}\frac{(m+v)!(n+w)!(2s2n)!s!}{(mv)!(sl)!(sm)!(sn)!(2s+1)!}


\times \ \sum_{t=p}^q (1)^t \frac{(l+u+t)!(m+nut)!}{t!(lut)!(mn+u+t)!(nwt)!}

This formula is to be used under the following assumptions:

the degrees are nonnegative integers l,m,n\ge0

all three orders are nonnegative integers u,v,w\ge 0

u is the largest of the three orders

the orders sum up u=v+w

the degrees obey m\ge n
Other quantities appearing in the formula are defined as

\ 2s = l+m+n

\ p = \max(0,\,nmu)

\ q = \min(m+nu,\,lu,\,nw)
The integral is zero unless

the sum of degrees is even so that s is an integer

the triangular condition is satisfied m+n\ge l \ge mn
Of course, Dong and Lemus (2002),^{[4]} succeeded to present true evaluation of integrals over multiple product of arbitrary order Associated Legendre Polynomials.
Generalization via hypergeometric functions
These functions may actually be defined for general complex parameters and argument:

P_{\lambda}^{\mu}(z) = \frac{1}{\Gamma(1\mu)} \left[\frac{1+z}{1z}\right]^{\mu/2} \,_2F_1 (\lambda, \lambda+1; 1\mu; \frac{1z}{2})
where \Gamma is the gamma function and _2F_1 is the hypergeometric function

\,_2F_1 (\alpha, \beta; \gamma; z) = \frac{\Gamma(\gamma)}{\Gamma(\alpha)\Gamma(\beta)} \sum_{n=0}^\infty\frac{\Gamma(n+\alpha)\Gamma(n+\beta)}{\Gamma(n+\gamma)\ n!}z^n,
They are called the Legendre functions when defined in this more general way. They satisfy the same differential equation as before:

(1z^2)\,y'' 2zy' + \left(\lambda[\lambda+1]  \frac{\mu^2}{1z^2}\right)\,y = 0.\,
Since this is a second order differential equation, it has a second solution, Q_\lambda^{\mu}(z), defined as:

Q_{\lambda}^{\mu}(z) = \frac{\sqrt{\pi}\ \Gamma(\lambda+\mu+1)}{2^{\lambda+1}\Gamma(\lambda+3/2)}\frac{1}{z^{\lambda+\mu+1}}(1z^2)^{\mu/2} \,_2F_1 \left(\frac{\lambda+\mu+1}{2}, \frac{\lambda+\mu+2}{2}; \lambda+\frac{3}{2}; \frac{1}{z^2}\right)
P_\lambda^{\mu}(z) and Q_\lambda^{\mu}(z) both obey the various recurrence formulas given previously.
Reparameterization in terms of angles
These functions are most useful when the argument is reparameterized in terms of angles, letting x = \cos\theta:

P_\ell^{m}(\cos\theta) = (1)^m (\sin \theta)^m\ \frac{d^m}{d(\cos\theta)^m}\left(P_\ell(\cos\theta)\right)\,
The first few polynomials, parameterized this way, are:

\begin{align} P_0^0(\cos\theta) & = 1 \\[8pt] P_1^0(\cos\theta) & = \cos\theta \\[8pt] P_1^1(\cos\theta) & = \sin\theta \\[8pt] P_2^0(\cos\theta) & = \tfrac{1}{2} (3\cos^2\theta1) \\[8pt] P_2^1(\cos\theta) & = 3\cos\theta\sin\theta \\[8pt] P_2^2(\cos\theta) & = 3\sin^2\theta \\[8pt] P_3^0(\cos\theta) & = \tfrac{1}{2} (5\cos^3\theta3\cos\theta) \\[8pt] P_3^1(\cos\theta) & = \tfrac{3}{2} (5\cos^2\theta1)\sin\theta \\[8pt] P_3^2(\cos\theta) & = 15\cos\theta\sin^2\theta \\[8pt] P_3^3(\cos\theta) & = 15\sin^3\theta \\[8pt] P_4^0(\cos\theta) & = \tfrac{1}{8} (35\cos^4\theta30\cos^2\theta+3) \\[8pt] P_4^1(\cos\theta) & =  \tfrac{5}{2} (7\cos^3\theta3\cos\theta)\sin\theta \\[8pt] P_4^2(\cos\theta) & = \tfrac{15}{2} (7\cos^2\theta1)\sin^2\theta \\[8pt] P_4^3(\cos\theta) & = 105\cos\theta\sin^3\theta \\[8pt] P_4^4(\cos\theta) & = 105\sin^4\theta \end{align}
For fixed m, P_\ell^m(\cos\theta) are orthogonal, parameterized by θ over [0, \pi], with weight \sin \theta:

\int_0^\pi P_k^{m}(\cos\theta) P_\ell^{m}(\cos\theta)\,\sin\theta\,d\theta = \frac{2 (\ell+m)!}{(2\ell+1)(\ellm)!}\ \delta _{k,\ell}
Also, for fixed ℓ:

\int_0^\pi P_\ell^{m}(\cos\theta) P_\ell^{n}(\cos\theta) \csc\theta\,d\theta = \begin{cases} 0 & \text{if } m\neq n \\ \frac{(\ell+m)!}{m(\ellm)!} & \text{if } m=n\neq0 \\ \infty & \text{if } m=n=0\end{cases}
In terms of θ, P_\ell^{m}(\cos \theta) are solutions of

\frac{d^{2}y}{d\theta^2} + \cot \theta \frac{dy}{d\theta} + \left[\lambda  \frac{m^2}{\sin^2\theta}\right]\,y = 0\,
More precisely, given an integer m\ge0, the above equation has nonsingular solutions only when \lambda = \ell(\ell+1)\, for ℓ an integer ≥ m, and those solutions are proportional to P_\ell^{m}(\cos \theta).
Applications in physics: spherical harmonics
In many occasions in physics, associated Legendre polynomials in terms of angles occur where spherical symmetry is involved. The colatitude angle in spherical coordinates is the angle \theta used above. The longitude angle, \phi, appears in a multiplying factor. Together, they make a set of functions called spherical harmonics. These functions express the symmetry of the twosphere under the action of the Lie group SO(3).
What makes these functions useful is that they are central to the solution of the equation \nabla^2\psi + \lambda\psi = 0 on the surface of a sphere. In spherical coordinates θ (colatitude) and φ (longitude), the Laplacian is

\nabla^2\psi = \frac{\partial^2\psi}{\partial\theta^2} + \cot \theta \frac{\partial \psi}{\partial \theta} + \csc^2 \theta\frac{\partial^2\psi}{\partial\phi^2}.
When the partial differential equation

\frac{\partial^2\psi}{\partial\theta^2} + \cot \theta \frac{\partial \psi}{\partial \theta} + \csc^2 \theta\frac{\partial^2\psi}{\partial\phi^2} + \lambda \psi = 0
is solved by the method of separation of variables, one gets a φdependent part \sin(m\phi) or \cos(m\phi) for integer m≥0, and an equation for the θdependent part

\frac{d^{2}y}{d\theta^2} + \cot \theta \frac{dy}{d\theta} + \left[\lambda  \frac{m^2}{\sin^2\theta}\right]\,y = 0\,
for which the solutions are P_\ell^{m}(\cos \theta) with \ell{\ge}m and \lambda = \ell(\ell+1).
Therefore, the equation

\nabla^2\psi + \lambda\psi = 0
has nonsingular separated solutions only when \lambda = \ell(\ell+1), and those solutions are proportional to

P_\ell^{m}(\cos \theta)\ \cos (m\phi)\ \ \ \ 0 \le m \le \ell
and

P_\ell^{m}(\cos \theta)\ \sin (m\phi)\ \ \ \ 0 < m \le \ell.
For each choice of ℓ, there are 2ℓ + 1 functions for the various values of m and choices of sine and cosine. They are all orthogonal in both ℓ and m when integrated over the surface of the sphere.
The solutions are usually written in terms of complex exponentials:

Y_{\ell, m}(\theta, \phi) = \sqrt{\frac{(2\ell+1)(\ellm)!}{4\pi(\ell+m)!}}\ P_\ell^{m}(\cos \theta)\ e^{im\phi}\qquad \ell \le m \le \ell.
The functions Y_{\ell, m}(\theta, \phi) are the spherical harmonics, and the quantity in the square root is a normalizing factor. Recalling the relation between the associated Legendre functions of positive and negative m, it is easily shown that the spherical harmonics satisfy the identity^{[5]}

Y_{\ell, m}^*(\theta, \phi) = (1)^m Y_{\ell, m}(\theta, \phi).
The spherical harmonic functions form a complete orthonormal set of functions in the sense of Fourier series. It should be noted that workers in the fields of geodesy, geomagnetism and spectral analysis use a different phase and normalization factor than given here (see spherical harmonics).
When a 3dimensional spherically symmetric partial differential equation is solved by the method of separation of variables in spherical coordinates, the part that remains after removal of the radial part is typically of the form

\nabla^2\psi(\theta, \phi) + \lambda\psi(\theta, \phi) = 0,
and hence the solutions are spherical harmonics.
Generalizations
The Legendre polynomials are closely related to hypergeometric series. In the form of spherical harmonics, they express the symmetry of the twosphere under the action of the Lie group SO(3). There are many other Lie groups besides SO(3), and an analogous generalization of the Legendre polynomials exist to express the symmetries of semisimple Lie groups and Riemannian symmetric spaces. Crudely speaking, one may define a Laplacian on symmetric spaces; the eigenfunctions of the Laplacian can be thought of as generalizations of the spherical harmonics to other settings.
See also
Notes and references

^ Courant & Hilbert 1953, V, §10.

^ Abramowitz, Milton; Stegun, Irene A., eds. (1965), "Chapter 8", .

^ From John C. Slater Quantum Theory of Atomic Structure, McGrawHill (New York, 1960), Volume I, page 309, which cites the original work of J. A. Gaunt, Philosophical Transactions of the Royal Society of London, A228:151 (1929)

^ Dong S.H., Lemus R., (2002), The overlap integral of three associated Legendre polynomials. Appl. Math. Lett. 15, 541546.

^ This identity can also be shown by relating the spherical harmonics to Wigner Dmatrices and use of the timereversal property of the latter. The relation between associated Legendre functions of ±m can then be proved from the complex conjugation identity of the spherical harmonics.

Arfken, G.B.; Weber, H.J. (2001), Mathematical methods for physicists, Academic Press, ; Section 12.5. (Uses a different sign convention.)

Belousov, S. L. (1962), Tables of normalized associated Legendre polynomials, Mathematical tables 18, Pergamon Press .

Condon, E. U.; Shortley, G. H. (1970), The Theory of Atomic Spectra, Cambridge, England: Cambridge University Press, ; Chapter 3.

.

Dunster, T. M. (2010), "Legendre and Related Functions", in

Edmonds, A.R. (1957), Angular Momentum in Quantum Mechanics, Princeton University Press, ; Chapter 2.

.

Koornwinder, Tom H.; Wong, Roderick S. C.; Koekoek, Roelof; Swarttouw, René F. (2010), "Orthogonal Polynomials", in

Schach, S. R. (1973) New Identities for Legendre Associated Functions of Integral Order and Degree , Society for Industrial and Applied Mathematics Journal on Mathematical Analysis, 1976, Vol. 7, No. 1 : pp. 59–69
External links

Associated Legendre polynomials in MathWorld

Legendre polynomials in MathWorld
This article was sourced from Creative Commons AttributionShareAlike License; additional terms may apply. World Heritage Encyclopedia content is assembled from numerous content providers, Open Access Publishing, and in compliance with The Fair Access to Science and Technology Research Act (FASTR), Wikimedia Foundation, Inc., Public Library of Science, The Encyclopedia of Life, Open Book Publishers (OBP), PubMed, U.S. National Library of Medicine, National Center for Biotechnology Information, U.S. National Library of Medicine, National Institutes of Health (NIH), U.S. Department of Health & Human Services, and USA.gov, which sources content from all federal, state, local, tribal, and territorial government publication portals (.gov, .mil, .edu). Funding for USA.gov and content contributors is made possible from the U.S. Congress, EGovernment Act of 2002.
Crowd sourced content that is contributed to World Heritage Encyclopedia is peer reviewed and edited by our editorial staff to ensure quality scholarly research articles.
By using this site, you agree to the Terms of Use and Privacy Policy. World Heritage Encyclopedia™ is a registered trademark of the World Public Library Association, a nonprofit organization.