logo

階段ポテンシャルに対するシュレーディンガー方程式の解 📂量子力学

階段ポテンシャルに対するシュレーディンガー方程式の解

概要

5D5500710.jpg

ポテンシャルが上の図のように階段関数であるとき、粒子がどのように運動するかを調べてみよう。ポテンシャルUU

U(x)={0x<0U0x>0 U(x) = \begin{cases} 0 & x<0 \\ U_{0} & x>0 \end{cases}

ポテンシャルがU(x)U(x)であるときの時間に無関係なシュレーディンガー方程式は

d2u(x)dx2+2m2[EU(x)]u(x)=0 \dfrac{d^2 u(x)}{dx^2}+\frac{2m}{\hbar ^2} \Big[ E-U(x) \Big]u(x)=0

解法1

E<0E<0

エネルギーがポテンシャルより小さい場合、解が存在しないため考慮する必要はない。

0<E<U00 < E < U_{0}

  • Part 2-1. x<0x<0

    この領域でのシュレーディンガー方程式は次のようになる。

    d2udx2+2m2Eu=0 \dfrac{d^2 u}{dx^2}+\frac{2m}{\hbar^2}Eu=0

    2m2E\frac{2m}{\hbar^2}Eが正なので、k2k^2で置き換えると

    d2udx2+k2u=0 \dfrac{d^2 u}{dx^2}+k^2u=0

    これは非常に簡単な2階常微分方程式である。微分方程式を解くと、その解は次のようになる。

    uI(x)=A+eikx+Aeikx u_{I}(x)=A_{+}e^{ikx} + A_{-}e^{-ikx}

    ここで、A+A_{+}AA_{-}は定数だ。ここで、u=A+eikxu=A_{+}e^{ikx}入射波AeikxA_{-}e^{-ikx}反射波である。

  • Part 2-2. x>0x>0

    この領域でのシュレーディンガー方程式は次の通り。

    d2udx2+2m2(EU0)u=0 \dfrac{d^2 u}{dx^2}+\frac{2m}{\hbar^2}(E-U_{0})u=0

    EU0<0E-U_{0}<0なので、2m2(EU0)=κ2\frac{2m}{\hbar^2}(E-U_{0})=-\kappa ^2で置き換えると

    d2udx2κ2u=0 \dfrac{d^2 u}{dx^2}-\kappa^2 u=0

    ここで、κ\kappaはギリシャ文字の「カパ」だ。kk(ケイ)とは違う文字だ。この微分方程式の解は次の通り

    uII(x)=B+eκx+Beκx u_{II}(x) = B_{+}e^{\kappa x}+B_{-}e^{-\kappa x}

    波動関数は二乗積分可能である必要があるので、limxuII(x)=0\lim \limits_{x \rightarrow \infty}u_{II}(x)=0でなければならない。ところが、limxB+eκx=\lim \limits_{x \rightarrow \infty} B_{+}e^{\kappa x}=\infty なので、B+=0B_{+}=0でなければならない。したがって、透過波utrans=Beκxu_{trans}=B_{-}e^{-\kappa x}

  • Part 2-3. 境界条件

    波動関数は滑らかだと仮定するので、x=0x=0では連続しており、波動関数の微分(勾配)もx=0x=0で連続している。したがって次の条件が得られる。

    uI(0)=uII(0)    A++A=BuI(0)=uII(0)    ikA+ikA=κB \begin{align*} u_{I}(0)=u_{II}(0) \quad \implies& \quad A_{+}+A_{-} = B_{-} \\ u_{I}^{\prime}(0)=u_{II}^{\prime}(0) \quad \implies& \quad ikA_{+}-ikA_{-} = -\kappa B_{-} \end{align*}

    二つの式を連立させると次のようになる。

    AA+=ik+κikκ \frac{A_{-}} {A_{+}} = \frac{ ik+\kappa}{ik -\kappa }

  • Part 2-4. 反射率透過率

    反射率と透過率を計算するために、入射波、反射波、透過波の確率流を求めてみると

jinc=kmA+2,jref=kmA2,jtrnas=0 j_{inc}=\frac{\hbar k}{m}|A_{+}|^2,\quad j_{ref}=\frac{\hbar k}{m}|A_{-}|^2,\quad j_{trnas}=0

したがって、反射率、透過率は

R=jrefjinc=AA+2=κ2+k2κ2+k2=1 R=\left| \frac{j_{ref}}{j_{inc}}\right|=\left| \frac{A_{-}} {A_{+}} \right|^2=\frac{\kappa^2+k^2}{\kappa^2 + k^2}=1

T=jtrnasjinc=0kmA+2=0 T=\left| \frac{ j_{trnas} }{j_{inc}}\right|=\frac{ 0}{\frac{\hbar k}{m}|A_{+}|^2}=0

上の結果により、反射のみが起こり、透過は起こらないという事実が分かる。しかし、厳密にはB0B_{-} \ne 0なので、x>0x>0で粒子が発見される確率が存在する。つまり、巨視的に見ると入射した波動は全て反射されたと見なせるが、微視的に見ると入射波の一部が古典的には不可能な領域に透過したことになる。これをトンネル効果tunnelingまたは量子トンネル効果quantum tunnel effectという。

U0<E U_{0} < E

  • Part 3-1. x<0x<0

    この領域での時間に無関係なシュレーディンガー方程式は

    d2udx2+2m2Eu=0 \dfrac{d^2 u}{dx^2}+\frac{2m}{\hbar^2}Eu=0

    2m2E\frac{2m}{\hbar^2}Eが正なので、k2k^2で置き換えると

    d2udx2+k2u=0 \dfrac{d^2 u}{dx^2}+k^2u=0

    上記の式は**Part 2-1.**で既に解いたので、解を知っている。

    uI(x)=A+eikx+Aeikx u_{I}(x)=A_{+}e^{ikx} + A_{-}e^{-ikx}

    入射波はuinc=A+eikxu_{inc}=A_{+}e^{ikx}、反射波はuref=Aeikxu_{ref}=A_{-}e^{-ikx}

  • Part 3-2. x>0x>0 この領域での時間に無関係なシュレーディンガー方程式は

    d2udx2+2m2(EU0)u=0 \dfrac{d^2 u}{dx^2}+\frac{2m}{\hbar^2}(E-U_{0})u=0

    EU0>0E-U_{0}>0なので、2m2(EU0)=κ2\frac{2m}{\hbar^2}(E-U_{0})=\kappa ^2で置き換えると

    d2udx2+κ2u=0 \dfrac{d^2 u}{dx^2}+\kappa ^2u=0

    **Part 3-1.**と同じ形なので解は

    uII(x)=B+eiκx+Beiκx u_{II}(x)=B_{+}e^{i \kappa x} + B_{-}e^{-i \kappa x}

    この時、x>0x>0領域において反射波がないのでB=0B_{-}=0だ。したがって

    uII(x)=B+eiκx u_{II}(x)=B_{+}e^{i \kappa x}

    これが透過波である。

  • Part 3-3. 境界条件

    波動関数は滑らかだと仮定するので、x=0x=0では連続しており、波動関数の微分(勾配)もx=0x=0で連続している。したがって

    uI(0)=uII(0)    A++A=B+ u_{I}(0)=u_{II}(0) \quad \implies \quad A_{+}+A_{-} = B_{+}

    uI(0)=uII(0)    ikA+ikA=iκB+ u_{I}^{\prime}(0)=u_{II}^{\prime}(0) \quad \implies \quad ikA_{+}-ikA_{-} = i\kappa B_{+}

    上記二つの式を連立させると次のようになる。

    AA+=kκk+κ \frac{A_{-}} {A_{+}} = \frac{ k-\kappa}{k +\kappa }

    B+A+=2kk+κ \frac{B_{+}} {A_{+}} = \frac{ 2k}{k +\kappa }

  • Part 3-4. 反射率透過率

    入射波、反射波、透過波の確率流は

    jinc=kmA+2,jref=kmA2,jtrnas=κmB+2 j_{inc}=\frac{\hbar k}{m}|A_{+}|^2,\quad j_{ref}=\frac{\hbar k}{m}|A_{-}|^2,\quad j_{trnas}=\frac{\hbar \kappa}{m}|B_{+}|^2

    したがって、反射率、透過率は

    R=jrefjinc=AA+2=(kκk+κ)2=k2kκ+κ2(k+κ)2 R=\left| \frac{j_{ref}}{j_{inc}}\right|=\left| \frac{A_{-}} {A_{+}} \right|^2=\left( \frac{k-\kappa}{ k + \kappa }\right) ^2= \frac{k^2 -k\kappa +\kappa^2}{ (k + \kappa )^2 }

    T=jtransjinc=κkB+A+2=κk(2kk+κ)2=2kκ(k+κ)2 T=\left| \frac{j_{trans}}{j_{inc}}\right|=\frac{\kappa}{k}\left| \frac{B_{+}} {A_{+}} \right|^2=\frac{\kappa}{k} \left( \frac{2k}{ k + \kappa }\right) ^2= \frac{2k \kappa }{ (k + \kappa )^2 }

    R+T=1R+T=1が成り立つことは簡単に確認できる。


  1. Stephen Gasiorowicz, 量子物理学(Quantum Physics, 西江大学物理学科共訳) (3rd Edition, 2005), p79-83 ↩︎