E. The Legendre equation

When solving the equation for the geometry in spherical coordinates, the following equation for the polar angle arises

ddz((1z2)dΘ(z)dz)\displaystyle \frac{\displaystyle d}{\displaystyle dz} \left( (1 - z^2) \cdot \frac{\displaystyle d \Theta(z)}{\displaystyle dz} \right)+[γ12m21z2]Θ(z)\displaystyle {} + \left[ \gamma_1^2 - \frac{\displaystyle m^2}{\displaystyle 1 - z^2} \right] \cdot \Theta(z)=0,\displaystyle {} = 0,1\displaystyle {} -1<z\displaystyle {} < z<1,\displaystyle {} < 1,
(E.1)

where Θ(z)\Theta(z) is the unknown function, mm is a nonnegative integer, γ12\gamma_1^2 is the separation constant. Let us expand the derivative in the first term

(1z2)d2Θ(z)dz2\displaystyle (1 - z^2) \cdot \frac{\displaystyle d^2 \Theta(z)}{\displaystyle dz^2}2zdΘ(z)dz\displaystyle {} - 2 \cdot z \cdot \frac{\displaystyle d \Theta(z)}{\displaystyle dz}+[γ12m21z2]Θ(z)\displaystyle {} + \left[ \gamma_1^2 - \frac{\displaystyle m^2}{\displaystyle 1 - z^2} \right] \cdot \Theta(z)=0.\displaystyle {} = 0.

Note that for mm=0{} = 0 the equation takes the form of the Legendre equation

(1z2)d2Θ(z)dz2\displaystyle (1 - z^2) \cdot \frac{\displaystyle d^2 \Theta(z)}{\displaystyle dz^2}2zdΘ(z)dz\displaystyle {} - 2 \cdot z \cdot \frac{\displaystyle d \Theta(z)}{\displaystyle dz}+γ12Θ(z)\displaystyle {} + \gamma_1^2 \cdot \Theta(z)=0,\displaystyle {} = 0,1\displaystyle {} -1<z\displaystyle {} < z<1.\displaystyle {} < 1.

This equation has solutions bounded on the interval 1{} -1<z{} < z<1{} < 1 only for the eigenvalues γ1k2\gamma_{1k}^2=k(k+1){} = k \cdot (k + 1), kk(0..){} \in (0..\infty); these solutions are given by the Rodrigues formula

Θk(z)\displaystyle \Theta_k(z) =Pk(z)\displaystyle {} = P_k(z)=12kk!dkdzk(z21)k,\displaystyle {} = \frac{\displaystyle 1}{\displaystyle 2^k \cdot k!} \cdot \frac{\displaystyle d^k}{\displaystyle dz^k} \left(z^2 - 1\right)^k,
(E.2)

where Pk(z)P_k(z) are the Legendre polynomials. We will look for solutions of equation (E.1) for the same eigenvalues γ1k2\gamma_{1k}^2=k(k+1){} = k \cdot (k + 1) — then it can be written as

(1z2)d2Θkm(z)dz2\displaystyle (1 - z^2) \cdot \frac{\displaystyle d^2 \Theta_{km}(z)}{\displaystyle dz^2}2zdΘkm(z)dz\displaystyle {} - 2 \cdot z \cdot \frac{\displaystyle d \Theta_{km}(z)}{\displaystyle dz}+[k(k+1)m21z2]\displaystyle {} + \left[ k \cdot (k + 1) - \frac{\displaystyle m^2}{\displaystyle 1 - z^2} \right]Θkm(z)\displaystyle {} \cdot \Theta_{km}(z)=0,\displaystyle {} = 0,k\displaystyle k(0..).\displaystyle {} \in (0..\infty).
(E.3)

Let us make the substitution Θkm(z)\Theta_{km}(z)=(1z2)m2Θ^km(z){} = \left( 1 - z^2 \right)^{\frac{\displaystyle m}{\displaystyle 2}} \cdot \widehat{\Theta}_{km}(z) into equation (E.3). The first derivative takes the form

