ホーム 例 スクリーンショット ユーザーマニュアル Blueskyロゴ YouTube ☰
OghmaNano 有機/ペロブスカイト太陽電池、OFET、OLEDをシミュレーション ダウンロード Quick Start guide

ランジュバン再結合

Illustration of diffusion-limited recombination showing electrons moving through a semiconductor towards a hole.
拡散律速型再結合とは、半導体内において電子と正孔が互いに出会う確率を指す。

ランジュバン再結合とは、有機半導体や非晶質シリコン、その他多くの薄膜材料といった無秩序な構造を持つ半導体における電荷キャリアの再結合を記述するために用いられる二分子反応型の再結合メカニズムである。

ランジュバンモデルにおいては、熱運動および静電気的引力の作用により電子と正孔が半導体中を移動することが再結合の原因となる。電子と正孔が十分に近づくと、両者の間に働くクーロン力が強くなり、急速に再結合が生じる。この過程は??に示されている。

ランジュバン再結合の本質的な特徴は、再結合速度が再結合現象そのものではなく、キャリアが材料中を移動して互いに出会う速度によって制限される点にある。このため、この過程は拡散律速型あるいは出会い律速型の再結合過程と呼ばれることも多い。

キャリアの運動はキャリア移動度によって制御されるため、ランジュバン再結合モデルによれば再結合率は電子および正孔の移動度に比例して変化するはずである。したがって、キャリア移動度が低い材料ではランジュバン再結合率も遅くなる傾向がある。

ランジュバン再結合の物理的描像

電子の海に囲まれた正孔を考えてみよう(??)。電子は熱エネルギーおよび局所的な電界の影響により半導体内を移動する。もし電子が正孔に十分近づくと、両キャリア間のクーロン引力が、それらを引き離そうとする熱エネルギーを上回ることになる。

ランジュバンは、クーロン相互作用エネルギーが熱エネルギーと等しくなる特徴的な捕捉半径\(r_c\)(??内の赤色円)を定義した。この半径は一般にランジュバン再結合半径と呼ばれる。

\[ E_{\mathrm{Coulomb}} = E_{\mathrm{thermal}} \]

クーロンポテンシャルエネルギーと熱エネルギーを等置すると次式が得られる:

\[ \frac{q^2}{4\pi \varepsilon r_c} = k_B T \]

式を整理すると、再結合半径が得られる:

\[ r_c = \frac{q^2}{4\pi \varepsilon k_B T} \]

この半径内では、クーロン力による引力が熱運動を上回るため、再結合が極めて起こりやすくなる。

誘電率が低下すると再結合半径は増大する。これは、有機半導体のような低誘電率材料においてはクーロン引力がより強く働くことを意味する。このことが、無秩序な構造を持つ有機電子材料においてランジュバン再結合が特に重要となる理由の一つである。

ランジュバン再結合率の導出

Carrier flux across the Langevin recombination sphere.
ホールを取り囲む再結合球の表面を電子がドリフトする。

ランジュバン再結合は、正孔を取り囲む再結合球内への電子フラックスを計算することで導出できる(??)。

電子のドリフト電流密度は次の通りである:

\[ j_n = q n \mu_n E \]

ここで\(n\)は電子密度、\(\mu_n\)は電子移動度、そして\(E\)は正孔によって生じる電界を表す。

正孔の周囲に生じるクーロン電場は次式で表される:

\[ E(r)=\frac{q}{4\pi \varepsilon r^2} \]

この電場をドリフト電流の式に代入すると、次式が得られる:

\[ j_n = q n \mu_n \frac{q}{4\pi \varepsilon r^2} \]

半径\(r\)の球面上を通過する電子流束の総量は、球の表面積\(4\pi r^2\)を乗じることで求められる。

\[ I_n = j_n 4\pi r^2 \]

したがって:

\[ I_n = \frac{q^2 \mu_n n}{\varepsilon} \]

正孔に対しても同様の式が存在する:

\[ I_p = \frac{q^2 \mu_p p}{\varepsilon} \]

