量子力学Ⅰ/3次元調和振動子 のバックアップ(No.14)

更新


前の単元 <<<              量子力学Ⅰ              >>> 次の単元

復習:$x,y,z$ 座標で解いた3次元調和振動子

バネ定数 K=m\omega^2 \omega は系の固有角振動数) の3次元調和振動子に対する時間を含まないシュレーディンガー方程式を x,y,z について変数分離して解けば、その解は \varphi(x,y,z)=X(x)Y(y)Z(z) の形となり、 X(x),Y(y),Z(z) はそれぞれ1次元調和振動子に対する解となることを 量子力学Ⅰ/調和振動子#y596d643 で学んだ。

すなわち r_0=\sqrt{\hbar/m\omega} およびエルミート多項式 H_n を使って、

  \varphi_{n_xn_yn_z}(x,y,z)=H_{n_x}(x/r_0)H_{n_y}(y/r_0)H_{n_z}(z/r_0)e^{-\frac{r^2}{2r_0^2}}

対応するエネルギーは N=n_x+n_y+n_z を用いて、

 &math( \varepsilon_{n_xn_yn_z}&=\hbar\omega(n_x+n_y+n_z+3/2)\\ &=\hbar\omega(N+3/2)\\ );

N=n_x+n_y+n_z       (n_xn_yn_z)      縮退度
0(0 0 0)1
1(1 0 0), (0 1 0), (0 0 1)3
2(2 0 0), (0 2 0), (0 0 2)6
(1 1 0), (1 0 1), (0 1 1)
3(3 0 0), (0 3 0), (0 0 3)10
(2 1 0), (0 2 1), (1 0 2)
(1 2 0), (0 1 2), (2 0 1)
(1 1 1)

球座標で解いた3次元調和振動子

一方、この問題は球対称なポテンシャル

  V(r,\theta,\phi)=\frac{1}{2}Kr^2

の中での運動であるから、球座標で展開して解けば

  \phi_{nlm}'(r,\theta,\phi)=R_n{}^l(r)Y_l{}^m(\theta,\phi)

の形の解が得られる。このとき、

  \varepsilon=\varepsilon_n{}^l

と表せ、エネルギーは m の値によらない。

状態lm縮退度
s001
p1-1, 0, +13
d2-2, -1, 0, 1, 25
f3-3, -2, -1, 0, 1, 2, 37

ただし、水素原子のところでも見たとおり、ポテンシャルの形状によっては m に対する縮退以外にも、異なる l を持つ状態が縮退することがある。

両者の関係は?

次に見るとおり、

  • N=0 の縮退度は 1 → l=0 (縮退度 1 の s 状態)
  • N=1 の縮退度は 3 → l=1 (縮退度 3 の p 状態)
  • N=2 の縮退度は 6 → l=2 (縮退度 5 の d 状態) と l=0 (縮退度 1 の s 状態) が縮退
  • N=3 の縮退度は 10 → l=3 (縮退度 7 の f 状態) と l=1 (縮退度 3 の p 状態) が縮退

にそれぞれ対応している。

演習:球座標による解の縮退

x,y,z 座標で変数分離して得られる解と、 r,\theta,\phi 座標で変数分離して得られる解との関係を明らかにしたい。

(1) 3次元調和振動子に対する動径方向の方程式は、 r/r_0=\xi と書き換え、 rR(r)=X(\xi) と置けば、

&math( X''(\xi)=\Big[\xi^2-\frac{2\varepsilon}{\hbar\omega}+\frac{l(l+1)}{\xi^2}\Big]X(\xi) );

となることを示せ。( l=0 のとき、この式は1次元調和振動子の式と一致する)

(2) \frac{d^2}{d\xi^2}\left(\xi^k e^{-\xi^2/2}\right)=\left(\xi^2-(2k+1)+\frac{(k-1)k}{\xi^2}\right)\xi^k e^{-\xi^2/2} を確かめよ。

