励起子ドメインチュートリアル（パート A）：3D 励起子ドメインシミュレーションを実行する

1. はじめに

励起子がどのように生成され、輸送され、解離するかを理解することは、 バルクヘテロ接合（BHJ）有機太陽電池の動作において中心的です。これらの系では、 光は主としてドナー相で吸収され、励起子として知られる束縛された電子–正孔対を生成します。 これらの励起子はドナー材料中を拡散し、ドナー–アクセプタ界面に到達しなければなりません。 界面では、それらは自由電子と正孔に解離する可能性があり、その確率は局所的な形態、 次元性、および関連する速度定数によって決まります。

したがって、励起子拡散、界面解離、および競合する損失過程のバランスは、 有効な光生成効率を決定するうえで決定的な役割を果たします。 しかし、明示的な励起子モデリングは、概念的困難さのためではなく、 パラメータの不確かさのために、歴史的には避けられてきました。励起子拡散長、寿命、 界面解離速度、放射および非放射減衰チャネル、ならびに消滅過程は、 与えられた材料系についてしばしば十分に制約されていませんでした。そのため、多くの実用的なモデリング ワークフローでは、これらの効果は単一のスカラー photon efficiency factor、\(\eta_{\mathrm{photon}}\) に吸収され、 基礎となる輸送と速度論を解像せずに、正味の geminate 損失を表していました。

この状況は現在変わりつつあります。実験的特性評価の進歩により、 BHJ 類似系における励起子寿命、拡散長、および損失チャネルの直接測定が提供され始めています。 これらのパラメータがよりよく制約されるにつれて、明示的な励起子ドメインモデリングは ますます有益になります。これは絶対的なデバイス効率を予測する手段としてではなく、 ドメインサイズ、次元性、および材料速度論がどのように共同で有効な 電荷生成収率を決定するかを探る方法として有益です。このチュートリアルで使用される 励起子ドメインモデルは、この最近の実験–モデリング研究の一部として開発され、 公表されました （たとえば Nature Materials 21, 55–61 (2022) を参照）。

したがって、このチュートリアルでは、理想化されてはいるものの完全に三次元のユニットセル形状を採用します： すなわち、アクセプタマトリクス中に埋め込まれた ドナードメインであり、最初はアクセプタボックス内のドナー球として表現されます。 形状は単純化されていますが、モデルは励起子生成、拡散、界面解離、および競合する損失過程を 明示的に扱います。これにより、実験的に測定されたパラメータがどのように相互作用するかを 調べ、それらの内部整合性を確認し、形態と速度論がどのように共同で三次元における 有効な 光生成効率を制御するかについて直感を育むための、 実用的かつ物理的に明快な枠組みとなります。

2. 支配方程式

このチュートリアルでは、シミュレーション領域全体で三次元 励起子拡散方程式を解くことにより、励起子輸送を明示的に扱います。励起子モデル が有効な場合、光学吸収は励起子集団へ直接供給され、それは次の 方程式に従って発展します：

\[ \frac{\partial X}{\partial t} = \nabla \cdot \left( D \nabla X \right) + G_{\mathrm{optical}} - k_{\mathrm{dis}} X - k_{\mathrm{FRET}} X - k_{\mathrm{PL}} X - \alpha X^2. \]

ここで \(X(\mathbf{r},t)\) は励起子密度（\(\mathrm{m^{-3}}\)）、 \(D\) は励起子拡散係数（\(\mathrm{m^2\,s^{-1}}\)）です。 ソース項 \(G_{\mathrm{optical}}\)（\(\mathrm{m^{-3}\,s^{-1}}\)）は、 光学吸収による励起子生成を表します。 残りの項は、集団から励起子を除去する競合過程を記述します：

