量子力学Ⅰ/球面調和関数 のバックアップ(No.7)

更新


目次

量子力学Ⅰ

球関数 $Y^m_l(\theta,\phi)$:角運動量の固有関数

\Theta の方程式

 &math(\Big[\sin\theta \frac{\PD}{\PD \theta} \Big(\sin\theta\frac{\PD}{\PD \theta}\Big)+l(l+1) \sin^2\theta\Big]\Theta(\theta)=m^2\Theta(\theta));

は、 l,m

  l=0,1,2,3,\dots

  m=-l,-(l-1),\dots,(l-1),l

zeta.png

の範囲の整数になるときのみ解を持ち、その固有関数はルジャンドルの陪関数を用いて表わすことができる。

  P_l^{|m|}(\zeta)=(1-\zeta^2)^{|m|/2}\frac{d^{|m|}}{d\zeta^{|m|}}P_l(\zeta)

ただし、 P_l(\zeta) ルジャンドルの多項式で、

  P_l(\zeta)=\frac{1}{\,2^l\,l!\,}\,\frac{d^l}{\,d\zeta^l\,}(\zeta^2-1)^l

によって与えられる。これらを用いた

 &math( Y_l^m(\theta,\phi)= \underbrace{(-1)^{(m+|m|)/2}\sqrt{\frac{2l+1}{2}\frac{(l-|m|)!}{(l+|m|)!}}P_l^{|m|}(\cos\theta)}_{\Theta(\theta)} \underbrace{\frac{1}{\sqrt{2\pi}}e^{im\phi}}_{\Phi(\phi)} );

は規格直交完全な固有関数となり、この関数を 球面調和関数 と呼ぶ。*1ここでは符号を \Theta に含めたが、符号を \Phi に含めても、両者で分け合っても、正規直交条件を満たすことはできる

Y_0^0=\frac{1}{2 \sqrt{\pi }}

Y_1^0=\frac{1}{2} \sqrt{\frac{3}{\pi }} \cos (\theta )

Y_1^{\pm 1}=\pm\frac{1}{2} \sqrt{\frac{3}{2 \pi }} e^{\pm i \phi } \sin (\theta )

Y_2^0=\frac{1}{4} \sqrt{\frac{5}{\pi }} \left(3 \cos ^2(\theta )-1\right)

Y_2^{\pm 1}=\pm\frac{1}{2} \sqrt{\frac{15}{2 \pi }} e^{\pm i \phi } \sin (\theta ) \cos (\theta )

Y_2^{\pm 2}=\frac{1}{4} \sqrt{\frac{15}{2 \pi }} e^{\pm 2 i \phi } \sin ^2(\theta )

Y_3^0=\frac{1}{4} \sqrt{\frac{7}{\pi }} \left(5 \cos ^3(\theta )-3 \cos (\theta )\right)

Y_3^{\pm 1}=\pm\frac{1}{8} \sqrt{\frac{21}{\pi }} e^{\pm i \phi } \sin (\theta ) \left(5 \cos ^2(\theta )-1\right)

Y_3^{\pm 2}=\frac{1}{4} \sqrt{\frac{105}{2 \pi }} e^{\pm 2 i \phi } \sin ^2(\theta ) \cos (\theta )

Y_3^{\pm 3}=\pm\frac{1}{8} \sqrt{\frac{35}{\pi }} e^{\pm 3 i \phi } \sin ^3(\theta )

・・・

特徴

  • (m-|m|)/2=\begin{cases}0&m>=0\\m&m<0\end{cases} より、 &math( (-1)^{(m-|m|)/2}=\begin{cases}
  1. 1\hspace{0.5cm}&m\,が偶数\\
  2. 1\hspace{0.5cm}&m\,が奇数でm>0\\
  • 1&m\,が奇数でm<0 \end{cases} );
  • \sin\theta \cos\theta l 次同次関数になっている ( 3\cos^2\theta-1=2\cos^2\theta-\sin^2\theta などとなることに注意せよ)

形状

\theta,\phi 方向別に原点から |Y_l^m(\theta,\phi)| の距離の点を結ぶ曲面をプロットした。 色が位相を表しており、黄色が1、青が-1に対応する。

File not found: "spherical-harmonic-y1.png" at page "量子力学Ⅰ/球面調和関数"[添付] 

  • \Phi は位相を回転させるだけで大きさを変えない
  • そのため、 Y_l^m Y_l^{-m} は位相のみが異なり、絶対値のプロットでは同じ形になる
  • 位相は \phi が一周する間に m 回だけ回転する
  • 位相は Y_l^m Y_l^{-m} とで xz 平面に対して対称になる
  • \Theta は隣り合う突出部で符号が反転することから、位相も反転する
  • Y_0^0 は球形
  • l>0 では Y_l^l はドーナツ型。 l が大きいほど扁平で、半径も大きい。 全角運動量が大きくなるため原点から遠ざかり、 l\sim l_z つまり l_x,l_y 成分がほぼゼロであるために扁平になると解釈できる。
  • \Theta に含まれる \cos\theta の次数を反映して、 z 軸方向に l-m+1 個の突出部が見られる。

より分かりやすい表示

|Y_l^m\pm Y_l^{-m}| をプロットすると、さらに球面調和関数の意味を理解しやすい。

