一次元箱形障壁のトンネル/メモ のバックアップソース(No.3)

• バックアップ一覧
• 差分 を表示
• 現在との差分 を表示
• バックアップ を表示
• 量子力学Ⅰ/一次元箱形障壁のトンネル/メモ へ行く。
  • 1 (2017-08-02 (水) 10:58:00)
  • 2 (2017-08-02 (水) 16:35:01)
  • 3 (2017-08-02 (水) 21:29:30)
  • 4 (2017-08-04 (金) 17:58:19)

* トンネル障壁再考 [#na9b6dec]

https://m-repo.lib.meiji.ac.jp/dspace/bitstream/10291/5021/1/kyouyoronshu_240_17.pdf

を参考にやってみる。

位置 &math(x_1); から始まる、
厚さ &math(a); 高さ &math(V); のトンネル障壁の手前で

　&math(\varphi(x)=Ae^{ikx}+Be^{-ikx});　&math((x<x_1));

直後で、

　&math(\varphi(x)=Ce^{ikx}+De^{-ikx});　&math((x_1+a<x));

とすると、障壁内での波動関数は

　&math(\kappa=\frac{1}{\hbar}\sqrt{2mV_0-\hbar^2k^2});

を用いて

　&math(\varphi(x)=Ee^{\kappa x}+Fe^{-\kappa x});　&math(x_1<x<x_1+a));

と表せることから、両端でなめらかに接続する条件を次の４つの式で表せる。

　(1) &math(\varphi(x_1)=Ae^{ikx_1}+Be^{-ikx_1}=Ee^{\kappa x_1}+Fe^{-\kappa x_1});