• 解離（\(k_{\mathrm{dis}} X\)）： 励起子がドナー–アクセプタ界面で自由電子と 正孔へ変換されること。これは最終的に電荷生成に寄与する生産的チャネルです。
• Förster 型移動（\(k_{\mathrm{FRET}} X\)）： 励起子が隣接サイトまたはドメインへ 非放射的にエネルギー移動すること。これは自由キャリアを直接生成せずに、 励起子エネルギーを再分配します。
• 放射損失（\(k_{\mathrm{PL}} X\)）： 励起子が光子放出 （フォトルミネセンス）によって再結合すること。これは固有の損失経路を表します。
• 消滅（\(\alpha X^2\)）： 高密度における二分子励起子–励起子消滅であり、 強励起や閉じ込められた形状で重要になります。

拡散係数 \(D\) は、励起子拡散長 \(L\)（m） と励起子寿命 \(\tau\)（s）を用いて、

\[ D = \frac{L^2}{\tau}. \]

このチュートリアルでは、励起子のみのシミュレーションを実行します：ソルバーは 三次元励起子密度 \(X(\mathbf{r})\) と関連する反応項を計算しますが、 電子および正孔輸送は解きません。したがって、解離項 \(k_{\mathrm{dis}} X\) は、 ドナー–アクセプタ界面で励起子が自由キャリアへどこでどの程度効率よく変換されるかを示す、 空間分解された潜在的な電荷生成速度として解釈されます。

3. 励起子ドメインシミュレーションを作成する

New simulation ウィンドウを開きます（Figure ??）。 これはメインウィンドウで New simulation をクリックすることでアクセスできます。 Exciton simulations をダブルクリックすると、励起子ライブラリが表示されます（Figure ??）。 Exciton domain をダブルクリックして、サンプルプロジェクトを開きます。

4. 形状とパラメータを確認する

サンプルを開くと、メインウィンドウには単純な三次元シーンが表示されます。すなわち、 アクセプタボックスの中に埋め込まれたドナー球です（Figure ??）。 緑色の光線は上方から入射する照明を示します。 励起子パラメータを確認するには、左側パネル （Device structure タブ内）の Electrical parameters をクリックします。これにより 電気パラメータエディタが開きます。 このウィンドウの下部までスクロールすると、 Excitons とラベル付けされたセクションがあり、 シーン内の各オブジェクト（周囲の層と埋め込み球）ごとに励起子固有パラメータが別々に表示されています。 これは Figures ?? および ?? に示されています。

Excitons セクションでは、拡散–反応方程式に現れるパラメータを定義します。 これには、励起子散乱長 \(L\) が含まれ、これは \(D = L^{2}/\tau\) によって拡散係数を定めます。また、励起子寿命 \(\tau\) は、 追加損失チャネルがない場合の輸送の全体的な時間スケールを制御します。放射再結合は フォトルミネセンス速度 \(k_{\mathrm{PL}}\) によって記述され、一方で非放射エネルギー移動過程は Förster 型速度 \(k_{\mathrm{FRET}}\) によって捉えられます。より高い励起子密度では、 二分子消滅は係数 \(\alpha\) によって含まれ、これは 支配方程式中の \(-\alpha X^{2}\) 項に対応します。最後に、解離速度 \(k_{\mathrm{dis}}\) は、励起子がドナー–アクセプタ界面でどれだけ効率よく電荷移動状態へ変換されるかを指定し、 このモデルにおける有効電荷生成効率を支配する主要パラメータです。

この例の重要な概念的特徴は、sphere（ドナー）と layer（アクセプタマトリクス）で異なる励起子物理を使用していることです。 上のスクリーンショットでは、ドナー球には非ゼロの \(k_{\mathrm{PL}}\)（損失チャネル）が含まれています。 周囲領域は非常に大きな \(k_{\mathrm{dis}}\) で設定されているため、励起子がその領域に到達すると（またはそこへ出ていくと）急速に解離します。 これにより、ドナー／アクセプタ境界は効率的なシンクとして振る舞い、 「中心で高く、境界近傍で枯渇する」という特徴的なプロファイルが生じます。

5. シミュレーションを実行する