spherical-harmonic-y2.png

  • (\Phi_m+\Phi_{-m})/2 \propto \cos m\phi , (\Phi_m-\Phi_{-m})/2i \propto \sin m\phi
  • そのため |Y_l^m\pm Y_l^{-m}| \phi の回転面内に 2m 本の突起を持つ
  • |Y_l^m+ Y_l^{-m}| |Y_l^m- Y_l^{-m}| は平面 x=y に対して鏡面対称
  • \theta 方向は実数関数なので位相は変化せず、大きさのみが変化する
  • \theta 方向には l-m+1 方向に分岐する

$z$ が特殊なわけではない

上のグラフを見るとあたかも z が特殊な方向であるかのように錯覚するがそんなことはない。

  \frac{1}{\sqrt{2}}\big(Y_1^{-1}(\theta,\phi)+Y_1^{1}(\theta,\phi)\big)

や、

  \frac{1}{\sqrt{2}}\big(Y_1^{-1}(\theta,\phi)-Y_1^{1}(\theta,\phi)\big)

は、 \big(Y_1^{0}(\theta,\phi) とそっくり同じ形で、それぞれ x,y 方向を向いた関数となる。

Y1-0z.jpg  Y1-0x.jpg  Y1-0y.jpg 

これらの関数は高校でも 2p_x, 2p_y, 2p_z などとして学んだ。

ある量子数 l に対して m の異なる 2l+1 個の固有関数が縮退している。 それらの任意の線形結合はすべて同じ固有値 l に属する固有関数となる。 z 軸を特殊な方向として取る球座標を使って変数分離したことにより、 これらの線形結合の中から、同時に \hat l_z の固有関数でもある物が Y_l^m として現れたのである。

\hat l_z の固有関数であるから z が特殊な方向になったというだけのことである。

同じ l に属する固有関数をすべて足し合わせてしまえば、

  \frac{1}{2l+1}\sum_{m=-1}^l|Y_l^m(\theta,\phi)|^2=\frac{1}{2\sqrt\pi}

のように球対称な定数関数が得られる。

参考資料

演習:$m$ に関する漸化式

  \Phi_m(\phi)=\frac{1}{\sqrt{2\pi}}e^{im\phi}

であるから、

  e^{\pm i\phi}\Phi_m(\phi)=\Phi_{m\pm 1}(\phi)

である。

しかし、 \Theta_l{}^m(\theta) をそのままにして \Phi_m(\phi) だけ m を増減させてしまうとつじつまが合わなくなってしまう。

Y_l{}^m(\theta,\phi)=\Theta_l{}^m(\theta)\Phi_m(\phi) Y_l{}^{m\pm 1}(\theta,\phi)=\Theta_l{}^{m\pm 1}(\theta)\Phi_{m\pm 1}(\phi) に変換するには、先に導入した演算子 \hat l_{\pm} が役に立つ。

具体的には次の等式が成り立つ。

 &math( \begin{cases} \hat l_+Y_l{}^{m}(\theta,\phi)=\hbar\sqrt{(l-m)(l+m+1)}Y_l{}^{m+1}\\ \hat l_-Y_l{}^{m}(\theta,\phi)=\hbar\sqrt{(l+m)(l-m+1)}Y_l{}^{m-1}\\ \end{cases} );

すなわち、 \hat l_\pm=\hat l_x\pm i\hat l_y は量子数 m を1だけ増やす/減らす演算子になっている。

(1) シュレーディンガー方程式を解くに当たり、 m -l\le m\le l の範囲に入らなければならないという制約があったが、 m=\pm l のとき、さらに1だけ増やそう/減らそうとすると何が起きるだろうか。

  Y_l{}^{\pm l}(\theta,\phi)\propto \sin^l\theta e^{\pm il\phi}

であることに注意して、 \hat l_\pm Y_l{}^{\pm l}=0 を導け。
ただし、  \hat l_\pm=\hat l_x\pm i\hat l_y=\hbar e^{\pm i\phi}\Big(\pm\frac{\PD}{\PD\theta}+\frac{i}{\tan\theta}\frac{\PD}{\PD\phi}\Big) である。

(2) Y_l{}^m \hat l_x,\hat l_y の固有関数でも、 \hat l_+,\hat l_- の固有関数でもないが、 \hat l_+\hat l_- \hat l_-\hat l_+ の固有関数になっていることを示し、 その固有値を求めよ。

(3) \hat l_x^2+\hat l_y^2&=\frac{1}{2}\left\{ \hat l_+\hat l_-+\hat l_-\hat l_+\right\} となることを確かめよ

(4) (2),(3) の結果を用いて \langle \hat l_x^2+\hat l_y^2\rangle=\hbar^2(l^2-m^2+l) を示せ。

●解答はこちら

解説

(4) より、 m=\pm l のときにも \langle \hat l_x^2+\hat l_y^2\rangle=\hbar^2 l となり x,y 方向の角運動量はゼロにならない。すなわち、

 &math( \bm l^2=\hbar^2 l(l+1)=\underbrace{\ \hbar^2l^2\ }_{l_z^2} + \underbrace{\ \hbar^2l\ }_{\langle l_x^2+l_y^2\rangle} );

と理解できるのである。


*1 ここでは符号を \Theta に含めたが、符号を \Phi に含めても、両者で分け合っても、正規直交条件を満たすことはできる

Counter: 89902 (from 2010/06/03), today: 12, yesterday: 0