ddz[(1z2)m2Θ^km(z)]\displaystyle \frac{\displaystyle d}{\displaystyle dz} \left[ (1 - z^2)^{\frac{\displaystyle m}{\displaystyle 2}} \cdot \widehat{\Theta}_{km}(z) \right] =(1z2)m2dΘ^km(z)dz\displaystyle {} = (1 - z^2)^{\frac{\displaystyle m}{\displaystyle 2}} \cdot \frac{\displaystyle d \widehat{\Theta}_{km}(z)}{\displaystyle dz}m\displaystyle {} - mz\displaystyle {} \cdot z(1z2)m21\displaystyle {} \cdot (1 - z^2)^{\frac{\displaystyle m}{\displaystyle 2} - 1}Θ^km(z).\displaystyle {} \cdot \widehat{\Theta}_{km}(z).

To compute the second derivative, we first compute

ddz[z1z2]\displaystyle \frac{\displaystyle d}{\displaystyle dz} \left[ \frac{\displaystyle z}{\displaystyle 1 - z^2} \right] =11z2\displaystyle {} = \frac{\displaystyle 1}{\displaystyle 1 - z^2}+2z2(1z2)2\displaystyle {} + \frac{\displaystyle 2 \cdot z^2}{\displaystyle (1 - z^2)^2}=1+z2(1z2)2.\displaystyle {} = \frac{\displaystyle 1 + z^2}{\displaystyle (1 - z^2)^2}.

Then the second derivative can be written as

d2dz2[(1z2)m2Θ^km(z)]\displaystyle \frac{\displaystyle d^2}{\displaystyle dz^2} \left[ (1 - z^2)^{\frac{\displaystyle m}{\displaystyle 2}} \cdot \widehat{\Theta}_{km}(z) \right] =(1z2)m2\displaystyle {} = (1 - z^2)^{\frac{\displaystyle m}{\displaystyle 2}}d2Θ^km(z)dz2\displaystyle {} \cdot \frac{\displaystyle d^2 \widehat{\Theta}_{km}(z)}{\displaystyle dz^2}m\displaystyle {} - mz\displaystyle {} \cdot z(1z2)m21\displaystyle {} \cdot (1 - z^2)^{\frac{\displaystyle m}{\displaystyle 2} - 1}dΘ^km(z)dz\displaystyle {} \cdot \frac{\displaystyle d \widehat{\Theta}_{km}(z)}{\displaystyle dz}mz1z2[(1z2)m2dΘ^km(z)dz\displaystyle {} - m \cdot \frac{\displaystyle z}{\displaystyle 1 - z^2} \cdot \biggl[ (1 - z^2)^{\frac{\displaystyle m}{\displaystyle 2}} \cdot \frac{\displaystyle d \widehat{\Theta}_{km}(z)}{\displaystyle dz}mz(1z2)m21Θ^km(z)]\displaystyle {} - m \cdot z \cdot (1 - z^2)^{\frac{\displaystyle m}{\displaystyle 2} - 1} \cdot \widehat{\Theta}_{km}(z) \biggr]m\displaystyle {} - m1+z2(1z2)2\displaystyle {} \cdot \frac{\displaystyle 1 + z^2}{\displaystyle (1 - z^2)^2}(1z2)m2\displaystyle {} \cdot (1 - z^2)^{\frac{\displaystyle m}{\displaystyle 2}}Θ^km(z).\displaystyle {} \cdot \widehat{\Theta}_{km}(z).

After collecting like terms it takes the form

d2dz2[(1z2)m2Θ^km(z)]\displaystyle \frac{\displaystyle d^2}{\displaystyle dz^2} \left[ (1 - z^2)^{\frac{\displaystyle m}{\displaystyle 2}} \cdot \widehat{\Theta}_{km}(z) \right] =(1z2)m2\displaystyle {} = (1 - z^2)^{\frac{\displaystyle m}{\displaystyle 2}}d2Θ^km(z)dz2\displaystyle {} \cdot \frac{\displaystyle d^2 \widehat{\Theta}_{km}(z)}{\displaystyle dz^2}2mz1z2\displaystyle {} - \frac{\displaystyle 2 \cdot m \cdot z}{\displaystyle 1 - z^2}(1z2)m2\displaystyle {} \cdot (1 - z^2)^{\frac{\displaystyle m}{\displaystyle 2}}dΘ^km(z)dz\displaystyle {} \cdot \frac{\displaystyle d \widehat{\Theta}_{km}(z)}{\displaystyle dz}+m\displaystyle {} + mmz21z2(1z2)2\displaystyle {} \cdot \frac{\displaystyle m \cdot z^2 - 1 - z^2}{\displaystyle (1 - z^2)^2}(1z2)m2\displaystyle {} \cdot (1 - z^2)^{\frac{\displaystyle m}{\displaystyle 2}}Θ^km(z).\displaystyle {} \cdot \widehat{\Theta}_{km}(z).

