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

更新


前の単元 <<<              量子力学Ⅰ              >>> 次の単元

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

電子のエネルギーはポテンシャルエネルギーと運動エネルギーの和である:

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

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

 (1)  \varepsilon \ge V(x)

であったが、量子力学では

 (2)  \varepsilon < V(x)

となる領域にも有限の確率密度を取り得る。

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

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

であるから、

(1) の領域では運動量 p および波数 k=p/\hbar は実数であり、

  e^{ikx}

は振動する解を与える。

(2) の領域では運動量 p および波数 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倍になるだけである。

流れの大きさ

確率密度の流れの大きさは、

  S=\mathrm{Re}\left[\psi^*\frac{\hat p}{m} \psi\right]

と表せるのであった。上記の \psi(x,t)\propto e^{i(kx-\omega t)} を代入すれば、

 &math( S&=\mathrm{Re}\left[\psi^*\frac{\hbar}{im}\frac{d}{dx} \psi\right]\\ &=\frac{\hbar k}{m} \mathrm{Re}\left[\psi^* \psi\right]\\ &=\frac{\hbar k}{m} );

すなわち、 \psi は単位時間あたり、( yz 面内の)単位面積あたり、 \hbar k/m 個の電子に相当する確率密度の流れを伴うことが分かる。 これは \psi 自身の振幅が1であり、単位面積と粒子の速度 \hbar k/m との積で表される体積に含まれる電子の数が体積自身に等しいことによる。

各領域における波動関数

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

  \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 入射電子 の頭文字)

これは、 x 軸の正方向への電子の確率密度の流れを表しており、 単位時間あたり、また、 yz 面内の単位面積あたり、電子 \hbar k/m 個が移動する大きさの流れである。 以下の計算はすべてに入射波の振幅に比例するから、後から計算結果に適当な係数をかけてやることで 実験と比較することができる。

反射波と透過波を

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

  \psi_T(x,t)=Te^{i(k_Tx-\omega_T t)}    ただし、 \varepsilon_T=\hbar\omega_T=\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

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

  |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: 46166 (from 2010/06/03), today: 4, yesterday: 0