量子力学Ⅰ/３次元調和振動子 の履歴(No.8)

３次元調和振動子†

原点からの距離に比例する大きさの、原点方向に引き戻す力を受ける粒子の運動を考える。

　

　

であるから、

　(1)

　(2)

と書ける。

以前、３次元調和振動子について学んだ際には、(1) の表式から について変数分離をして、１次元の調和振動子の積として解けることを見た。

　

ただし、 はそれぞれエルミート多項式 を使って、

　

の形に書けるから、

　

ただし、 と表される。このとき、

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

ただし、 と置いた。

一方、３次元調和振動子は上の の形を見ても明らかな通り球対称な中心力中における運動であるから、 これを球座標で展開して解くことにより、

　

と表せるはずである。ただしこのとき、

　

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

$x,y,z$ 座標形で求めた解と球座標解との対応†

３次元調和振動子を 座標で説いた結果得られた解をエネルギーの低いものから並べれば、

などとなる。

すなわち、 のそれぞれに対して 独立な解が 個存在することになる。

一方、同じ系を球座標で解けばその角度依存性は先に見たとおり球面調和関数で表され、

のように、 の解は縮退がなく、 の解は３重に、 は５重に縮退するはずである。 また、水素原子のところでも見たとおり、ポテンシャルの形状によっては に対する縮退以外にも、異なる を持つ状態が縮退することがある。

これはどういうことだろう？

• 縮退のない は s 状態であるべき
  • もし p 状態なら縮退した３つの独立な解があるはず！
• ３重に縮退した は３重に縮退した p 状態であるか、あるいは３つの s 状態が縮退した状態である（今の場合前者が正解）
• ６重に縮退した は縮退のない s 状態と５重に縮退した d 状態が縮退した状態、あるいは３重に縮退した p 状態が２つ縮退した状態である（今の場合前者が正解）
• １０重に縮退した は３重に縮退した p 状態と７重に縮退した f 状態が縮退した状態、あるいは・・・

にそれぞれ対応していることが期待され、以下に見るように確かにそのようになっている。

$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}} );

は に依存しない。

一方、 も に依存していないので、

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

と書ける。

すなわち の解は s 状態であり、 に対して であることが分かる。

$n=1$ の解†

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

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

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

ここで、

　&math( \begin{cases} x=r\sin\theta\cos\phi\\ y=r\sin\theta\sin\phi\\ z=r\cos\theta \end{cases} );

でることに注意せよ。

一方、

　

　

であったから、

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

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

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

すなわち、 に対応する はすべて、

　

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

同じ固有値に属する固有関数の線形結合は、やはり同じ固有値に属する固有関数になるから、 に対する の固有関数の線形結合として表されるこれらの関数は、 やはり の固有関数になる。すなわち p 状態である。

はそのまま の に対する固有関数でもある。

はそのままでは の固有関数ではないが、

　

　

とすることにより のそれぞれ の固有関数となる。

物理的には、 方向に振動する については 軸周りの角運動量はゼロであり、 方向に振動する についてもそれら単独ではやはり 軸周りの角運動量の期待値はゼロである。

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

ここでは をかけて足しており、 これは物理的には 方向の振動に対して 方向の振動の位相が だけずれており、なおかつ両方向の振幅が等しいことに対応する。 すなわち、 が 的な振動をするとき が 的に振動することになり、 これはすなわち 軸の周りの左回り／右回りの円運動に他ならない。

$n=2$ および $n=3$ の解†

上では、

　&math(\begin{array}{lll} r (Y_1{}^{1}-Y_1{}^{-1})&=r\phi_{px}&\propto x\\ r i(Y_1{}^{1}+Y_1{}^{-1})&=r\phi_{py}&\propto y\\ r Y_1{}^0&=r\phi_{pz}&\propto z\\ \end{array} );

を用いて式変形を行ったが、これと同様に、

　&math( \begin{array}{lll} r^2\, i(Y_2{}^{2}-Y_2{}^{-2}) &=r^2\varphi_{dxy} &\propto xy\\ r^2\, i(Y_2{}^{1}+Y_2{}^{-1})&=r^2\varphi_{dyz} &\propto yz\\ r^2\, (Y_2{}^{1}-Y_2{}^{-1}) &=r^2\varphi_{dzx} &\propto zx\\ r^2\, (Y_2{}^{2}+Y_2{}^{-2})&=r^2\varphi_{d(x^2-y^2)} &\propto x^2-y^2\\ r^2\, Y_2{}^{0} &=r^2\varphi_{d(3z^2-r^2)}&\propto 3z^2-r^2 \end{array} );

などの関係を用いることにより、

• ６重に縮退した は s 状態と d 状態が１つずつ縮退した状態

• １０重に縮退した は p 状態と f 状態の縮退した状態

になっていることを確かめられる。具体的には、

の解は、

　 　と、 　

の線形結合で、

の解は、

　 　と、 　

の線形結合で、それぞれ表せる。詳しい計算はこちら

演習：３次元調和振動子の動径方向の方程式†

上で出てきた動径方向の関数が、対応する の値に対して３次元調和振動子の方程式の解になっていることと、 その固有値が与えられた から予想されるエネルギー値に等しいことを確かめたい。

(1) ３次元調和振動子に対する動径方向の方程式は、 と書き換え、 と置けば、

&math(

• \frac{X''(\xi)}{X(\xi)}+\xi^2+\frac{l(l+1)}{\xi^2}=\frac{2\varepsilon}{\hbar\omega}\\ );

となることを示せ。

(2) となることを確かめよ。

(3) が に対する固有関数となることを示せ。

(4) が に対する固有関数となることを示せ。

前の単元 ＜＜＜ 　　　　　　　　　　　　　量子力学Ⅰ　　　　　　　　　　　　　 ＞＞＞ 次の単元

質問・コメント†