球面調和関数⓪:ルジャンドリアンの固有値問題


 微分演算子であるルジャンドリアン\Lambda(\theta,\phi) の固有値問題である

    \begin{eqnarray*} &&\Lambda(\theta,\phi)Y(\theta,\phi)=-l(l+1)Y(\theta,\phi)\\ &&\Lambda(\theta,\phi)\equiv\frac{1}{\sin \theta}\frac{\partial}{\partial \theta}\left(\sin \theta \frac{\partial}{\partial \theta}\right) +\frac{1}{\sin ^2\theta}\frac{\partial^2}{\partial \phi^2} \end{eqnarray*}

を解いて、固有関数を求めていく。ここで固有値は -l(l+1) としている。 あとにわかるように、この固有関数 Y(\theta,\phi) は球面調和関数になっている。


導入

かんたんな導入をおこなう。


なぜ -l(l+1) と置くか?

 一般的に演算子 \hat{A} に対する固有値方程式は

    \begin{eqnarray*} \hat{A}f=af \end{eqnarray*}

の形で、固有値 a, 固有関数 f を求める問題である。

\Lambda(\theta,\phi) の固有値を -l(l+1) と置くのは、\Lambda(\theta,\phi) に対する固有値問題を解く時に、ルジャンドルの微分方程式の形が現れるようにするためである。


ルジャンドルの微分方程式/ルジャンドル多項式

ルジャンドルの微分方程式について簡単に書いておく。微分方程式の形は

    \begin{eqnarray*} \frac{d}{dx}\left[ (1-x^2)\frac{dy}{dx} \right] +l(l+1)y=0 \end{eqnarray*}

である。

【微分方程式】ルジャンドル多項式の導出(ルジャンドルの微分方程式)

 この方程式を満たす y(x) の一つに、以下のルジャンドル多項式がある。

    \begin{eqnarray*} P_l(x)\equiv y(x)&=& \sum_{k=0}^{k=\left[\frac{l}{2}\right]} (-1)^k \frac{(2l-2k)!}{2^l k! (l-k)!(l-2k)!}x^{l-2k}\\ &=& \frac{(2l)!}{2^l (l!)^2} \left\{ x^l - \frac{l(l-1)}{2(2l-1)}x^{l-2}+ \frac{l(l-1)(l-2)(l-3)}{2\cdot 4 \cdot (2l-1)(2l-3)}x^{l-4}\cdots\right\} \end{eqnarray*}

ここで、[\frac{n}{2}]\frac{n}{2} を超えない最大の整数を表す。 いくつかの項を書いておくと

    \begin{eqnarray*} P_0(x)&=&1\\ P_1(x)&=&x\\ P_2(x)&=&\frac{1}{2}(3x^2-1)\\ P_3(x)&=&\frac{1}{2}(5x^3-3x)\\ P_4(x)&=&\frac{1}{8}(35x^4-30x^2+3)\\ P_5(x)&=&\frac{1}{8}(63x^5-70x^3+15x)\\ P_6(x)&=&\frac{1}{16}(231x^6-315x^4+105x^2-5)\\ \end{eqnarray*}

となる。


固有値・固有関数を求める

もう一度、解くべき式を書いておく。

    \begin{eqnarray*} &&\Lambda(\theta,\phi)Y(\theta,\phi)=-l(l+1)Y(\theta,\phi)\quad\cdots\quad (*)\\ &&\Lambda(\theta,\phi)\equiv\frac{1}{\sin \theta}\frac{\partial}{\partial \theta}\left(\sin \theta \frac{\partial}{\partial \theta}\right) +\frac{1}{\sin ^2\theta}\frac{\partial^2}{\partial \phi^2} \end{eqnarray*}


2変数関数を分離する

多変数関数の偏微分方程式を解く時の定石である、求める関数を複数の関数に分ける手法を用いる。 ここでは、\theta\phi を分離する。

    \begin{eqnarray*} Y(\theta,\phi)=f(\theta)g(\phi) \end{eqnarray*}

と置き、式(*)に代入する。

    \begin{eqnarray*} &&\frac{1}{\sin \theta}\frac{\partial}{\partial \theta}\left(\sin \theta \frac{\partial}{\partial \theta}\right) f(\theta)g(\phi)+\frac{1}{\sin ^2\theta}\frac{\partial^2}{\partial \phi^2}f(\theta)g(\phi)=-l(l+1)f(\theta)g(\phi)\\ \Leftrightarrow\quad &&\frac{\textcolor{blue}{g(\phi)}}{\sin \theta}\frac{\partial}{\partial \theta}\left(\sin \theta \frac{\partial}{\partial \theta}\right) f(\theta)+\frac{\textcolor{blue}{f(\theta)}}{\sin ^2\theta}\frac{\partial^2}{\partial \phi^2}g(\phi)=-l(l+1)\textcolor{blue}{f(\theta)g(\phi)}\\ \Leftrightarrow\quad &&\frac{1}{\textcolor{blue}{f(\theta)}\sin \theta}\frac{\partial}{\partial \theta}\left(\sin \theta \frac{\partial}{\partial \theta}\right) f(\theta)+\frac{1}{\textcolor{blue}{g(\phi)}\textcolor{red}{\sin ^2\theta}}\frac{\partial^2}{\partial \phi^2}g(\phi)=-l(l+1)\\ \Leftrightarrow\quad &&\left[\frac{\sin \theta}{\textcolor{blue}{f(\theta)}}\frac{\partial}{\partial \theta}\left(\sin \theta \frac{\partial}{\partial \theta}\right)+l(l+1)\sin ^2\theta\right] +\frac{1}{\textcolor{blue}{g(\phi)}}\frac{\partial^2}{\partial \phi^2}g(\phi)=0 \end{eqnarray*}

