量子力学Ⅰ/箱の中の自由粒子 の履歴(No.9)

概要†

シュレーディンガー方程式を解くという問題は数学的に難しいため、 特に解析的な関数として解が得られる例は非常に限られている（それ以外は数値的に、あるいは近似的に解くことになる）。

この節と次の節では最も簡単な２つの例について時間に依存しないシュレーディンガー方程式を解き、 定常的な解を求めながら、波動関数の一般的性質 について学ぶ。

時間に依らないシュレーディンガー方程式：

　&math( \left(-\frac{\hbar^2}{2m}\nabla^2+V(\bm r)\right)\psi(\bm r)=E\psi(\bm r) );

演習：１次元の箱の中の自由粒子†

を正の定数として、 の領域に閉じ込められた電子の定常状態を考える。

このような状況は、上記の範囲内で 、範囲外で と仮定することで実現され、井戸型ポテンシャルの問題とも呼ばれる。

(1) 箱の中の領域での、時間に依存しない一次元シュレーディンガー方程式を書け。 （一次元なので 、 とすればよい。 また、上記の通り箱の中では である。）

(2) の形の波動関数が (1) の解となることを確かめ、 の値を定めよ。 この式は２つの任意パラメータを含むから、(1)の２次微分方程式に対する一般解である。

(3) 箱の外では となる。その理由を答えよ。

(4) 波動関数の連続性より、箱の内側でも壁面において でなければならない。 この境界条件を満たすために に課される条件を 任意の自然数を表わす変数 を用いて書け。 *1通常は の１次微分 も連続となるが、ここではポテンシャルの無限大の不連続性を反映して微係数は不連続になる

以下、 で指定される を 、 対応する波動関数を 、 対応するエネルギー固有値を と書く。 このように離散化した固有値や固有関数を指定する指標 は 量子数と呼ばれる。

(5) を求めよ。

(6) を求めよ。ただし、規格化すること。

(7) に対応する、 時間に依存するシュレーディンガー方程式の解 を求めよ。

非定常状態の解†

非定常状態の例として、 における波動関数を

　

と置いてみる。右図は を表わしており、 中心から左に寄った位置にピークを持つ確率分布になることが分かる。 係数 は規格化因子で、 は第１種変形ベッセル関数である。

の中に現れるそれぞれの の時間発展は上で見たとおり

　

であり、シュレーディンガー方程式は線型であるから、上の初期条件で与えられる波動関数の時間発展は

　

で与えられる。

の時間変化をグラフに示せば、

右は縦軸を時間軸として表示したもので、確率密度を色で示した。 赤線は 座標の期待値を表わす。

　

電子が２つの壁にはね返り、振動する様子が見て取れる。

このように、シュレーディンガー方程式を満たす関数群 の一次結合で初期状態を表わせるならば、その時間発展は容易に求められる。

上記のように得られた解は において振動を続けるが、 もしこれが現実の電子の運動であれば、荷電粒子の加速度運動は光子の放出を伴うため、 徐々にエネルギーを失い最終的には時間に依存しない状態に落ち着くはずである。 そのような現象について議論するには相対論を取り入れた波動方程式を解かなければならない。 上で学んだシュレーディンガー方程式は相対論的効果の少ない条件でのみ成り立つ近似的な方程式である。 *2特殊相対論の効果を入れた波動方程式は知られており、電子の運動による光子の放出などを記述できる。一方、一般相対論と整合性を持つような、重力まで扱える量子力学の方程式は未だ確立されていない。

１次元の箱の中の自由粒子（有限ポテンシャル）†

次に、有限のポテンシャル障壁 に挟まれた一次元の箱の中の自由粒子を考える。 粒子の持つエネルギーは であり、古典的には一切箱の外には出られない。

無限のポテンシャル障壁では壁面で となることが境界条件となったが、
有限のポテンシャル障壁では壁面から外へ波動関数がしみ出すため境界条件が変化する。
この様子を見ていこう。

箱の中の解は上と同様に として、

　

箱の外での時間に依らないシュレーディンガー方程式は障壁高さを とすれば、

　

　

