関連するルジャンドル微分方程式と多項式
📂微分方程式関連するルジャンドル微分方程式と多項式
定義
下の微分方程式を、関連ルジャンドル微分方程式という。
or(1−x2)dx2d2y−2xdxdy+[+l(l+1)−1−x2m2]y=dxd[(1−x2)y′]+[l(l+1)−1−x2m2]y= 0 0
関連ルジャンドル微分方程式の解をPlm(x)として表記し、これを関連ルジャンドル多項式や一般化されたルジャンドル多項式という。
Plm(x)=(1−x2)2∣m∣dx∣m∣d∣m∣Pl(x)=(1−x2)2∣m∣dx∣m∣d∣m∣[2ll!1dxldl(x2−1)l]
ここで、Pl(x)はルジャンドル多項式だ。 mの符号によって区分される場合、
Plm(x)=(1−x2)2m2ll!1dxl+mdl+m(x2−1)l
Pl−m=(−1)m(l+m)!(l−m)!Plm(x)
関連ルジャンドル多項式は、球座標系のラプラス方程式を解く際に出現する。ここで、定数l、mは、量子力学で量子数と関連している。
解
m=0の場合は、ルジャンドル微分方程式だ。この場合の解を基に、m=0の場合の解も見つけることができる。まず、関連ルジャンドル微分方程式の解は定数l、mによって決まるので、以下のように表記しよう。
y=Plm(x)
これを(1)に代入して整理すると以下のようになる。
dxd[(1−x2)dxdPlm(x)]+[l(l+1)−1−x2m2]Plm(x)=0
そして、解が以下の形であると仮定しよう。
Plm(x)=(1−x2)2∣m∣u(x)
xを一度微分すると
dxdPlm(x)=−∣m∣x(1−x2)2∣m∣−1u(x)+(1−x2)2∣m∣u′(x)
これを(2)の最初の項に代入して整理すると以下のようになる。
dxd[(1−x2)dxdPlm(x)]=== dxd[−∣m∣x(1−x2)2∣m∣u(x)+(1−x2)2∣m∣+1u′(x)] −∣m∣(1−x2)2∣m∣u(x)+∣m∣2x2(1−x2)2∣m∣−1u(x)−∣m∣x(1−x2)2∣m∣u′(x)−(∣m∣+2)x(1−x2)2∣m∣u′(x)+(1−x2)2∣m∣+1u′′(x) (1−x2)2∣m∣+1u′′(x)−2(∣m∣+1)(1−x2)2∣m∣u′(x)−[∣m∣(∣m∣+1)x2−∣m∣](1−x2)2∣m∣−1u(x)
両辺に(1−x2)∣m∣/21を掛けると次のようになる。
=(1−x2)∣m∣/21dxd[(1−x2)dxdPlm(x)] (1−x2)u′′(x)−2(∣m∣+1)xu′(x)−[∣m∣(∣m∣+1)x2−∣m∣](1−x2)−1u(x)
したがって、(2)の両辺に(1−x2)∣m∣/21を掛けると
(1−x2)u′′(x)−2(∣m∣+1)xu′(x)−(1−x2∣m∣(∣m∣+1)x2−∣m∣+l(l+1)−1−x2m2)u(x)=0
u(x)の係数を整理すると以下のようになる。
====1−x2∣m∣(∣m∣+1)x2−∣m∣+l(l+1)−1−x2m2 1−x2∣m∣(∣m∣+1)x2−∣m∣+l(l+1)(1−x2)−m2 1−x2−m2(1−x2)−∣m∣(1−x2)+l(l+1)(1−x2) l(l+1)−m2−∣m∣ l(l+1)−∣m∣(∣m∣+1)
従って、(3)は以下のような形で整理される。
(1−x2)dx2d2u−2(∣m∣+1)xdxdu+[l(l+1)−∣m∣(∣m∣+1)]u=0
m=0であれば、実際にルジャンドル微分方程式になる。従って、∣m∣=0の場合の解はPl0(x)=Pl(x)である。今一度、(4)をxに関して微分してみよう。係数を整理すると以下の式を得る。
(1−x2)dx3d3u−2[(∣m∣+1)+1]xdx2d2u+[l(l+1)−(∣m∣+1)(∣m∣+2)]dxdu=0
再び、(5)をxに関して微分して係数を整理すると以下の式を得る。
(1−x2)dx4d4u−2[(∣m∣+2)+1]xdx3d3u+[l(l+1)−(∣m∣+2)(∣m∣+3)]dx2d2u=0
これらの式をよく見ると、(4)をuに、∣m∣を∣m∣+1に置き換えたときに(5)を得ることができることがわかる。(5)で同様の方法で置換すると、(6)を得る。これにより、∣m∣=0の場合の解はPl(x)、∣m∣=1の場合の解はdxdPl(x)、∣m∣=2の場合の解はdx2d2Pl(x)であることがわかる。したがって、これを一般化すると次のようになる。
u(x)=dx∣m∣d∣m∣Pl(x)
したがって、関連ルジャンドル多項式は以下のようになる。
Plm(x)=(1−x2)2∣m∣dx∣m∣d∣m∣Pl(x)=(1−x2)2∣m∣dx∣m∣d∣m∣[2ll!1dxldl(x2−1)l]
■