ここでの計算の注意は以下のとおり。

  • 2行目:微分に関係ない関数は左へ通り抜ける
  • 3行目:両辺を f(\theta)g(\phi) で割った
  • 4行目:両辺に \sin ^2\theta をかけた

φについての関数を求める

最後の式の左辺は、

  • 第一項[\cdots]が \theta のみの関数
  • 第二項が \phi のみの関数

となっている。その和が0(右辺)なので、第一項と第二項はそれぞれ定数 \pm m^2 である。 あとのため、定数を m^2 と置いた。

つまり、

    \begin{eqnarray*} &&\frac{\sin \theta}{f(\theta)}\frac{\partial}{\partial \theta}\left(\sin \theta \frac{\partial}{\partial \theta}\right)+l(l+1)\sin ^2\theta =m^2\\ &&\frac{1}{g(\phi)}\frac{\partial^2}{\partial \phi^2}g(\phi)=-m^2 \end{eqnarray*}

となる。2つ目の式から

    \begin{eqnarray*} \frac{\partial^2 g(\phi)}{\partial \phi^2}&=&-m^2 \phi\\ g(\phi)&=& g_0 \exp(im\phi)\quad(m=0,\pm1,\pm2,\cdots)\quad\blacksquare \end{eqnarray*}

ここで物理的な観点から、周期性 g(\phi+2\pi)=g(\phi) を満たすため、m は整数になる。 規格化定数 g_0=\frac{1}{\sqrt{2\pi}} は後で求める。


θについての関数を求める

つぎに、\theta について、

    \begin{eqnarray*} &&\frac{\sin \theta}{f(\theta)}\frac{\partial}{\partial \theta}\left(\sin \theta \frac{\partial}{\partial \theta}\right)+l(l+1)\sin ^2\theta =m^2\\ \Leftrightarrow\quad &&\sin \theta\frac{\partial}{\partial \theta}\left(\sin \theta \frac{\partial}{\partial \theta}\right)+\left[l(l+1)\sin ^2\theta-m^2\right]f(\theta) =0 \end{eqnarray*}

ここで、z=\cos\theta;\,\sin \theta=\sqrt{1-z^2} と置く。これは、極座標変換と同じである(図)。 また、f(\theta)\to f(z) と置き換えておく。

チェーンルールより

    \begin{eqnarray*} \frac{d}{d\theta}&=&\frac{dz}{d\theta}\cdot\frac{d}{dt}\\ &=&-\sin \theta\frac{d}{dt}\\ &=&-\sqrt{1-z^2}\frac{d}{dt} \end{eqnarray*}

を利用する(また、1変数なので微分は \partial\to d としておく)。

    \begin{eqnarray*} &&\sqrt{1-z^2}\left( -\sqrt{1-z^2}\frac{d}{dz} \right)\left(-(1-z^2)\frac{d}{dz}\right)f(z) +\left[ l(l+1)(1-z^2) -m^2\right]f(z)=0\\ \Leftrightarrow\quad&& -(1-z^2)\frac{d}{dz}\left( -(1-z^2)\frac{df(z)}{dz} \right)+ \left[l(l+1)(1-z^2)-m^2\right]f(z)=0\\ \Leftrightarrow\quad&& \frac{d}{dz}\left[(1-z^2)\frac{df(z)}{dz}\right] -\left[l(l+1)-\frac{m^2}{1-z^2}\right]f(z)=0 \end{eqnarray*}


m=0のとき

m=0,\pm1,\pm2,\cdots であったが、m=0 とすると、

    \begin{eqnarray*} \frac{d}{dz}\left[(1-z^2)\frac{df(z)}{dz}\right] -l(l+1)f(z)=0 \end{eqnarray*}

ルジャンドルの微分方程式そのものである。したがって、はじめに見たように

    \begin{eqnarray*} f(z)=P_l(z)= \sum_{k=0}^{k=\left[\frac{l}{2}\right]} (-1)^k \frac{(2l-2k)!}{2^l k! (l-k)!(l-2k)!}x^{l-2k} \end{eqnarray*}

となる。最初に -l(l+1) と固有値を決めておいたことで、ルジャンドルの微分方程式が現れたのである。


m≠0のとき

m\neq 0 のときはルジャンドルの微分方程式ではないが、その解はルジャンドルの陪多項式 P_l^{m}(z) になる。 このとき、微分方程式は l,m に依存していたため、多項式解が l,m に依存している。 また、m については m^2 の形で微分方程式に現れているため、P_l^{|m|}(z) と表す。