より、

　

ただし、 が一般解となる。 よりこれは２つの指数関数の重ね合わせだから、 で となるためには で 、 で でなければならない。

すなわち、''自身のエネルギーを越える障壁中に染み出した粒子の波動関数は 距離と共に指数関数的に減衰する''。 減衰定数は となる。

および で波動関数 とその１次微分 がどちらも連続であるという条件の下、 を決定すると 箱の中心を基準として 的な解と、 的な解の２種類が現れる。

それぞれに対して に対する条件は

　 ただし

および

　 ただし

となる。(→ 詳しい導出過程)

と をグラフに表わすと、

縦軸の の値によって には一般に有限個の解が得られ、 それらが離散化した波数、ひいては離散化したエネルギー固有値を与える。 グラフから が大きくなるほど解の個数が増えていくことが分かる。

のときは がすべて解になる。 これらの解は の場合に求めた に一致する。

例えば においては は５つの解を持ち、その形状は次のようになる。

対比のため、 の場合も合わせて示した。 波動関数は、そのエネルギー期待値だけ上方にオフセットして表示してある。

ポテンシャルが無限の場合と比べると、 ポテンシャルが有限の場合には障壁から左右に染み出す分だけ箱内部の の値が小さくなり、 エネルギーも低下する。しかし両者で波動関数の特徴は一致しており、 特に障壁エネルギーよりもずっと小さなエネルギーに対応する波動関数では差が小さい。

上記の計算は例えば金属中に閉じ込められた電子のモデルと考えることができる。 金属中には結晶の周期でイオン核が存在するにも関わらず、 電子は自由電子に運動できることを固体物理学で学ぶ。 金属柱に比べて真空中は電子の感じるエネルギーが高いため、 金属の端面は上で見たポテンシャル障壁と見なすことができる。

金属物体の端面は必ずしも急峻でもないし、平面的でもないため、 そのようなポテンシャル形状に対する境界条件は非常に複雑な物となる。 ただ上でも見たとおり、一般に「境界条件」はシュレーディンガー方程式の解を定量的に変えるのみで、 定性的には変化させしない。そこで、不定形の金属の境界条件を箱形であると仮定したり、 あるいは周期的であると仮定したり（箱の一方の端がもう一方と繋がっていることを仮定）して 問題を解いたとしても、端面近傍や極端に波数の長い解の性質を除いて得られる結果は現実の物と 変わらない。そのようにして問題を易しくして解くことが良く行われる。

３次元の箱の中の自由粒子†

を正の定数として、 の領域に閉じ込められた電子の定常状態を考える。

　

のように変数分離できることを仮定すれば、

　&math( &-\frac{\hbar^2}{2m}\left(

 \frac{\PD^2}{\PD x^2}+\frac{\PD^2}{\PD y^2}+\frac{\PD^2}{\PD z^2}\right)X(x)Y(y)Z(z)\\

&=-\frac{\hbar^2}{2m}\left[

 \left(\frac{\PD^2}{\PD x^2}X(x)\right)Y(y)Z(z) + 
 X(x)\left(\frac{\PD^2}{\PD y^2}Y(y)\right)Z(z) +
 X(x)Y(y)\left(\frac{\PD^2}{\PD z^2}Z(z)\right)

\right]\\ &=EX(x)Y(y)Z(z) );

　&math( \left(\frac{\PD^2}{\PD x^2}X(x)\right)\frac{1}{X(x)} + \left(\frac{\PD^2}{\PD y^2}Y(y)\right)\frac{1}{Y(y)} + \left(\frac{\PD^2}{\PD z^2}Z(z)\right)\frac{1}{Z(z)} =\frac{-2mE}{\hbar^2} );

左辺の各項はそれぞれ のみの関数であり、右辺は定数である。 任意の に対してこの式が成り立つためには、左辺の各項が定数でなければならない。

すなわち、

