大面積有機太陽電池シミュレーション (PM6:Y6): 3Dコンタクトネットワークから照射デバイスへ
1. はじめに
このチュートリアルでは、大面積デバイスシミュレーションのシリーズを継続し、小型の実験室セルを超えて実用スケールのジオメトリへ移行したデバイスに焦点を当てます。 これは、マニュアル内の既存の2つのワークフローの中間に位置します: 大面積コンタクトの設計（アクティブデバイスなし） と 複雑な3Dペロブスカイトモジュール です。 この例はその2つの橋渡しを行い、以前に設計した コンタクト構造 にダイオードと光吸収を導入しますが、完全な 相互接続モジュール までは進みません。
ここでは、印刷された六角形金属コンタクトを持つPM6:Y6型アーキテクチャに基づく単一の大面積有機太陽電池をシミュレーションします。 目的は、小面積では良好に動作するデバイスが、スケールアップされたときにどのように振る舞うかを理解することです。 そこでは、横方向電流、電圧降下、およびコンタクト抵抗が支配的になり始めます。 このチュートリアルではOPV材料を使用しますが、ワークフロー自体は一般的です: 光学n/kデータを差し替え、部品値を調整することで、同じモデルはあなたの研究デバイスを表現できます。
このワークフローを自分のデバイスへ適用するには、2組の入力を与えるだけです: 層の抵抗率（通常は文献または簡単な試験構造から取得）、 および理想因子、光電流、飽和電流などの少数のダイオードパラメータです。 これらを設定すれば、このモデルによって実験室スケールのデバイスを仮想的にスケールアップし、 それが大面積の実環境ジオメトリへ移されたときに、どの損失機構が効率を制限するかを特定できます。
ここでは完全な3Dドリフト–拡散ソルバは使用しませんが、OghmaNanoにはそれも用意されています。 大面積デバイスでは、完全に解像された3Dドリフト–拡散はすぐに計算的に非現実的になり、 支配的な制限が各層を通るキャリア輸送ではなく、横方向輸送とコンタクト抵抗から生じる場合には しばしば不要です。 その代わりに、OghmaNanoはジオメトリを考慮した3D回路モデルを構築し、 横方向輸送とコンタクトを抵抗要素で表現し、光起電力挙動を埋め込みダイオードを用いて捉えます。
電気問題は回路形式で解かれますが、光学は厳密に扱われます。 フォトジェネレーションは、三次元ジオメトリ全体にわたる一次元スライスで評価されるtransfer-matrix光学ソルバを用いて計算され、 薄膜干渉と吸収効果を保持しつつ、回路ソルバが大面積電流収集を処理します。
2. シミュレーションの作成
メインウィンドウで New simulation をクリックします。シミュレーションカテゴリの一覧が表示されます（?? を参照）。次に Large area PM6:Y6 solar cell をダブルクリックします （?? を参照）。
シミュレーションを作成すると、メインウィンドウが開き、デバイスジオメトリが三次元で表示されます （?? を参照）。 上部コンタクトは Ag と表示された六角形（ハニカム）金属メッシュです。 この銀層は、高導電性の電流収集グリッドとして機能します。 メッシュとして実装されているため、表面の大部分は開いたままであり、光がコンタクトを通過できます。銀メッシュの下には PEDOT:PSS があり、横方向の電流拡散層として機能します。大面積コンタクトのチュートリアルと同様に、電荷はまずポリマー層内で横方向に収集され、 その後、効率的な抽出のために低抵抗の銀グリッドへ移されます。
光活性層は PM6:Y6 です。 電気的には、この層は光学生成を有効にしたダイオード要素として表現され、 追加の抵抗成分により、ソルバはその有限な導電性を考慮します。スタックの最下部には、デバイス全面積にわたって広がる連続的な Ag コンタクトがあり、電気接続を完結させます。
メインウィンドウに見える緑色のバーは浮いているように見えるかもしれませんが、 それは単に抽出境界条件の視覚的なインジケータです。 これは、ソルバが回路ネットワークから電流を抽出するために使用する電気コンタクトを適用する領域を示しています。
3. エピタキシャル層スタックとコンタクトの確認
このシミュレーション（および以前の大面積コンタクトの例）は、層状のエピタキシャル構造を用いて定義されています。 デバイスは層エディタを用いて薄膜層のスタックとして構築され、これにより各層が次の層のきれいな上に配置され、 ギャップがなく、厚さが一貫して扱われます。 このアプローチは、垂直スタックが固定され、横方向挙動が回路モデルによって導入される平面大面積デバイスに適しています。 これは、後のモジュールスタイルのシミュレーション（ 相互接続モジュールチュートリアル を参照）とは異なり、 そこではフィンガー、バスバー、およびより複雑なレイアウトが完全な3Dジオメトリを用いて明示的に定義されます。
構造を確認するには、Layer editor を開いて層スタックを表示し （?? を参照）、 次にデバイス構造エディタで Contacts パネルを開きます （?? を参照）。 ジオメトリ的には、この構成は単純ですが重要です。下部コンタクトはデバイス全面積に広がり、 一方で上部コンタクトははるかに小さく、電流が抽出される領域を定義しており、 これは3D表示に示された緑色の抽出マーカと一致しています。 このコンタクトジオメトリを理解することは不可欠であり、それが大面積デバイス挙動を支配する横方向電流経路を設定します。
4. メッシュの確認
Circuit diagram タブへ切り替え、左下の refresh アイコン （リサイクル風のボタン）をクリックします。この操作により、定義されたジオメトリと層スタックから直接、デバイスの回路表現が構築されます。生成された回路メッシュを ?? に示します。 メッシュ内の各リンクは、抵抗やダイオードなどの回路要素に対応し、その色はそれが由来する層または材料領域を示します。
表示を回転させることで、上部および下部コンタクト領域を明確に識別できます。 ?? および ?? では、 青い点 が回路メッシュから電流が抽出されるノードを示しています。 これらのノードはコンタクト境界条件を定義し、回路ソルバが デバイスから電流を収集する位置を表します。
5. 電気パラメータと光学生成
デバイス構造パネルから Electrical parameters を開きます。このエディタは、大面積回路モデルで使用される抵抗およびダイオード特性を制御します。 設定例を ?? および ?? に示します。PM6:Y6 層では、ideality factor n および reverse-bias (saturation) current I0 を含む、ダイオード要素固有のパラメータが表示されます。これらの値を自分の実験曲線から得た値に置き換えると、自分のデバイスがどのようにスケールアップするかを見ることができます。また、optical charge-carrier generation を有効にすることもできます。光学生成が有効な場合、transfer-matrix光学ソルバは光生成項を与え、それが電気モデルへ photocurrent Iph として照射ダイオード方程式に入ります:
I(V) = I0 [ exp( qV / (n k T) ) − 1 ] − Iph
このチュートリアルでは、このトグルは実質的に Iph 項を制御します。光学生成を無効にすると、回路は暗ダイオードネットワークとして振る舞います。 有効にすると、照明がダイオード要素を通じて回路へ電流を注入し、Iph はtransfer-matrix光学から得られる吸収光子フラックスによって決定されます。
💡 実装上の詳細: 活性層は、バルク輸送および直列抵抗を捉えるために、その厚さ方向に複数の抵抗要素へ離散化されることがあります。 ただし、回路モデルには常に単一のダイオード層が含まれ、通常は活性領域の下部に配置されます。 これは、ダイオードがデバイス全体のJV応答とフォトジェネレーションを表現し、 一方で抵抗要素が材料内部の電圧降下と電流拡散を記述することを反映しています。
6. シミュレーションの実行
シミュレーションを開始するには、Run simulation（メインツールバーの青い三角形）をクリックするか、F9 を押します。 ソルバはただちに実行を開始し、出力ターミナルはリアルタイムで更新されます （?? を参照）。
出力の初期に現れる色付きの行は光学計算に対応しています。 この段階では、transfer-matrixソルバがデバイス全体にわたるフォトジェネレーションを スライスごとに評価し、空間的に分解された生成プロファイルを構築しています。 これが完了すると、ソルバは電気問題へ移行し、回路ネットワークの解法を開始します。
シミュレーション実行中のターミナル出力には注意を払う価値があります。 この情報は診断ノイズではありません: シミュレーションが物理的に妥当に振る舞っているかどうかについて有用な洞察を与えます。 示されたスナップショットでは、およそ −9.48 × 102 A m−2 が一方のコンタクトから流れ出ており、 他方からは −9.48 × 101 A m−2 が流れ出ていることが分かります。 同時に、報告されているソルバ残差 F = 1.7 × 10−14 は、電流連続性が極めて高い精度で満たされていることを示しており、 最初のJV点は約 17 ms で解かれています。
2つのコンタクト間の電流密度の差は完全に予想されるものです。 下部コンタクトはデバイス全面積に広がる一方、上部コンタクトは小さな領域しか占めません。 その結果、同じ総電流が上部コンタクトでははるかに小さい面積を通じて抽出されるため、 そこでは有意に高い電流密度になります。 このようにターミナル出力を読むことは、シミュレーションのジオメトリと境界条件が物理的に一貫していることを確認するための、迅速で効果的なサニティチェックを提供します。 残差誤差が1のオーダーまで増大したり、実行時間が異常に長くなったりする場合には、 通常それは壊れたメッシュや切断されたコンタクトのような問題を示しています。シミュレーションが完了したら、Output タブを開きます （?? を参照）。 利用可能な出力は、 前のチュートリアル でより詳しく説明され、 次のチュートリアル でさらに詳しく扱われます。 ここでは、コンタクトごとに分解されたJV曲線のみに焦点を当てます。
各電気コンタクトには、それぞれ専用のJVファイル jv_contact0.csv および jv_contact1.csv があり、 その特定のコンタクトから流れ出る電流を記録します。 電流密度は異なりますが、コンタクト面積を考慮すると、抽出される総電流は同じです。 jv_contact0.csv を開くと、以下に示すJV曲線が得られます。
7. 材料パラメータの編集
材料特性はオブジェクトエディタを通じて直接変更できます。 3D表示内の任意の層を右クリックして Edit object を選択すると、 オブジェクトエディタが開きます （?? および ?? を参照）。
ここから、オブジェクトに割り当てられた 光学材料 を変更したり、 その幾何学的形状を修正したり（たとえばハニカムメッシュのようなパターン化コンタクトを使用する場合）、 色のような視覚属性を調整したりできます。 このエピタキシャルスタイルのシミュレーションでは、絶対的なオブジェクト位置は主に層スタックによって定義されるため、 並進や回転のようなパラメータは通常ほとんど影響しません。 最も重要なのは、光学材料の割り当てとオブジェクトジオメトリです。
💡 試してみてください: PM6:Y6 層を右クリックし、その光学材料をペロブスカイトに置き換えてください。 その後、シミュレーションを再実行し、JV曲線がどのように変化するかを観察してください。 必要であれば、理想因子、 飽和電流、および光電流のようなダイオードパラメータも、ペロブスカイト文献で報告されている値を用いて更新できます。