球対称井戸型ポテンシャル のバックアップソース(No.15)

更新

[[前の単元 <<<>量子力学Ⅰ/水素原子]]
             [[量子力学Ⅰ]]             
[[>>> 次の単元>量子力学Ⅰ/3次元調和振動子]]~

* 目次 [#gad10ece]

#contents
&katex();

* 球形の箱の中の粒子 [#h9e0cd19]

3D の箱型ポテンシャル中の粒子について考える。

 &math(
V(r)=\begin{cases}
0&(r<=a)\\
\infty&(r>a)\\
\end{cases}
);

この場合には、&math(\chi(r)=rR(r)); を考えるよりも &math(R(r)); をそのまま扱った方が都合がよい。

 &math(\rho=\sqrt{\frac{2m\varepsilon}{\hbar^2}}r);

と置くことにより、箱の内部の方程式は

 &math(
\frac{d^2R}{d\rho^2}+\frac{2}{\rho}\frac{dR}{d\rho}+\left\{1-\frac{l(l+1)}{\rho^2}\right\}R=0
);

となる。&math(r\to\infty); での極限形は、

 &math(
\frac{d^2R}{d\rho^2}+\frac{2}{\rho}\frac{dR}{d\rho}+R=0
);

を満たすはずで、

 &math(
\Big(\frac{\sin\rho}\rho\Big)'
=-\frac{\sin\rho}{\rho^2}+\frac{\cos\rho}\rho
);

 &math(
\begin{aligned}
\Big(\frac{\sin\rho}\rho\Big)^{\prime\prime}
&=\frac{2\sin\rho}{\rho^3}-\frac{2\cos\rho}{\rho^2}-\frac{\sin\rho}\rho\\
&=\frac2\rho\Big(\frac{\sin\rho}\rho\Big)'-\frac{\sin\rho}\rho
\end{aligned}
);

を参考にすると、&math(R\to A\frac{\sin\rho}\rho+B\frac{\cos\rho}\rho); 
となればよいことが分かる。
実際この形は &math(l=0); の解になっている。

そこで &math(l\ne 0); では、

 &math(\Big(\frac1{\rho}\ \text{の多項式}\Big)\frac{\cos\rho}{\rho}+\Big(\frac1{\rho}\ \text{の多項式}\Big)\frac{\sin\rho}{\rho});

の形の解が存在することが予想される。

特にこのうち原点で発散しない解は''球ベッセル関数'' &math(j_l(\rho)); と呼ばれる(発散するものは''球ノイマン関数'' &math(n_l(\rho)); と呼ばれる)。

 &math(j_l(\rho)=(-\rho)^l\left(\frac{1}{\rho}\frac{d}{d\rho}\right)^l\frac{\sin\rho}{\rho},\hspace{1cm}n_l(\rho)=-(-\rho)^l\left(\frac{1}{\rho}\frac{d}{d\rho}\right)^l\frac{\cos\rho}{\rho});

 &math(j_0(\rho)=\frac{\sin\rho}{\rho});

 &math(j_1(\rho)=\frac1\rho \frac{\sin\rho}{\rho}-\frac{\cos\rho}{\rho});

 &math(j_2(\rho)=\left(\frac{3}{\rho^2}-1\right)\frac{\sin\rho}\rho-\frac{3}{\rho}\frac{\cos\rho}\rho);

 &math(j_3(\rho)=\left(-\frac{6}{\rho}+\frac{15}{\rho^3}\right) \frac{\sin\rho}\rho+\left(1-\frac{15}{\rho^2}\right)\frac{\cos\rho}\rho);

 ...

[[詳しい導出はこちら>@量子力学Ⅰ/球対称井戸型ポテンシャル/メモ#q4e20ed3]]

** 特徴 [#mf1aed18]

&math(|\rho j_l(\rho)|^2); をプロットした。

&attachref(SphericalBesselJ2.png,,66%);

- &math(|\rho j_l(\rho)|^2); は &math(\rho); の大きいところで、
-- &math(l=2k+0); なら &math(\sin^2\rho);
-- &math(l=2k+1); なら &math(\cos^2\rho);
- &math(\sin); や &math(\cos); の周期性を反映して &math(j_l(\rho)=0); を満たす根を無限個持つ
- &math(\sin); の項と &math(\cos); の項がうまく打消し合い、原点で発散しないようになっている
-- &math(j_0(\rho)=1);
-- &math(l\ge 1); では &math(j_l(\rho)=0);
-- 任意の &math(l\ge 0); に対して &math(\rho j_l(\rho)=0);
-- 実際、テーラー展開してみると原点付近で &math(\rho j_l(\rho)\propto\rho^{l+1}); になっている。
- &math(\rho); の小さいところでは &math(l); が大きいほどゆっくり振動する
-- &math(|\rho j_{l+2}|^2); は &math(|\rho j_l|^2); に比べて振動回数が1回少なくなる

** 境界条件 [#g9f280dc]

今まで学んできたとおり、微分方程式の一般解に境界条件を要求することでエネルギーが量子化するのであった。
この問題における境界条件は 

- &math(r=0); で発散しない(さもないと原点においてシュレーディンガー方程式が満たされない)
- &math(r=a); でゼロとなる

である。

上記微分方程式の一般解は、原点で発散しない球ベッセル関数 &math(j_l(\rho)); と、原点で発散する球ノイマン関数 &math(n_l(\rho)); の線形結合で表されるが、

 &math(R(r)=Aj_l(\rho)+Bn_l(\rho));

&math(B=0); と置いたことにより、すでに原点での境界条件は満たしている。

次に &math(r=a); でゼロになることを要求するのであるが、
上記の「箱の内部の方程式」には見た目上エネルギー固有値が現れず、動かせるパラメータがないように見える。よく見てみると、
&math(r); と &math(\rho); との変換式にエネルギーが押し込められてしまっているのだ。

そこで &math(r); と &math(\rho); との比を変えることで境界条件を満たすことを考えると、

 &math(j_l\Big(\sqrt{\frac{2m\varepsilon}{\hbar^2}}a\Big)=0);

により &math(\varepsilon); が決定されることになる。すなわち、&math(n); 番目の根を &math(\rho_n{}^l); として、

 &math(\sqrt{\frac{2m\varepsilon_n{}^l}{\hbar^2}}a=\rho_n{}^l);

が成り立つようエネルギー &math(\varepsilon_n{}^l); を決めれば良いことになるが、
このとき &math(r\in[0,a]); に対する &math(R(r)); は、
&math(\rho\in[0,\rho_n{}^l]); に相当するから、グラフは次のようになる。
ただし、見やすいように最大値で規格化した。横軸は &math(r/a); である。

&math(\varphi(r));~
&attachref(SphericalBesselJScaled.png);

&math(|r\varphi(r)|^2);~
&attachref(SphericalBesselJScaled3.png);

&math(rR(r)); と一次元井戸型ポテンシャルの解との類似性に注意せよ。
- &math(l=0); については一次元井戸型ポテンシャルの解と完全に一致する
- &math(l>0); については原点付近の存在確率が下がり、外側に寄っていく

** エネルギー固有値 [#h0cab4d2]

エネルギー固有値は、&math(n); 番目の根を &math(\rho_n{}^l); として

 &math(\varepsilon_n{}^l=\frac{\hbar^2}{2m}\left(\frac{\rho_n{}^l}{a}\right)^2\propto\big(\rho_n{}^l\big)^2);

のように、根の位置の2乗に比例する。

&math(n); と &math(l); との大小関係に制約はないから、任意の &math(n\ge 1); に対して
任意の &math(l\ge 0); が対応する。

水素の時には $l<n$ だったのとの違いは何だろうか?

ここでの &math(n); は動径方向の運動量を表す量子数で、
[[量子力学Ⅰ/水素原子#z46f54fd]] における &math(n'); に相当するから、
&math(l+n); としたものが水素原子の時の &math(n); と同等である。

水素の時に要求された $l<n$ という条件は、ここでは $l<l+n$ に相当して、それはまあここでも当然成り立っているわけである。

実際に値を入れてみると、

 &math(\varepsilon_1{}^0<\varepsilon_1{}^1<\varepsilon_1{}^2<\varepsilon_2{}^0<\varepsilon_1{}^3<\varepsilon_2{}^1<\dots);

となる。

#ref(spherical-box-energies.svg,right,around,400x271);

|&math(l);|&math(n);|&math(l+n);|対応|&math((\rho_n{}^l)^2);|
|0|1| 1|1s|9.86959|
|1|1| 2|2p|20.1907|
|2|1| 3|3d|33.2175|
|0|2| 2|2s|39.4785|
|3|1| 4|4f|48.8312|
|1|2| 3|3p|59.6795|
|4|1| 5|5g|66.9543|
|2|2| 4|4d|82.7192|
|0|3| 3|3s|88.8265|
|3|2| 4|4f|108.516|
|1|3| 4|4p|118.9|

エネルギー準位の大小関係を視覚化するため、横軸に &math(l); を、
縦軸に &math((\rho_n{}^l)^2); を取り、&math(n=1,2,3,\dots); に対応する値をプロットした。

ここから、&math(n); に比べて &math(l); に対するエネルギーの増加が少ないことが分かる。

クーロンポテンシャルでは角運動量が増加して軌道が外に寄ると、
その分、ポテンシャルエネルギーも増えたが、
箱型ポテンシャルではポテンシャルが一定であるために角運動量の増加は純粋な運動エネルギーの増加としてしか寄与しない。これがクーロンポテンシャルに比べて &math(l); に対するエネルギー増加が小さい原因となっている。

その結果、$n+l$ が等しい準位も大きく異なるエネルギーを持つ(クーロンポテンシャルに対して見られた縮退は解けている)。

** 箱の外のポテンシャルが有限の場合 [#s9bc8294]

箱の外のポテンシャルが有限の場合に、箱の外にも波動関数が漏れ出す。

 &math(V(r)=\begin{cases}
0&(r\le a)\\
V_0&(a<r)\\
\end{cases});

箱の中も外も、シュレーディンガー方程式の形はエネルギー固有値が異なるだけ(中では &math(\varepsilon);、外では &math(\varepsilon-V_0); となる)なので、
どちらも球ベッセル関数 &math(j_l(\rho)); と球ノイマン関数 &math(n_j(\rho)); の線形結合が一般解を与える。

ただし、箱の外では &math(\rho=\sqrt{\frac{2m(\varepsilon-V_0)}{\hbar^2}}\,r); の平方根の中味が負になるため、&math(\rho); は純虚数となる。これは一次元運動で見た、運動エネルギーが負になり、波数が虚数になる状況に対応する。

箱の中については &math(r); と &math(\rho); の関係も、原点における境界条件も変らないため、上と同様に球ベッセル関数が箱の中の波動関数を与える。ただし &math(r=a); においては &math(R(a)=0); となるのではなく、箱の外の波動関数と滑らかに繋がることが条件となるためエネルギー固有値は変化する(実際には低下する)。

箱の外の波動関数は原点における発散があるものでも構わず、すなわち球ノイマン関数の係数はゼロでなくても構わない。一方で、無限遠までの積分で値が収束しなければならないため &math(R(\infty)=0); でなければならない。

ところが箱の外では &math(\rho); が純虚数となるから、

 &math(\cos\rho=\frac{e^{i\rho}+e^{-i\rho}}{2}\propto\cosh|\rho|=\frac{e^{|\rho|}+e^{-|\rho|}}{2}); 

や 

 &math(\sin\rho=\frac{e^{i\rho}-e^{-i\rho}}{2i}\propto\cosh|\rho|=\frac{e^{|\rho|}-e^{-|\rho|}}{2}); 

が成り立ち、&math(r\to\infty); において &math(j_l(\rho),n_l(\rho)); はどちらも指数関数的に発散してしまう。

そこで &math(j_l(\rho)); と &math(n_j(\rho)); との線形結合により、
それぞれ第1種および第2種球ハンケル関数と呼ばれる

 &math(
h^{(1)}(\rho)=j_l(\rho)+i\,n_l(\rho)=-i(-\rho)^l\Big(\frac1d\frac d{d\rho}\Big)^le^{i\rho}\ \ \propto \ \ e^{i\rho} \ \sim\ e^{-|\rho|}
);

 &math(
h^{(2)}(\rho)=j_l(\rho)-i\,n_l(\rho)=i(-\rho)^l\Big(\frac1\rho\frac d{d\rho}\Big)^le^{-i\rho}\ \ \ \propto \ \ e^{-i\rho} \ \sim\ e^{+|\rho|}
);

を作ると、このうち &math(r\to\infty); すなわち &math(\rho\to \infty i); において &math(h^{(1)}(\rho)); はゼロに収束する一方、&math(h^{(2)}(\rho)); は発散してしまう(ただし &math(r_0=\frac{\hbar}{\sqrt{2m(V_0-\varepsilon)}});)。
すなわち、外部の解は第1種球ハンケル関数 &math(h^{(1)}(\rho)); で書かれる。

内側の球ベッセル関数と外側の第1種ハンケル関数とが &math(r=a); で滑らかに接続する条件

 &math(Aj_l(a-0)=Bh^{(1)}(a+0)\hspace{2cm}\text{(連続)});

 &math(A{j_l}'(a-0)=B{h^{(1)}}'(a+0)\hspace{1.8cm}\text{(滑らか)});

により、両者の振幅比 &math(A/B); とエネルギーが決まる。

1次元箱型ポテンシャルの問題で見たのと同様に、
&math(V_0); が小さくなり閉じ込めが弱くなると、同じ &math(n,l); に対するエネルギーは低下する。
このときポテンシャルエネルギーの期待値はむしろ増加するから、
エネルギーの低下は運動エネルギーの低下によるものである。

~
[[前の単元 <<<>量子力学Ⅰ/水素原子]]
             [[量子力学Ⅰ]]             
[[>>> 次の単元>量子力学Ⅰ/3次元調和振動子]]~

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

#article_kcaptcha

Counter: 20705 (from 2010/06/03), today: 1, yesterday: 0