(3) X(\xi)=\xi^{l+1} e^{-\frac{\xi^2}{2}} \varepsilon=\left(l+\frac{3}{2}\right)\hbar\omega に対する固有関数となることを示せ。

(4) X(\xi)=\left(\xi^{l+1}-\frac{2}{2l+3}\xi^{l+3}\right)e^{-\xi^2/2} \varepsilon=\left(l+\frac{7}{2}\right)\hbar\omega に対する固有関数となることを示せ。

(5) (発展) X_n{}^l(\xi)=\left(\sum_{k=0}^n c_k\xi^{2k}\right)\xi^{l+1}e^{-\xi^2/2} \varepsilon=\left(l+2n+\frac{3}{2}\right)\hbar\omega に対する固有関数となるよう c_k を定めよ。


● 解答はこちら

解説

(3) で l=0 とすれば R(r)\propto e^{-\xi^2/2} となって、 N=0 の解に相当する。

(3) で l=1 とすれば R(r)\propto \xi e^{-\xi^2/2} となって、 N=1 の解に相当する。

(3) で l=2 とした R(r)\propto \xi^2 e^{-\xi^2/2} と、 (4) で l=0 とした R(r)\propto \left(\xi-\frac{2}{3}\xi^3\right)e^{-\xi^2/2} とは同じエネルギーを持ち(縮退しており)、 N=2 の解に相当する。

(3) で l=3 とした R(r)\propto \xi^3 e^{-\xi^2/2} と、 (4) で l=1 とした R(r)\propto \left(\xi^2-\frac{2}{5}\xi^4\right)e^{-\xi^2/2} とは同じエネルギーを持ち(縮退しており)、 N=3 の解に相当する。

(5) のように、一般には

  X_n{}^l(\xi)=(\xi の 2n 次多項式)\ \xi^{l+1} e^{-\xi^2/2}

の形で \varepsilon=(l+2n+3/2)\hbar\omega に対する固有関数を作れるから (n=0,1,2,\dots) N=l+2n となるように l,n を選べば N と同じエネルギーを与える。

たとえば N=4 なら以下の通り15重に縮退する。

  • n=0,\ l=4 の 9 個
  • n=1,\ l=2 の 5 個
  • n=2,\ l=0 の 1 個

グラフ形状

spherical_harmonic_oscillation.svg 

球対称な箱形ポテンシャルと同様に、

  • n が増加すると振動回数が増える → r 方向の運動量が大きくなる
  • l が増加すると角運動量が増える → 遠心力で外側へ寄る

エネルギーは (l+2n+3/2)\hbar\omega なので、 n の増加は l の増加の2倍、エネルギーを増やす。

n の増加は調和振動子の振幅の増加に当たり、 ポテンシャルエネルギーが大きく増加するためだ。

$N=0$ の解

 &math( \varphi_{000}(\bm r) &=\left(\frac{4\pi}{r_0}\right)^{3/4}e^{-\frac{r^2}{2r_0^2}} &\propto e^{-\frac{r^2}{2r_0^2}} );

\theta,\phi に依存しない。

一方、 Y_0{}^0=1/\sqrt{4\pi} \theta,\phi に依存しないので、

 &math( \varphi_{000}(\bm r) &\propto e^{-\frac{r^2}{r_0^2}}\,Y_0{}^0\\ );

と書ける。

すなわち n=0 の解は s 状態であり、 \varphi_{000} に対して \hat l^2=0,\hbar l_z=0 であることが分かる。

$N=1$ の解

 &math( \varphi_{100}(\bm r)\propto xe^{-\frac{r^2}{2r_0^2}} =r\sin\theta\cos\phi e^{-\frac{r^2}{2r_0^2}} );

 &math( \varphi_{010}(\bm r)\propto ye^{-\frac{r^2}{2r_0^2}} =r\sin\theta\sin\phi e^{-\frac{r^2}{2r_0^2}} );

 &math( \varphi_{001}(\bm r)\propto ze^{-\frac{r^2}{2r_0^2}} =r\cos\theta e^{-\frac{r^2}{2r_0^2}} );

