量子力学Ⅰ/球面調和関数 のバックアップの現在との差分(No.10)
更新- バックアップ一覧
- 差分 を表示
- ソース を表示
- バックアップ を表示
- 量子力学Ⅰ/球面調和関数 へ行く。
- 追加された行はこの色です。
- 削除された行はこの色です。
[[前の単元 <<<>量子力学Ⅰ/球座標を用いた変数分離]] [[量子力学Ⅰ]] [[>>> 次の単元>量子力学Ⅰ/水素原子]]~ * 目次 [#d96730f5] [[量子力学Ⅰ]] #contents &mathjax(); * 球面調和関数 $Y^m_l(\theta,\phi)$:角運動量の固有関数 [#s564caed] * 球関数 $Y^m_l(\theta,\phi)$:角運動量の固有関数 [#s564caed] 全角運動量の二乗と、$z$ 軸周り角運動量との同時固有関数となる球面調和関数(球関数)の性質について学ぶ。 &math(\Theta); の方程式 &katex(); &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)); 中心力に対する時間を含まないシュレーディンガー方程式を変数分離した際の $Y(\theta,\phi)$ に対する方程式 は、&math(l,m); が $$\begin{aligned} \hat\Lambda Y(\theta,\phi)=-l(l+1)Y(\theta,\phi) \end{aligned}$$ &math(l=0,1,2,3,\dots); は、$Y(\theta,\phi)=\Theta(\theta)\Phi(\phi)$ と分離して、さらに &math(m=-l,-(l-1),\dots,(l-1),l); $$\begin{aligned} \frac{\partial}{\partial\phi}\Phi=im\Phi \end{aligned}$$ の範囲の整数になるときのみ解を持ち、その固有関数は''ルジャンドルの陪関数''を用いて表わすことができる。 を満たすことを仮定すると、 &math(P_l^{|m|}(\zeta)=(1-\zeta^2)^{|m|/2}\frac{d^{|m|}}{d\zeta^{|m|}}P_l(\zeta)); $$\begin{aligned}\Phi(\phi)=\Phi_0e^{im\phi}\hspace{1cm}\text{ただし}\ \ \ \ m=0,\pm 1,\pm 2,\dots \end{aligned}$$ #ref(zeta.png,right,around,33%); ただし、&math(P_l(\zeta)); は''ルジャンドルの多項式''で、 および、$\Theta$ に関する方程式 &math(P_l(\zeta)=\frac{1}{\,2^l\,l!\,}\,\frac{d^l}{\,d\zeta^l\,}(\zeta^2-1)^l); $$\Big[-\frac{1}{\sin\theta} \frac{\partial}{\partial \theta} \Big(\sin\theta\frac{\partial}{\partial \theta}\Big)+\frac{m^2}{\sin^2\theta}\Big]\Theta(\theta)=l(l+1)\Theta(\theta)$$ によって与えられる。これらを用いた を与えることを [[前項>量子力学Ⅰ/球座標を用いた変数分離#mfb957f5]] で見た。この方程式は、$l$ が $l\ge|m|$ を満たす 整数であるときのみ解を持ち、 &math( $$\begin{aligned} 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)} ); \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)} \end{aligned}$$ は規格直交完全な固有関数となり、この関数を ''球面調和関数'' と呼ぶ。((ここでは符号を &math(\Theta); に含めたが、符号を &math(\Phi); に含めても、両者で分け合っても、正規直交条件を満たすことはできる)) が、 $$\begin{aligned} \int_0^\pi \Theta_l^m(\theta)^*\Theta_{l'}^m(\theta)\,\underbrace{\sin\theta\rule[-6pt]{0pt}{0pt}}_{\text{重み}} d\theta=\delta_{ll'} \end{aligned}$$ $$\begin{aligned} \int_0^{2\pi}\Phi_m(\phi)^*\Phi_{m'}(\phi) d\phi=\delta_{mm'} \end{aligned}$$ を満たす正規直交完全な固有関数となる。この関数を ''球面調和関数'' と呼ぶ。((ここでは符号を $\Theta$ に含めたが、符号を $\Phi$ に含めても、両者で分け合っても、正規直交条件を満たすことはできる)) ここで、$P_l^{|m|}(\zeta)$ は''ルジャンドルの&ruby(ばい){陪};関数''で、 #ref(zeta.png,right,around,25%); $$P_l^{|m|}(\zeta)=\underbrace{(1-\zeta^2)^{|m|/2}\rule[-10pt]{0pt}{0pt}}_{\sin^{|m|}\theta}\,\underbrace{\frac{d^{|m|}}{d\zeta^{|m|}}P_l(\zeta)\rule[-10pt]{0pt}{0pt}}_{\cos\theta\text{\,の}\,l-|m|\,\text{次}}\hspace{1cm}\text{ただし\ }\zeta=\cos\theta$$ $P_l(\zeta)$ は''ルジャンドルの多項式''で、 $$P_l(\zeta)=\frac{1}{\,2^l\,l!\,}\,\underbrace{\frac{d^l}{\,d\zeta^l\,}\underbrace{(\zeta^2-1)^l\rule[-10pt]{0pt}{0pt}}_{\zeta\,\text{の}\,2l\,\text{次}}}_{\zeta\,\text{の}\,l\,\text{次}}$$ によって与えられる。(ここでは微分方程式の解き方には踏み込まない) #multicolumns &math(Y_0^0=\frac{1}{2 \sqrt{\pi }}); $Y_0^0=\frac{1}{2 \sqrt{\pi }}$ &math(Y_1^0=\frac{1}{2} \sqrt{\frac{3}{\pi }} \cos (\theta )); $Y_1^0=\frac{1}{2} \sqrt{\frac{3}{\pi }} \cos (\theta )$ &math(Y_1^{\pm 1}=\pm\frac{1}{2} \sqrt{\frac{3}{2 \pi }} e^{\pm i \phi } \sin (\theta )); $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)); $Y_2^0=\frac{1}{4} \sqrt{\frac{5}{\pi }} \left(3 \cos ^2(\theta )-1\right)$ &math(Y_2^{\pm 1}=\pm\frac{1}{2} \sqrt{\frac{15}{2 \pi }} e^{\pm i \phi } \sin (\theta ) \cos (\theta )); $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 )); $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)); $Y_3^0=\frac{1}{4} \sqrt{\frac{7}{\pi }} \left(5 \cos ^3(\theta )-3 \cos (\theta )\right)$ &math(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 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 )); $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}=\pm\frac{1}{8} \sqrt{\frac{35}{\pi }} e^{\pm 3 i \phi } \sin ^3(\theta )); $Y_3^{\pm 3}=\mp\frac{1}{8} \sqrt{\frac{35}{\pi }} e^{\pm 3 i \phi } \sin ^3(\theta )$ ・・・ #multicolumns(end) 正規直交完全かつ $\phi$ での微分に対して固有関数という使い勝手の良さから、普通の周期関数をフーリエ級数展開するのと同じように、任意の $\theta,\phi$ の関数を球面調和関数で展開して処理するようなことが行われるため、これは応用数学的にも重要な関数である。 ルジャンドル多項式は、$\{1,x,x^2,x^3,\dots\}$ を、内積 $(f,g)=\int_{-1}^1 f^*(x)g(x)\,dx$ の下でシュミットの直交化により正規直交化して得られる関数系であった。→ [[線形代数II/関数空間#obfa5335]] $m=0$ において、$\Theta_l^0(\theta)$ はルジャンドル多項式に $\cos\theta$ を代入した形になるが、 $x=\cos\theta$ の変数変換により $dx=-\sin\theta\,d\theta$ となり、 $$\begin{aligned} (\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^*(x)P_{l'}(x)\,dx \end{aligned}$$ となることから、$\Theta_l^0(\theta)$ の直交性はルジャンドル多項式の直交性から直接導かれる。 同様に、$\Theta_l^m(\theta)$ の直交性もルジャンドル陪関数の直交性から導かれる。 ** 特徴 [#t568ee09] - &math((m-|m|)/2=\begin{cases}0&m>=0\\m&m<0\end{cases}); より、 &math( (-1)^{(m-|m|)/2}=\begin{cases} +1\hspace{0.5cm}&m\,が偶数\\ +1\hspace{0.5cm}&m\,が奇数でm>0\\ -1&m\,が奇数でm<0 \end{cases} ); - &math(\sin\theta); と &math(\cos\theta); の &math(l); 次同次関数になっている (&math(3\cos^2\theta-1=2\cos^2\theta-\sin^2\theta); などとなることに注意せよ) $$Y_l^m(\theta,\phi) \propto \underbrace{(\sin\theta)^{|m|}\rule[-6.6pt]{0pt}{0pt}}_{0,\pi\text{に}|m|\text{重根}}\underbrace{\Big[\cos\theta\,\text{の}\,l-|m|\,\text{次多項式}\Big]}_{0\,<\,\theta\,<\,\pi\,\text{に}\,l-|m|\,\text{個の根を持つ}}e^{im\phi}$$ の形になっている。 - $\phi$ 依存性は $e^{im\phi}$ なので、$0\le\phi\le 2\pi$ で $|m|$ 回振動する - 絶対値 $|Y_l^m|$ は $\phi$ に依存しない - $\Theta_l^m$ は -- $\theta=0,\pi$ に $|m|$ 重根を持つ -- $0<\theta<\pi$ に $l-m$ 個の単根を持つ -- つまり、$0\le\theta\le\pi$ に $l+1-m$ 個の山を持つ 「振動が多い = 運動量が大きい」であることを考えると、 $\Phi$ が $m$ 回、$\Theta$ が $l+1-m$ 回振動し、合わせて $l+1$ になる 全角運動量 $l$ を $z$ 軸周り成分とその他とに分けたことに対応する。 ** 形状 [#ac710070] &math(\theta,\phi); 方向別に原点から &math(|Y_l^m(\theta,\phi)|); の距離の点を結ぶ曲面をプロットした。 色は位相を表しており、黄色が1、青が-1に対応する。 $|\Theta_l^m(\theta)|$ の形状は次のようになる。 &attachref(spherical-harmonic-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); &ref(量子力学Ⅰ/球面調和関数/メモ/theta2.svg); - &math(\Theta); は実数関数である - &math(\Phi); は位相を回転させるだけで大きさを変えない - そのため絶対値のプロットは &math(z); 軸を中心とする回転体となり、また、&math(Y_l^m); と &math(Y_l^{-m}); は同じ形になる - 位相は &math(\phi); が一周する間に &math(m); 回だけ回転する - 位相は &math(Y_l^m); と &math(Y_l^{-m}); とで &math(xz); 平面に対して対称になる - &math(\theta); 方向は隣り合う突出部で符号が反転することから、位相も反転する - &math(Y_0^0); は球形 - &math(l>0); では &math(Y_l^l); はドーナツ型。&math(l); が大きいほど扁平で、半径も大きい。 全角運動量が大きくなるため原点から遠ざかり、&math(l\sim l_z); つまり &math(l_x,l_y); 成分がほぼゼロであるために扁平になると解釈できる。 - &math(\theta); 方向には &math(l-m+1); 個の突出部が見られる ここでは縦軸に $\theta$ を取り、右方向に $|\Theta_l^m(\theta)|$ を取った。 ** より分かりやすい表示 [#a4d982f1] - グラフは左から、$l-m=0,1,2,3,4$ に相当する。$l-m$ が増えると山の数が増える = $x,y$軸周りの角運動量が大きくなる。 - 1つのプロットの中では $m=0,1,2,3,4$ のグラフを重ねて示した。 - $m$ が増えると $\theta=0,\pi$ 付近の立ち上がりが遅くなり、山が中央に寄る。 - これは $l_z$ が大きくなることで遠心力により分布が扁平になったことに対応する。 &math((Y_l^m\pm Y_l^{-m})/\sqrt{2}); をプロットすると、さらに球面調和関数の意味を理解しやすい。 上と同じものを動径プロットにしたのが次のグラフである。 (原点から $\theta$ の向きに $|\Theta_l^m(\theta)|$ だけ離れた点を結んだ曲線) - &math((\Phi_m+\Phi_{-m})/2 \propto \cos m\phi);, &math((\Phi_m-\Phi_{-m})/2i \propto \sin m\phi); - したがって、&math(\Psi_l^m=\begin{cases} \Big(Y_l^{|m|}-(-1)^{m}Y_l^{-|m|}\Big)/i\sqrt 2\hspace{5mm}&(m<0)\\ \hspace{1cm}Y_l^0&(m=0)\\ \Big(Y_l^{|m|}+(-1)^{m}Y_l^{-|m|}\Big)/\sqrt 2&(m>0) \end{cases}); と置けば、すべて実数関数となる - &math(\cos m\phi); と &math(\sin m\phi); も互いに直交するため、&math(\{\Psi_l^m\}); は正規直交系をなす &ref(量子力学Ⅰ/球面調和関数/メモ/theta1.svg); &attachref(spherical-harmonic-y2.png,,66%); さらに $\phi$ 依存性も合わせて表示するため、 $\theta,\phi$ 方向別に原点から $|Y_l^m(\theta,\phi)|$ の距離の点を結ぶ曲面をプロットしたのが次の図である。 色は複素数としての位相 $\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); - $\Theta$ は実数関数である - $\Phi$ は位相を回転させるだけで大きさを変えない - そのため絶対値のプロットは $z$ 軸を中心とする回転体となり、また、$Y_l^m$ と $Y_l^{-m}$ は同じ形になる - 位相 ($\Phi$ で決まる) は -- $\phi$ が一周する間に $m$ 回だけ回転する。 -- $Y_l^m$ と $Y_l^{-m}$ とは $xz$ 平面に対して対称になる - 形状 ($\Theta$ で決まる) は -- $\theta$ 方向は隣り合う突出部で符号が反転することから、位相も反転する -- $l-m$ が等しいもの(上図で縦に並んだもの)について形状が似ている --- $\theta$ 方向に $l-m+1$ 個の突出部が見られる。 --- これらは $l^2-l_z^2=l_x^2+l_y^2$ が等しい状態である -- $l-m$ が等しいもの同士では、$m$ が大きくなるほど扁平になり、半径も大きくなる。$l_z$ が大きくなるため遠心力により引き延ばされている。 ** $x,y,z$ との関係 [#zb464047] $Y_l^m$ は位相成分 $e^{im\phi}$ を含むため、$m\ne 0$ に対しては必然的に複素関数となる。 一方、$Y_l^m\pm Y_l^{-m}$ を作れば、 $$\begin{aligned}Y_l^m\pm Y_l^{-m}\propto e^{im\phi}\pm e^{-im\phi}\propto \sin m\phi,\cos m\phi\end{aligned}$$ であり(ただし、複合の順序は両辺で必ずしも一致しない)、係数を適当に選ぶと全体を実数関数にできる。 以下に例を挙げる: *** $p\ (l=1)$ 状態 [#p90ee679] $$\begin{aligned}\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} \end{aligned}$$ であるから、これらの軌道は $p_x,p_y,p_z$ と呼ばれる。 *** $d\ (l=2)$ 状態 [#n1c41846] $$\begin{aligned} \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} \end{aligned}$$ であるから、これらの軌道は $d_{xy},d_{yz},d_{zx},d_{x^2-y^2},d_{3z^2-r^2}$ などと呼ばれる。 同様に、 $$\begin{aligned}Y_1{}^{1}\propto (x-iy)/r\end{aligned}$$ のように $Y_l^m$ を $x,y,z$ で表すことも可能である。 ** 実関数表示 [#a4d982f1] 以下に $(Y_l^m\pm Y_l^{-m})/\sqrt{2}$ をプロットした。($m=0$ についてはそのまま $Y_l^0$ をプロットした) $\cos m\phi$ と $\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); は &math(\phi); の回転面内に &math(2m); 本の突起を持つ - &math(\Psi_l^m); と &math(\Psi_l^{-m}); は半周期だけ &math(\phi); の位相がずれている - &math(\theta); 方向には &math(l-m+1); 方向に分岐する - 隣り合う突出部では波動関数の符号が反転している(節をまたぐと符号が反転する理由は下図参照) この表示から、$\phi$ 方向に $m$ 個、$\theta$ 方向に $l-m+1$ 個の突出部を持つことが明確に読み取れる。 #ref(wave-function-polarity.png,,75%); 実際にこれらをキレイな形の実関数とするためには符号に注意する必要があり、正確な表示は次のようになる。 ** $z$ が特殊なわけではない [#oe386d3b] $$\begin{aligned} \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} \end{aligned}$$ 上のグラフを見るとあたかも &math(z); が特殊な方向であるかのように錯覚するがそんなことはない。 これらの関数は、$l_z=\hbar m$ の状態と $l_z=-\hbar m$ の状態とを等しい重みで重ね合わせたものであるので、もはや $l_z$ の固有状態ではなく、$l_z$ を測定すれば等確率で $\pm\hbar m$ のどちらかの値を取る。したがって、$\overline{l_z}=0$ である。 &math(\Psi_1{}^1=\frac{1}{\sqrt{2}}\big(Y_1^{-1}(\theta,\phi)+Y_1^{1}(\theta,\phi)\big)); 隣り合う突出部で波動関数の符号が反転している理由は下図を参照せよ。 &math(\Psi_1{}^{-1}=\frac{1}{\sqrt{2}}\big(Y_1^{-1}(\theta,\phi)-Y_1^{1}(\theta,\phi)\big)); #ref(wave-function-polarity.png,,75%); は、&math(Y_1^{0}(\theta,\phi)); とそっくり同じ形で、それぞれ &math(x,y); 方向を向いた関数となる。 ** $z$ の特殊性 [#oe386d3b] これらの関数は高校でも &math(p_x, p_y, p_z); 軌道として学んだ。 球面調和関数は $\hat l^2$ および $\hat l_z$ の同時固有関数となるように取られているため、 表式の上では $z$ が特殊な方向となっている。しかし、 &attachref(Y1-0z.jpg,,33%); $$\begin{aligned} \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} \end{aligned}$$ のようにして作った $p_x, p_y$ が $p_z$ と異なる方向を向いているものの、まったく同じ形になることからも分かるとおり、 中心力に対するシュレーディンガー方程式自体は当然等方的であり、その解も $z$ 方向を特別な方向としているわけではない。 &attachref(Y1-0z.jpg,,33%,ogp); &attachref(Y1-0x.jpg,,33%); &attachref(Y1-0y.jpg,,33%); ある量子数 &math(l); に対して &math(m); の異なる &math(2l+1); 個の固有関数が縮退している。 それらの任意の線形結合はすべて同じ固有値 &math(l); に属する固有関数となる。 &math(z); 軸を特殊な方向として取る球座標を使って変数分離したことにより、 これらの線形結合の中から、同時に &math(\hat l_z); の固有関数でもある物が &math(Y_l^m); として現れたのである。 実際、ある特定の $l$ に対する $2l+1$ 個の関数 $\{Y_l{}^{m}\}\ (m=-l,-l+1,\dots,l-1,l)$ を用いることで、$\hat l^2=\hbar^2l(l+1)$ となる「任意の関数」をこれらの線形結合(重ね合わせ)で表現できる。 &math(\hat l_z); の固有関数であるから &math(z); が特殊な方向になったというだけのことである。 言い換えると、$\hat l^2$ の固有値 $\hbar^2l(l+1)$ の固有空間は $2l+1$ 次元であり、 $\{Y_l{}^{m}\}\ (m=-l,-l+1,\dots,l-1,l)$ はその空間に取った正規直交完全系である。 同じ &math(l); に属する固有関数をすべて足し合わせてしまえば、 次のように球対称な定数関数が得られる。 この正規直交完全系のそれぞれが $\hat l_z$ の固有関数であるように選んだために $z$ が特殊な方向になったのである。 &math(\sum_{m=-1}^l|Y_l^m(\theta,\phi)|^2=\frac{1}{2\sqrt\pi}(2l+1)); このことを示すように、同じ $l$ に属する $2l+1$ 個の固有関数の確率密度をすべて足し合わせてしまえば、次のように球対称な定数関数が得られる。これは、仮想的に全ての $\{Y_l{}^{m}\}\ (m=-l,-l+1,\dots,l-1,l)$ に1つずつ電子が入るなら、その電荷密度は空間的に等方的に分布することを示している。 ** 参考資料 [#q0452a7a] $$\begin{aligned}\sum_{m=-1}^l|Y_l^m(\theta,\phi)|^2=\frac{2l+1}{4\pi}\end{aligned}$$ - http://www.sccj.net/CSSJ/jcs/v3n1/a5/textj.html * 演習:$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)); $\displaystyle\begin{aligned}\Phi_m(\phi)=\frac{1}{\sqrt{2\pi}}e^{im\phi}\end{aligned}$ より、 $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)); 一方、$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)$ に変換するには、先に導入した演算子 &math(\hat l_\pm=\hat l_x\pm i\hat l_y); $$\begin{aligned}\hat l_\pm=\hat l_x\pm i\hat l_y\end{aligned}$$ が役に立つ。 具体的には &math(m); に対する次の漸化式が成り立つ。 具体的には $m$ に対する次の漸化式が成り立つ。 &math( $$\begin{aligned} \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} ); \end{aligned}$$ すなわち、&math(\hat l_\pm); は量子数 &math(m); を1だけ増やす/減らす演算子になっている。 すなわち、$\hat l_\pm$ は量子数 $m$ を1だけ増やす/減らす演算子になっている。 (1) &math(Y_l{}^m); の &math(m); は &math(-l\le m\le l); の範囲に入らなければならなかった。 &math(m=\pm l); のとき、さらに1だけ &math(m); を増やそう/減らそうとすると何が起きるか? $Y_l{}^m$ の $m$ は $-l\le m\le l$ の範囲に入らなければならなかった。 $m=\pm l$ のとき、さらに1だけ $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_+); の固有関数になっている。 (2) $Y_l{}^m$ は $\hat l_x,\hat l_y$ の固有関数でも、$\hat l_+,\hat l_-$ の固有関数でもないが、$\hat l_+\hat l_-$ や $\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); となることを確かめよ (3) $\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); を示せ。 (4) (2),(3) の結果を用いて $\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(m=\pm l); のときにも &math(\big( \hat l_x^2+\hat l_y^2\big)Y_l{}^l=\hbar^2 l\,Y_l{}^l); となり &math(x,y); 方向の角運動量はゼロにならず、 次のように解釈できる。 (1) から分かるように、 通常、固有関数に生成・消滅演算子を適用すると、量子数を1だけ増やしたり、 減らしたりした固有関数が得られるのであるが、 量子数の取りうる範囲が限られているとき、その範囲を超えて量子数を増減させようとすると、 演算結果はゼロになってしまい、有効な固有関数を生成しない。 これは生成・消滅演算子が一般的に持つ性質である。 &math( \bm l^2Y_l{}^l=\hbar^2 l(l+1)Y_l{}^l=\big(\underbrace{\ \hbar^2l^2\ }_{l_z^2} + \underbrace{\ \hbar^2l\ }_{l_x^2+l_y^2}\big)Y_l{}^l ); 調和振動子では下降演算子と上昇演算子の合成 $\hat a^\dagger\hat a$ は数演算子 $\hat n$ として働き、 波動関数から量子数を取り出すのに使えたが、(2) から分かるように $\hat l_+\hat l_-$ は数演算子にはならない。ボゾンやフェルミオンで $\hat a^\dagger\hat a$ が数演算子 $\hat n$ として働くのはボゾンやフェルミオンの生成・消滅演算子が持つ特異な交換関係 $\hat a\hat a^\dagger\pm\hat a^\dagger\hat a=1$ に依存していることを知っておくように。 (4) より、$Y_l{}^m$ は $\hat l^2$ および $\hat l_z$ の固有関数であると共に、 $\hat l_x^2+\hat l_y^2$ の固有関数でもある。(ただし、$\hat l_x^2$ や $\hat l_y^2$ の固有関数ではない) このことは、次のように書き表してみると明らかである。 $$\begin{aligned} \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 \end{aligned}$$ すなわち、上記 (4) の結果は $\hat l^2Y_l{}^m=\hbar^2l(l+1)Y_l{}^m$ と $\hat l_zY_l{}^m=\hbar mY_l{}^m$ とから直接導けるものである。 この関係を図形的に理解しよう。 たとえば $l=5$ に対応する f 状態では $\hat l^2=\hbar^2 5(5+1)$ であるから、 対応する角運動量ベクトル $\bm l$ は半径 $|\bm l|=\hbar \sqrt{30}$ の球面上にある。 この状態に対して適当な方向に $z$ 軸を取ると、角運動量のその方向の成分 $l_z$ は 全角運動量 $|\bm l|^2$ と同時に確定値を取ることができて、その値は $m=-l,\dots,l$ に対応する $2l+1$ 個の値を取りうる。 例えば $l=5,m=5$ の状態に対しては $\hat l^2=30\hbar^2, \hat l_z=5 \hbar$ が確定し、 さらに $\hat l_x^2+\hat l_y^2=\hat l^2-\hat l_z^2=30\hbar^2-(5\hbar)^2=5\hbar^2$ が確定するものの、 $\hat l_z$ と非可換な $\hat l_x$ や $\hat l_y$ 自体は、 $\langle l_x\rangle=\langle l_y\rangle=0$ のように期待値がゼロであることはわかるものの、 ${\sigma_{l_x}}^2+{\sigma_{l_y}}^2=\hat l_x^2+\hat l_y^2=5\hbar^2$ であるから、 それぞれの値は完全に不定となる。 すなわち、角運動量ベクトル $\bm l$ は左図に示した紫の円錐と球との交線上にあることになる。 &attachref(spherical-harmonics-angular-momentum.png,,25%); 他の $m$ についてもそのような円錐を考えることができ、その方向は中図に示したようになる。 $l$ は全角運動量を決める量子数であるため、''(軌道)角運動量量子数'' あるいは ''軌道量子数''と呼ばれる。また歴史的経緯により ''方位量子数'' と呼ばれることもある。 一方、$m$ はこれが原子内での電子の運動を表す場合には $z$ 方向への磁化を決める量であるため、''(軌道)磁気量子数''と呼ばれる(荷電粒子である電子の磁気モーメントは角運動量に比例するため)。 量子数の名前に(軌道)と付いているのは、スピン(角運動量)量子数やスピン磁気量子数と区別する際に用いるためだ。 ** 方位量子数という名前について [#ne699023] 「方位量子数」は英語で Azimuthal Quantum Number といい、Azimuth という単語を方位と訳したのがこの名前だ。 通常 Azimuth とは角度 $\phi$ のことを表すので、この名前はむしろ $l_z$ に相当する量子数である $m$ にこそふさわしい感じがして、紛らわしいのだが、どうやらこの命名は歴史的経緯によるものであるらしい。 この量子数を導入したアルノルト・ゾンマーフェルトも始めは Azimuth を $z$ 軸周りの角度の意味で使っていたのであるが、後に「(原子核の周りを惑星のように一定の軌道を描いて回転する電子の)軌道面上で測った角度」の意味で使い始め、 結果的に Azimuthal Quantum Number が「軌道面に垂直な軸周りの角運動量」すなわち「全軌道角運動量」 を表す量子数を指すようになったらしい。((https://hsm.stackexchange.com/questions/7107/why-is-the-azimuthal-quantum-number-so-named を鵜呑みにした解説なのだけれど、、、信憑性は高そう?)) 非常に紛らわしいので、特に国内では「軌道量子数」「(軌道)角運動量量子数」との呼び方が一般的になりつつある。 ~ [[前の単元 <<<>量子力学Ⅰ/球座標を用いた変数分離]] [[量子力学Ⅰ]] [[>>> 次の単元>量子力学Ⅰ/水素原子]]~ * 質問・コメント [#se65bd6b] #article_kcaptcha **質問 [#s00f1787] >[[こいけ]] (&timetag(2022-10-21T06:58:15+09:00, 2022-10-21 (金) 15:58:15);)~ ~ 球面調和関数Y(θ,φ):角運動量の固有関数~ の項の最後から4行目が~ dx=sinθdθではなくdx=-sinθdθではないでしょうか?~ // - おっしゃる通りでした。修正いたしました。 -- [[武内(管理人)]] &new{2022-10-25 (火) 11:27:25}; #comment_kcaptcha **球面調和関数の実数化 [#w6411b6f] >[[質問]] (&timetag(2020-11-01T05:18:39+09:00, 2020-11-01 (日) 14:18:39);)~ ~ なぜ、確率密度の分布を考える際に、球面調和関数の2乗ではなく、実数化した球面調和関数の2乗を考えるのでしょうか?~ 言い換えると、球面調和関数の実数化は何を目的として、行っているのでしょうか?~ ~ 理解が浅く、的を得ていない質問であるかもしれませんが、よろしくお願いします~ // - 球面調和関数の虚数版と実数版との関係は、フーリエ級数の虚数版 $e^{\pm ikx}$ と実数版 $\sin kx,\cos kx$ との関係に似ていると考えるといいと思います。複素関数を念頭に置く場合には複素版がふさわしいけれど、主に実数関数を考える場合には実数系がふさわしい。また、実数系は $x,y,z,r$ で書き換えることのできるため、そういう意味でも重畳されますね。化学の時間に $p_x$ 軌道や $d_{xy}$ 軌道などというのが出てきたと思いますが、それらが実数版の関数系と対応しています。 -- [[武内(管理人)]] &new{2020-11-02 (月) 11:50:08}; #comment_kcaptcha
Counter: 92943 (from 2010/06/03),
today: 5,
yesterday: 0