　&math( &\left(\frac{\PD^2}{\PD x^2}X(x)\right)\frac{1}{X(x)} = \frac{-2mE_x}{\hbar^2}\\ &\left(\frac{\PD^2}{\PD y^2}Y(y)\right)\frac{1}{Y(y)} = \frac{-2mE_y}{\hbar^2}\\ &\left(\frac{\PD^2}{\PD z^2}Z(z)\right)\frac{1}{Z(z)} = \frac{-2mE_z}{\hbar^2}\\ &E_x+E_y+E_z=E );

に対する方程式は１次元の箱形ポテンシャルの問題に帰着して、

　

　

等の解を得る。 の解は量子数 により指定できて、

　

　

となる。

例えば電子( ) を に閉じ込めれば、ゼロ点エネルギーは となる。

次の準位は である。

このように異なる量子数に対応する波動関数のエネルギーが等しいとき、 それらの準位は縮退していると言う。この様子を示したのが下図左である。

ではこのうちいくつかの縮退が解けて、準位の分裂が生じる。 としたときのエネルギー準位と、分裂前の縮退した準位との関係を下図右に示した。

１次元の調和振動子†

調和振動子のポテンシャルは であるから、時間に依存しないシュレーディンガー方程式は

　&math( \left(-\frac{\hbar^2}{2m}\frac{d^2}{dx^2}+\frac{k}{2}x^2\right)\psi(x)=E\psi(x) );

このような方程式を解く場合には、変数を無次元化するのが常套手段である。 すなわち、長さの次元を持つ自由変数 を変数変換して、無次元の量 で記述する。 ここでは、

　 ,

と置くと良い。ただし、 は古典論から得られる角振動数である。 すると与式は、

　&math( \left(-\frac{d^2}{d\xi^2}+\xi^2-\lambda\right)\psi(\xi)=0 );

となる。 の大きなところでは となるから、 そこでは は近似的に次の方程式を満たす。

　&math( \frac{d^2}{d\xi^2}\psi(\xi)=\xi^2\psi(\xi) );

ここから予想されるのは、

　&math( \psi(\xi)=H(\xi)e^{\pm\xi^2/2} );

という解の形である。系が 付近に束縛されていることから、 複号は負を取る。

　&math( &-\frac{d^2}{d\xi^2}\big[H(\xi)e^{-\xi^2/2}\big]+\xi^2H(\xi)e^{-\xi^2/2}-\lambda H(\xi)e^{-\xi^2/2}\\ &=-\frac{d}{d\xi}\big[H'(\xi)e^{-\xi^2/2}-\xi H(\xi)e^{-\xi^2/2}\big]+\xi^2H(\xi)e^{-\xi^2/2}-\lambda H(\xi)e^{-\xi^2/2}\\ &=-H''(\xi)e^{-\xi^2/2}+2\xi H'(\xi)e^{-\xi^2/2}+H(\xi)e^{-\xi^2/2}-\lambda H(\xi)e^{-\xi^2/2}\\ &=0\\ );

両辺を で割れば、

　&math( H''(\xi)=2\xi H'(\xi)+(1-\lambda) H(\xi) );

を得る。 と置いて代入すれば、

　&math( \sum_{l=0}^\infty l(l-1)c_l\xi^{l-2}=2\xi \sum_{l=0}^\infty l c_l\xi^{l-1}+(1-\lambda) \sum_{l=0}^\infty c_l\xi^l );

より において、

　

　

を得る。この式によれば、 を適当に決めると が、 を適当に決めると が、 それぞれすべて決まることになる。

あるいは あるいは が成立すれば、 それより大きな に対して がゼロになるが、 そうでない限り がゼロになることはない。

がゼロにならない場合、 において

　

が成り立つ。これは

　

とした時の係数の比と同じであり、このようになっていては が でゼロに収束するという境界条件を満たさない。

すなわち、 あるいは のどちらかがゼロであり、 もう一方と同じ偶奇性（パリティ）を持つある において

&math(\lambda=2l+1); が成立することが境界条件から要求され、

その結果 となる項は有限個となる。

• のとき ,
• のとき ,
• のとき ,
• のとき ,
• のとき ,
• ・・・

ここで現れた多項式 はエルミートの多項式と呼ばれる。

３次元の調和振動子†