ルジャンドルの倍多項式 P_l^{|m|} はルジャンドルの多項式 P_l(z) によって、

    \begin{eqnarray*} P_l^{|m|}(z)=(1-z^2)^{\frac{|m|}{2}}\frac{d^{|m}}{dz^{|m|}}P_l(z)\quad(m=-l,-l+1,\cdots,l-1,l) \end{eqnarray*}

となることが知られている。m|m|\leq l を満たす 2l+1 個の整数であり、l に依存するため m_l と書く。 この m_l は磁気量子数と呼ばれる。

以上より、z=\cos\theta であったので、f(\theta)=f_0 P_l^{|m_l|}(\cos\theta) となる。


規格化定数

以上の結果をまとめて、

    \begin{eqnarray*} f(\theta)&=&f_0 P_l^{|m_l|}(\cos\theta)\\ g(\phi)&=&g_0 \exp(im_l\phi) \end{eqnarray*}

となる。次にそれぞれの規格化定数を求める。

g(\phi) については簡単で、

    \begin{eqnarray*} \int_0^{2\pi}|g(\phi)|^2\,d\phi&=&1\\ |g_0|^2 \int_0^{2\pi} \exp(+im_l\phi)\exp(-im_l\phi)\,d\phi&=&1\\ g_0&=&\frac{1}{\sqrt{2\pi}}\quad\blacksquare \end{eqnarray*}

f(\theta) については、極座標における \thetaの積分に \sin \theta が入ることに注意して、

    \begin{eqnarray*} &&\int_0^{\pi} \left|f_0 P_l^{|m_l|}(\cos\theta)\right|^2\sin \theta d\theta=1\\ &&\int_{-1}^{1} \left|f_0 P_l^{|m_l|}(z)\right|^2\,dz=1\quad(\because z=\cos\theta)\\ &&|f_0|^2 \int_{-1}^{1} \left[P_l^{|m_l|}(z)\right]^2\,dz=1\\ \end{eqnarray*}

となる。このルジャンドル陪多項式の規格化については煩雑なので結果のみを示す。

    \begin{eqnarray*} f_0=(-1)^{\frac{(m_l+|m_l|)}{2}}\sqrt{\frac{(2l+1)(l-|m_l|)!}{2(l+|m_l|)!}} \end{eqnarray*}


球面調和関数を得る

規格化定数を求めることができたので球面調和関数 Y_{l,m_l}(\theta,\phi) を得る。

    \begin{eqnarray*} Y_{l,m_l}(\theta,\phi) &=&f(\theta)g(\theta)\\\\ &=&(-1)^{\frac{(m_l+|m_l|)}{2}}\sqrt{\frac{(2l+1)(l-|m_l|)!}{2(l+|m_l|)!}} P_l^{|m_l|}(\cos\theta)\,\frac{1}{\sqrt{2\pi}}\exp(im_l\phi)\quad\blacksquare \end{eqnarray*}

具体的な値は以下のとおり。

l=0;\, m=0:

    \begin{eqnarray*} Y_{0,0} =\frac{1}{\sqrt{4\pi}} \end{eqnarray*}

l=1;\,m=0,\pm 1:

    \begin{eqnarray*} Y_{1,0} &=&\sqrt{\frac{3}{4\pi}}\cos\theta\\ Y_{1,\pm1} &=&\sqrt{\frac{3}{8\pi}}\sin \theta e^{\pm i\phi} \end{eqnarray*}

l=2;\,m=0,\pm 1,\pm 2:

    \begin{eqnarray*} Y_{2,0} &=&{\sqrt\frac{5}{16\pi}}(3\cos^2\theta-1)\\ Y_{2,\pm1} &=&\sqrt{\frac{15}{8\pi}}\sin \theta\cos\theta e^{\pm i\phi}\\ Y_{2,\pm2} &=&\sqrt{\frac{15}{32\pi}}\sin ^2\theta e^{\pm 2i\phi} \end{eqnarray*}

l=3;\,m=0,\pm 1,\pm 2,\pm 3:

    \begin{eqnarray*} Y_{3,0} &=&{\sqrt\frac{7}{16\pi}}(5\cos^3\theta-3\cos\theta)\\ Y_{3,\pm1} &=&\sqrt{\frac{21}{64t\pi}}(5\cos^2\theta-1)\sin \theta e^{\pm i\phi}\\ Y_{3,\pm2} &=&\sqrt{\frac{105}{32\pi}}\sin ^2\theta\cos\theta e^{\pm 2i\phi}\\ Y_{3,\pm3} &=&\sqrt{\frac{35}{64\pi}}\sin ^3\theta e^{\pm 3i\phi} \end{eqnarray*}


まとめ

m_l=0 のときはルジャンドルの微分方程式(mに依存しない)が現れ、その解がルジャンドル多項式 P_l(z) になる。 m_l\neq0 のときは mに依存する微分方程式が現れ、その解がルジャンドル陪多項式 P_l^{|m_l|}(z) になる。



コメントを残す

メールアドレスが公開されることはありません。