球対称井戸型ポテンシャル のバックアップ差分(No.14)
更新- バックアップ一覧
- 現在との差分 を表示
- ソース を表示
- バックアップ を表示
- 量子力学Ⅰ/球対称井戸型ポテンシャル へ行く。
- 追加された行はこの色です。
- 削除された行はこの色です。
[[前の単元 <<<>量子力学Ⅰ/水素原子]] [[量子力学Ⅰ]] [[>>> 次の単元>量子力学Ⅰ/3次元調和振動子]]~ * 目次 [#gad10ece] #contents * 球形の箱の中の粒子 [#h9e0cd19] 3D の箱型ポテンシャル中の粒子について考える。 &math( V(r)=\begin{cases} 0&(r<=a)\\ V_0&(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(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}); &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=0); で発散しない(さもないと原点においてシュレーディンガー方程式が満たされない) - &math(r=a); でゼロとなる である。 上記微分方程式の一般解は球ベッセル関数 &math(j_l(\rho)); と球ノイマン関数 &math(n_l(\rho)); の和で表されるが、 上記微分方程式の一般解は、原点で発散しない球ベッセル関数 &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(\varepsilon); が決定されることになる。すなわち、&math(n); 番目の根を &math(\rho_n{}^l); として、 このときグラフは次のようになる。 見やすいように最大値で規格化した。横軸は &math(r/a); である。 &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(|\varphi(r)|^2);~ &attachref(SphericalBesselJScaled2.png); &math(|r\varphi(r)|^2);~ &attachref(SphericalBesselJScaled3.png); &math(rR(r)); と一次元井戸型ポテンシャルの解との類似性に注意せよ。 - &math(l=0); については一次元井戸型ポテンシャルの解と完全に一致する - &math(l>0); については原点付近の存在確率が下がるが、原点から遠いところではやはり一次元の解と近くなる - &math(l>0); については原点付近の存在確率が下がる ** エネルギー固有値 [#h0cab4d2] &math(n); 番目の根を &math(\rho_n{}^l); とすれば、 エネルギー固有値は、&math(n); 番目の根を &math(\rho_n{}^l); として &math(\sqrt{\frac{2m\varepsilon}{\hbar^2}}a=\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(\varepsilon_n{}^l=\frac{\hbar^2}{2m}\left(\frac{\rho_n{}^l}{a}\right)^2); として、エネルギーの大きさは根の位置の2乗に比例する &math(\varepsilon_n{}^l\propto\big(\rho_n{}^l\big)^2);。 実際に値を入れてみると、 &math(\varepsilon_1{}^0<\varepsilon_1{}^1<\varepsilon_1{}^2<\varepsilon_2{}^0<\varepsilon_1{}^3<\varepsilon_2{}^1<\dots); となる。 &math(n); と &math(l); との大小関係に制約はないから、任意の &math(n\ge 0); に対して 任意の &math(l\ge 0); が対応する。 ここでの &math(n); は [[量子力学Ⅰ/水素原子#z46f54fd]] における &math(n'); に相当するから、 &math(l+n); が水素原子の時の &math(n); と同等である。 #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); に対するエネルギー増加が小さい原因となっている。 ** 箱の外のポテンシャルが有限の場合 [#s9bc8294] 箱の外のポテンシャルが有限の場合にも、&math(r=a); における位相が少しずれるものの、 外部の解(第1種球ハンケル関数となる)となめらかに接続する条件から &math(\varepsilon); が決定される。 箱の外のポテンシャルが有限の場合に、箱の外にも波動関数が漏れ出す。 閉じ込めが弱くなると、同じ &math(n,l); に対するエネルギーは低下する。 &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} ); &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} ); を作ると、このうち &math(r\to\infty); において &math(h^{(1)}(ir/r_0)\propto e^{-r/r_0}/r); はゼロに収束してくれる一方、&math(h^{(2)}(ir/r_0)\propto e^{r/r_0}/r); は発散してしまう(ただし &math(r_0=\frac{\hbar}{\sqrt{2m(V_0-\varepsilon)}});)。 すなわち、外部の解は第1種球ハンケル関数 &math(h^{(1)}(\rho)); で書かれる。 両者が &math(r=a); で滑らかに接続する条件(&math(R(a-0)=R(a+0),R'(a-0)=R'(a+0));) によりエネルギーが決まる。 &math(V_0); が小さくなり閉じ込めが弱くなると、同じ &math(n,l); に対するエネルギーは低下する。 このときポテンシャルエネルギーの期待値はむしろ増加するから、 エネルギーの低下は運動エネルギーの低下によるものである。 これは1次元箱型ポテンシャルの問題で見たのと同じ現象である。 ~ [[前の単元 <<<>量子力学Ⅰ/水素原子]] [[量子力学Ⅰ]] [[>>> 次の単元>量子力学Ⅰ/3次元調和振動子]]~ * 質問・コメント [#z49eb7a9] #article_kcaptcha
Counter: 20792 (from 2010/06/03),
today: 1,
yesterday: 0