一次元箱形障壁のトンネル のバックアップ差分(No.8)

更新


  • 追加された行はこの色です。
  • 削除された行はこの色です。
[[前の単元 <<<>量子力学Ⅰ/3次元調和振動子]]
             [[量子力学Ⅰ]]             
[[>>> 次の単元>量子力学Ⅰ/一次元の散乱現象]]~
#contents

* 電子のエネルギーとポテンシャルエネルギー [#ha47b60e]

電子のエネルギーはポテンシャルエネルギーと運動エネルギーの和であるから、次式が成り立つ。

 &math(\varepsilon = V(x) + \frac{p^2}{2m});

このため古典力学においては、常に

(1) (電子のエネルギー) > (ポテンシャルエネルギー)
 (1) (電子のエネルギー) > (ポテンシャルエネルギー)

が成り立った。

一方、量子力学では

(2) (電子のエネルギー) < (ポテンシャルエネルギー)
 (2) (電子のエネルギー) < (ポテンシャルエネルギー)

となる領域にも有限の確率密度を取り得ることを見てきた。

 &math(p^2=\hbar^2k^2=2m\{\varepsilon-V(r)\});

と書きなおせば分かるとおり、

 &math(k=\pm\frac{1}{\hbar}\sqrt{2m\{\varepsilon-V(r)\}});

であり、(1) では √ 内部が正であるが、(2) では負である。つまり、

(1) では運動エネルギーは正である。これは数学的には運動量が実数であるためである。

(2) では運動エネルギーは負である。すなわち運動量は虚数になる。

これらに対応して、(1) に対応する位置においては波数 &math(k=p/\hbar); は実数であり、

 &math(e^{ikx});

は振動する解を与える。
&math(k); はポテンシャルエネルギーの関数であるから、場所によって波長も異なる。

一方、(2) に対応する位置においては波数 &math(k=p/\hbar); は虚数であるから、

 &math(k=i\kappa);

と置けば、

 &math(e^{ikx}=e^{-\kappa x});

となり、指数関数的に減衰する/増加する解を与える。

次は調和振動子に対する波動関数を図示したものである(左が &math(\varphi(x)); 
右は &math(|\varphi(x)|^2);)。
上記の関係をこの図にあてはめて理解せよ。

&attachref(量子力学Ⅰ/調和振動子/harmonic2.png,,50%);  
&attachref(量子力学Ⅰ/調和振動子/harmonic1.png,,50%);

~
特徴:
- 二次曲線は &math(y=Kx^2/2); であり、古典的な調和振動子ではこの外には出られない
- 量子力学的な解は外側にも少しはみ出している
- &math(\varepsilon>V); の領域では振動する
- &math(\varepsilon<V); では指数関数的に減衰する
- &math(V(x)); に比べて &math(\varepsilon); が大きいほど波長が短く=波数が大きくなる
- &math(V(x)); に比べて &math(\varepsilon); が小さいほど早く減衰する

* トンネル現象 [#xa22c51e]

量子力学においては、電子が自身のエネルギーよりも高いポテンシャル中にも存在できるため、
「トンネル現象」と呼ばれる量子力学に特有の現象が生じる。

#ref(no-tunnel.png,right,around,25%);
図のように、電子が自身のエネルギーよりも高いエネルギー障壁 

 &math(
V(x)=\begin{cases}
0                  &(x<0,a\le x)\\
V_0\ \ >\varepsilon  &(0<x\le a)\\
\end{cases}
);

へ左から入射する場合を考えよう。

古典論ではこのような障壁は完全弾性障壁とみなせ、
入射方向と逆方向に入射時と同じ大きさの運動量を持って跳ね返される。

古典論でも電子はエネルギー障壁に跳ね返される(反射する)が、
確率密度の一部は自身よりも高いエネルギーをもつ障壁の中へ染み込む。

染み込んだ障壁中では確率密度が距離とともに指数関数的に減衰するが、
障壁の右端でも完全には零とならず、
その成分はエネルギー障壁を通り抜けて進む電子を表す確率密度となる。

すなわち量子論においては、電子が ''自身のエネルギーよりも高いエネルギー障壁を通り抜けて進む'' 
確率が存在する。

あたかもエネルギー障壁にトンネルを空けてその中を電子が通るかのようであるという意味で、
この現象は「トンネル現象」と呼ばれる。

現実の問題としては、わずかな間隙を挟んで平面的な金属電極が向かい合わされている状況が、
上図の状況に相当する。
電子にとって真空部分のポテンシャルエネルギーは金属内部のポテンシャルエネルギーよりも高く、
その部分がエネルギー障壁となる。
電極間距離が 1 nm 程度まで近づけば、計測可能な程度の「トンネル電流」が計測されることになる。

以下ではこの現象をシュレーディンガー方程式から理解しよう。

* 確率密度の流れを伴う「定常状態」 [#m88a76fc]

古典論には「定常状態」とはすべての粒子が静止している状態、
あるいは一定の運動を繰り返している状態を指す言葉であった。

量子論では確率密度の流れを伴う「定常状態」を考えることができる。

 &math(\psi(x,t)\propto e^{i(kx-\omega t)});

は、&math(x); 軸の正方向へ確定した運動量 &math(p_0=\hbar k_0); 
を持って進む電子を表す波動関数であった。

電子は静止しておらず、一定の運動を繰り返しているわけでもないが、
確率密度の空間分布は時間とともに変化しない。

このような状態を考える上で問題となるのは、この波動関数を規格化できないことである。
現実にはこのような確率密度の流れがあれば左遠方の電子の確率密度が減少し、
右遠方の確率密度が増加しなければならず、そのような電子の供給源、
吸収先が存在しないことが規格化できないという問題の根底にある。

それでも、上記のような波動関数を考えることには物理的な意味がある。

先に学んだように電子の波数が &math(k_0); に確定しておらず、
&math(k=k_0\pm\Delta k); の範囲にあるような波束を考えれば、
そのような波束は空間的に &math(\Delta x=\hbar/2\Delta k); 
程度かそれ以上に広がることを学んだ。

&math(\Delta k); が十分小さい場合、&math(\Delta x); は十分大きく、
そのような波束の中央部分ではほぼ &math(\psi(x,t)\propto e^{i(kx-\omega t)});
と同じ状態が実現する。

そのような「準定常状態」は上記のような「定常状態」の解とほとんど同じであろうというのがその理由である。

シュレーディンガー方程式は線形であるから、
「準定常状態」において生じる物理現象を議論する目的には
波束全体の規格化によって決まる波動関数の振幅を確定する必要はない。
入射波の振幅が2倍になれば、反射波や透過波の振幅も2倍になるだけである。

* 各領域における波動関数 [#n41f28fb]

そこで以下では、左から入射する電子の波動関数を

 &math(\psi_I(x,t)=e^{i(k_Ix-\omega_I t)});  
ただし、&math(\varepsilon_I=\hbar\omega_I=\frac{\hbar^2k_I^2}{2m}); 

と置く。(I は incident electron 入射電子 の頭文字)

また、反射波と透過波を

 &math(\psi_R(x,t)=Re^{i(-k_Rx-\omega_R t)});  
ただし、&math(\varepsilon_R=R\hbar\omega_0=\frac{\hbar^2k_R^2}{2m});

 &math(\psi_T(x,t)=Te^{i(k_Tx-\omega_T t)});  
ただし、&math(\varepsilon_T=T\hbar\omega_0=\frac{\hbar^2k_T^2}{2m});

と置く。(R は reflected electron 反射電子、T は transmitted electron 透過電子 の頭文字)

さらに、障壁内部をトンネルする電子を、

 &math(\psi_B(x,t)=B_-e^{-\kappa_B x-i\omega_B t}+B_+e^{\kappa_B x-i\omega_B t});  
ただし、&math(\varepsilon_B=\hbar\omega_B=V_0-\frac{\hbar^2\kappa_B^2}{2m});

と置く。(B は electron in barrier 障壁中の電子 の頭文字)

#ref(tunneling-waves.png,right,around,25%);

&math(k_I); が与えられたとして、
&math(R,k_R,T,k_T,B_1,B_2,\kappa_t); を境界条件から決定するのがここでの問題である。

&math(\varepsilon_I,\omega_I,\varepsilon_R,\omega_R,\varepsilon_T,\omega_T,\varepsilon_B,\omega_B);
は &math(k_I,k_R,k_T,\kappa_B); より上記の分散関係を用いて容易に求まることに注意せよ。

* 波動関数の連続性 [#zcdcd945]

上記の問題では障壁端面での波動関数の連続性が境界条件を与える。

シュレーディンガー方程式は &math(x); に対する2階微分を含んでいるから、
&math(V(x)); がなめらかである限り、波動関数は
「&math(x); で2階微分可能」になる。

ただし、上記の箱型障壁の端点ように &math(V(x)); が ''不連続に変化する点'' では、
「&math(x); に対する1階微分が連続」
ではあるが、「2階微分は不可能」になる。
そもそもその点では &math(V(x)); の値が決まらないので、
シュレーディンガー方程式自体が意味を持たないことに注意せよ。

#ref(potential-step.png,right,around,25%);

波動関数の空間微分 &math(\varphi'(x)); が連続であることは、
&math(V(x)); が &math(a<x<b); の範囲で連続に、
しかし急峻に変化する状況を考えれば証明できる。

この区間でシュレーディンガー方程式を積分すれば、

 &math(
\int_a^b\left(-\frac{\hbar^2}{2m}\frac{d^2}{dx^2}\varphi(x)+V(x)\varphi(x)-\varepsilon\varphi(x)\right)\,dx=0
);

 &math(
\left[\frac{\hbar^2}{2m}\frac{d}{dx}\varphi(x)\right]_a^b
=\int_a^b\Big(V(x)-\varepsilon\Big)\varphi(x)\,dx
);

 &math(
\frac{\hbar^2}{2m}\Big(\varphi'(b)-\varphi'(a)\Big)
&=\overline{\Big(V(x)-\varepsilon\Big)\varphi(x)}\int_a^b\,dx\\
&=\overline{\Big(V(x)-\varepsilon\Big)\varphi(x)}\cdot (b-a)
);

ただし、&math(\overline{\Big(V(x)-\varepsilon\Big)\varphi(x)}); は区間 &math([a,b]);
における &math(\Big(V(x)-\varepsilon\Big)\varphi(x)); の平均値である。

&math(V(x)); が非常に急峻に変化しており、&math(a); と &math(b); を十分に近く取れる場合には
&math(b-a=0); とみなせるから、

 &math(\varphi'(b)-\varphi'(a)=0);

すなわち、&math(V(x)); が不連続に変化する点を挟んで、&math(\varphi'(x)); が連続であることが示される。

ただし、無限大の深さをもつ緯度型ポテンシャルの時のように不連続点の片側で &math(V(x)); 
が &math(\pm\infty); となる場合には、
「&math(\Big(V(x)-\varepsilon\Big)\varphi(x)); の平均値」が有限とならないため、
&math(\varphi'(b)-\varphi'(a)); の値も不定となる。

すなわち、そのような点では &math(\varphi'(x)); は不連続になりうる。

* 境界条件 [#nd265f03]

障壁の左端と右端とで波動関数が

- &math(t); に対して連続
- &math(x); に対してその1回微分が連続

となる条件から、上記のすべての未知変数を決定できる。

波動関数は、

障壁の左側: &math(\psi_L(x,t)=\psi_I(x,t)+\psi_R(x,t));~
障壁の内部: &math(\psi_M(x,t)=\psi_B(x,t));~
障壁の右側: &math(\psi_R(x,t)=\psi_T(x,t));

として、

左端: &math(\psi_L(0,t)=\psi_M(0,t),\ \psi_L'(0,t)=\psi_M'(0,t));~
右端: &math(\psi_M(a,t)=\psi_R(a,t),\ \psi_M'(a,t)=\psi_R'(a,t));

となる。具体的には、

 (1) &math(
e^{i(k_I\cdot 0-\omega_I t)}+Re^{i(-k_R\cdot 0-\omega_R t)}=
B_-e^{-\kappa_B \cdot 0-i\omega_B t}+B_+e^{\kappa_B \cdot 0-i\omega_B t}
);

 (2) &math(
ik_Ie^{i(k_I\cdot 0-\omega_I t)}-ik_R Re^{i(-k_R\cdot 0-\omega_R t)}=
-\kappa_B B_-e^{-\kappa_B \cdot 0-i\omega_B t}+\kappa_B B_+e^{\kappa_B \cdot 0-i\omega_B t}
);

 (3) &math(
B_-e^{-\kappa_B  a-i\omega_B t}+B_+e^{\kappa_B  a-i\omega_B t}=
Te^{i(k_T a-\omega_T t)}
);

 (4) &math(
-\kappa_B B_-e^{-\kappa_B  a-i\omega_B t}+\kappa_B B_+e^{\kappa_B  a-i\omega_B t}=
ik_T Te^{i(k_T a-\omega_T t)}
);

となる。

(1) より、 &math(
1+Re^{-i(\omega_R-\omega_I) t}=\big(B_-+B_+)e^{-i(\omega_B-\omega_I) t}
);

ここから &math(\omega_R=\omega_B=\omega_I); でなければならないことが理解できる。
同様に、(3) より &math(\omega_T=\omega_B); であることが得られるから、
すなわち、&math(\varepsilon_R=\varepsilon_B=\varepsilon_T=\varepsilon_I); となり、
すべての箇所で電子のエネルギーは等しい。

ここから、&math(k_R=k_T=k_I); (&math(=k); と置く) および &math(\kappa_B=\sqrt{\frac{2mV_0}{\hbar^2}-k_I^2}); (&math(=\kappa); と置く) を得る。

このとき &math(\kappa/k=\lambda); と置けば、

 (1)' &math(1+R=B_-+B_+);

 (2)' &math(1 - R=-i\lambda \big(- B_-+B_+));

 (3)' &math(B_-e^{-\kappa  a}+B_+e^{\kappa  a}=Te^{ik a});

 (4)' &math(i\lambda B_-e^{-\kappa  a}-i\lambda B_+e^{\kappa  a}= Te^{ik a});

(1)' と (2)' より、&math(2=(1+i\lambda)B_-+(1-i\lambda)B_+);

(3)' と (4)' より、&math((1-i\lambda)e^{-\kappa a}B_-+(1+i\lambda)e^{\kappa a}B_+=0);

したがって、

 &math(B_-=-\frac{1+i\lambda}{1-i\lambda}e^{2\kappa a}B_+);

 &math(2=\left[-\frac{1+i\lambda}{1-i\lambda}e^{2\kappa a}+(1-i\lambda)\right]B_+);

 &math(
B_+&=\frac{2(1-i\lambda)e^{-\kappa a}}{(1-i\lambda)^2e^{-\kappa a}-(1+i\lambda)^2e^{\kappa a}}\\
&=\frac{-(1-i\lambda)e^{-\kappa a}}{2i\lambda\cosh\kappa a+(1-\lambda^2)\sinh\kappa a}
=\frac{-(1-i\lambda)e^{-\kappa a}}{X}\\);

 &math(
B_-=\frac{(1+i\lambda)e^{\kappa a}}{X}
);

より、

 &math(
R&=B_-+B_+-1\\
&=\frac{2i\lambda\cosh\kappa a+2\sinh\kappa a}{X}-1\\
&=\frac{(1+\lambda^2)\sinh\kappa a}{X}
);

 &math(
T&=\frac{(1+i\lambda)-(1-i\lambda)}{X}=\frac{i2\lambda}{X}
);

ここから、

 &math(
|X|^2&=4\lambda^2\cosh^2\kappa a+(1-\lambda^2)^2\sinh\kappa a\\
&=4\lambda^2(1+\sinh^2\kappa a)+(1-\lambda^2)^2\sinh\kappa a\\
&=4\lambda^2+\underbrace{(1+\lambda^2)^2\sinh^2\kappa a}_{4\lambda^2Y^2}\\
&=4\lambda^2(1+Y^2)
);

 &math(
|R|^2&=\frac{Y^2}{1+Y^2}\\
);

 &math(
|T|^2=\frac{1}{1+Y^2}
);

となり、&math(|R|^2+|T|^2=1); を確かめられる。

これは入射電子の振幅が 1 であることと対応している。
反射した量と透過した量を加えると入射した量に等しくなるのは期待通りと言える。

 &math(
Y^2&=\left(\frac{1+\lambda^2}{2\lambda}\right)^2\sinh^2\kappa a\\
&=\frac{k^2+\kappa^2}{4k^2\kappa^2}\sinh^2\kappa a\\
&=\frac{V_0^2}{4\epsilon(V_0-\epsilon)}\sinh^2\kappa a\\
);

* 波形 [#ra980378]

波動関数の実部をグラフにした結果は下記の通りになる。
ここでは &math(\hbar=1,\frac{\hbar^2}{2m}=1,k_I=1); とした上で、
&math(\epsilon=1,V_0=1.1,a=5); と置いた。

&attachref(tunnel.gif); &qr(http://dora.bk.tsukuba.ac.jp/~takeuchi/?plugin=attach&refer=%E9%87%8F%E5%AD%90%E5%8A%9B%E5%AD%A6%E2%85%A0%2F%E4%B8%80%E6%AC%A1%E5%85%83%E7%AE%B1%E5%BD%A2%E9%9A%9C%E5%A3%81%E3%81%AE%E3%83%88%E3%83%B3%E3%83%8D%E3%83%AB&openfile=tunnel.gif);

- &math(x<0); では進行波(赤)と反射波(青)とが干渉し、定在波が立っている
- 障壁内部と外部の波動関数はなめらかにつながる
- 障壁内部では振幅が急激に減衰する
- 透過波は入射波と同じ波数、同じ周期を持つが、振幅が減少し位相がずれている

波動関数の振幅の二乗は以下のようになる。障壁の左側では進行波と反射波の干渉を反映して振幅が波打つ。
障壁内では振幅が急速に減少し、右端で残った成分が透過波となる。

&attachref(tunnel.png);

&math(1/\kappa); よりも厚い障壁に対しては、透過率は障壁厚さに対して指数関数的に減少する。

 &math(|T|^2=\frac{1}{1+Y^2}\simeq \frac{1}{Y^2}\propto \sinh^2\kappa a\simeq (\exp 2\kappa a)/2);

&attachref(transmission-vs-width.png);

~
[[前の単元 <<<>量子力学Ⅰ/3次元調和振動子]]
             [[量子力学Ⅰ]]             
[[>>> 次の単元>量子力学Ⅰ/一次元の散乱現象]]~

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

#article_kcaptcha


Counter: 46166 (from 2010/06/03), today: 4, yesterday: 0