青い Run simulation 三角形をクリックしてソルバーを開始します。Terminal 出力 がメインウィンドウに表示されます。これは Figure ?? に示されています。 最初の行では、シミュレーション設定と形状に関する一般情報が報告されます。 重要な診断は励起子ソルバーが動作し始めると現れ、 Exciton f() = … , steps = … という形式を持ちます。この行は、 三次元励起子ソルバーの進行状況を報告します。 量 f() は、現在の励起子密度場が 自己無撞着性からどれだけ離れているかを測る残差であり、ソルバーが反復を進めるにつれて この値は減少するはずです。 併記される steps カウンタは反復番号を示します。 図示された例では、残差は初期値から step 2 で 2.97 × 104 に低下し、その後ソルバーは 方程式が解かれたと報告して終了します。

一般的な現代のノートパソコンでは、このサンプルはおよそ 5–10 秒で完了するはずです。 実行時間が数分に及ぶ場合は、通常、メッシュが過度に細かいか、 あるいは形状やパラメータが変更されてメッシュ点数が大幅に 増加していることを示します。その場合は、先へ進む前に メッシュ設定を見直すか、前の設定手順をたどり直す価値があります。

6. 出力をプロットする

実行が終了したら、Output タブへ切り替えます（Figure ??）。 このチュートリアルで最も重要な項目は、 空間分解された励起子結果を含む exciton_output/ ディレクトリと、 全体の生成および損失統計を要約したファイル exciton_sim_info.json です。これは後で Part C で使用します。 exciton_output/ をダブルクリックすると、Figure ?? に示される内容が表示されます。

光学モデルによって生成された励起子生成率を見るには G.csv を開いてください（Figure ??）。 3D プロットウィンドウでは、マウスでシーンを回転でき、 下部の Z および Y スライダーを使って体積をスライスできます。 これは、形状、照明方向、およびメッシュが 意図した通りに動作しているかを確認する最も速い方法であることがよくあります。

次に、exciton.csv を開いてください（Figure ??）。 これは定常状態の励起子密度を示します。 界面シンクへの拡散の特徴的なシグネチャが観測されるはずです： すなわち、ドナー球内部では励起子密度が高く、 解離が強いドナー–アクセプタ境界に向かって枯渇します。 この枯渇が表示角度のアーティファクトではなく、 界面形状に従っていることをスライススライダーで確認してください。

最後に、Gn.csv を開いてください（Figure ??）。 これは励起子 解離から生じる空間分解された電子生成率を示します。 この例では、ドナー–アクセプタ界面近傍に局在した 特徴的なリング（またはシェル）の生成が見えるはずです。 Z および Y スライダーを用いて スライス面を体積内で移動し、この解離プロファイルが位置に応じてどのように変化するかを調べてください。

7. exciton_output/ にあるファイルとは？

exciton_output/ ディレクトリには、2 種類の異なるファイルが含まれます。 第 1 のグループは、励起子ソルバーによって生成される 計算された場 であり、 励起子密度や解離に由来する生成率などが含まれます。 第 2 のグループは パラメータマップ です。すなわち、 GUI で定義された励起子パラメータの三次元コピーです。これらは主として 整合性確認およびデバッグの補助として書き出され、 意図した材料パラメータが空間内で正しく適用されているかを検証できるようになっています。Table 1 および 2 は、この 2 つのカテゴリを要約しています。すべての量は SI 単位で報告されます。 （対応するサブモデルが無効である場合、一部のファイルは書き出されないことがあります。）

Table 2 にあるパラメータマップファイルは、通常の意味でのソルバー出力ではありません。 それらは単に GUI の Electrical → Excitons セクションで定義された励起子パラメータを、 解場と同じ三次元メッシュ上に書き出したものです。 その目的は、特に複雑な形状で作業する際に、パラメータが ドメイン内の正しいオブジェクトと領域に割り当てられていることを 透過的に確認する方法を提供することです。

8.1 exciton_sim_info.json を解釈する（解離効率）