Let us substitute everything into equation (E.3), cancelling along the way (1z2)m2(1 - z^2)^{\frac{\displaystyle m}{\displaystyle 2}}:

(1z2)[d2Θ^km(z)dz2\displaystyle (1 - z^2) \cdot \biggl[ \frac{\displaystyle d^2 \widehat{\Theta}_{km}(z)}{\displaystyle dz^2}2mz1z2dΘ^km(z)dz\displaystyle {} - \frac{\displaystyle 2 \cdot m \cdot z}{\displaystyle 1 - z^2} \cdot \frac{\displaystyle d \widehat{\Theta}_{km}(z)}{\displaystyle dz}+mmz21z2(1z2)2Θ^km(z)]\displaystyle {} + m \cdot \frac{\displaystyle m \cdot z^2 - 1 - z^2}{\displaystyle (1 - z^2)^2} \cdot \widehat{\Theta}_{km}(z) \biggr]2z[dΘ^km(z)dz\displaystyle {} - 2 \cdot z \cdot \biggl[ \frac{\displaystyle d \widehat{\Theta}_{km}(z)}{\displaystyle dz}mz1z2Θ^km(z)]\displaystyle {} - \frac{\displaystyle m \cdot z}{\displaystyle 1 - z^2} \cdot \widehat{\Theta}_{km}(z) \biggr]+[k(k+1)m21z2]\displaystyle {} + \left[ k \cdot (k + 1) - \frac{\displaystyle m^2}{\displaystyle 1 - z^2} \right]Θ^km(z)\displaystyle {} \cdot \widehat{\Theta}_{km}(z)=0.\displaystyle {} = 0.

Let us group the coefficients of the derivatives:

ζ1d2Θ^km(z)dz2\displaystyle \zeta_1 \cdot \frac{\displaystyle d^2 \widehat{\Theta}_{km}(z)}{\displaystyle dz^2}+ζ2dΘ^km(z)dz\displaystyle {} + \zeta_2 \cdot \frac{\displaystyle d \widehat{\Theta}_{km}(z)}{\displaystyle dz}+ζ3Θ^km(z)\displaystyle {} + \zeta_3 \cdot \widehat{\Theta}_{km}(z)=0,\displaystyle {} = 0,
ζ1\displaystyle \zeta_1 =1\displaystyle {} = 1z2,\displaystyle {} - z^2,
ζ2\displaystyle \zeta_2 =2mz\displaystyle {} = -2 \cdot m \cdot z2z\displaystyle {} - 2 \cdot z=2(m+1)z,\displaystyle {} = -2 \cdot (m + 1) \cdot z,
ζ3\displaystyle \zeta_3 =mmz21z21z2\displaystyle {} = m \cdot \frac{\displaystyle m \cdot z^2 - 1 - z^2}{\displaystyle 1 - z^2}+2mz21z2\displaystyle {} + \frac{\displaystyle 2 \cdot m \cdot z^2}{\displaystyle 1 - z^2}+k(k+1)\displaystyle {} + k \cdot (k + 1)m21z2\displaystyle {} - \frac{\displaystyle m^2}{\displaystyle 1 - z^2}=m2z2mmz2+2mz2m21z2\displaystyle {} = \frac{\displaystyle m^2 \cdot z^2 - m - m \cdot z^2 + 2 \cdot m \cdot z^2 - m^2}{\displaystyle 1 - z^2}+k(k+1)\displaystyle {} + k \cdot (k + 1)=m2(1z2)+m(1z2)1z2\displaystyle {} = - \frac{\displaystyle m^2 \cdot (1 - z^2) + m \cdot (1 - z^2)}{\displaystyle 1 - z^2}+k(k+1)\displaystyle {} + k \cdot (k + 1)=m2\displaystyle {} = - m^2m\displaystyle {} - m+k(k+1)\displaystyle {} + k \cdot (k + 1)=(km)(k+m+1),\displaystyle {} = (k - m) \cdot (k + m + 1),

