モット・ガーニー方程式の導出

空間電荷制限電流（SCLC）は、半導体デバイス物理学において最も重要な輸送機構の一つである。これは、電極から注入されるキャリアの密度が、半導体内部に既に存在する平衡状態のキャリア密度をはるかに上回る場合に生じる。このような条件下では、電流は材料本来の導電性によって制限されるのではなく、注入された電荷自体が生成する電界によって制限されるのである。

SCLC測定法は、有機半導体やOLED、ペロブスカイトデバイス、および無秩序な構造を持つ材料において、キャリア移動度の算出やトラップ準位の解析に広く用いられている。実験的には、SCLCの測定には通常ホール専用または電子専用の構造が用いられ、この場合両電極からは同一種類のキャリアのみが優先的に注入される。これにより電子・ホールの再結合が抑制され、単一種類のキャリアのみによる輸送が支配的となる。

トラップが存在しない場合の空間電荷制限電流を記述する最もよく知られた解析式がモット・ガーニー方程式であり、この式によれば電流密度は印加電圧の二乗に比例することになる：

$$ J = \frac{9}{8}\varepsilon\mu\frac{V^2}{L^3} $$

ここで\(J\)は電流密度、\(\varepsilon\)は材料の誘電率、\(\mu\)はキャリア移動度、\(V\)は印加電圧、そして\(L\)はデバイスの厚みを表します。

この式自体は簡潔であるものの、その背後にある物理的仮定が極めて重要である。導出過程では、トラップ準位が存在しない輸送状態、単一種類のキャリアのみの存在、オーム性注入条件、均一なキャリア移動度、拡散電流の無視、および再結合現象の無視が仮定されている。実際のSCLC測定結果を解釈する際には、こうした仮定を正しく理解しておくことが不可欠である。

空間電荷制限電流の物理的起源

通常のオーム性伝導においては、半導体内のキャリア密度はほぼ一定であり、電流密度は電圧に比例して増加する。この状態ではデバイス内部の電界はほぼ一様であり、注入されるキャリア密度も比較的少ない。

しかし、オーム性注入が効率的に起こる場合、低導電性半導体内部に注入された電荷が次第に蓄積される。この蓄積電荷はポアソン方程式に従って内部電界を変化させる。その結果、電界は位置に依存するようになり、電流は注入された電荷自体の静電相互作用によって制限されることになる。

この現象は??に示されている。オーム性伝導とは異なり、SCLC領域における内部電界は均一ではない。実際には、注入された電荷が静電ポテンシャルを継続的に変化させるため、電界強度はデバイス内の距離に応じて増大していく。

出発点となる方程式

モット・ガーニー方程式の導出は、ポアソン方程式およびドリフト電流方程式という2つの基本方程式から始まる。

ポアソンの方程式は、電界の勾配と電荷密度との関係を示す式である：

$$ \frac{dE}{dx} = \frac{\rho}{\varepsilon} $$

ここで \(E\) は電界強度、\(\rho\) は電荷密度、\(\varepsilon\) は半導体の誘電率を表す。

単一種類のキャリアが存在しその密度が\(n\)であると仮定すると、電荷密度は次のように表される：

$$ \rho = qn $$

ここで\(q\)は素電荷を表す。

ドリフト電流密度は次式で与えられる：

$$ J = q\mu nE $$

ここで\(\mu\)はキャリア移動度を表す。

モット・ガーニー式の導出においては拡散電流は無視される。この近似は、電界強度が十分に強くドリフト輸送が拡散輸送よりも支配的となる場合に成り立つ。

完全なドリフト拡散方程式は次の通りである：

$$ J = q\mu nE + qD\frac{dn}{dx} $$

しかし空間電荷制限電流条件下では、拡散項はドリフト項に比べて小さいと仮定される：

$$ qD\frac{dn}{dx} \ll q\mu nE $$

モット・ガーニー方程式の導出

まずポアソンの方程式とドリフト電流方程式を組み合わせる。次に代入を行う。

$$ n = \frac{\varepsilon}{q}\frac{dE}{dx} $$

ドリフト電流の式に代入すると：

$$ J = q\mu nE $$

結果として：

$$ J = q\mu \left( \frac{\varepsilon}{q}\frac{dE}{dx} \right) E $$

これを簡略化すると：

$$ J = \varepsilon\mu E\frac{dE}{dx} $$

式の整理：

$$ E\,dE = \frac{J}{\varepsilon\mu}\,dx $$

デバイス全体にわたって積分すると：