三次元プロットは、励起子がどこで 生成され、失われ、解離するかを可視化するうえで非常に有用ですが、 それだけでは簡潔な数値要約を与えません。 ファイル exciton_sim_info.json はこの役割を果たし、シミュレーション中の各オブジェクト （周囲の layer、埋め込まれた sphere、および追加した任意の形状）ごとの、 空間積分された 生成率および損失率を報告します。

この情報を見るには、メインの Output タブを開いて exciton_sim_info.json をダブルクリックします。このファイルには、 各オブジェクトについて、全励起子生成率と、各競合する 損失または解離チャネルに対応する全速度を列挙したネストされた JSON 構造が含まれています。 これらの量により、材料 パラメータ、形状、またはモルフォロジーの変化が、電荷生成効率の変化へどのように変換されるかを定量化できます。

最上位キー segment0 は、単にレポートの開始を示します。 その内部で、フィールド name は囲っているオブジェクト （ここでは layer）を識別し、ネストされた shape0 ブロックは埋め込まれた オブジェクト（sphere）に対応します。各ブロックは、そのオブジェクトの 生成、解離、および損失過程への総寄与を報告します。

Table 2 は、このファイルで遭遇するフィールドと、その解釈方法を要約しています。

三次元プロットと exciton_sim_info.json の数値要約を合わせることで、シミュレーションに対する補完的な見方が得られます： プロットは過程がどこで起こるかを示し、JSON ファイルは 各過程がどれだけ寄与するかを定量化します。これにより、 パラメータや形状の変化が励起子 解離および有効電荷生成にどう影響するかを系統的に解析するための 有用なツールになります。

9. 全体をまとめる

三次元結果を総合すると、首尾一貫した物理像が得られます。 光学吸収は、主としてドナードメイン内部で励起子を生成します。 これらの励起子はドナー内部を拡散し、ドナー–アクセプタ界面近傍で枯渇し、 周囲のアクセプタ領域に到達すると効率よく解離します。 空間プロットは、これらの過程がどこで起こるかを示し、 要約レポート exciton_sim_info.json は各過程がどれだけ寄与するかを定量化します。

この例では、要約レポートは空間領域ごとに基礎物理を明瞭に分離しており、 生成、損失、および解離過程を透過的に解釈できます。

• Sphere（ドナー）： 光学吸収の大部分、したがって励起子生成の大部分はドナー球内部で起こり、 これは shape0/sphere ブロック内の大きな tot_G 値に反映されています。 ドナー内部では、この例における主要な競合損失チャネルは放射減衰であり、 tot_Rk_pl によって定量化されます。 tot_Rk_pl と tot_G の比は、 励起子がドナー–アクセプタ界面に到達する前に放射損失がどれほど強く 励起子集団を抑制するかの直接的尺度を与えます。
• Layer（アクセプタマトリクス）： 周囲の layer はこの例では弱くしか吸収しないため、 励起子生成（tot_G）への寄与は小さいです。 その主たる役割は、励起子の変換にあります。layer は強い 解離チャネルで設定されているため、layer ブロックでは tot_Rk_dis が支配的な項になります。 これは空間プロットに見られる界面解離と直接対応しており、 Gn.csv ではドナー–アクセプタ境界に集中した リング状またはシェル状の電子生成領域が示されます。

最終的な指標 dis_eff は、この系全体の解離効率 （ここではおよそ \(89\%\)）を報告します。 物理的には、これはドナーで生成された励起子の大部分が界面へ拡散して 解離し、比較的少ない割合だけが放射減衰に失われることを示します。 Förster 移動や励起子–励起子消滅のような他のチャネルは、この 特定の例では非活性です。

言い換えると、このシミュレーションは 界面解離が励起子減衰より優勢である 領域で動作しており、 その結果として高い有効光子–自由キャリア変換収率が得られます。 この三次元可視化と定量レポートの組み合わせにより、 形状と速度パラメータがどのように共同して 励起子ドメインモデルにおける電荷生成効率を決定するかについて、 コンパクトでありながら物理的に完全な像が得られます。