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

更新


量子力学Ⅰ

電子のエネルギーとポテンシャルエネルギー

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

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

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

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

が成り立った。

一方、量子力学では

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

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

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

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

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

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

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

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

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

  e^{ikx}

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

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

  k=i\kappa

と置けば、

  e^{ikx}=e^{-\kappa x}

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

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

harmonic2.png   harmonic1.png


特徴:

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

トンネル現象

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

no-tunnel.png

図のように、電子が自身のエネルギーよりも高いエネルギー障壁

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

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

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

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

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

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

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

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

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

確率密度の流れを伴う「定常状態」

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

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

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

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

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

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

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

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

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

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

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

各領域における波動関数

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

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

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

また、反射波と透過波を

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

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

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

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

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

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

tunneling-waves.png

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

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

波動関数の連続性

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

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

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

potential-step.png

波動関数の空間微分 \varphi'(x) が連続であることは、 V(x) 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) );

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

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

  \varphi'(b)-\varphi'(a)=0

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

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

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

境界条件

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

  • t に対して連続
  • x に対してその1回微分が連続

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

波動関数は、

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

として、

左端:  \psi_L(0,t)=\psi_M(0,t),\ \psi_L'(0,t)=\psi_M'(0,t)
右端:  \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} );

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

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

このとき \kappa/k=\lambda と置けば、

 (1)'  1+R=B_-+B_+

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

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

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

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

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

したがって、

  B_-=-\frac{1+i\lambda}{1-i\lambda}e^{2\kappa a}B_+

  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} );

となり、 |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\\ );

グラフ

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

tunnel.gif

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

波動関数の絶対値の二乗は以下のようになる。

tunnel.png

波長よりも厚い障壁に対しては、透過率は障壁厚さに対して指数関数的に減少する。

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

transmission-vs-width.png

質問・コメント





Counter: 46173 (from 2010/06/03), today: 11, yesterday: 0