thus, after the substitution, equation (E.3) takes the form

(1z2)d2Θ^km(z)dz2\displaystyle (1 - z^2) \cdot \frac{\displaystyle d^2 \widehat{\Theta}_{km}(z)}{\displaystyle dz^2}2(m+1)zdΘ^km(z)dz\displaystyle {} - 2 \cdot (m + 1) \cdot z \cdot \frac{\displaystyle d \widehat{\Theta}_{km}(z)}{\displaystyle dz}+(km)(k+m+1)Θ^km(z)\displaystyle {} + (k - m) \cdot (k + m + 1) \cdot \widehat{\Theta}_{km}(z)=0.\displaystyle {} = 0.
(E.4)

Let us take equation (E.3), set mm=0{} = 0 and take into account that in this case its solutions are the Legendre polynomials Pk(z)P_k(z) — we obtain the following equation

(1z2)d2Pk(z)dz2\displaystyle (1 - z^2) \cdot \frac{\displaystyle d^2 P_k(z)}{\displaystyle dz^2}2zdPk(z)dz\displaystyle {} - 2 \cdot z \cdot \frac{\displaystyle d P_k(z)}{\displaystyle dz}+k(k+1)Pk(z)\displaystyle {} + k \cdot (k + 1) \cdot P_k(z)=0,\displaystyle {} = 0,
(E.5)

which we differentiate mm times using the Leibniz formula, defined as follows

dmdzm(u(z)v(z))\displaystyle \frac{\displaystyle d^m}{\displaystyle dz^m} \left( u(z) \cdot v(z) \right) =i=0m(mi)\displaystyle {} = \sum_{i=0}^m \binom{m}{i}dmidzmiu(z)\displaystyle {} \cdot \frac{\displaystyle d^{m-i}}{\displaystyle dz^{m-i}} u(z)didziv(z),\displaystyle {} \cdot \frac{\displaystyle d^i}{\displaystyle dz^i} v(z),

where (mi)\binom{m}{i}=m!i!(mi)!{} = \frac{\displaystyle m!}{\displaystyle i! \cdot (m - i)!} is the binomial coefficient.

Let us start with the first term of equation (E.5), taking v(z)v(z)=1{} = 1z2{} - z^2

d0dz0v\displaystyle \frac{\displaystyle d^0}{\displaystyle dz^0} v =1\displaystyle {} = 1z2;d1dz1v\displaystyle {} - z^2; \quad \frac{\displaystyle d^1}{\displaystyle dz^1} v=2z;d2dz2v\displaystyle {} = -2 \cdot z; \quad \frac{\displaystyle d^2}{\displaystyle dz^2} v=2;didziv\displaystyle {} = -2; \quad \frac{\displaystyle d^i}{\displaystyle dz^i} v=0,\displaystyle {} = 0,i\displaystyle i(3..),\displaystyle {} \in (3..\infty),

now let us take v(z)v(z)=2z{} = -2 \cdot z, we obtain

d0dz0v\displaystyle \frac{\displaystyle d^0}{\displaystyle dz^0} v =2z;d1dz1v\displaystyle {} = -2 \cdot z; \quad \frac{\displaystyle d^1}{\displaystyle dz^1} v=2;didziv\displaystyle {} = -2; \quad \frac{\displaystyle d^i}{\displaystyle dz^i} v=0,\displaystyle {} = 0,i\displaystyle i(2..),\displaystyle {} \in (2..\infty),

for the case v(z)v(z)=k(k+1){} = k \cdot (k + 1) everything is trivial

d0dz0v\displaystyle \frac{\displaystyle d^0}{\displaystyle dz^0} v =k(k+1);didziv\displaystyle {} = k \cdot (k + 1); \quad \frac{\displaystyle d^i}{\displaystyle dz^i} v=0,\displaystyle {} = 0,i\displaystyle i(1..).\displaystyle {} \in (1..\infty).

The binomial coefficients are respectively equal to