全再結合率は、電子と正孔が互いに出会う確率に比例しており、これによりランジュバン再結合の式が得られる。

\[ R_{\mathrm{L}} = k_{\mathrm{L}} np \]

ここで、ランジュバン再結合の前因子は次の通りである:

\[ k_{\mathrm{L}} = \frac{q}{\varepsilon}\left(\mu_n+\mu_p\right) \]

多くの有機半導体においては、一方のキャリアの移動度が支配的であるため、再結合率は概ね移動度の遅い方のキャリアによって制御される。

無秩序な半導体におけるランジュバン再結合

Illustration of recombination in disordered semiconductors.
無秩序な構造を持つ半導体においては、キャリアの輸送は局在化した状態間のホッピング運動によって行われるため、再結合現象はキャリア移動度に強く依存する。

ランジュバン再結合は、再結合現象とキャリア移動度とを自然に結び付けるため、無秩序な半導体におけるドリフト拡散シミュレーションで広く用いられている。

有機半導体のような材料においては、電荷輸送は結晶シリコンの場合のような広域なブロッホ状態を介した輸送ではなく、局在化した状態間での熱活性化に基づくホッピングによって行われるのが一般的である。そのため、キャリア移動度はエネルギー的無秩序度、温度、およびキャリア密度に強く依存することになる。

ランジュバン再結合はキャリア移動度に直接依存するため、再結合率もまた移動度の値に強く左右される。このような挙動は、有機太陽電池や有機LED材料において実験的に頻繁に観測される。

ドリフト拡散デバイスシミュレーションにおいては、ランジュバン再結合率は通常、二分子反応型の再結合項として直接組み込まれる。

\[ R_{\mathrm{L}} = \frac{q}{\varepsilon} \left( \mu_n+\mu_p \right) \left( np-n_i^2 \right) \]

再結合率を次の形式で表す \(np-n_i^2\) これにより、熱平衡状態において純粋な再結合率がゼロになることが保証される。

ランジュバンモデルの仮定事項

標準的なランジュバン再結合モデルには、いくつかの重要な仮定が含まれている:

これらの仮定は無秩序な半導体においては概ね成り立つが、相分離が顕著である系や、キャリアが強くトラップされる場合、キャリアの局在化が生じる場合、あるいは界面における再結合障壁が存在する場合には成立しなくなる可能性がある。

低減型ランジュバン再結合

Experimental evidence showing reduced Langevin recombination.
実験測定の結果によれば、実際の再結合率は古典的なランジュバン理論から予測される値よりも低くなる場合が多い。

ランジュバン理論は無秩序な半導体における再結合現象の多くを説明することに成功しているが、実験で観測される再結合率は古典的なランジュバン方程式から予測される値よりもしばしば大幅に低くなっている。

この現象は一般に「低減型ランジュバン再結合」と呼ばれる。実際には、再結合の前係数は次のように表記されることが多い:

\[ k = \gamma k_{\mathrm{L}} \]

ここで \(\gamma\) は通常1より小さい値をとる減少係数である。

いくつかの物理的メカニズムが、ランジュバン再結合の抑制を引き起こす可能性がある:

有機太陽電池デバイスにおいては、実験で得られた電流電圧特性曲線や過渡光電圧測定結果、さらには電荷抽出実験の結果を再現するため、理論値よりも低いランジュバン再結合率がしばしば必要とされる。

次に何をすべきか

実際のデバイスシミュレーションにおいてランジュバン再結合がどのように現れるかを確認するには、有機太陽電池の再結合チュートリアルへ進んでください。そこではドリフト拡散シミュレーションにおいてランジュバン再結合とショックレー・リード・ホール再結合が直接比較されています。

トラップを介した再結合メカニズムについては、以下をご参照ください ショックレー・リード・ホール再結合(SRH再結合)無秩序な材料中における局在状態の役割については、こちらをご覧ください。 なぜ無秩序な半導体モデルにおいてトラップ状態が必要なのか.

ドリフト拡散方程式の導出法およびランジュバン再結合との組み合わせについては ドリフト拡散理論 を参照されたい。