　(2) &math(\varphi'(x_1)=ikAe^{ikx_1}-ikBe^{-ikx_1}=\kappa Ee^{\kappa x_1}-\kappa Fe^{-\kappa x_1});

　(3) &math(\varphi(x_1+a)=Ee^{\kappa (x_1+a)}+Fe^{-\kappa (x_1+a)}=Ce^{ik(x_1+a)}+De^{-ik(x_1+a)});

　(4) &math(\varphi'(x_1+a)=\kappa Ee^{\kappa (x_1+a)}-\kappa Fe^{-\kappa (x_1+a)}
                       =ikCe^{ik(x_1+a)}-ikDe^{-ik(x_1+a)});

この４つから &math(E,F); を消去し、&math(C,D); を &math(A,B); で表す。

(2) に &math(\lambda=\kappa/k); を導入すると、

　&math(iAe^{ikx_1}-iBe^{-ikx_1}=\lambda Ee^{\kappa x_1}-\lambda Fe^{-\kappa x_1});

(1) と合わせて、

　&math(
\begin{pmatrix}
1&1\\
i&-i\\
\end{pmatrix}
\begin{pmatrix}Ae^{ikx_1}\\Be^{-ikx_1}\end{pmatrix}
=
\begin{pmatrix}
1&1\\
\lambda&-\lambda\\
\end{pmatrix}
\begin{pmatrix}Ee^{\kappa x_1}\\Fe^{-\kappa x_1}\end{pmatrix}
=
\begin{pmatrix}
1&1\\
\lambda&-\lambda\\
\end{pmatrix}
\begin{pmatrix}e^{\kappa x_1}&\\&e^{-\kappa x_1}\end{pmatrix}
\begin{pmatrix}E\\F\end{pmatrix}
);

(4) に &math(\lambda=\kappa/k); を導入すると、

　&math(\lambda Ee^{\kappa (x_1+a)}-\lambda Fe^{-\kappa (x_1+a)}
=iCe^{ik(x_1+a)}-iDe^{-ik(x_1+a)});

(3) と合わせて、

　&math(
\begin{pmatrix}
1&1\\
\lambda&-\lambda\\
\end{pmatrix}
\begin{pmatrix}
Ee^{\kappa (x_1+a)}\\Fe^{-\kappa (x_1+a)}
\end{pmatrix}=
\begin{pmatrix}
1&1\\
i&-i\\
\end{pmatrix}
\begin{pmatrix}
e^{\kappa (x_1+a)}&\\&e^{-\kappa (x_1+a)}
\end{pmatrix}
\begin{pmatrix}E\\F\end{pmatrix}=
\begin{pmatrix}
1&1\\
\lambda&-\lambda\\
\end{pmatrix}
\begin{pmatrix}
Ce^{ik(x_1+a)}\\De^{-ik(x_1+a)}\\
\end{pmatrix}
);

したがって、

　&math(
\begin{pmatrix}E\\F\end{pmatrix}&=
\begin{pmatrix}e^{\kappa x_1}&\\&e^{-\kappa x_1}\end{pmatrix}^{-1}
\begin{pmatrix}
1&1\\
\lambda&-\lambda\\
\end{pmatrix}^{-1}
\begin{pmatrix}
1&1\\
i&-i\\
\end{pmatrix}
\begin{pmatrix}Ae^{ikx_1}\\Be^{-ikx_1}\end{pmatrix}
\\&=
\begin{pmatrix}
e^{\kappa (x_1+a)}&\\&e^{-\kappa (x_1+a)}
\end{pmatrix}^{-1}
\begin{pmatrix}
1&1\\
\lambda&-\lambda\\
\end{pmatrix}^{-1}
\begin{pmatrix}
1&1\\
i&-i\\
\end{pmatrix}
\begin{pmatrix}
Ce^{ik(x_1+a)}\\De^{-ik(x_1+a)}\\
\end{pmatrix}
);

ここから、

　&math(
&\begin{pmatrix}
Ce^{ik(x_1+a)}\\De^{-ik(x_1+a)}\\
\end{pmatrix}

\\&=
\begin{pmatrix}
1&1\\
i&-i\\
\end{pmatrix}^{-1}
\begin{pmatrix}
1&1\\
\lambda&-\lambda\\
\end{pmatrix}
\begin{pmatrix}
e^{\kappa (x_1+a)}&\\&e^{-\kappa (x_1+a)}
\end{pmatrix}
\begin{pmatrix}e^{\kappa x_1}&\\&e^{-\kappa x_1}\end{pmatrix}^{-1}
\begin{pmatrix}
1&1\\
\lambda&-\lambda\\
\end{pmatrix}^{-1}
\begin{pmatrix}
1&1\\
i&-i\\
\end{pmatrix}
\begin{pmatrix}Ae^{ikx_1}\\Be^{-ikx_1}\end{pmatrix}

\\&=
\frac{1}{4i\lambda}
\begin{pmatrix}
-i&-1\\
-i&1\\
\end{pmatrix}
\begin{pmatrix}
1&1\\
\lambda&-\lambda\\
\end{pmatrix}
\begin{pmatrix}
e^{\kappa a}&\\&e^{-\kappa a}
\end{pmatrix}
\begin{pmatrix}
-\lambda&-1\\
-\lambda&1\\
\end{pmatrix}
\begin{pmatrix}
1&1\\
i&-i\\
\end{pmatrix}
\begin{pmatrix}Ae^{ikx_1}\\Be^{-ikx_1}\end{pmatrix}

\\&=
\frac{1}{4i\lambda}
\begin{pmatrix}
-(i+\lambda)e^{\kappa a}&-(i-\lambda)e^{-\kappa a}\\
-(i-\lambda)e^{\kappa a}&-(i+\lambda)e^{-\kappa a}\\
\end{pmatrix}
\begin{pmatrix}
-(i+\lambda)&i-\lambda\\
i-\lambda&-(i+\lambda)\\
\end{pmatrix}
\begin{pmatrix}Ae^{ikx_1}\\Be^{-ikx_1}\end{pmatrix}

\\&=
\frac{1}{4i\lambda}
\begin{pmatrix}
(i+\lambda)^2e^{\kappa a}-(i-\lambda)^2e^{-\kappa a}&
(1+\lambda^2)e^{\kappa a}-(1+\lambda^2)e^{-\kappa a}\\
-(1+\lambda^2)e^{\kappa a}+(1+\lambda^2)e^{-\kappa a}&
-(i+\lambda)^2e^{\kappa a}+(i-\lambda)^2e^{-\kappa a}\\
\end{pmatrix}
\begin{pmatrix}Ae^{ikx_1}\\Be^{-ikx_1}\end{pmatrix}
);

ここで、

　&math(
\alpha&=\frac{1}{4i\lambda}\Big[(i+\lambda)^2e^{\kappa a}-(i-\lambda)^2e^{-\kappa a}\Big]
=\frac{-\overline X}{2i\lambda}
);

　&math(
Y=\frac{1}{4\lambda}\big[(1+\lambda^2)e^{\kappa a}-(1+\lambda^2)e^{-\kappa a}\big]
=\frac{1+\lambda^2}{2\lambda}\sinh \kappa a
);

と置けば、

　&math(
\begin{pmatrix}Ce^{ik(x_1+a)}\\De^{-ik(x_1+a)}\end{pmatrix}=
\begin{pmatrix}
\alpha&-iY\\
iY&\overline\alpha
\end{pmatrix}
\begin{pmatrix}Ae^{ikx_1}\\Be^{-ikx_1}\end{pmatrix}\\
);

を得る。

以前の結果から、&math(|X|^2=4\lambda^2(1+Y^2)); すなわち、

　&math(|\alpha|^2=1+Y^2);

　&math(\left|\frac{\alpha}{\sqrt{1+Y^2}}\right|^2=1);

である。そこで、

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

　&math(\frac{\alpha}{\sqrt{1+Y^2}}=e^{i\theta});

と置けば、

　&math(
\begin{pmatrix}Ce^{ik(x_1+a)}\\De^{-ik(x_1+a)}\end{pmatrix}=\frac{1}{\sqrt{1-r^2}}
\begin{pmatrix}
e^{i\theta}&-ir\\
ir&e^{-i\theta}
\end{pmatrix}
\begin{pmatrix}Ae^{ikx_1}\\Be^{-ikx_1}\end{pmatrix}\\
);

トンネル問題を考え、&math(A=1, D=0); の下にこれを解けば

　&math(B=R=r\,e^{i2kx_1+i(\theta-\pi/2)});　→ 反射波

　&math(
Ce^{ik(x_1+a)}&=\frac{1}{\sqrt{1-r^2}}\left(
e^{ikx_1+i\theta}-r^2\,e^{ikx_1+i\theta}
\right)\\
&=\frac{1}{\sqrt{1-r^2}}e^{ikx_1+i\theta}\left(
1-r^2
\right)\\
&=\sqrt{1-r^2}\,e^{ikx_1+i\theta}
);

　&math(C=T=\sqrt{1-r^2}\,e^{-ikL+i\theta});　→ 透過波

となる。

&math(r=|R|); であり、&math(\theta); は反射波の位相の進みを表す。


これだけだと苦労した意味があまりないので、
トンネル障壁が複数ある場合について考える。

* ２重トンネル [#u0e0b560]

厚さ &math(a); の障壁が &math(b); だけ間を置いて連続して２つあるとする。

障壁前を &math(A,B);、障壁の間を &math(C,D);、障壁の後を &math(E,F); で表せば、

　&math(
&\begin{pmatrix}Ce^{ika}\\De^{-ika}\end{pmatrix}=\frac{1}{\sqrt{1-r^2}}
\begin{pmatrix}
e^{i\theta}&-ir\\
ir&e^{-i\theta}
\end{pmatrix}
\begin{pmatrix}A\\B\end{pmatrix}\\
);

　&math(
&\begin{pmatrix}Ee^{ik(2a+b)}\\Fe^{-ik(2a+b)}\end{pmatrix}=\frac{1}{\sqrt{1-r^2}}
\begin{pmatrix}
e^{i\theta}&-ir\\
ir&e^{-i\theta}
\end{pmatrix}
\begin{pmatrix}Ce^{ik(a+b)}\\De^{-ik(a+b)}\end{pmatrix}\\
);

より、

　&math(
\begin{pmatrix}Ee^{ik(2a+b)}\\Fe^{-ik(2a+b)}\end{pmatrix}
&=
\frac{1}{1-r^2}
\begin{pmatrix}
e^{i\theta}&-ir\\
ir&e^{-i\theta}
\end{pmatrix}
\begin{pmatrix}e^{ikb}&\\&e^{-ikb}\end{pmatrix}
\begin{pmatrix}
e^{i\theta}&-ir\\
ir&e^{-i\theta}
\end{pmatrix}
\begin{pmatrix}A\\B\end{pmatrix}\\
&=
\frac{1}{1-r^2}
\begin{pmatrix}
e^{ikb+i\theta}&-ire^{-ikb}\\
ire^{ikb}&e^{-ikb-i\theta}
\end{pmatrix}
\begin{pmatrix}
e^{i\theta}&-ir\\
ir&e^{-i\theta}
\end{pmatrix}
\begin{pmatrix}A\\B\end{pmatrix}\\
&=
\frac{1}{1-r^2}
\begin{pmatrix}
e^{ikb+2i\theta}+r^2e^{-ikb}&
-ire^{ikb+i\theta}-ire^{-ikb-i\theta}\\
ire^{ikb+i\theta}+ire^{-ikb-i\theta}&
r^2e^{ikb}+e^{-ikb-2i\theta}
\end{pmatrix}
\begin{pmatrix}A\\B\end{pmatrix}\\
);

&math(A=1,F=0); と置くと、

　&math(
|B|^2
&=r^2\left|\frac{e^{ikb+i\theta}+e^{-ikb-i\theta}}{r^2e^{ikb+i\theta}+e^{-ikb-i\theta}}\right|^2\\
&=\frac{4r^2\cos^2(kb+\theta)}{(1+r^2)^2\cos^2(kb+\theta)+(1-r^2)^2\sin^2(kb+\theta)}\\
);

より、透過率は

　&math(
|E|^2=1-|B|^2
&=\frac{(1+r^2)^2\cos^2(kb+\theta)+(1-r^2)^2\sin^2(kb+\theta)-4r^2\cos^2(kb+\theta)}{(1+r^2)^2\cos^2(kb+\theta)+(1-r^2)^2\sin^2(kb+\theta)}\\
&=\frac{(1-r^2)^2}{(1+r^2)^2\cos^2(kb+\theta)+(1-r^2)^2\sin^2(kb+\theta)}\\
&=\frac{1}{\left(\frac{1+r^2}{1-r^2}\right)^2\cos^2(kb+\theta)+\sin^2(kb+\theta)}\\
&=\frac{1}{\left[\left(\frac{1+r^2}{1-r^2}\right)^2-1\right]\cos^2(kb+\theta)+1}\\
);

となる。

この値は反射率が高く &math(r\sim 1); である場合にも、&math(\cos^2(kb+\theta)=0); 
すなわち &math(kb+\theta=\pi(n+1/2)); 
の条件で透過率が 100% になることを示している。これが共鳴トンネルと呼ばれる現象である。

このとき、２つの障壁に挟まれた井戸の中での確率分布は、

　&math(
&\begin{pmatrix}Ce^{ika}\\De^{-ika}\end{pmatrix}=\frac{1}{\sqrt{1-r^2}}
\begin{pmatrix}
e^{i\theta}&-ir\\
ir&e^{-i\theta}
\end{pmatrix}
\begin{pmatrix}A\\0\end{pmatrix}\\
);

より、

　&math(
&|\varphi(x)|^2=|Ce^{ikx}+De^{-ikx}|^2\\
&=\frac{1}{1-r^2}\Big[1-r^2
-ire^{i2k(x-a)+i\theta}
+ire^{-i2k(x-a)-i\theta}
\Big]\\
&=1+\frac{2r}{1-r^2}\sin\{2k(x-a)+\theta\}\\
);

このとき、

　&math(
|\varphi(a)|^2=1+\frac{2r}{1-r^2}\sin\{\theta\}
);

　&math(
|\varphi(a+b)|^2
&=1+\frac{2r}{1-r^2}\sin(2kb+\theta)\\
&=1+\frac{2r}{1-r^2}\sin\{2(kb+\theta)-\theta\}\\
&=1+\frac{2r}{1-r^2}\sin\{2\pi(n+1/2)-\theta\}\\
&=1+\frac{2r}{1-r^2}\sin\theta\\
);

となり、位相 &math(\theta); だけ染み出すものの、
ほぼ井戸の両端でゼロになっている。

すなわち、共鳴トンネルは入射波エネルギーが
井戸の準定常状態のエネルギーと一致することが条件となっている。

&math(k); が &math(n); 番目の共振波長を &math(k_n); に近いと近似してみると、

　&math(
|E|^2
&=\frac{1}{\left[\left(\frac{1+r^2}{1-r^2}\right)^2-1\right]\sin^2(b(k-k_n))+1}\\
&\sim\frac{1}{\left(\frac{1-r^2}{2br}\right)^{-2}(k-k_n)^2+1}\\
);

のようにローレンツ型の共振特性が得られる。

このときの半値幅は、&math(
\Delta k=\frac{1-r^2}{br}
); となる。

反射率が高く &math(r\sim 1); のとき、

　&math(\Delta k\sim\frac{1}{b}(障壁１つあたりの透過率)^2);

また、

　&math(
k-k_n&=\sqrt{2m\varepsilon/\hbar^2}-k_n\\
&=k_n\left(\sqrt{1+(\varepsilon-\varepsilon_n)/\varepsilon_n}-1\right)\\
&=k_n\big\{1+(\varepsilon-\varepsilon_n)/2\varepsilon_n-1\big\}\\
&=\frac{k_n}{2\varepsilon_n}(\varepsilon-\varepsilon_n)\\
);

より、

　&math(
|E|^2
&\sim\frac{1}{\left(\frac{1-r^2}{2r}\frac{2\varepsilon_n}{k_nb}\right)^{-2}(\varepsilon-\varepsilon_n)^2+1}\\
);

として、エネルギー表示においてもやはりローレンツ型の表示が可能である。

このときのエネルギー幅は、&math(
\Delta E
=\frac{1-r^2}{br}\frac{2\varepsilon_n}{k_n}
=\frac{1-r^2}{br}\sqrt{\frac{2\hbar^2\varepsilon_n}{m}}
);
で、反射率が高いときには &math((1-r^2)/r); はほぼ透過率と等しくなるから、

　&math(\Delta E
=\frac{2\varepsilon_n}{k_nb}(障壁１つあたりの透過率)^2
=\sqrt{\frac{2\hbar^2\varepsilon_n}{b^2m}}(障壁１つあたりの透過率)^2
);

となる。

これらのローレンツ型の式は、あくまで &math(k); や &math(\varepsilon); 
が共振点から離れていないことを前提としており、&math(r); の &math(k); 依存性や、
&math(\sin(bk-bk_n)\sim b(k-k_n)); などの１次近似の範囲内でのみ成立することに注意が必要である。

　&math(
(障壁１つあたりの透過率)=\frac{4\epsilon(V_0-\epsilon)}{V_0^2}\exp\left[-2a\sqrt{\frac{2m}{\hbar^2}(V_0-\varepsilon)}\right]
);

なので、エネルギーが高くなればエネルギー幅は広くなる。

* 非対称２重トンネル [#s83e4299]

　&math(
\begin{pmatrix}Ee^{ik(a+b+a')}\\Fe^{-ik(a+b+a')}\end{pmatrix}
&=
\frac{1}{\sqrt{1-r^2}\sqrt{1-r'^2}}
\begin{pmatrix}
e^{i\theta}&-ir\\
ir&e^{-i\theta}
\end{pmatrix}
\begin{pmatrix}e^{ikb}&\\&e^{-ikb}\end{pmatrix}
\begin{pmatrix}
e^{i\theta'}&-ir'\\
ir'&e^{-i\theta'}
\end{pmatrix}
\begin{pmatrix}A\\B\end{pmatrix}\\
&=
\frac{1}{\sqrt{1-r^2}\sqrt{1-r'^2}}
\begin{pmatrix}
e^{ikb+i\theta}&-ire^{-ikb}\\
ire^{ikb}&e^{-ikb-i\theta}
\end{pmatrix}
\begin{pmatrix}
e^{i\theta'}&-ir'\\
ir'&e^{-i\theta'}
\end{pmatrix}
\begin{pmatrix}A\\B\end{pmatrix}\\
&=
\frac{1}{\sqrt{1-r^2}\sqrt{1-r'^2}}
\begin{pmatrix}
e^{ikb+i(\theta+\theta')}+rr'e^{-ikb}&
-ir'e^{ikb+i\theta}-ire^{-ikb-i\theta'}\\
ire^{ikb+i\theta'}+ir'e^{-ikb-i\theta}&
rr'e^{ikb}+e^{-ikb-i(\theta+\theta')}
\end{pmatrix}
\begin{pmatrix}A\\B\end{pmatrix}\\
);

&math(A=1,F=0); と置いて、

　&math(
|B|^2
&=\left|\frac{ire^{ikb+i\theta'}+ir'e^{-ikb-i\theta}}{rr'e^{ikb}+e^{-ikb-i(\theta+\theta')}}
\right|^2\\
&=\left|\frac{re^{ikb+i(\theta+\theta')/2}+r'e^{-ikb-i(\theta+\theta')/2}}
{rr'e^{ikb+i(\theta+\theta')/2}+e^{-ikb-i(\theta+\theta')/2}}
\right|^2\\
&=\frac{(r+r')^2\cos^2\{kb+(\theta+\theta')/2\}+(r-r')^2\sin^2\{kb+(\theta+\theta')/2\}}
{(rr'+1)^2\cos^2\{kb+(\theta+\theta')/2\}+(rr'-1)^2\sin^2\{kb+(\theta+\theta')/2\}}\\
&=\frac{\{(r+r')^2-(r-r')^2\}\cos^2\{kb+(\theta+\theta')/2\}+(r-r')^2}
{\{(rr'+1)^2-(rr'-1)^2\}\cos^2\{kb+(\theta+\theta')/2\}+(rr'-1)^2}\\
&=\frac{4rr'\cos^2\{kb+(\theta+\theta')/2\}+(r-r')^2}
{4rr'\cos^2\{kb+(\theta+\theta')/2\}+(1-rr')^2}\\
);

　&math(
|E|^2=1-|B|^2
&=\frac{(1-rr')^2-(r-r')^2}
{4rr'\cos^2\{kb+(\theta+\theta')/2\}+(1-rr')^2}\\
&=\frac{1-2rr'+r^2r'^2-r^2+2rr'-r'^2}
{4rr'\cos^2\{kb+(\theta+\theta')/2\}+(1-rr')^2}\\
&=\frac{(1-r^2)(1-r'^2)}{(1-rr')^2}\frac{1}
{\left\{\frac{1-rr'}{2\sqrt{rr'}}\right\}^{-2}\cos^2\{kb+(\theta+\theta')/2\}+1}\\
);

この関数は &math(kb+(\theta+\theta')/2=\pi(n+1/2)); の時に最大値

　&math(\frac{(1-r^2)(1-r'^2)}{(1-rr')^2});

を取る。このとき、

　&math(
|\varphi(a)|^2=1+\frac{2r}{1-r^2}\sin\{\theta\}
);

　&math(
|\varphi(a+b)|^2
&=1+\frac{2r}{1-r^2}\sin\theta'\\
);

であり、ピーク幅は、

　&math(
\frac{1-rr'}{b\sqrt{rr'}}
);

となる。両側の反射率の相乗平均を取って &math(r_\mathrm{mean}=\sqrt{rr'}); とすると、

　&math(
\Delta k=\frac{1-r_\mathrm{mean}^2}{br_\mathrm{mean}}
);

となり、対称な場合の結果と一致する。

対称２重トンネルの結果が、

　&math(
\Delta k=\frac{1}{b}(壁１枚の透過率)^2
);

だったのに、

　&math(
\Delta k\ne\frac{1}{b}(壁１の透過率)(壁２の透過率)
);

なのは興味深い。

むしろ、対称な時を

　&math(
\Delta k=\frac{1}{b}\Big[1-(壁１枚の反射率)^2\Big]
);

と書き、非対称なときに

　&math(
\Delta k=\frac{1}{b}\Big[1-(壁１の反射率)(壁２の反射率)\Big]
);

と書けることが本質であるということか？

* 左右対称３重トンネル [#d6769b17]

　&math(
&\begin{pmatrix}Ge^{ik(2a'+2b+a)}\\He^{-ik(2a'+2b+a)}\end{pmatrix}\\
&=
\frac{1}{\sqrt{1-r^2}(1-r'^2)}
\begin{pmatrix}
e^{i\theta'}&-ir'\\
ir'&e^{-i\theta'}
\end{pmatrix}
\begin{pmatrix}e^{ikb}&\\&e^{-ikb}\end{pmatrix}
\begin{pmatrix}
e^{i\theta}&-ir\\
ir&e^{-i\theta}
\end{pmatrix}
\begin{pmatrix}e^{ikb}&\\&e^{-ikb}\end{pmatrix}
\begin{pmatrix}
e^{i\theta'}&-ir'\\
ir'&e^{-i\theta'}
\end{pmatrix}
\begin{pmatrix}A\\B\end{pmatrix}\\
&=
\frac{1}{\sqrt{1-r^2}(1-r'^2)}
\begin{pmatrix}
e^{i\theta'}&-ir'\\
ir'&e^{-i\theta'}
\end{pmatrix}
\begin{pmatrix}e^{ikb}&\\&e^{-ikb}\end{pmatrix}
\begin{pmatrix}
e^{ikb+i(\theta+\theta')}+rr'e^{-ikb}&
-ir'e^{ikb+i\theta}-ire^{-ikb-i\theta'}\\
ire^{ikb+i\theta'}+ir'e^{-ikb-i\theta}&
rr'e^{ikb}+e^{-ikb-i(\theta+\theta')}
\end{pmatrix}
\begin{pmatrix}A\\B\end{pmatrix}\\
&=
\frac{1}{\sqrt{1-r^2}(1-r'^2)}
\begin{pmatrix}
e^{i\theta'}&-ir'\\
ir'&e^{-i\theta'}
\end{pmatrix}
\begin{pmatrix}
e^{2ikb+i(\theta+\theta')}+rr'&
-ir'e^{2ikb+i\theta}-ire^{-i\theta'}\\
ire^{i\theta'}+ir'e^{-2ikb-i\theta}&
rr'+e^{-2ikb-i(\theta+\theta')}
\end{pmatrix}
\begin{pmatrix}A\\B\end{pmatrix}\\
&=
\frac{1}{\sqrt{\dots}(\dots)}
\begin{pmatrix}
e^{2ikb+i(\theta+2\theta')}+r'^2e^{-2ikb-i\theta}+2rr'e^{i\theta'}&
-ir'e^{2ikb+i(\theta+\theta')}-ir(1+r'^2)
-ir'e^{-2ikb-i(\theta+\theta')}\\
ir'e^{2ikb+i(\theta+\theta')}+ir(1+r'^2)+ir'e^{-2ikb-i(\theta+\theta')}&
e^{-2ikb-i(\theta+2\theta')}+r'^2e^{2ikb+i\theta}+2rr'e^{-i\theta'}
\end{pmatrix}
\begin{pmatrix}A\\B\end{pmatrix}\\
);

&math(A=1,H=0); と置いて、

　&math(
|B|^2
&=\left|\frac{ir'e^{2ikb+i(\theta+\theta')}+ir(1+r'^2)+ir'e^{-2ikb-i(\theta+\theta')}}
{e^{-2ikb-i(\theta+2\theta')}+r'^2e^{2ikb+i\theta}+2rr'e^{-i\theta'}}\right|^2\\
&=\frac{\big\{2r'\cos(2kb+\theta+\theta')+r(1+r'^2)\big\}^2}
{\left|e^{-2ikb-i(\theta+\theta')}+r'^2e^{2ikb+i(\theta+\theta')}+2rr'\right|^2}\\
&=\frac{\big\{2r'\cos(2kb+\theta+\theta')+r(1+r'^2)\big\}^2}
{\left|(1+r'^2)\cos(2kb+\theta+\theta')+2rr'+i(1-r'^2)\sin(kb+\theta+\theta')\right|^2}\\
&=\frac{\big\{2r'\cos(2kb+\theta+\theta')+r(1+r'^2)\big\}^2}
{\big\{(1+r'^2)\cos(2kb+\theta+\theta')+2rr'\big\}^2+(1-r'^2)^2\sin^2(kb+\theta+\theta')}\\
&=\frac{\big\{2r'\cos(2kb+\theta+\theta')+r(1+r'^2)\big\}^2}
{\big\{(1+r'^2)\cos(2kb+\theta+\theta')+2rr'\big\}^2+(1-r'^2)^2-(1-r'^2)^2\cos^2(kb+\theta+\theta')}\\
);

　&math(
(分母)
=&(1+r'^2)^2\cos^2(2kb+\theta+\theta')+4rr'(1+r'^2)\cos(2kb+\theta+\theta')+4r^2r'^2+\\
&\ \ \ \ (1-r'^2)^2-(1-r'^2)^2\cos^2(kb+\theta+\theta')\\
=&4r'^2\cos^2(2kb+\theta+\theta')+4rr'(1+r'^2)\cos(2kb+\theta+\theta')+4r^2r'^2+(1-r'^2)^2\\
=&\big\{2r'^2\cos^2(2kb+\theta+\theta')+r(1+r'^2)\big\}^2-r^2(1+r'^2)^2+4r^2r'^2+(1-r'^2)^2\\
=&\big\{2r'^2\cos^2(2kb+\theta+\theta')+r(1+r'^2)\big\}^2-r^2(1-r'^2)^2+(1-r'^2)^2\\
=&\big\{2r'^2\cos^2(2kb+\theta+\theta')+r(1+r'^2)\big\}^2+(1-r^2)(1-r'^2)^2\\
);

　&math(
|E'|^2=1-|B|^2
&=\frac{(1-r^2)(1-r'^2)^2}
{\big\{2r'\cos(2kb+\theta+\theta')+r(1+r'^2)\big\}^2+(1-r^2)(1-r'^2)^2}\\
&=\frac{1}
{\frac{4r'^2}{(1-r^2)(1-r'^2)^2}\big\{\cos(2kb+\theta+\theta')+\frac{r(1+r'^2)}{2r'}\big\}^2+1}\\
);

透過確率が１になるのは

　&math(
\cos(2kb+\theta+\theta')=-\frac{r(1+r'^2)}{2r'}
);

&math(\frac{r(1+r'^2)}{2r'}<1); は &math(r<\frac{2r'}{1+r'^2}); だから、
両側に比べてそこそこ中央が通りやすければ解がある。

その周辺に於いて、

　&math(
|E'|^2
&=\frac{1}
{\frac{4r'^2}{(1-r^2)(1-r'^2)^2}4b^2\left[1-\left(\frac{r(1+r'^2)}{2r'}\right)^2\right](k-k_n)^2+1}\\
&=\frac{1}
{4b^2\frac{4r'^2-r^2(1+r'^2)^2}{(1-r^2)(1-r'^2)^2}(k-k_n)^2+1}\\
&=\frac{1}
{4b^2\frac{(2r'-r-rr'^2)(2r'+r+rr'^2)}{(1-r^2)(1-r'^2)^2}(k-k_n)^2+1}\\
);

だから、ピーク幅は

　&math(
\Delta k=\frac{\sqrt{(1-r^2)}(1-r'^2)}{b\sqrt{(2r'-r-rr'^2)(2r'+r+rr'^2)}}
);

となる。

PukiWiki 1.5.0 Copyright © 2001-2006 PukiWiki Developers Team. License is GPL.
Based on "PukiWiki" 1.3 by yu-ji. Powered by PHP 7.4.33. HTML convert time: 0.002 sec.
Also based on Bootstrap, Start Bootstrap - SB Admin 2, highlight.js, Google fonts, Google recaptcha, Google chart, Bitly.