励起子シミュレーションチュートリアル
1. はじめに
このチュートリアルでは、有機太陽電池（OPV）デバイスにおける励起子ダイナミクスをどのようにシミュレーションするかを調べます。 有機半導体で光子が吸収されると、励起子 と呼ばれる束縛電子–正孔対が生成されます。 OPVデバイスでは、これらの励起子はドナー–アクセプタ界面までドリフトまたは拡散し、 そこで自由電荷キャリアへ解離しなければなりません。 この過程は、励起子輸送と最終的な電子および正孔への変換の両方を記述する 励起子解離方程式 を用いてモデル化できます。 OghmaNanoでは、これらの励起子ダイナミクスは標準的なdrift–diffusion方程式の上に導入され、 完全な光物理像を捉えます。励起子モデリングが必要となる場合 （および省略できる場合）についてのより詳しい説明は、 励起子とgeminate再結合 に示されています。
2. 開始方法
ファイルリボンの New simulation タブから、 New simulation ウィンドウを開くためにクリックしてください （Figure 1a を参照）。 このウィンドウには、多くの事前構成されたサンプルカテゴリがあります。 Exciton simulations フォルダをダブルクリックして、 利用可能な励起子関連サンプルのサブメニューを開いてください （Figure 1b を参照）。 このチュートリアルでは、Exciton device テンプレートを使用します。 これは、OPV構造における励起子生成、輸送、解離をシミュレーションするための シンプルな出発点を提供します。
3. シミュレーションを実行して励起子出力を確認する
デバイスを開くと、メインウィンドウ内にスタックが表示されます。 シミュレーションを実行する前に、Electrical タブへ移動し、 Exciton solver ボタンが押し込まれた状態になっていることを確認してください — 通常これはデフォルトで有効ですが、確認する価値があります （??）。 これにより、モデル内で励起子ダイナミクスが有効であることが保証されます。 これが完了したら、File メニューへ戻り、Run simulation ボタンを押して モデルを実行してください（??）。
シミュレーションが終了したら、
Output タブで結果を確認できます
（??）。
ここには、励起子関連出力を含む、ディスクに書き込まれたすべてのファイルが一覧表示されます。
jv.csv をダブルクリックすると、JV曲線をプロットできます
（??）。
このプロットは通常のデバイスシミュレーションから得られる標準的なJV曲線のように見えますが、
この場合は Exciton solver が有効であるため、
基礎となる物理には励起子ダイナミクスが含まれています。
jv.csv からのJV曲線。
標準的なJV曲線のように見えますが、Exciton solver が有効になっています。
出力ディレクトリ内の exciton_output フォルダをダブルクリックすると、励起子ソルバーの詳細な結果にアクセスできます（?? を参照）。このディレクトリには、位置の関数としての定数や、励起子ダイナミクスから導出された計算量など、すべてのソルバー出力が含まれています。たとえば、exciton.csv をダブルクリックすると、デバイス膜厚方向にわたる励起子分布のプロットが得られます（?? を参照）。同様に、Gn.csv は位置の関数として電子生成率を与え、Gp.csv は正孔生成率を示します（?? を参照）。単純な1Dデバイスの場合、これら2つの生成率ファイルは実質的に同一です。
exciton_output フォルダの内容。
4. 励起子輸送方程式とパラメータ
デバイス内の励起子分布は、次の励起子輸送方程式によって支配されます:
\[ \frac{\partial X}{\partial t} = \nabla \!\cdot \!\big(D\,\nabla X\big) + G_{\mathrm{optical}} - k_{\mathrm{dis}}\,X - k_{\mathrm{FRET}}\,X - k_{\mathrm{PL}}\,X - \alpha\,X^{2} \]
ここで \(X(\mathbf{r},t)\) は励起子密度（m\(^{-3}\)）、 \(D\) は励起子拡散係数（m\(^2\)s\(^{-1}\)）、 \(G_{\mathrm{optical}}\) は光学モデルによって与えられる局所励起子生成率 （吸収光子に比例）です。 \(k_{\mathrm{dis}}\) は自由電荷への解離速度、 \(k_{\mathrm{FRET}}\) はFörster共鳴エネルギー移動速度、 \(k_{\mathrm{PL}}\) は放射減衰速度、 \(\alpha\) は励起子–励起子消滅係数（m\(^3\)s\(^{-1}\)）です。 このモデルが有効なとき、電子drift–diffusion方程式は \(G = k_{\mathrm{dis}}\,X\) と定義された生成項を使用し、 設定に応じて界面領域に制限される場合があります。
これらすべてのパラメータは、メインメニューからアクセスできる Electrical parameter editor で表示および変更できます （?? を参照）。 励起子固有のフィールドは Excitons 見出しの下にまとめられています。以下の表に要約を示します。
| パラメータ | 意味 | 単位 |
|---|---|---|
| Scattering length | 散乱前の励起子の有効拡散長。 | m |
| Lifetime | 励起子が減衰または解離する前に存続する平均時間。 | s |
| kPL | 放射減衰速度（フォトルミネセンス）。 | s⁻¹ |
| kFRET | Förster共鳴エネルギー移動速度。 | s⁻¹ |
| kα | 励起子–励起子消滅係数。 | m³ s⁻¹ |
| kdis | 励起子を自由電荷へ変換する解離速度定数。 | s⁻¹ |
5. ソルバーがシミュレーションプロセスにどのように組み込まれるか
励起子ソルバーがオフの場合（?? の上部）、transfer-matrix光学は キャリア生成率 を計算し、それをdrift–diffusionソルバーへ直接渡します。励起子ソルバーがオンの場合（?? の下部）、光学モデルは代わりに 励起子生成率 を励起子ソルバーへ供給します。励起子ソルバーはその後、Electrical parameter editor で設定されたパラメータを用いて、その集団を拡散、移動（FRET）、解離、放射減衰、および消滅させながら時間発展させ、drift–diffusionソルバー用の最終的な キャリア生成率 を出力します。要するに、励起子ソルバーは光学と電気輸送の間に挿入され、光学設定やdrift–diffusion方程式を変更することなく励起子物理をモデル化できます。
励起子ソルバーを有効にしたシミュレーションにおける処理順序は、 ?? に示されています。 まず、光学ソルバーが実行され、デバイス全体をスライスごとに横切って光子吸収プロファイルを計算します。 次に、励起子ソルバーが実行され、生成された励起子が伝播、移動、または解離し、 ソルバーが収束するまで処理されます。通常は数十ステップ以内です。 最後に、励起子生成率がデバイス深さと波長の両方の関数として調べられます。 これはtransfer matrixモデルによって計算されます。 これらの出力をまとめると、光学モデルから励起子ソルバー、そして最終的にdrift–diffusion方程式へとつながるパイプラインが構成されます。
💡 タスク: これらの簡単な編集を試して励起子モデルを調べてください（ヒント: 効果を見るには1〜2桁変更してください。）:
- Electrical parameter editor で exciton lifetime を変更し、シミュレーションを再実行する。
- dissociation rate (kdis) を下げ、キャリア生成への影響を観察する。
- より強い励起子–励起子相互作用をシミュレートするために、exciton–exciton annihilation coefficient (α) を増加させる。
- 励起子ソルバーを enabled にした場合と disabled にした場合のJV曲線を比較する。
✅ 期待される結果
- 励起子寿命が長いほど励起子密度が高くなり、キャリア生成が増加します。
- 解離速度が低いと、自由キャリアへ変換される励起子が少なくなるため、光電流は低下します。
- 高い消滅速度は励起子密度プロファイルを広げる、または抑制し、効率を低下させます。
- 励起子ソルバーを有効にしたJV曲線では、特にトラッピングや消滅が顕著な場合、標準的なdrift–diffusionモデルと比較してJSC が低下することがあります。