数学 において、関連するルジャンドル多項式は、 一般ルジャンドル微分方程式 の標準解である。
( 1 − x 2 ) d 2 d x 2 P ℓ m ( x ) − 2 x d d x P ℓ m ( x ) + [ ℓ ( ℓ + 1 ) − m 2 1 − x 2 ] P ℓ m ( x ) = 0 、 {\displaystyle \left(1-x^{2}\right){\frac {d^{2}}{dx^{2}}}P_{\ell }^{m}(x)-2x{\frac {d}{dx}}P_{\ell }^{m}(x)+\left[\ell (\ell +1)-{\frac {m^{2}}{1-x^{2}}}\right]P_{\ell }^{m}(x)=0,}
または同等に
d d x [ ( 1 − x 2 ) d d x P ℓ m ( x ) ] + [ ℓ ( ℓ + 1 ) − m 2 1 − x 2 ] P ℓ m ( x ) = 0 、 {\displaystyle {\frac {d}{dx}}\left[\left(1-x^{2}\right){\frac {d}{dx}}P_{\ell }^{m}(x)\right]+\left[\ell (\ell +1)-{\frac {m^{2}}{1-x^{2}}}\right]P_{\ell }^{m}(x)=0,}
ここで、添え字ℓ とm (整数) は、それぞれ関連ルジャンドル多項式の次数と次数と呼ばれます。この方程式は、ℓ とm が0 ≤ m ≤ ℓ を 満たす整数である場合、または自明に同等の負の値を持つ場合に限り、 [ −1, 1 ] 上で非特異な非ゼロ解を持ちます。さらにm が偶数の場合、この関数は多項式に なります。m がゼロで ℓ が整数の場合、 これらの関数は ルジャンドル多項式 と同一です。一般に、ℓ とmが整数の場合、正則解は、 m が奇数の場合には多項式 ではないにもかかわらず、「関連ルジャンドル多項式」と呼ばれることがあります。ℓ とm が任意の実数または複素数である関数の完全な一般クラスは、ルジャンドル関数 です。この場合、パラメータは通常ギリシャ文字でラベル付けされます。
ルジャンドル常微分方程式は、 物理学を はじめとする様々な技術分野で頻繁に登場します。特に、球面座標系 におけるラプラス方程式 (および関連する偏微分方程式)を解く際に現れます。関連ルジャンドル多項式は、 球面調和関数 の定義において重要な役割を果たします。
非負整数パラメータℓ およびmの定義 これらの関数は次のように表記されます。P ℓ m ( x ) {\displaystyle P_{\ell }^{m}(x)} ここで、上付き文字はP の次数を示し、べき乗ではありません。最も直接的な定義は、通常のルジャンドル多項式 ( m ≥ 0)の導関数によるものです。
P ℓ m ( x ) = ( − 1 ) m ( 1 − x 2 ) m / 2 d m d x m ( P ℓ ( x ) ) 、 {\displaystyle P_{\ell }^{m}(x)=(-1)^{m}(1-x^{2})^{m/2}{\frac {d^{m}}{dx^{m}}}\left(P_{\ell }(x)\right),}
この式の(−1) m 因子は、コンドン–ショートリー位相として知られています。一部の著者はこれを省略しています。この方程式で記述される関数が、指定されたパラメータ ℓ とm の値を持つ一般ルジャンドル微分方程式を満たすことは、P ℓ のルジャンドル方程式をm 回微分することによってわかります。[ 1 ] ( 1 − x 2 ) d 2 d x 2 P ℓ ( x ) − 2 x d d x P ℓ ( x ) + ℓ ( ℓ + 1 ) P ℓ ( x ) = 0. \displaystyle \left(1-x^{2}\right){\frac {d^{2}}{dx^{2}}}P_{\ell }(x)-2x{\frac {d}{dx}}P_{\ell }(x)+\ell (\ell +1)P_{\ell }(x)=0.}
さらに、ロドリゲスの公式 によれば、 P ℓ ( x ) = 1 2 ℓ ℓ ! d ℓ d x ℓ [ ( x 2 − 1 ) ℓ ] 、 {\displaystyle P_{\ell }(x)={\frac {1}{2^{\ell }\,\ell !}}\ {\frac {d^{\ell }}{dx^{\ell }}}\left[(x^{2}-1)^{\ell }\right],} P m ℓ は 次の形式で表すことができます。 P ℓ m ( x ) = ( − 1 ) m 2 ℓ ℓ ! ( 1 − x 2 ) m / 2 d ℓ + m d x ℓ + m ( x 2 − 1 ) ℓ 。 {\displaystyle P_{\ell }^{m}(x)={\frac {(-1)^{m}}{2^{\ell }\ell !}}(1-x^{2})^{m/2}\ {\frac {d^{\ell +m}}{dx^{\ell +m}}}(x^{2}-1)^{\ell }.}
この式により、 m の範囲を− ℓ ≤ m ≤ ℓ に拡張できます。この式に± mを代入して得られる P ℓ ± m の定義は比例関係にあります。実際、左辺と右辺の同じべき乗の係数を等しくすると、 d ℓ − m d x ℓ − m ( x 2 − 1 ) ℓ = c l m ( 1 − x 2 ) m d ℓ + m d x ℓ + m ( x 2 − 1 ) ℓ 、 {\displaystyle {\frac {d^{\ell -m}}{dx^{\ell -m}}}(x^{2}-1)^{\ell }=c_{lm}(1-x^{2})^{m}{\frac {d^{\ell +m}}{dx^{\ell +m}}}(x^{2}-1)^{\ell },} すると比例定数は次のようになる。 c l m = ( − 1 ) m ( ℓ − m ) ! ( ℓ + m ) ! 、 {\displaystyle c_{lm}=(-1)^{m}{\frac {(\ell -m)!}{(\ell +m)!}},} となることによって P ℓ − m ( x ) = ( − 1 ) m ( ℓ − m ) ! ( ℓ + m ) ! P ℓ m ( x ) 。 {\displaystyle P_{\ell }^{-m}(x)=(-1)^{m}{\frac {(\ell -m)!}{(\ell +m)!}}P_{\ell }^{m}(x).}
別表記法 文献では以下の代替表記も使用されています: [ 2 ] P ℓ m ( x ) = ( − 1 ) m P ℓ m ( x ) {\displaystyle P_{\ell m}(x)=(-1)^{m}P_{\ell }^{m}(x)}
ルジャンドル多項式 に関する記事で示された明示的な形式から始める
P l ( x ) = 2 l ∑ k = 0 l x k ( l k ) ( ( l + k − 1 ) / 2 l ) {\displaystyle P_{l}(x)=2^{l}\sum _{k=0}^{l}x^{k}{\binom {l}{k}}{\binom {(l+k-1)/2}{l}}}
標準ルールでは、m {\displaystyle m} べき乗の-階微分
P l m ( x ) = ( − 1 ) m ⋅ 2 l ⋅ ( 1 − x 2 ) m / 2 ⋅ ∑ k = m l k ! ( k − m ) ! ⋅ x k − m ⋅ ( l k ) ( l + k − 1 2 l ) ${\displaystyle P_{l}^{m}(x)=(-1)^{m}\cdot 2^{l}\cdot (1-x^{2})^{m/2}\cdot \sum _{k=m}^{l}{\frac {k!}{(km)!}}\cdot x^{km}\cdot {\binom {l}{k}}{\binom {\frac {l+k-1}{2}}{l}}}$
単純な単項式と二項係数の一般化された形式 を用いて。和は実質的に以下の項のみに及ぶ。l − k {\displaystyle lk} は偶数です。奇数の場合はl − k {\displaystyle lk} 二項係数( ( l + k − 1 ) / 2 l ) {\displaystyle {\binom {(l+k-1)/2}{l}}} ゼロです。
ドーハ [ 3 ] の結果をまとめると 、導関数をルジャンドル多項式に展開すると係数が定義されます。τ \displaystyle \tau }
d m d x m P l ( x ) = ∑ t = 0 ⌊ ( l − m ) / 2 ⌋ τ l 、 m 、 t P l − m − 2 t ( x ) 、 {\displaystyle {\frac {d^{m}}{dx^{m}}}P_{l}(x)=\sum _{t=0}^{\lfloor (lm)/2\rfloor }\tau _{l,m,t}P_{lm-2t}(x),}
どこ
τ l 、 m 、 t = ϵ l − t l − m − 2 t + 1 / 2 2 l − 2 t + 1 ( 2 m ) ! 2 m m ! ( 2 l − 2 t + 1 2 m ) m m + t ( m + t t ) 1 ( l − t m ) 、 {\displaystyle \tau _{l,m,t}=\epsilon _{l-t}{\frac {l-m-2t+1/2}{2l-2t+1}}{\frac {(2m)!}{2^{m}m!}}{\binom {2l-2t+1}{2m}}{\frac {m}{m+t}}{\binom {m+t}{t}}{\frac {1}{\binom {l-t}{m}}},}
そしてどこで
ϵ q ≡ { 1 、 q = 0 ; 2 、 q ≥ 1 {\displaystyle \epsilon _{q}\equiv {\begin{cases}1,&q=0;\\2,&q\geq 1\end{cases}}}
これはノイマン因子です。
直交性 関連ルジャンドル多項式は一般に互いに直交しません。例えば、P 1 1 {\displaystyle P_{1}^{1}} は直交しないP 2 2 {\displaystyle P_{2}^{2}} ただし、一部の部分集合は直交しています。0 ≤ m ≤ ℓと仮定すると、固定された m に対して直交条件を満たします。
∫ − 1 1 P k m P ℓ m d x = 2 ( ℓ + m ) ! ( 2 ℓ + 1 ) ( ℓ − m ) ! δ k 、 ℓ {\displaystyle \int _{-1}^{1}P_{k}^{m}P_{\ell }^{m}dx={\frac {2(\ell +m)!}{(2\ell +1)(\ell -m)!}}\ \delta _{k,\ell }}
ここで、δ k , ℓ はクロネッカーのデルタ です。
また、固定されたℓ に対して直交条件を満たす。
∫ − 1 1 P ℓ m P ℓ n 1 − x 2 d x = { 0 もし m ≠ n ( ℓ + m ) ! m ( ℓ − m ) ! もし m = n ≠ 0 ∞ もし m = n = 0 {\displaystyle \int _{-1}^{1}{\frac {P_{\ell }^{m}P_{\ell }^{n}}{1-x^{2}}}dx={\begin{cases}0&{\text{if }}m\neq n\\{\frac {(\ell +m)!}{m(\ell -m)!}}&{\text{if }}m=n\neq 0\\\infty &{\text{if }}m=n=0\end{cases}}}
負のm および/または負のℓ 微分方程式は、 m の符号の変化に対して明らかに不変である。
負のm に対する関数は、正のm に対する関数に比例することが上記で示されました。 P ℓ − m = ( − 1 ) m ( ℓ − m ) ! ( ℓ + m ) ! P ℓ m {\displaystyle P_{\ell }^{-m}=(-1)^{m}{\frac {(\ell -m)!}{(\ell +m)!}}P_{\ell }^{m}}
(これはロドリゲスの公式の定義に基づくものである。この定義により、様々な漸化式は正または負のm に対して適用可能となる。) もし | m | > ℓ それから P ℓ m = 0. {\displaystyle {\text{If}}\quad |m|>\ell \,\quad {\text{then}}\quad P_{\ell }^{m}=0.\,}
微分方程式は、 ℓから −ℓ − 1 への変化に対しても不変であり、負のℓ に対する関数は次のように定義される。
P − ℓ m = P ℓ − 1 m 、 ( ℓ = 1 、 2 、 … ) 。 {\displaystyle P_{-\ell }^{m}=P_{\ell -1}^{m},\ (\ell =1,\,2,\,\dots ).}
パリティ 定義から、関連ルジャンドル関数は偶関数か奇関数かが確認できる。
P ℓ m ( − x ) = ( − 1 ) ℓ − m P ℓ m ( x ) {\displaystyle P_{\ell }^{m}(-x)=(-1)^{\ell -m}P_{\ell }^{m}(x)}
最初のいくつかの関連するルジャンドル関数 m = 0の場合のルジャンドル関連関数m = 1の場合のルジャンドル関連関数m = 2の場合のルジャンドル関連関数負の値のm に対するものも含め、最初のいくつかの関連ルジャンドル関数は次のとおりです。
P 0 0 ( x ) = 1 {\displaystyle P_{0}^{0}(x)=1}
P 1 − 1 ( x ) = − 1 2 P 1 1 ( x ) P 1 0 ( x ) = x P 1 1 ( x ) = − ( 1 − x 2 ) 1 / 2 {\displaystyle {\begin{aligned}P_{1}^{-1}(x)&=-{\tfrac {1}{2}}P_{1}^{1}(x)\\P_{1}^{0}(x)&=x\\P_{1}^{1}(x)&=-(1-x^{2})^{1/2}\end{aligned}}}
P 2 − 2 ( x ) = 1 24 P 2 2 ( x ) P 2 − 1 ( x ) = − 1 6 P 2 1 ( x ) P 2 0 ( x ) = 1 2 ( 3 x 2 − 1 ) P 2 1 ( x ) = − 3 x ( 1 − x 2 ) 1 / 2 P 2 2 ( x ) = 3 ( 1 − x 2 ) {\displaystyle {\begin{aligned}P_{2}^{-2}(x)&={\tfrac {1}{24}}P_{2}^{2}(x)\\P_{2}^{-1}(x)&=-{\tfrac {1}{6}}P_{2}^{1}(x)\\P_{2}^{0}(x)&={\tfrac {1}{2}}(3x^{2}-1)\\P_{2}^{1}(x)&=-3x(1-x^{2})^{1/2}\\P_{2}^{2}(x)&=3(1-x^{2})\end{aligned}}}
P 3 − 3 ( x ) = − 1 720 P 3 3 ( x ) P 3 − 2 ( x ) = 1 120 P 3 2 ( x ) P 3 − 1 ( x ) = − 1 12 P 3 1 ( x ) P 3 0 ( x ) = 1 2 ( 5 x 3 − 3 x ) P 3 1 ( x ) = 3 2 ( 1 − 5 x 2 ) ( 1 − x 2 ) 1 / 2 P 3 2 ( x ) = 15 x ( 1 − x 2 ) P 3 3 ( x ) = − 15 ( 1 − x 2 ) 3 / 2 {\displaystyle {\begin{aligned}P_{3}^{-3}(x)&=-{\tfrac {1}{720}}P_{3}^{3}(x)\\P_{3}^{-2}(x)&={\tfrac {1}{120}}P_{3}^{2}(x)\\P_{3}^{-1}(x)&=-{\tfrac {1}{12}}P_{3}^{1}(x)\\P_{3}^{0}(x)&={\tfrac {1}{2}}(5x^{3}-3x)\\P_{3}^{1}(x)&={\tfrac {3}{2}}(1-5x^{2})(1-x^{2})^{1/2}\\P_{3}^{2}(x)&=15x(1-x^{2})\\P_{3}^{3}(x)&=-15(1-x^{2})^{3/2}\end{aligned}}}
P 4 − 4 ( x ) = 1 40320 P 4 4 ( x ) P 4 − 3 ( x ) = − 1 5040 P 4 3 ( x ) P 4 − 2 ( x ) = 1 360 P 4 2 ( x ) P 4 − 1 ( x ) = − 1 20 P 4 1 ( x ) P 4 0 ( x ) = 1 8 ( 35 x 4 − 30 x 2 + 3 ) P 4 1 ( x ) = − 5 2 ( 7 x 3 − 3 x ) ( 1 − x 2 ) 1 / 2 P 4 2 ( x ) = 15 2 ( 7 x 2 − 1 ) ( 1 − x 2 ) P 4 3 ( x ) = − 105 x ( 1 − x 2 ) 3 / 2 P 4 4 ( x ) = 105 ( 1 − x 2 ) 2 {\displaystyle {\begin{aligned}P_{4}^{-4}(x)&={\tfrac {1}{40320}}P_{4}^{4}(x)\\P_{4}^{-3}(x)&=-{\tfrac {1}{5040}}P_{4}^{3}(x)\\P_{4}^{-2}(x)&={\tfrac {1}{360}}P_{4}^{2}(x)\\P_{4}^{-1}(x)&=-{\tfrac {1}{20}}P_{4}^{1}(x)\\P_{4}^{0}(x)&={\tfrac {1}{8}}(35x^{4}-30x^{2}+3)\\P_{4}^{1}(x)&=-{\tfrac {5}{2}}(7x^{3}-3x)(1-x^{2})^{1/2}\\P_{4}^{2}(x)&={\tfrac {15}{2}}(7x^{2}-1)(1-x^{2})\\P_{4}^{3}(x)&=-105x(1-x^{2})^{3/2}\\P_{4}^{4}(x)&=105(1-x^{2})^{2}\end{aligned}}}
これらの関数には、いくつかの再帰特性があります。
( ℓ − m + 1 ) P ℓ + 1 m ( x ) = ( 2 ℓ + 1 ) x P ℓ m ( x ) − ( ℓ + m ) P ℓ − 1 m ( x ) {\displaystyle (\ell -m+1)P_{\ell +1}^{m}(x)=(2\ell +1)xP_{\ell }^{m}(x)-(\ell +m)P_{\ell -1}^{m}(x)}
2 m x P ℓ m ( x ) = − 1 − x 2 [ P ℓ m + 1 ( x ) + ( ℓ + m ) ( ℓ − m + 1 ) P ℓ m − 1 ( x ) ] {\displaystyle 2mxP_{\ell }^{m}(x)=-{\sqrt {1-x^{2}}}\left[P_{\ell }^{m+1}(x)+(\ell +m)(\ell -m+1)P_{\ell }^{m-1}(x)\right]}
1 1 − x 2 P ℓ m ( x ) = − 1 2 m [ P ℓ − 1 m + 1 ( x ) + ( ℓ + m − 1 ) ( ℓ + m ) P ℓ − 1 m − 1 ( x ) ] {\displaystyle {\frac {1}{\sqrt {1-x^{2}}}}P_{\ell }^{m}(x)={\frac {-1}{2m}}\left[P_{\ell -1}^{m+1}(x)+(\ell +m-1)(\ell +m)P_{\ell -1}^{m-1}(x)\right]}
1 1 − x 2 P ℓ m ( x ) = − 1 2 m [ P ℓ + 1 m + 1 ( x ) + ( ℓ − m + 1 ) ( ℓ − m + 2 ) P ℓ + 1 m − 1 ( x ) ] {\displaystyle {\frac {1}{\sqrt {1-x^{2}}}}P_{\ell }^{m}(x)={\frac {-1}{2m}}\left[P_{\ell +1}^{m+1}(x)+(\ell -m+1)(\ell -m+2)P_{\ell +1}^{m-1}(x)\right]}
1 − x 2 P ℓ m ( x ) = 1 2 ℓ + 1 [ ( ℓ − m + 1 ) ( ℓ − m + 2 ) P ℓ + 1 m − 1 ( x ) − ( ℓ + m − 1 ) ( ℓ + m ) P ℓ − 1 m − 1 ( x ) ] {\displaystyle {\sqrt {1-x^{2}}}P_{\ell }^{m}(x)={\frac {1}{2\ell +1}}\left[(\ell -m+1)(\ell -m+2)P_{\ell +1}^{m-1}(x)-(\ell +m-1)(\ell +m)P_{\ell -1}^{m-1}(x)\right]}
1 − x 2 P ℓ m ( x ) = − 1 2 ℓ + 1 [ P ℓ + 1 m + 1 ( x ) − P ℓ − 1 m + 1 ( x ) ] {\displaystyle {\sqrt {1-x^{2}}}P_{\ell }^{m}(x)={\frac {-1}{2\ell +1}}\left[P_{\ell +1}^{m+1}(x)-P_{\ell -1}^{m+1}(x)\right]}
1 − x 2 P ℓ m + 1 ( x ) = ( ℓ − m ) x P ℓ m ( x ) − ( ℓ + m ) P ℓ − 1 m ( x ) {\displaystyle {\sqrt {1-x^{2}}}P_{\ell }^{m+1}(x)=(\ell -m)xP_{\ell }^{m}(x)-(\ell +m)P_{\ell -1}^{m}(x)}
1 − x 2 P ℓ m + 1 ( x ) = ( ℓ − m + 1 ) P ℓ + 1 m ( x ) − ( ℓ + m + 1 ) x P ℓ m ( x ) {\displaystyle {\sqrt {1-x^{2}}}P_{\ell }^{m+1}(x)=(\ell -m+1)P_{\ell +1}^{m}(x)-(\ell +m+1)xP_{\ell }^{m}(x)}
1 − x 2 d d x P ℓ m ( x ) = 1 2 [ ( ℓ + m ) ( ℓ − m + 1 ) P ℓ m − 1 ( x ) − P ℓ m + 1 ( x ) ] {\displaystyle {\sqrt {1-x^{2}}}{\frac {d}{dx}}{P_{\ell }^{m}}(x)={\frac {1}{2}}\left[(\ell +m)(\ell -m+1)P_{\ell }^{m-1}(x)-P_{\ell }^{m+1}(x)\right]}
( 1 − x 2 ) d d x P ℓ m ( x ) = 1 2 ℓ + 1 [ ( ℓ + 1 ) ( ℓ + m ) P ℓ − 1 m ( x ) − ℓ ( ℓ − m + 1 ) P ℓ + 1 m ( x ) ] {\displaystyle (1-x^{2}){\frac {d}{dx}}{P_{\ell }^{m}}(x)={\frac {1}{2\ell +1}}\left[(\ell +1)(\ell +m)P_{\ell -1}^{m}(x)-\ell (\ell -m+1)P_{\ell +1}^{m}(x)\right]}
( x 2 − 1 ) d d x P ℓ m ( x ) = ℓ x P ℓ m ( x ) − ( ℓ + m ) P ℓ − 1 m ( x ) {\displaystyle (x^{2}-1){\frac {d}{dx}}{P_{\ell }^{m}}(x)={\ell }xP_{\ell }^{m}(x)-(\ell +m)P_{\ell -1}^{m}(x)}
( x 2 − 1 ) d d x P ℓ m ( x ) = − ( ℓ + 1 ) x P ℓ m ( x ) + ( ℓ − m + 1 ) P ℓ + 1 m ( x ) {\displaystyle (x^{2}-1){\frac {d}{dx}}{P_{\ell }^{m}}(x)=-(\ell +1)xP_{\ell }^{m}(x)+(\ell -m+1)P_{\ell +1}^{m}(x)}
( x 2 − 1 ) d d x P ℓ m ( x ) = 1 − x 2 P ℓ m + 1 ( x ) + m x P ℓ m ( x ) {\displaystyle (x^{2}-1){\frac {d}{dx}}{P_{\ell }^{m}}(x)={\sqrt {1-x^{2}}}P_{\ell }^{m+1}(x)+mxP_{\ell }^{m}(x)}
( x 2 − 1 ) d d x P ℓ m ( x ) = − ( ℓ + m ) ( ℓ − m + 1 ) 1 − x 2 P ℓ m − 1 ( x ) − m x P ℓ m ( x ) {\displaystyle (x^{2}-1){\frac {d}{dx}}{P_{\ell }^{m}}(x)=-(\ell +m)(\ell -m+1){\sqrt {1-x^{2}}}P_{\ell }^{m-1}(x)-mxP_{\ell }^{m}(x)}
( ℓ − m − 1 ) ( ℓ − m ) P ℓ m ( x ) = − P ℓ m + 2 ( x ) + P ℓ − 2 m + 2 ( x ) + ( ℓ + m ) ( ℓ + m − 1 ) P ℓ − 2 m ( x ) {\displaystyle (\ell -m-1)(\ell -m)P_{\ell }^{m}(x)=-P_{\ell }^{m+2}(x)+P_{\ell -2}^{m+2}(x)+(\ell +m)(\ell +m-1)P_{\ell -2}^{m}(x)}
有用な識別子(最初の再帰の初期値):
P ℓ + 1 ℓ + 1 ( x ) = − ( 2 ℓ + 1 ) 1 − x 2 P ℓ ℓ ( x ) {\displaystyle P_{\ell +1}^{\ell +1}(x)=-(2\ell +1){\sqrt {1-x^{2}}}P_{\ell }^{\ell }(x)} P ℓ ℓ ( x ) = ( − 1 ) ℓ ( 2 ℓ − 1 ) ! ! ( 1 − x 2 ) ( ℓ / 2 ) {\displaystyle P_{\ell }^{\ell }(x)=(-1)^{\ell }(2\ell -1)!!(1-x^{2})^{(\ell /2)}} P ℓ + 1 ℓ ( x ) = x ( 2 ℓ + 1 ) P ℓ ℓ ( x ) {\displaystyle P_{\ell +1}^{\ell }(x)=x(2\ell +1)P_{\ell }^{\ell }(x)}
!! は二重階乗 です。
3 つの関連するルジャンドル多項式の積 (次数が以下に示すように一致する) の積の積の積は、ルジャンドル多項式の積をルジャンドル多項式に関して線形な級数に展開する際に必要な要素です。たとえば、クーロン演算子 の行列要素が必要となるハートリー・フォック 型の原子計算を行う場合、これは必要になります。このために、ガウントの公式[ 4 ] [ 5 ]があります。 1 2 ∫ − 1 1 P l u ( x ) P m v ( x ) P n w ( x ) d x = ( − 1 ) s − m − w ( m + v ) ! ( n + w ) ! ( 2 s − 2 n ) ! s ! ( m − v ) ! ( s − l ) ! ( s − m ) ! ( s − n ) ! ( 2 s + 1 ) ! × ∑ t = p q ( − 1 ) t ( l + u + t ) ! ( m + n − u − t ) ! t ! ( l − u − t ) ! ( m − n + u + t ) ! ( n − w − t ) ! {\displaystyle {\begin{aligned}{\frac {1}{2}}\int _{-1}^{1}P_{l}^{u}(x)P_{m}^{v}(x)P_{n}^{w}(x)dx={}&{}(-1)^{s-m-w}{\frac {(m+v)!(n+w)!(2s-2n)!s!}{(m-v)!(s-l)!(s-m)!(s-n)!(2s+1)!}}\\&{}\times \ \sum _{t=p}^{q}(-1)^{t}{\frac {(l+u+t)!(m+n-u-t)!}{t!(l-u-t)!(m-n+u+t)!(n-w-t)!}}\end{aligned}}} この式は、以下の前提条件の下で使用してください。
次数は非負の整数であるl 、 m 、 n ≥ 0 {\displaystyle l,m,n\geq 0} 3つの順序はすべて非負の整数です。u 、 v 、 w ≥ 0 {\displaystyle u,v,w\geq 0} u {\displaystyle u} 3つの目の中で最大である命令をまとめるとu = v + w {\displaystyle u=v+w} 度数は従うm ≥ n {\displaystyle m\geq n} 式に現れるその他の量は次のように定義される。 2 s = l + m + n {\displaystyle 2s=l+m+n} p = 最大 ( 0 、 n − m − u ) {\displaystyle p=\max(0,\,n-m-u)} q = ミニ ( m + n − u 、 l − u 、 n − w ) {\displaystyle q=\min(m+n-u,\,l-u,\,n-w)}
積分は、
度数の合計は偶数なので、s {\displaystyle s} 整数です 三角形条件が満たされているm + n ≥ l ≥ m − n {\displaystyle m+n\geq l\geq m-n} DongとLemus(2002)[ 6 ] は、この公式の導出を任意の数の関連ルジャンドル多項式の積の積分に一般化した。
超幾何関数による一般化 これらの関数は、実際には一般的な複素パラメータと引数に対して定義される可能性があります。[ 7 ]
P λ μ ( z ) = 1 Γ ( 1 − μ ) [ 1 + z 1 − z ] μ / 2 2 F 1 ( − λ 、 λ + 1 ; 1 − μ ; 1 − z 2 ) {\displaystyle P_{\lambda }^{\mu }(z)={\frac {1}{\Gamma (1-\mu )}}\left[{\frac {1+z}{1-z}}\right]^{\mu /2}\,_{2}F_{1}(-\lambda ,\lambda +1;1-\mu ;{\frac {1-z}{2}})}
どこΓ {\displaystyle \Gamma } はガンマ関数 であり、2 F 1 {\displaystyle _{2}F_{1}} は超幾何関数 である
2 F 1 ( α 、 β ; γ ; z ) = Γ ( γ ) Γ ( α ) Γ ( β ) ∑ n = 0 ∞ Γ ( n + α ) Γ ( n + β ) Γ ( n + γ ) n ! z n 、 {\displaystyle \,_{2}F_{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},}
このようにより一般的な方法で定義されると、それらはルジャンドル関数 と呼ばれます。それらは以前と同じ微分方程式を満たします。
( 1 − z 2 ) y 」 − 2 z y ′ + ( λ [ λ + 1 ] − μ 2 1 − z 2 ) y = 0. {\displaystyle (1-z^{2})\,y''-2zy'+\left(\lambda [\lambda +1]-{\frac {\mu ^{2}}{1-z^{2}}}\right)\,y=0.\,}
これは2階微分方程式なので、2つ目の解があります。Q λ μ ( z ) {\displaystyle Q_{\lambda }^{\mu }(z)} 定義:
Q λ μ ( z ) = π Γ ( λ + μ + 1 ) 2 λ + 1 Γ ( λ + 3 / 2 ) 1 z λ + μ + 1 ( 1 − z 2 ) μ / 2 2 F 1 ( λ + μ + 1 2 、 λ + μ + 2 2 ; λ + 3 2 ; 1 z 2 ) {\displaystyle Q_{\lambda }^{\mu }(z)={\frac {{\sqrt {\pi }}\ \Gamma (\lambda +\mu +1)}{2^{\lambda +1}\Gamma (\lambda +3/2)}}{\frac {1}{z^{\lambda +\mu +1}}}(1-z^{2})^{\mu /2}\,_{2}F_{1}\left({\frac {\lambda +\mu +1}{2}},{\frac {\lambda +\mu +2}{2}};\lambda +{\frac {3}{2}};{\frac {1}{z^{2}}}\right)}
P λ μ ( z ) {\displaystyle P_{\lambda }^{\mu }(z)} そしてQ λ μ ( z ) {\displaystyle Q_{\lambda }^{\mu }(z)} どちらも、先に述べた様々な漸化式に従う。
角度による再パラメータ化 これらの関数は、引数が角度に関して再パラメータ化されている場合に最も役立ちます。x = コス θ {\displaystyle x=\cos \theta } :
P ℓ m ( コス θ ) = ( − 1 ) m ( 罪 θ ) m d m d ( コス θ ) m ( P ℓ ( コス θ ) ) {\displaystyle P_{\ell }^{m}(\cos \theta )=(-1)^{m}(\sin \theta )^{m}\ {\frac {d^{m}}{d(\cos \theta )^{m}}}\left(P_{\ell }(\cos \theta )\right)}
関係性を利用する( 1 − x 2 ) 1 / 2 = 罪 θ {\displaystyle (1-x^{2})^{1/2}=\sin \theta } 上記のリストから 、次のようにパラメータ化された最初のいくつかの多項式が得られます。
P 0 0 ( コス θ ) = 1 P 1 0 ( コス θ ) = コス θ P 1 1 ( コス θ ) = − 罪 θ P 2 0 ( コス θ ) = 1 2 ( 3 コス 2 θ − 1 ) P 2 1 ( コス θ ) = − 3 コス θ 罪 θ P 2 2 ( コス θ ) = 3 罪 2 θ P 3 0 ( コス θ ) = 1 2 ( 5 コス 3 θ − 3 コス θ ) P 3 1 ( コス θ ) = − 3 2 ( 5 コス 2 θ − 1 ) 罪 θ P 3 2 ( コス θ ) = 15 コス θ 罪 2 θ P 3 3 ( コス θ ) = − 15 罪 3 θ P 4 0 ( コス θ ) = 1 8 ( 35 コス 4 θ − 30 コス 2 θ + 3 ) P 4 1 ( コス θ ) = − 5 2 ( 7 コス 3 θ − 3 コス θ ) 罪 θ P 4 2 ( コス θ ) = 15 2 ( 7 コス 2 θ − 1 ) 罪 2 θ P 4 3 ( コス θ ) = − 105 コス θ 罪 3 θ P 4 4 ( コス θ ) = 105 罪 4 θ {\displaystyle {\begin{aligned}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}\theta -1)\\[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}\theta -3\cos \theta )\\[8pt]P_{3}^{1}(\cos \theta )&=-{\tfrac {3}{2}}(5\cos ^{2}\theta -1)\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}\theta -30\cos ^{2}\theta +3)\\[8pt]P_{4}^{1}(\cos \theta )&=-{\tfrac {5}{2}}(7\cos ^{3}\theta -3\cos \theta )\sin \theta \\[8pt]P_{4}^{2}(\cos \theta )&={\tfrac {15}{2}}(7\cos ^{2}\theta -1)\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{aligned}}}
上記の直交関係は、この定式化では次のようになる。固定されたm に対して、P ℓ m ( コス θ ) {\displaystyle P_{\ell }^{m}(\cos \theta )} は直交しており、θでパラメータ化されている。[ 0 、 π ] {\displaystyle [0,\pi ]} 重量罪 θ {\displaystyle \sin \theta } :
∫ 0 π P k m ( コス θ ) P ℓ m ( コス θ ) 罪 θ d θ = 2 ( ℓ + m ) ! ( 2 ℓ + 1 ) ( ℓ − m ) ! δ k 、 ℓ {\displaystyle \int _{0}^{\pi }P_{k}^{m}(\cos \theta )P_{\ell }^{m}(\cos \theta )\,\sin \theta \,d\theta ={\frac {2(\ell +m)!}{(2\ell +1)(\ell -m)!}}\ \delta _{k,\ell }}
また、固定されたℓ の場合:
∫ 0 π P ℓ m ( コス θ ) P ℓ n ( コス θ ) csc θ d θ = { 0 もし m ≠ n ( ℓ + m ) ! m ( ℓ − m ) ! もし m = n ≠ 0 ∞ もし m = n = 0 {\displaystyle \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(\ell -m)!}}&{\text{if }}m=n\neq 0\\\infty &{\text{if }}m=n=0\end{cases}}}
θに関して言えば、P ℓ m ( コス θ ) {\displaystyle P_{\ell }^{m}(\cos \theta )} は
d 2 y d θ 2 + ベビーベッド θ d y d θ + [ λ − m 2 罪 2 θ ] y = 0 {\displaystyle {\frac {d^{2}y}{d\theta ^{2}}}+\cot \theta {\frac {dy}{d\theta }}+\left[\lambda -{\frac {m^{2}}{\sin ^{2}\theta }}\right]\,y=0\,}
より正確には、整数mが与えられた場合 ≥ {\displaystyle \geq } 0の場合、上記の方程式は、次の場合にのみ非特異解を持ちます。λ = ℓ ( ℓ + 1 ) {\displaystyle \lambda =\ell (\ell +1)\,} ℓ は 整数でm 以上であり、これらの解は P ℓ m ( コス θ ) {\displaystyle P_{\ell }^{m}(\cos \theta )} 。
物理学における応用例:球面調和関数 物理学では、球 対称 が関係する多くの場面で、角度に関するルジャンドル多項式が現れる。球 座標 における余緯度 は、角度である。θ {\displaystyle \theta } 上記で使用。経度角、ϕ {\displaystyle \phi } は乗数として現れます。これらを合わせると、球面調和関数と呼ばれる関数群になります。これらの関数は 、リー群 SO(3)の作用下における2次元球面 の対称性を表します。
これらの関数が有用な理由は、方程式の解を求める上で中心的な役割を果たすからである。 ∇ 2 ψ + λ ψ = 0 {\displaystyle \nabla ^{2}\psi +\lambda \psi =0} 球面上において。球面座標θ(余緯度)とφ(経度)では、ラプラシアン は
∇ 2 ψ = ∂ 2 ψ ∂ θ 2 + ベビーベッド θ ∂ ψ ∂ θ + csc 2 θ ∂ 2 ψ ∂ ϕ 2 。 {\displaystyle \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}}}.}
偏微分方程式
∂ 2 ψ ∂ θ 2 + ベビーベッド θ ∂ ψ ∂ θ + csc 2 θ ∂ 2 ψ ∂ ϕ 2 + λ ψ = 0 {\displaystyle {\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}
変数分離 法によって解くと、φに依存する部分が得られる。罪 ( m ϕ ) {\displaystyle \sin(m\phi )} またはコス ( m ϕ ) {\displaystyle \cos(m\phi )} 整数m≥0の場合、およびθ依存部分の式
d 2 y d θ 2 + ベビーベッド θ d y d θ + [ λ − m 2 罪 2 θ ] y = 0 {\displaystyle {\frac {d^{2}y}{d\theta ^{2}}}+\cot \theta {\frac {dy}{d\theta }}+\left[\lambda -{\frac {m^{2}}{\sin ^{2}\theta }}\right]\,y=0\,}
その解決策はP ℓ m ( コス θ ) {\displaystyle P_{\ell }^{m}(\cos \theta )} とℓ ≥ m {\displaystyle \ell {\geq }m} そしてλ = ℓ ( ℓ + 1 ) {\displaystyle \lambda =\ell (\ell +1)} 。
したがって、方程式
∇ 2 ψ + λ ψ = 0 {\displaystyle \nabla ^{2}\psi +\lambda \psi =0}
非特異分離解を持つのは、λ = ℓ ( ℓ + 1 ) {\displaystyle \lambda =\ell (\ell +1)} 、そしてそれらの解は
P ℓ m ( コス θ ) コス ( m ϕ ) 0 ≤ m ≤ ℓ {\displaystyle P_{\ell }^{m}(\cos \theta )\ \cos(m\phi )\ \ \ \ 0\leq m\leq \ell }
そして
P ℓ m ( コス θ ) 罪 ( m ϕ ) 0 < m ≤ ℓ 。 {\displaystyle P_{\ell }^{m}(\cos \theta )\ \sin(m\phi )\ \ \ \ 0<m\leq \ell .}
ℓ のそれぞれの選択に対して、 m の様々な値と正弦関数および余弦関数の選択に対応する関数が2ℓ + 1 個存在する。これらの関数はすべて、球面上で積分すると、ℓ とm の 両方に関して直交する。
解は通常、複素指数関数 を用いて表される。
Y ℓ 、 m ( θ 、 ϕ ) = ( 2 ℓ + 1 ) ( ℓ − m ) ! 4 π ( ℓ + m ) ! P ℓ m ( コス θ ) e 私 m ϕ − ℓ ≤ m ≤ ℓ 。 {\displaystyle Y_{\ell ,m}(\theta ,\phi )={\sqrt {\frac {(2\ell +1)(\ell -m)!}{4\pi (\ell +m)!}}}\ P_{\ell }^{m}(\cos \theta )\ e^{im\phi }\qquad -\ell \leq m\leq \ell .} 機能Y ℓ 、 m ( θ 、 ϕ ) {\displaystyle Y_{\ell ,m}(\theta ,\phi )} は球面調和 関数であり、平方根の中の量は正規化係数です。正と負のm の関連ルジャンドル関数の関係を思い出すと、球面調和関数が恒等式[ 8 ] を満たすことが容易に示されます。
Y ℓ 、 m * ( θ 、 ϕ ) = ( − 1 ) m Y ℓ 、 − m ( θ 、 ϕ ) 。 {\displaystyle Y_{\ell ,m}^{*}(\theta ,\phi )=(-1)^{m}Y_{\ell ,-m}(\theta ,\phi ).}
球面調和関数は、フーリエ級数 の意味で完全な正規直交関数系を形成します。測地学、地磁気学、スペクトル解析の分野の研究者は、ここで示されているものとは異なる位相と正規化係数を使用します(球面調和関数 を参照)。
球面座標系における変数分離法を用いて3次元球対称偏微分方程式を解くと、半径方向の項を取り除いた後に残る部分は、通常、次の形式となる。
∇ 2 ψ ( θ 、 ϕ ) + λ ψ ( θ 、 ϕ ) = 0 、 {\displaystyle \nabla ^{2}\psi (\theta ,\phi )+\lambda \psi (\theta ,\phi )=0,}
したがって、解は球面調和関数となる。
注釈と参考文献 ↑ クーラント& ヒルベルト 1953 、V、§10 。↑ アブラモウィッツ、ミルトン ;ステガン、アイリーン・アン 編(1983年)[1964年6月]。「第8章」。『数式、グラフ、数表付き数学関数ハンドブック』 。応用数学シリーズ。第 55巻(第10版(1972年12月)の修正を加えた第9版再版;初版 )。ワシントンDC;ニューヨーク:米国商務省国立標準局;ドーバー出版。332 ページ。ISBN 978-0-486-61272-0 . LCCN 64-60036 . MR 0167642 . LCCN 65-12253 . ↑ Doha, EH (1991). "超球面多項式の微分展開と導関数の係数". Computers & Mathematics with Applications . 21 (2): 115– 122. doi : 10.1016/0898-1221(91)90089-M . ISSN 0898-1221 . ↑ ジョン・C・スレーター著『原子構造の量子論』 、マグロウヒル社(ニューヨーク、1960年)、第1巻、309ページ より引用。同書は、JA・ゴーントの原著論文『ロンドン王立協会の哲学的紀要』 、A228:151(1929年) ↑ Xu, Yu-Lin (1996). "Gaunt係数の高速評価". Math. Comp . 65 (216): 1601– 1612. Bibcode : 1996MaCom..65.1601X . doi : 10.1090/S0025-5718-96-00774-0 . ↑ Dong SH、Lemus R.、(2002)、「3 つの関連ルジャンドル多項式の重なり積分」、Appl. Math. Lett. 15、541-546。 ↑ Mavromatis, HA; Alassar, RS (1999). "3 つの関連ルジャンドル多項式の積分の一般化された公式". Appl. Math. Lett . 12 (3): 101– 105. doi : 10.1016/S0893-9659(98)00180-3 . ↑ この恒等式は、球面調和関数をウィグナーD行列 に関連付け、後者の時間反転特性を利用することによっても示すことができる。そして、± m の関連ルジャンドル関数間の関係は、球面調和関数の複素共役恒等式から証明できる。 ↑ Campos, LMBC; Cunha, FSRP (2012). "On hyperspherical Legendre polynomials and higher dimensional multipole expansions" (PDF) . J. Inequal. Spec. Func . 3 (3). Arfken, GB; Weber, HJ (2001), Mathematical methods for physicists , Academic Press, ISBN 978-0-12-059825-0 セクション12.5。(異なる符号規則を使用する。)Belousov, SL (1962)、「正規化された関連ルジャンドル多項式の表」 、数学表、第 18巻、パーガモン・プレス 。コンドン、EU;ショートリー、GH(1970)、『原子スペクトルの理論』 、ケンブリッジ、イングランド:ケンブリッジ大学出版局、OCLC 5388084 第3章クーラント、リチャード ;ヒルベルト、デイヴィッド (1953)、『数理物理学の方法』第1巻 、ニューヨーク:インターサイエンス・パブリッシャーズ社 。Dunster, TM (2010)、「ルジャンドル関数と関連関数」、Olver, Frank WJ 、Lozier, Daniel M.、Boisvert, Ronald F.、Clark, Charles W. (編)、NIST 数学関数ハンドブック 、ケンブリッジ大学出版局、ISBN 978-0-521-19225-5 MR 2723248 。 エドモンズ、AR(1957)、『量子力学における角運動量』 、プリンストン大学出版局、ISBN 978-0-691-07912-7 ; 第 2 章。Gaunt, JA (1929). "IV. ヘリウムの三重項" . Phil. Trans. Royal Soc. A . 228 ( 659– 669): 151– 196. doi : 10.1098/rsta.1929.0004 . ヒルデブランド、FB (1976)、『応用のための高度な微積分』 、プレンティスホール、ISBN 978-0-13-011189-0 。Koornwinder, Tom H.; Wong, Roderick SC; Koekoek, Roelof; Swarttouw, René F. (2010)、「直交多項式」、Olver, Frank WJ ; Lozier, Daniel M.; Boisvert, Ronald F.; Clark, Charles W. (編)、NIST 数学関数ハンドブック 、ケンブリッジ大学出版局、ISBN 978-0-521-19225-5 MR 2723248 。 Schach, SR (1976). "整数位数と次数を持つルジャンドル関連関数の新しい恒等式". SIAM J. Math. Anal . 7 (1): 59–69 . doi : 10.1137/0507007 .
外部リンク MathWorldにおける関連ルジャンドル多項式 MathWorldにおけるルジャンドル多項式 DLMFにおけるルジャンドル関数および関連関数