Legendre differential equation arises from the study of axisymmetric solutions of Laplace's equation in R3. Specifically, let the spherical coordinate be configured as
x=ρsinθcosϕ,y=ρsinθsinϕ,z=ρcosθ.
Then Laplace's equation Δu=0 becomes
ρ−2∂ρ(ρ2uρ)+ρ−2csc2θuϕϕ+ρ−2cscθ∂θ(sinθuθ)=0.(1)
Let u be symmetric about the z-axis, so it does not depend on ϕ, which means we should investigate (1) by the separation of variables u=f(ρ)g(θ):
gρ−2∂ρ(ρ2fρ)+fρ−2cscθ∂θ(sinθgθ)=0.
Rearranging the terms gives
f∂ρ(ρ2fρ)=−gcscθ∂θ(sinθgθ).(2)
Because the left-hand side does not depend on θ and the right-hand side does not depend on ρ, both sides of (2) should be equal to some constant λ. As a result, we convert Laplace's equation into two ordinary differential equations:
ρ2fρρ+2ρfρ=λf,(3)
−cscθ∂θ(sinθgθ)=λg.(4)
The radial equation (3) is a Cauchy-Euler equation, whose general solution is given by f=Aρm+Bρ−m−1, where m∈[0,+∞) is chosen to satisfy λ=m(m+1). Because the solution needs to be regular at the origin, we must have B=0, so f=Aρm.
More interesting is the polar equation (4). By making the substitution μ=cosθ, we have ∂θ=−sinθ∂μ, so (4) becomes
Lg:=−dμd[(1−μ2)dμdg]=λg,(5)
which is known as the Legendre differential equation. Because g(θ) needs to be 2π-periodic and continuous on [0,2π], g(μ) needs to be continuous on [−1,1]. In this article, we prove that under these conditions, the admissible eigenvalues λ of the operator L in (5) are of the form ℓ(ℓ+1) with ℓ∈Z.
Application of the Frobenius method
Let an be a sequence of numbers defined on non-negative integers and assume g(μ) to be in the form of
g(μ)=n≥0∑anμn+α.(6)
Then we have
dμdg=n≥0∑(n+α)anμn+α−1=αa0μα−1+(α+1)a1μα+n≥0∑(n+α+2)an+2μn+α+1.
Consequently, (5) becomes
α(α−1)a0μα−2+α(α+1)a1μα−1+n≥0∑[(n+α+1)(n+α+2)an+2−(n+α+1)(n+α)an+λ]μn+α=0.
Comparing coefficients gives α=0, so we obtain a recursive relation for an:
an+2=(n+1)(n+2)n(n+1)−λan.(7)
Asymptotic analysis of an
Let ν∈{0,1}, so substituting n=2k+ν−2 into (7) gives
a2k+ν={1−2k+ν2+O(k21)}a2(k−1)+ν={1−k1+O(k21)}a2(k−1)+ν.
Because log(1+z)=z+O(z2) for small z, there exists some k0 such that whenever k>k0,
1−k1+O(k21)=exp{−k1+O(k21)}.
Consequently, we obtain
a2k+ν=exp{−k0<r≤k∑r1+O(k0<r≤k∑r21)}a2k0+ν=exp{−logk+O(1)}a2k0+ν.
This indicates that if a2k0+ν=0, then there exists some constant A>0 such that
ke−A<a2k0+νa2k+ν<keA.(8)
Setting α=0 in (6), we have
g(μ)=k≥0∑a2kμ2k+k≥0∑a2k+1μ2k+1.(9)
Because a2k+ν is a multiple of aν, we conclude that the odd and even components in (9) are exactly the fundamental solutions of (5), and we can adjust the values of a0 and a1 to represent all general solutions.
The spectrum of the Legendre differential operator
According to (8) and the comparison test, g(μ) is regular when −1<μ<1. When μ=±1, we see that the fundamental solution
k≥0∑a2k+νμ2k+ν
will diverge if a2k+ν is nonzero for all sufficiently large k. Consequently, for (5) to have meaningful solutions, at least of the sequences a2k,a2k+1 needs to vanish at all large k.
From the properties of the recursive relation (7), this is possible if and only if λ=ℓ(ℓ+1) for some integer ℓ. In this situation, the meaningful solutions are exactly the constant multiples of one fundamental solution satisfying ν≡ℓ(mod2).
Legendre polynomials
Let λ=ℓ(ℓ+1) for some non-negative integer ℓ and ν∈{0,1} be such that ℓ−ν is even, so the sequence a2k+ν becomes zero after 2k+ν>ℓ, so it follows from (7) that for 0≤m≤(ℓ−ν)/2, there is
aℓ−2m=(ℓ−2m)(ℓ−2m+2)−ℓ(ℓ+1)(ℓ−2m+1)(ℓ−2m+2)aℓ−2m+2=(ℓ−2m)(−2m)+(−2m)(ℓ+1)(ℓ−2m+1)(ℓ−2m+2)aℓ−2m+2=(2ℓ+1−2m)(−2m)(ℓ−2m+1)(ℓ−2m+2)aℓ−2m+2=(2m)!!(−1)m⋅(2ℓ+1−2m)(2ℓ+3−2m)⋯(2ℓ−1)(ℓ−2m+1)(ℓ−2m+2)⋯(ℓ−1)ℓaℓ=(2m)!!(−1)m⋅(ℓ−2m)!(2ℓ−2m−1)!!⋅(2ℓ−1)!!ℓ!aℓ=(2m)!!(−1)m⋅(ℓ−2m)!(2ℓ−2m)!!(2ℓ−2m)!⋅(2ℓ−1)!!ℓ!aℓ=2ℓm!(−1)m⋅(ℓ−m)!(ℓ−2m)!(2ℓ−2m)!⋅(2ℓ−1)!!ℓ!aℓ.
Setting aℓ=(2ℓ−1)!!/ℓ! gives the Legendre polynomial Pℓ(μ) of order ℓ:
Pℓ(μ)=0≤m≤ℓ/2∑2ℓm!(−1)m(ℓ−m)!(ℓ−2m)!(2ℓ−2m)!μℓ−2m.(10)
Notice that
(ℓ−2m)!(2ℓ−2m)!=dμℓdℓ(μ2ℓ−2m),
so we have
Pℓ(μ)=dμℓdℓ0≤m≤ℓ/2∑2ℓm!(ℓ−m)!(−1)mμ2ℓ−2m=2ℓℓ!1dμℓdℓm=0∑ℓ(mℓ)(−1)m(μ2)ℓ−m.
Finally, applying the binomial theorem, we obtain a neat expression for Pℓ(μ):
Pℓ(μ)=2ℓℓ!1dμℓdℓ[(μ2−1)ℓ].(11)
This is known as Rodrigues' formula for Legendre polynomials.
Conclusion
In this article, we began our discussion by seeking axisymmetric solutions from Laplace's equation. By the separation of variables, we introduced Legendre's differential equation. Subsequently, we used the Frobenius method and some asymptotic analysis to demonstrate that the eigenvalues λ of the Legendre differential operator L are exactly in the form of ℓ(ℓ+1) for ℓ∈Z. Finally, we derived expressions for the solutions to Legendre's differential equation when λ=ℓ(ℓ+1).