logo

関連するルジャンドル微分方程式と多項式 📂微分方程式

関連するルジャンドル微分方程式と多項式

定義1

下の微分方程式を、関連ルジャンドル微分方程式という。

(1x2)d2ydx22xdydx+[+l(l+1)m21x2]y= 0orddx[(1x2)y]+[l(l+1)m21x2]y= 0 \begin{equation} \begin{aligned} &&(1-x^{2})\frac{ d^{2}y }{ dx^{2} }-2x \frac{dy}{dx}+\left[ +l(l+1)-\frac{m^{2}}{1-x^{2}} \right]y =&\ 0 \\ \mathrm{or} && \frac{ d }{ dx } \left[ (1-x^{2})y^{\prime} \right] +\left[ l(l+1)-\frac{m^{2}}{1-x^{2}} \right]y =&\ 0 \end{aligned} \label{1} \end{equation}

関連ルジャンドル微分方程式の解をPlm(x)P_{l}^{m}(x)として表記し、これを関連ルジャンドル多項式一般化されたルジャンドル多項式という。

Plm(x)=(1x2)m2dmdxmPl(x)=(1x2)m2dmdxm[12ll!dldxl(x21)l] \begin{align*} P_{l}^{m}(x)&= (1-x ^{2})^{\frac{|m|}{2}} \frac{ d^{|m|} }{ dx^{|m|} } P_{l}(x) \\ &=(1-x ^{2})^{\frac{|m|}{2}} \frac{ d^{|m|} }{ dx^{|m|} }\left[ \dfrac{1}{2^{l} l!} \dfrac{d^{l}}{dx^{l}}(x^2-1)^{l}\right] \end{align*}

ここで、Pl(x)P_{l}(x)ルジャンドル多項式だ。 mmの符号によって区分される場合

Plm(x)=(1x2)m212ll!dl+mdxl+m(x21)l P_{l}^{m}(x) = (1-x ^{2})^{\frac{m}{2}} \dfrac{1}{2^{l} l!} \dfrac{d^{l+m}}{dx^{l+m}}(x^2-1)^{l}

Plm=(1)m(lm)!(l+m)!Plm(x) P_{l}^{-m}=(-1)^{m}\frac{(l-m)!}{(l+m)!}P_{l}^{m}(x)

関連ルジャンドル多項式は、球座標系のラプラス方程式を解く際に出現する。ここで、定数llmmは、量子力学で量子数と関連している。

m=0m=0の場合は、ルジャンドル微分方程式だ。この場合の解を基に、m0m\ne 0の場合の解も見つけることができる。まず、関連ルジャンドル微分方程式の解は定数llmmによって決まるので、以下のように表記しよう。

y=Plm(x) y=P_{l}^{m}(x)

これを(1)\eqref{1}に代入して整理すると以下のようになる。

ddx[(1x2)dPlm(x)dx]+[l(l+1)m21x2]Plm(x)=0 \begin{equation} \frac{ d }{ dx }\left[ (1-x^{2})\frac{ d P_{l}^{m}(x)}{ dx } \right]+\left[ l(l+1)-\frac{m^{2}}{1-x^{2}} \right]P_{l}^{m}(x)=0 \label{2} \end{equation}

そして、解が以下の形であると仮定しよう。

Plm(x)=(1x2)m2u(x) P_{l}^{m}(x)=(1-x^{2})^{\frac{|m|}{2}}u(x)

xxを一度微分すると

dPlm(x)dx=mx(1x2)m21u(x)+(1x2)m2u(x) \frac{ d P_{l}^{m}(x)}{ d x }=-|m|x(1-x^{2})^{\frac{|m|}{2}-1}u(x)+(1-x^{2})^{\frac{|m|}{2}}u^{\prime}(x)

これを(2)\eqref{2}の最初の項に代入して整理すると以下のようになる。

ddx[(1x2)dPlm(x)dx]= ddx[mx(1x2)m2u(x)+(1x2)m2+1u(x)]= m(1x2)m2u(x)+m2x2(1x2)m21u(x)mx(1x2)m2u(x)(m+2)x(1x2)m2u(x)+(1x2)m2+1u(x)= (1x2)m2+1u(x)2(m+1)(1x2)m2u(x)[m(m+1)x2m](1x2)m21u(x) \begin{align*} \frac{ d }{ dx }\left[ (1-x^{2})\frac{ d P_{l}^{m}(x)}{ dx } \right] =&\ \frac{ d }{ dx }\left[ -|m|x(1-x^{2})^{\frac{|m|}{2}}u(x)+(1-x^{2})^{\frac{|m|}{2}+1}u^{\prime}(x) \right] \\ =&\ -|m|(1-x^{2})^{\frac{|m|}{2}}u(x)+|m|^{2}x^{2}(1-x^{2})^{\frac{|m|}{2}-1}u(x) \\ & -|m|x(1-x^{2})^{\frac{|m|}{2}}u^{\prime}(x)-(|m|+2)x(1-x^{2})^{\frac{|m|}{2}}u^{\prime}(x) \\ & +(1-x^{2})^{\frac{|m|}{2}+1}u^{\prime \prime}(x) \\ =&\ (1-x^{2})^{\frac{|m|}{2}+1}u^{\prime \prime}(x)-2(|m|+1)(1-x^{2})^{\frac{|m|}{2}}u^{\prime}(x) \\ & -[|m|(|m|+1)x^{2}-|m|] (1-x^{2})^{\frac{|m|}{2}-1}u(x) \end{align*}

両辺に1(1x2)m/2\dfrac{1}{(1-x^{2})^{|m|/2}}を掛けると次のようになる。

