量子力学Ⅰ/球面調和関数 の変更点

更新


[[前の単元 <<<>量子力学Ⅰ/球座標を用いた変数分離]]
             [[量子力学Ⅰ]]             
[[>>> 次の単元>量子力学Ⅰ/水素原子]]~

* 目次 [#d96730f5]

[[量子力学Ⅰ]]
#contents

* 球面調和関数 $Y^m_l(\theta,\phi)$:角運動量の固有関数 [#s564caed]

球面調和関数(球関数)が全角運動量および &math(z); 方向の角運動量の固有関数となることと、その性質について学ぶ。

&mathjax();

&math(Y(\theta,\phi)); に対する方程式は、
&math(Y(\theta,\phi)); に対する方程式

 &math(
\hat\Lambda Y(\theta,\phi)=-l(l+1)Y(\theta,\phi)
);

&math(Y(\theta,\phi)=\Theta(\theta)\Phi(\phi)); と分離でき、さらに
は、&math(Y(\theta,\phi)=\Theta(\theta)\Phi(\phi)); と分離して、さらに

 &math(
\frac{\partial}{\partial\phi}\Phi=im\Phi
);

を満たすことを仮定すると、

 &math(\Phi(\phi)=\Phi_0e^{im\phi}); ただし &math(m=0,\pm 1,\pm 2,\dots);

および、&math(\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));

を与えることを前項で見た。この方程式は、&math(l); が &math(l\ge|m|); を満たす
整数であるときのみ解を持ち、その固有関数は''ルジャンドルの陪関数''を用いて表わせる。

#ref(zeta.png,right,around,25%);
 &math(P_l^{|m|}(\zeta)=(1-\zeta^2)^{|m|/2}\frac{d^{|m|}}{d\zeta^{|m|}}P_l(\zeta));

ただし、&math(P_l(\zeta)); は''ルジャンドルの多項式''で、

 &math(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{\rule[-15pt]{0pt}{0pt}(-1)^{(m+|m|)/2}\sqrt{\frac{2l+1}{2}\frac{(l-|m|)!}{(l+|m|)!}}P_l^{|m|}(\cos\theta)}_{\displaystyle \Theta_l{}^m(\theta)}
\underbrace{\rule[-15pt]{0pt}{0pt}\frac{1}{\sqrt{2\pi}}e^{im\phi}}_{\displaystyle\Phi_m(\phi)}
);

は、

 &math(
\int_0^\pi \Theta_l^m(\theta)^*\Theta_{l'}^m(\theta)\,\sin\theta d\theta=\delta_{ll'}
);

 &math(
\int_0^{2\pi}\Phi_m(\phi)^*\Phi_{m'}(\phi) d\phi=\delta_{mm'}
);

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

#multicolumns
&math(Y_0^0=\frac{1}{2 \sqrt{\pi }});

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

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

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

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

&math(Y_2^{\pm 2}=\frac{1}{4} \sqrt{\frac{15}{2 \pi }} e^{\pm 2 i \phi } \sin ^2(\theta ));
#multicolumns
&math(Y_3^0=\frac{1}{4} \sqrt{\frac{7}{\pi }} \left(5 \cos ^3(\theta )-3 \cos (\theta )\right));

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

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

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

・・・
#multicolumns(end)

ルジャンドル多項式は、&math(\set{1,x,x^2,x^3,\dots}); を、内積 &math((f,g)=\int_{-1}^1 f^*(x)g(x)\,dx); の下で正規直交化して得られる関数系であった。→ [[線形代数II/関数空間#obfa5335]]

&math(m=0); において、&math(\Theta(\theta)); はルジャンドル多項式で表されるが、

 &math(
(\Theta_l^0,\Theta_{l'}^0)\propto\int_0^\pi P_l^*(\cos\theta)P_{l'}(\cos\theta)\sin\theta\,d\theta
&=\int_1^{-1} P_l^*(\cos\theta)P_{l'}(\cos\theta)\,d(\cos\theta)
);

となることから、&math(\Theta(\theta)); の直交性はルジャンドル多項式の直交性から直接導かれることに注意せよ。
** 特徴 [#t568ee09]

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

** 形状 [#ac710070]

&math(\theta,\phi); 方向別に原点から &math(|Y_l^m(\theta,\phi)|); の距離の点を結ぶ曲面をプロットした。
色は複素数としての位相 &math(\mathrm{arg}\,Y_l^m(\theta,\phi)); を表しており、黄色が+1、青が-1に対応する。

&attachref(spherical-harmonics-y1.png,,50%); &qr(http://dora.bk.tsukuba.ac.jp/~takeuchi/?plugin=attach&refer=%E9%87%8F%E5%AD%90%E5%8A%9B%E5%AD%A6%E2%85%A0%2F%E8%A7%92%E9%81%8B%E5%8B%95%E9%87%8F%E3%81%AE%E5%9B%BA%E6%9C%89%E9%96%A2%E6%95%B0&openfile=spherical-harmonic-y1.png);

- &math(\Theta); は実数関数である
- &math(\Phi); は位相を回転させるだけで大きさを変えない
- そのため絶対値のプロットは &math(z); 軸を中心とする回転体となり、また、&math(Y_l^m); と &math(Y_l^{-m}); は同じ形になる
- 位相 (&math(\Phi); で決まる) は
-- &math(\phi); が一周する間に &math(m); 回だけ回転する。回転数=振動数なので、振動数が多いことはそのまま、&math(z); 周りの角運動量が大きいことに対応する。
-- &math(Y_l^m); と &math(Y_l^{-m}); とは &math(xz); 平面に対して対称になる
- 形状 (&math(\Theta); で決まる) は
-- &math(\theta); 方向は隣り合う突出部で符号が反転することから、位相も反転する
-- &math(l-m); が等しいもの(上図で縦に並んだもの)について形状が似ている
--- 後に見るとおり、これらは互いに &math(l^2-l_z^2=l_x^2+l_y^2); が等しい状態である
--- &math(l-m); &math(\theta); 方向には &math(l-m+1); 個の突出部が見られる。突出部の数=振動数なので、振動数が高いことはそのまま、&math(x); 軸や &math(y); 軸周りの角運動量が大きいことに対応する。
-- 同じ縦列では、&math(m); が大きくなるほど扁平になり、半径も大きくなる。&math(l_z); が大きくなるため遠心力により引き延ばされている。

** $x,y,z$ との関係 [#zb464047]

&math(Y_l^m); は位相成分 &math(e^{im\phi}); を含むため、&math(m\ne 0); に対しては必然的に複素関数となる。

一方、&math(Y_l^m\pm Y_l^{-m}); を作れば、

  &math(Y_l^m\pm Y_l^{-m}\propto e^{im\phi}\pm e^{im\phi}\propto \sin m\phi,\cos m\phi);

であり(ただし、複合の順序は両辺で必ずしも一致しない)、係数を適当に選ぶと全体を実数関数にできる。

以下に例を挙げる:

*** $p\ (l=1)$ 状態 [#p90ee679]

 &math(\begin{array}{lll}
-(Y_1{}^{1}-Y_1{}^{-1}) &\propto \sin\theta\cos\phi&\propto x/r\\
i(Y_1{}^{1}+Y_1{}^{-1})&\propto \sin\theta\sin\phi&\propto y/r\\
Y_1{}^0                &\propto \cos\theta&\propto z/r\\
\end{array}
);

であるから、これらの軌道は &math(p_x,p_y,p_z); と呼ばれる。

*** $d\ (l=2)$ 状態 [#n1c41846]

 &math(
\begin{array}{lll}
-i(Y_2{}^{2}-Y_2{}^{-2})&\propto \sin^2\theta\sin 2\phi&\propto xy/r^2\\
i(Y_2{}^{1}+Y_2{}^{-1})&\propto \sin\theta\cos\theta\sin\phi&\propto yz/r^2\\
-(Y_2{}^{1}-Y_2{}^{-1}) &\propto \sin\theta\cos\theta\cos\phi&\propto zx/r^2\\
(Y_2{}^{2}+Y_2{}^{-2}) &\propto \sin^2\theta\cos 2\phi&\propto (x^2-y^2)/r^2\\
Y_2{}^{0}              &\propto 3\cos^2\theta-1&\propto (3z^2-r^2)/r^2
\end{array}
);

であるから、これらの軌道は &math(d_{xy},d_{yz},d_{zx},d_{x^2-y^2},d_{3z^2-r^2});
などと呼ばれる。

同様に、

 &math(Y_1{}^{1}\propto (x-iy)/r); 

などとも表せる。

** 実関数表示 [#a4d982f1]

以下に &math((Y_l^m\pm Y_l^{-m})/\sqrt{2}); をプロットした。(&math(m=0); についてはそのまま &math(Y_l^0); をプロットした)

&math(\cos m\phi); と &math(\sin m\phi); は互いに直交するため、これらの関数も正規直交系をなす

&attachref(spherical-harmonic-y2.png,,33%);
&qr(http://dora.bk.tsukuba.ac.jp/~takeuchi/?plugin=attach&refer=%E9%87%8F%E5%AD%90%E5%8A%9B%E5%AD%A6%E2%85%A0%2F%E8%A7%92%E9%81%8B%E5%8B%95%E9%87%8F%E3%81%AE%E5%9B%BA%E6%9C%89%E9%96%A2%E6%95%B0&openfile=spherical-harmonic-y2.png);

実際にこれらをキレイな形の実関数とするためには符号に注意する必要があり、正確な表示は次のようになる。

 &math(
\Psi_l^m=\begin{cases}
i\Big\{(-1)^{|m|}Y_l^{|m|}-Y_l^{-|m|}\Big\}/\sqrt 2\hspace{4mm}&(m<0)\hspace{5mm}\propto\ \sin m\phi\\
\hspace{1cm}Y_l^0&(m=0)\hspace{4mm}\propto \ \ \ 1\\
\Big\{(-1)^{|m|}Y_l^{|m|}+Y_l^{-|m|}\Big\}/\sqrt 2&(m>0)\hspace{4mm}\propto\ \cos m\phi\end{cases}
);

これらの関数は、&math(l_z=\hbar m); の状態と &math(l_z=-\hbar m); の状態とを等しい重みで重ね合わせたものであるので、&math(\overline l_z=0); となる。当然これらは &math(l_z); の固有状態ではなく、&math(l_z); の値は等確率で &math(\pm\hbar m); のどちらかの値を取ることになる。

隣り合う突出部では波動関数の符号が反転しているが、節をまたぐと符号が反転する理由は下図を参照せよ。

#ref(wave-function-polarity.png,,75%);

** $z$ の特殊性 [#oe386d3b]

球面調和関数は &math(\hat l^2); および &math(\hat l_z); の同時固有関数となるように取られているため、
表式の上では &math(z); が特殊な方向となっている。しかし、

 &math(
\begin{array}{lll}
p_x:\ \Psi_1{}^1\hspace{2.3mm}=(-Y_1{}^{1}+Y_1{}^{-1})/\sqrt{2} &\propto \sin\theta\cos\phi&\propto x/r\\
p_y:\ \Psi_1{}^{-1}=i(Y_1{}^{1}+Y_1{}^{-1})/\sqrt{2}&\propto \sin\theta\sin\phi&\propto y/r\\
p_z:\ \Psi_1{}^0\hspace{2.3mm}=Y_1{}^0                &\propto \cos\theta&\propto z/r\\
\end{array}
);

のようにして作った &math(p_x, p_y); が &math(p_z); 
と異なる方向を向いているものの、まったく同じ形になることからも分かるとおり、
中心力に対するシュレーディンガー方程式自体は当然等方的であり、その解も &math(z); 方向を特別な方向としているわけではない。

&attachref(Y1-0z.jpg,,33%,ogp); 
&attachref(Y1-0x.jpg,,33%); 
&attachref(Y1-0y.jpg,,33%); 

実際、ある特定の &math(l); に対する &math(2l+1); 個の関数 &math(\set{Y_l{}^{m}}\ (m=-l,-l+1,\dots,l-1,l)); を用いることで、&math(\hat l^2=\hbar^2l(l+1)); となる任意の関数をこれらの線形結合(重ね合わせ)で表現できる。

言い換えると、&math(\hat l^2); の固有値 &math(\hbar^2l(l+1)); の固有空間は &math(2l+1); 次元であり、
&math(\set{Y_l{}^{m}}\ (m=-l,-l+1,\dots,l-1,l)); はその空間に取った正規直交完全系である。

この正規直交完全系のそれぞれが &math(\hat l_z); の固有関数であるように選んだために &math(z); 
が特殊な方向になったのである。

このことを示すように、同じ &math(l); に属する &math(2l+1); 個の固有関数の確率密度をすべて足し合わせてしまえば、次のように球対称な定数関数が得られる。これは、仮想的に全ての &math(\set{Y_l{}^{m}}\ (m=-l,-l+1,\dots,l-1,l)); に1つずつ電子が入るなら、その電荷密度は空間的に等方的に分布することを示している。

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

* 演習:$m$ に関する漸化式 [#de06da6a]

 &math(\Phi_m(\phi)=\frac{1}{\sqrt{2\pi}}e^{im\phi}); より、
&math(e^{\pm i\phi}\Phi_m(\phi)=\Phi_{m\pm 1}(\phi));

一方、&math(Y_l{}^m(\theta,\phi)=\Theta_l{}^m(\theta)\Phi_m(\phi)); を 
&math(Y_l{}^{m\pm 1}(\theta,\phi)=\Theta_l{}^{m\pm 1}(\theta)\Phi_{m\pm 1}(\phi)); 
に変換するには、先に導入した演算子 

 &math(\hat l_\pm=\hat l_x\pm i\hat l_y);

が役に立つ。
具体的には &math(m); に対する次の漸化式が成り立つ。

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

すなわち、&math(\hat l_\pm); は量子数 &math(m); を1だけ増やす/減らす演算子になっている。

(1)
&math(Y_l{}^m); の &math(m); は &math(-l\le m\le l); の範囲に入らなければならなかった。
&math(m=\pm l); のとき、さらに1だけ &math(m); を増やそう/減らそうとすると何が起きるか?
上記の漸化式を元に確かめよ。

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

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

(4) (2),(3) の結果を用いて &math(\left(\hat l_x^2+\hat l_y^2\right)Y_l{}^m=\hbar^2(l^2-m^2+l)Y_l{}^m); を示せ。

[[●解答はこちら>@量子力学Ⅰ/球面調和関数/メモ]]

** 解説 [#b6480998]

(4) より、&math(Y_l{}^m); は &math(\hbar l^2); および &math(\hat l_z); の固有関数であると共に、
&math(\hat l_x^2+\hat l_y^2); の固有関数でもある。(ただし、&math(\hat l_x^2); や 
&math(\hat l_y^2); の固有関数ではない)

このことは、次のように書き表してみると明らかである。

 &math(
\underbrace{\hat l^2}_{\hbar^2l(l+1)}Y_l{}^m=
(\underbrace{\hat l_x^2+\hat l_y^2}_{\hbar^2(l^2+l-m^2)}+
\underbrace{\hat l_z^2}_{\hbar^2m^2})Y_l{}^m
);

すなわち、上記 (4) の結果は &math(\hat l^2Y_l{}^m=\hbar^2l(l+1)Y_l{}^m); と
&math(\hat l_zY_l{}^m=\hbar mY_l{}^m); とから直接導けるものである。

この関係を図形的に理解しよう。

たとえば &math(l=5); に対応する f 状態では &math(\hat l^2=\hbar^2 5(5+1)); であるから、
対応する角運動量ベクトル &math(\bm l); は半径 &math(|\bm l|=\hbar \sqrt{30}); の球面上にある。

この状態に対してある方向の角運動量成分 &math(l_z); は &math(|\bm l|^2); 
と同時に確定値を取ることができて、その値は &math(m=-l,\dots,l); に対応する 
&math(2l+1); 個の値を取りうる。

例えば &math(l=5,m=5); の状態に対しては &math(\hat l^2=\hbar^2 l(l+1), \hat l_z=5 \hbar); であり、
球面調和関数の一般系から &math(\langle l_x\rangle=\langle l_y\rangle=0); を導ける。
このとき、&math({\sigma_{l_x}}^2+{\sigma_{l_y}}^2=\hat l_x^2+\hat l_y^2=\hat l^2-\hat l_z^2=\hbar^2 5\cdot 6-\hbar^2 5^2)=5\hbar^2); であるから、&math(\hat l_x,\hat l_y); についてはその二乗の和が分かるだけであり、それぞれの値は完全に不定である。

すなわち、&math(l); は左図に示した紫の円錐上にあることになる。

&attachref(spherical-harmonics-angular-momentum.png,,25%);

他の &math(m); についてもそのような円錐を考えることができ、その方向は中図に示したようになる。 

これは角運動量ベクトルの成分を確定した方向(ここでは &math(z); 軸)と、
実際に角運動量ベクトルが向く方向との間の角度、
すなわち方位が &math(2l+1); 通りに決まると言ってもよい。
そこで &math(l); は''方位量子数''あるいは''軌道(角運動量)量子数''と呼ばれる。

一方、&math(m); はある方向への磁化を決める量であるため、(軌道)''磁気量子数''と呼ばれる(荷電粒子である電子の磁気モーメントは角運動量に比例するため)。

** 方位量子数という名前について [#ne699023]

「方位量子数」は英語で Azimuthal Quantum Number といい、Azimuth の語を方位と訳したのがこの名前だ。

通常 Azimuth とは角度 &math(\phi); のことを表すので、この名前はむしろ &math(l_z); 
にこそふさわしい感じがして、紛らわしいのだが、どうやらこの命名は歴史的経緯によるものであるらしい。

この量子数を導入したアルノルト・ゾンマーフェルトも始めは Azimuth を &math(z); 
軸周りの角度の意味で使っていたのであるが、後に「軌道面上で測った角度」の意味で使い始め、
結果的に Azimuthal Quantum Number が「軌道面に垂直な軸周りの角運動量」すなわち「全軌道角運動量」
を表す量子数を指すようになったらしい。((https://hsm.stackexchange.com/questions/7107/why-is-the-azimuthal-quantum-number-so-named を鵜呑みにした解説なのだけれど、、、信憑性は高そう?)) 

非常に紛らわしいので、特に国内では「軌道(角運動量)量子数」との呼び方が一般的になりつつある。


~
[[前の単元 <<<>量子力学Ⅰ/球座標を用いた変数分離]]
             [[量子力学Ⅰ]]             
[[>>> 次の単元>量子力学Ⅰ/水素原子]]~


* 質問・コメント [#se65bd6b]

#article_kcaptcha

Counter: 12606 (from 2010/06/03), today: 17, yesterday: 0