渡邉孝信（早稲田大学・電子物理システム学科）

ねぇおかしいでしょ1/2

半導体デバイスの授業では、キャリアの基本的な輸送機構として、

ドリフト電流　と　拡散電流

の概念をまず学びます。まだ非平衡統計力学を学んでいない場合がほとんどでしょうから、ドリフトや拡散の概念を説明する際は、高校物理で登場する気体分子運動論が用いられます。熱運動する気体分子に電子をなぞらえ、平均速度 $\bar{v}$ 、平均自由行程 $\bar{l}$ 、平均自由時間$\bar{t}$ を使って、ドリフト移動度 $\mu$ と拡散係数 $D$ を表してみせるのです。

　その際に困るのは、この初歩的に導出される拡散係数 $D$ の式が本によってまちまちで、しかもイマイチ釈然としない説明が多いことです。一方、ドリフト移動度 $\mu$ の式はほぼ１通りに落ち着いています。そこで、金科玉条のアインシュタインの関係式

$$\begin{array}{}\text{(1.1)}& D=\frac{\mu }{q}{k}_{B}T\end{array}$$

を使って（ $q$は素電荷、 $k_B$はボルツマン定数、 $T$は絶対温度）、ドリフト移動度 $\mu$ から拡散係数$D$ を決定することで、そうした混乱を回避している教科書も多数あります。というよりむしろ、拡散係数$D$ は式$\text{(1.1)}$で定めるべきでしょう。

　そうは言っても、拡散というキャリア輸送の微視的描像をつかむには、気体分子運動論に基づく説明に触れておくこともたいへん重要です。

　表１に、筆者が知っている範囲ですが、代表的な教科書や専門書に書かれている拡散係数の式をまとめておきました。どうしてこうもいろいろあるのでしょうか？

それは、平均速度 $\bar{v}$ 、平均自由行程 $\bar{l}$、平均自由時間 $\bar{t}$といった量の「平均」の取り方に、いろいろな方法や考え方があるからです。

表１　様々な拡散係数の式。 $\bar{v}$ :平均速度（熱運動速度）、 $\bar{l}$：平均自由行程、$\bar{t}$：平均自由時間 

まず注意すべきは、考えている系の次元が必ずしも同じではないということです。３次元空間における拡散流を考える際、その流れの方向に対して斜影をとるため、$\bar{v}$と$\bar{l}$はそれぞれ$1/\sqrt{3}$倍されます。３次元モデルの多くに1/3や1/6という係数がついているのはそのためです。終始１次元で議論している本では、1/3倍しないままにしています。

また、平均熱速度 には「速さの平均」や「根二乗平均速度」など、複数の定義があることも念頭に置いておく必要があります。ただしこの差はそれほど大きくありません。マクスウェル－ボルツマン分布を仮定した場合

$$\begin{array}{}\text{(1.2)}& \mathrm{平}\mathrm{均}\mathrm{速}\mathrm{さ}：\overline{|\mathit{v}|}=\sqrt{\frac{8k_{B}T}{\pi m}}\end{array}$$

$$\begin{array}{}\text{(1.3)}& \mathrm{根}\mathrm{二}\mathrm{乗}\mathrm{平}\mathrm{均}\mathrm{速}\mathrm{度}：\sqrt{\overline{{\mathit{v}}^2}}=\sqrt{\frac{3k_{B}T}{m}}\end{array}$$

となり、根二乗平均速度は速さの平均の0.92倍と、若干ですが小さめになります。

表1の文献[9]のMeyerの本は、気体分子運動論に基づいて拡散係数を定式化した最初期の仕事です。Meyerの式では、式$\text{(1.2)}$の速さの平均$\overline{|\mathit{v}|}$が$\bar{v}$として採用されていて、平均自由行程も$\bar{l}=\overline{|\mathit{v}|}\bar{t}$とおいています。これらを根二乗平均速度$\sqrt{\overline{{\mathit{v}}^2}}$に置き換えると、$\overline{|\mathit{v}|}=\sqrt{8/(3\pi )}\sqrt{\overline{{\mathit{v}}^2}}$、$\bar{l}=\sqrt{8/(3\pi )}\sqrt{\overline{{\mathit{v}}^2}}\bar{t}$となるので、

$$D=\frac{\pi }{8}\bar{v}\bar{l}=\frac{\pi }{8}{\overline{|\mathit{v}|}}^2\bar{t}=\frac{\pi }{8}\frac{8}{3\pi }\overline{{\mathit{v}}^2}\bar{t}=\frac{1}{3}\overline{{\mathit{v}}^2}\bar{t}$$

となり、$D=\bar{v}\bar{l}/3$としている他の文献と一致します。

　もう一つ、拡散係数$D$の違いの要因となっているのは、平均自由時間$\bar{t}$と平均自由行程$\bar{l}$の「平均」の取り方です。1/2倍の係数がかかったりかからなかったりするのは、この平均操作の仕方の違いに起因します（注：パウリ[3]とChapman[11]はまた別の理由で1/2の係数をつけています）。気体分子運動論では、気体分子同士の衝突イベントと次の衝突イベントの間の時間の平均が$\bar{t}$、衝突イベントから次の衝突イベントが起きるまでに進める距離の平均を$\bar{l}$、としています。これらの統計平均を計算する際、どのような確率分布を仮定するかによって、平均値（期待値）が変わってしまうのです。このシリーズ記事では、この確率分布の取り扱い方に焦点を当てて、拡散係数の式の違いを議論していきます。