$$ \int_0^{E} E\,dE = \frac{J}{\varepsilon\mu} \int_0^x dx $$

結果として：

$$ \frac{E^2}{2} = \frac{Jx}{\varepsilon\mu} $$

したがって：

$$ E(x) = \sqrt{ \frac{2Jx}{\varepsilon\mu} } $$

デバイス両端の電圧は、電界を積分することで求められる：

$$ V = \int_0^L E(x)\,dx $$

電界の式を代入すると：

$$ V = \int_0^L \sqrt{ \frac{2Jx}{\varepsilon\mu} } \,dx $$

定数を積分記号の外に出す：

$$ V = \sqrt{ \frac{2J}{\varepsilon\mu} } \int_0^L x^{1/2}\,dx $$

積分を実行すると：

$$ V = \sqrt{ \frac{2J}{\varepsilon\mu} } \cdot \frac{2}{3}L^{3/2} $$

両辺を二乗すると：

$$ V^2 = \frac{8JL^3}{9\varepsilon\mu} $$

最後に、電流密度について整理すると、モット・ガーニー方程式が得られる。

$$ J = \frac{9}{8} \varepsilon\mu \frac{V^2}{L^3} $$

この導出の背後にある物理的仮定

この導出過程は数学的には単純に見えるものの、いくつかの重要な物理的仮定に基づいている。

トラップ準位が存在しない輸送状態

この導出過程では、注入されたすべてのキャリアが移動性を維持し、電気輸送に直接寄与すると仮定されている。しかし実際の無秩序な半導体においては、キャリアはしばしば局在化したトラップ準位に捕捉される。このようなトラップ現象によりキャリア分布が変化し、電流の電圧依存性にも影響が生じる。

単一キャリア輸送

モット・ガーニー方程式では、輸送過程において単一種類のキャリアのみが主要な役割を果たすと仮定されている。電子と正孔の両方を含むバイポーラ型デバイスにおいては、再結合現象や電荷中性条件の影響により、単純な空間電荷制限電流モデルの導出が成立しなくなる。

オーム性接触

この導出では、キャリアが半導体中へ自由に注入可能であると仮定している。もし接触面にショットキーバリアまたは注入障壁が存在する場合、電流はデバイス内部の空間電荷ではなく接触部からの注入現象によって制限されることになる。

均一なキャリア移動度

この導出では、キャリア移動度が電界強度やキャリア密度に依存しないことを前提としている。だが多くの有機半導体においては、移動度は電界強度や物質の無秩序性、さらにはトラップ準位の占有状態に強く依存する。

拡散電流は無視できる程度である

この導出過程においては拡散輸送は無視されている。この近似は電界が弱い場合や平衡状態に近い条件下では精度が低下する。

モット・ガーニー方程式が適用できなくなる場合

実際の半導体デバイスでは、理想的なモット・ガーニー挙動からの逸脱が頻繁に見られる。その最も一般的な原因の一つがトラップ準位の存在である。トラップが支配的な条件下における空間電荷制限電流では、電流はしばしばべき乗則に従う。

$$ J \propto V^m $$

ここで、\(m > 2\)となる。

電界依存性を持つキャリア移動度もまた、電流と電圧の関係に大きな影響を与える。多くの無秩序な有機半導体ではポール・フレンケル効果が見られ、この現象において移動度は電界強度に対して指数関数的に増加する。

極めて高い電圧下では、直列抵抗や自己発熱の影響も無視できなくなる。同様に、電子と正孔が同時に存在する場合、再結合現象によって電界分布が変化し、単一キャリアという仮定が成り立たなくなる可能性がある。

これらの理由から、実際の半導体デバイスを正確に記述するためには数値的なドリフト拡散シミュレーションがしばしば必要となる。

ドリフト拡散シミュレーションとの関連性

モット・ガーニー方程式は、空間電荷制限輸送現象に関する解析的な知見を与えるため極めて有用である。だが実際の半導体デバイスには、トラップ準位の分布や電界依存性を持つキャリア移動度、非オーム性接触、再結合過程、さらには空間的に変化する電界がしばしば存在し、これらは単純な解析モデルでは記述できない。

OghmaNanoをはじめとするドリフト拡散シミュレーションツールは、ポアソン方程式および連続の式を数値的に解くことで、上記のような非理想的な効果を明示的に取り入れることができる。これにより、不規則性、トラップ準位、移動度の勾配、さらには複雑な多層構造を有する実際のSCLCデバイスを解析することが可能となる。