HOME > natural science Laboratory > コンピュータ・シミュレーション講座 > 仮想物理実験室

【量子力学再入門7】
箱型ポテンシャル障壁によるトンネル効果

文責:遠藤 理平 (2018年5月 4日) カテゴリ:仮想物理実験室(325)計算物理学(165)

量子力学とは原子や電子といったナノスケールにおける粒子の運動を記述する物理学の一分野です。 メートルスケールで成り立つ「ニュートンの運動方程式」に対して、量子力学では「シュレディンガー方程式」と呼ばれる基礎方程式が成り立つと考えられています。本項では量子力学の基本的な問題を数値的に解くことで、量子力学の理解を深めることが目的とします。7回目のテーマは「箱型ポテンシャル障壁によるトンネル効果」です(波束の箱型ポテンシャル障壁への衝突は次回行います)。

箱型ポテンシャル障壁によるトンネル効果とは

古典力学で対象とする質点の運動において、物体が持っている力学的エネルギー(運動エネルギー+ポテンシャルエネルギー)よりも大きなポテンシャル障壁が存在する場合、外からの力が加わらない限りその障壁を越えることはありません。しかしながら、波の場合には禁止された領域に染み出すことは様々な物理系で実験的にも理論的にも認められています。最も有名な例は光が境界面で全反射する際に反対側に染み出すエバネッセント光です。この現象は、波動現象一般で生じるため、波動方程式であるシュレディンガー方程式の解にも存在します。
トンネル効果とは量子力学的粒子が持つ運動エネルギーよりも大きなポテンシャル障壁が存在する系にて、ポテンシャル障壁を突き抜ける現象を指します。

ポテンシャル障壁中の波動関数

ポテンシャル障壁中の波動関数は前項で示したとおり、

\psi_{\rm II}(x,t) = \left[A^{(+)}_{\rm II}e^{ik_{\rm II}x} + A^{(-)}_{\rm II}e^{-ik_{\rm II}x}\right] e^{-i\omega t}

です。量子力学的粒子の運動エネルギーがポテンシャル障壁の大きさよりも小さい場合、量子力学的粒子の波数が純虚数となります。つまり、 k_{\rm II}= i\kappa_{\rm II}とおいて、波動関数は

\psi_{\rm II}(x,t) = \left[A^{(+)}_{\rm II}e^{-\kappa_{\rm II}x} + A^{(-)}_{\rm II}e^{\kappa_{\rm II}x}\right] e^{-i\omega t}

となります。この表式はxとともに指数関数的に減少する項と増大する項の和で表されることを意味します。そのため、障壁の端からの距離に対して指数関数的に減衰するわけではありません。ここでは示しませんが、ポテンシャル障壁の高さが高いほど、幅が広いほど A^{(-)}_{\rm II} は小さくなり、実質的には第1項目のみが寄与します。

シミュレーション結果

箱型ポテンシャル障壁でトンネル効果が生じている単一平面波の様子です。左より入射してます。

20180503.gif

計算パラメータ

ポテンシャルの高さ: V2 = 20.0 * eV
平面波のエネルギー:E0 = 10.0 * eV;
空間スケール 1E-13[m];//横軸の値


考察と次の課題

・境界面でもなめらかに接続する様子がわかります。
・絶対値が時間に依存しないのは定常解であるためです。
・次は波束を箱型ポテンシャル障壁へ衝突させます。

プログラムソース(C++)

http://www.natural-science.or.jp/files/physics/QuantumPhysics7.cpp
※VisualC++、GCC(MinGW)で動作確認しています。



▲このページのトップNPO法人 natural science トップ

関連記事

仮想物理実験室







計算物理学

▲このページのトップNPO法人 natural science トップ




Warning: mysqli_connect(): (28000/1045): Access denied for user 'xsvx1015071_ri'@'sv102.xserver.jp' (using password: YES) in /home/xsvx1015071/include/natural-science/include_counter-d.php on line 8
MySQL DBとの接続に失敗しました