1(1x2)m/2ddx[(1x2)dPlm(x)dx]= (1x2)u(x)2(m+1)xu(x)[m(m+1)x2m](1x2)1u(x) \begin{align*} &\frac{1}{(1-x^{2})^{|m|/2}}\frac{ d }{ dx }\left[ (1-x^{2})\frac{ d P_{l}^{m}(x)}{ dx } \right] \\ =&\ (1-x^{2})u^{\prime \prime}(x)-2(|m|+1)xu^{\prime}(x) -[|m|(|m|+1)x^{2}-|m|] (1-x^{2})^{-1}u(x) \end{align*}

したがって、(2)\eqref{2}の両辺に1(1x2)m/2\dfrac{1}{(1-x^{2})^{|m|/2}}を掛けると

(1x2)u(x)2(m+1)xu(x)(m(m+1)x2m1x2+l(l+1)m21x2)u(x)=0 \begin{equation} \begin{aligned} &(1-x^{2})u^{\prime \prime}(x)-2(|m|+1)xu^{\prime}(x) \\ &-\left( \frac{|m|(|m|+1)x^{2}-|m|}{1-x^{2}}+l(l+1)-\frac{m^{2}}{1-x^{2}}\right)u(x)=0 \end{aligned} \label{1} \end{equation}

u(x)u(x)の係数を整理すると以下のようになる。

m(m+1)x2m1x2+l(l+1)m21x2= m(m+1)x2m+l(l+1)(1x2)m21x2= m2(1x2)m(1x2)+l(l+1)(1x2)1x2= l(l+1)m2m= l(l+1)m(m+1) \begin{align*} &\frac{|m|(|m|+1)x^{2}-|m|}{1-x^{2}}+l(l+1)-\frac{m^{2}}{1-x^{2}} \\ =&\ \frac{|m|(|m|+1)x^{2}-|m|+l(l+1)(1-x^{2})-m^{2}}{1-x^{2}} \\ =&\ \frac{-m^{2}(1-x^{2})-|m|(1-x^{2})+l(l+1)(1-x^{2})}{1-x^{2}} \\ =&\ l(l+1)-m^{2}-|m| \\ =&\ l(l+1)-|m|(|m|+1) \end{align*}

従って、(3)\eqref{3}は以下のような形で整理される。

(1x2)d2udx22(m+1)xdudx+[l(l+1)m(m+1)]u=0 \begin{equation} (1-x^{2})\frac{ d^{2} u }{ d x^{2} }-2(|m|+1)x\frac{ d u}{ dx }+[l(l+1)-|m|(|m|+1)]u=0 \label{4} \end{equation}

m=0m=0であれば、実際にルジャンドル微分方程式になる。従って、m=0|m|=0の場合の解はPl0(x)=Pl(x)P_{l}^{0}(x)=P_{l}(x)である。今一度、(4)(4)xxに関して微分してみよう。係数を整理すると以下の式を得る。

(1x2)d3udx32[(m+1)+1]xd2udx2+[l(l+1)(m+1)(m+2)]dudx=0 \begin{equation} (1-x^{2}) \frac{ d^{3} u }{ d x^{3} } -2[(|m|+1)+1]x\frac{ d^{2} u}{ dx^{2} }+[l(l+1)-(|m|+1)(|m|+2)]\frac{ d u}{ d x}=0 \label{5} \end{equation}

再び、(5)\eqref{5}xxに関して微分して係数を整理すると以下の式を得る。

(1x2)d4udx42[(m+2)+1]xd3udx3+[l(l+1)(m+2)(m+3)]d2udx2=0 \begin{equation} (1-x^{2}) \frac{ d^{4} u }{ d x^{4} } -2[(|m|+2)+1]x\frac{ d^{3} u}{ dx^{3} }+[l(l+1)-(|m|+2)(|m|+3)]\frac{ d^{2} u}{ d x^{2}}=0 \label{6} \end{equation}

これらの式をよく見ると、(4)\eqref{4}uuに、m|m|m+1|m|+1に置き換えたときに(5)\eqref{5}を得ることができることがわかる。(5)\eqref{5}で同様の方法で置換すると、(6)\eqref{6}を得る。これにより、m=0|m|=0の場合の解はPl(x)P_{l}(x)m=1|m|=1の場合の解はddxPl(x)\dfrac{ d }{ d x }P_{l}(x)m=2|m|=2の場合の解はd2dx2Pl(x)\dfrac{d^{2}}{dx^{2}}P_{l}(x)であることがわかる。したがって、これを一般化すると次のようになる。

u(x)=dmdxmPl(x) u(x)=\frac{d^{|m|}}{dx^{{|m|}}}P_{l}(x)

したがって、関連ルジャンドル多項式は以下のようになる。

Plm(x)=(1x2)m2dmdxmPl(x)=(1x2)m2dmdxm[12ll!dldxl(x21)l] \begin{align*} P_{l}^{m}(x)&= (1-x ^{2})^{\frac{|m|}{2}} \frac{ d^{|m|} }{ dx^{|m|} } P_{l}(x) \\ &=(1-x ^{2})^{\frac{|m|}{2}} \frac{ d^{|m|} }{ dx^{|m|} }\left[ \dfrac{1}{2^{l} l!} \dfrac{d^{l}}{dx^{l}}(x^2-1)^{l}\right] \end{align*}


  1. Mary L. Boas, 数理物理学(Mathematical Methods in the Physical Sciences, 最中峻輔 訳) (3rd Edition, 2008), p597-598 ↩︎