(m0)\displaystyle \binom{m}{0} =1;(m1)\displaystyle {} = 1; \quad \binom{m}{1}=m;(m2)\displaystyle {} = m; \quad \binom{m}{2}=m(m1)2.\displaystyle {} = \frac{\displaystyle m \cdot (m - 1)}{\displaystyle 2}.

The result of differentiating the first term of equation (E.5) equals

(1z2)dm+2dzm+2Pk(z)\displaystyle (1 - z^2) \cdot \frac{\displaystyle d^{m+2}}{\displaystyle dz^{m+2}} P_k(z)2zmdm+1dzm+1Pk(z)\displaystyle {} - 2 \cdot z \cdot m \cdot \frac{\displaystyle d^{m+1}}{\displaystyle dz^{m+1}} P_k(z)m(m1)dmdzmPk(z),\displaystyle {} - m \cdot (m - 1) \cdot \frac{\displaystyle d^m}{\displaystyle dz^m} P_k(z),

for the second term

2zdm+1dzm+1Pk(z)\displaystyle {} -2 \cdot z \cdot \frac{\displaystyle d^{m+1}}{\displaystyle dz^{m+1}} P_k(z)2mdmdzmPk(z),\displaystyle {} - 2 \cdot m \cdot \frac{\displaystyle d^m}{\displaystyle dz^m} P_k(z),

for the third term

k(k+1)dmdzmPk(z),k \cdot (k + 1) \cdot \frac{\displaystyle d^m}{\displaystyle dz^m} P_k(z),

let us add all three terms

(1z2)dm+2dzm+2Pk(z)\displaystyle (1 - z^2) \cdot \frac{\displaystyle d^{m+2}}{\displaystyle dz^{m+2}} P_k(z)2\displaystyle {} - 2(m+1)\displaystyle {} \cdot (m + 1)z\displaystyle {} \cdot zdm+1dzm+1Pk(z)\displaystyle {} \cdot \frac{\displaystyle d^{m+1}}{\displaystyle dz^{m+1}} P_k(z)+(km)\displaystyle {} + (k-m)(k+m+1)\displaystyle {} \cdot (k + m + 1)dmdzmPk(z).\displaystyle {} \cdot \frac{\displaystyle d^m}{\displaystyle dz^m} P_k(z).

The left-hand side of equation (E.5) is identically zero, hence so is its mm-fold derivative:

(1z2)dm+2dzm+2Pk(z)\displaystyle (1 - z^2) \cdot \frac{\displaystyle d^{m+2}}{\displaystyle dz^{m+2}} P_k(z)2\displaystyle {} - 2(m+1)\displaystyle {} \cdot (m + 1)z\displaystyle {} \cdot zdm+1dzm+1Pk(z)\displaystyle {} \cdot \frac{\displaystyle d^{m+1}}{\displaystyle dz^{m+1}} P_k(z)+(km)\displaystyle {} + (k-m)(k+m+1)\displaystyle {} \cdot (k + m + 1)dmdzmPk(z)\displaystyle {} \cdot \frac{\displaystyle d^m}{\displaystyle dz^m} P_k(z)=0.\displaystyle {} = 0.

The resulting equation coincides with equation (E.4), whence Θ^km(z)\widehat{\Theta}_{km}(z)=dmdzmPk(z){} = \frac{\displaystyle d^m}{\displaystyle dz^m} P_k(z), which means the solution of equation (E.1) has the form

Θkm(z)\displaystyle \Theta_{km}(z) =Pk(m)(z)\displaystyle {} = P_k^{(m)}(z)=(1z2)m2dmdzmPk(z),\displaystyle {} = (1 - z^2)^{\frac{\displaystyle m}{\displaystyle 2}} \cdot \frac{\displaystyle d^m}{\displaystyle dz^m} P_k(z),
(E.6)

where γ1k2\gamma_{1k}^2=k(k+1){} = k \cdot (k + 1) are the eigenvalues, Pk(m)(z)P_k^{(m)}(z) are the associated Legendre polynomials. For mm>k{} > k the derivative dmdzmPk(z)\frac{\displaystyle d^m}{\displaystyle dz^m} P_k(z) vanishes, so nontrivial solutions exist only for kkm{} \ge m.