一方、

  Y_1{}^0\propto \cos\theta

  Y_1{}^{\pm 1}\propto \pm e^{i\phi}\sin\theta

より、

 &math( \varphi_{100}(\bm r)\propto re^{-\frac{r^2}{2r_0^2}}\,(Y_1{}^1-Y_1{}^{-1}) );

 &math( \varphi_{010}(\bm r)\propto re^{-\frac{r^2}{2r_0^2}}\,(Y_1{}^1+Y_1{}^{-1}) );

 &math( \varphi_{001}(\bm r)\propto re^{-\frac{r^2}{2r_0^2}}\,Y_1{}^0 );

すなわち、 N=1 に対応する \varphi_{n_xn_yn_z} はすべて、

  \frac{r}{r_0}e^{-\frac{r^2}{2r_0^2}}\,Y_1{}^m(\theta,\phi)

の形の関数の線形結合で表されることが分かる。

同じ固有値に属する固有関数の線形結合はやはり固有関数であるため、 これらの \varphi に対して、

  \hat H\varphi=(1+3/2)\hbar\omega\varphi=\frac{5}{2}\hbar\omega\varphi

  \hat{\bm l}^2\varphi=1(1+1)\hbar\omega\varphi=2\hbar\omega\varphi

が成り立つ。 \hat{\bm l}^2=2\hbar より、これらは p 状態である。

一方、 \hat l_z に対しては、

  \hat l_z\varphi_{001}=0\hbar \varphi_{001}=0

であるものの、 \varphi_{100},\varphi_{010} \hat l_z の固有関数ではない。 これら2つについて l_z を測定すれば、 1/2 の確率で \hbar が、 1/2 の確率で -\hbar が、観測されることになる。

\varphi_{100},\varphi_{010} はそのままでは \hat l_z の固有関数ではないが、

  \varphi_{100}+i\varphi_{010}\propto Y_1{}^1

  \varphi_{100}-i\varphi_{010}\propto Y_1{}^{-1}

とすることにより \hat l_z のそれぞれ m=1,-1 の固有関数となる。

物理的に見れば、 z 方向に振動する \varphi_{001} については z 軸周りの角運動量はゼロであり、 x,y 方向に振動する \varphi_{100},\varphi_{010} についてもそれら単独ではやはり z 軸周りの角運動量の期待値はゼロである。

しかし、 x 方向に振動しつつ同時に y 方向にも振動する場合、 その位相によっては z 軸周りの角運動量が生じることになる。

ここでは \pm i=e^{\pm \pi/2} をかけて足しており、 物理的には x 方向の振動に対して y 方向の振動の位相が \pm \pi/2 だけずれていることに対応する。 すなわち、 x \sin\omega t 的な振動をするとき y \pm\cos\omega t 的に振動することになり、 これはすなわち z 軸の周りの左回り/右回りの円運動に他ならない。

$N=2$ および $N=3$ の解

N=2 の解のうち、

  \varphi_{200},\ \varphi_{020},\ \varphi_{002} l=0 の解と l=2 の解の線形結合で、

  \varphi_{110},\ \varphi_{011},\ \varphi_{101} l=2 の解のみで

それぞれ表せる。すなわち後者は d 状態に対応するが、前者は s 状態と d 状態の混合、 すなわち \hat{\bm l}^2 の固有関数ではなく、 全角運動量を測定すれば 1/3 の確率で l=0 を、2/3 の確率で l=2 を取る。

N=3 の解のうち、

  \varphi_{111} は純粋な f 状態であるが、

 他は、p 状態と f 状態の混合となっている。

詳しい計算はこちら


前の単元 <<<              量子力学Ⅰ              >>> 次の単元

質問・コメント





Counter: 58043 (from 2010/06/03), today: 4, yesterday: 0