転送行列法（TMM）を用いた有機太陽電池の光学シミュレーション
OghmaNanoにおいて転送行列法を用いて薄膜型有機太陽電池のシミュレーションを行う方法を学びましょう。本チュートリアルでは、光学的干渉現象や吸収過程、多層構造におけるキャビティ効果が、光子の吸収、電荷キャリアの生成、および有機太陽光発電（OPV）デバイスで得られる最大光電流値にどのように影響するかを解説します。
1. はじめに
このチュートリアルでは、転送行列法（TMM）を用いて薄膜型有機太陽電池（OPV）の光学的挙動をシミュレーションします。このシミュレーションにより、光電場分布、吸収された光子密度、キャリア生成分布、ならびにデバイスの反射・透過スペクトルが算出されます。OghmaNanoにおける転送行列法ソルバーの実装に関する詳細情報はこちらに、またこの手法の完全な導出過程についてはこちらに記載されています。
光学的シミュレーションは、太陽電池デバイスの性能を初歩的に把握するための最も信頼性の高い手段の一つである。なぜなら、必要となる材料パラメータが比較的少ないからだ。各層の屈折率\(n\)および光吸収係数\(\alpha\)のみが必要であり、これらは文献で広く公開されており、通常は正確に測定も可能である。これに対し電気的なドリフト拡散モデルでは、移動度、再結合率、トラップ密度、界面特性といったパラメータは動作条件に大きく依存し、一意に決定することがしばしば困難である。
有機太陽電池デバイスは、活性層の厚みが通常数百ナノメートル程度であるため、光干渉の影響を強く受ける。各材料界面での反射により多層構造内に定在波パターンが形成され、それに伴って光電場や光生成分布が波長、層厚、および光学定数に応じて大きく変化する。このような光学的モデルを用いることで、デバイス内部で光子がどこで吸収され、電荷キャリアがどこで生成されるかを直接的に把握することが可能となる。
転送行列法では、各層内においてマクスウェル方程式を解析的に解き、界面を跨いで進行する前進波および後退波を整合させる。この手法は仮定の数が比較的少なく時間発展の計算も行わないため、計算効率に優れ、FDTDのような全波解析手法よりもはるかに高速である。このような特性から、転送行列法によるシミュレーションは、電気輸送や再結合モデルを導入する前段階として、薄膜太陽電池の最適化において広く用いられている。
この光学モデルでは、達成可能な短絡電流密度の上限値も算出される\(J_{SC}\)。多層構造体が光を効率的に吸収できない場合、たとえ電気的輸送が理想的であっても高性能なデバイスは実現できない。光学最適化に関する具体例としては、光学フィルター、反射防止コーティング、ファブリ・ペロー共振器、ペロブスカイト太陽電池、およびOLEDの光取り出し特性に関するチュートリアルが挙げられる。
OghmaNanoにおいては、転送行列法ソルバーによって算出された光生成分布を電気的なドリフト拡散モデルに直接結合させることが可能であり、これにより光吸収、キャリア生成、輸送、再結合、および抽出の各過程を自己整合的にシミュレーションできる。
2. シミュレーションの作成
??に示される新規シミュレーションウィンドウを開きます。次に「有機太陽電池」カテゴリを選択し、??に示されるPM6:Y6_E6_0hrsのサンプルを選択してください。
シミュレーションをディスク上の適切なディレクトリに保存してください。シミュレーションが読み込まれると、??に示すメインインターフェースが表示されます。
このシミュレーションでは、透明電極、輸送層、PM6:Y6活性層、および金属製裏面電極からなる完全な多層構造の有機太陽電池モデルが用いられている。デバイスの厚みが光の波長と同程度であるため、干渉効果が内部の光場分布に大きく影響を与える。
3. デバイス構造の確認
左側のツールバーからレイヤーエディターを開くと、??に示されている多層構造が表示されます。
このデバイスは、インジウム錫酸化物（ITO）製の透明電極、ZnOからなる電子輸送層、PM6:Y6からなる活性層、MoO\(_x\)からなる正孔輸送層、および銀製の裏面電極で構成されている。
各層には実験的に測定された光学定数 \(n\) および \(k\) が設定されており、これにより波長依存性の吸収現象や干渉効果を正確にシミュレーションすることが可能となる。
薄膜型有機太陽電池においては、建設的干渉の影響により光電場はしばしば活性層内部で最大値に達する。一方で、破壊的干渉により特定領域における光電場が抑制されることもある。こうした干渉効果は、電荷キャリアの生成分布が光電場分布と直接結びついているため、最終的な電流電圧特性に強く影響を与える。
4. 転送行列ソルバーの起動
??に示される光学タブを開いてください。このリボンには転送行列の制御機能、光学メッシュの設定項目、光学検出器、および光源エディタが含まれています。
「転送行列」ボタンをクリックすると、??に示す転送行列シミュレーションエディタが開きます。光学的シミュレーションを実行するには「光学的シミュレーションを実行」ボタンを押してください。実行中には??に示すターミナル出力が表示されます。このソルバーは、多層構造を持つ有機薄膜太陽電池における光子密度分布、吸収された光子密度、生成率分布、反射スペクトル、透過スペクトルを算出します。
転送行列法では各層内における光学的問題が解析的に解かれるため、多数の波長や複数の吸光層を有するデバイスであっても計算速度は極めて速い。
5. 光子密度と吸収された光子
シミュレーションが完了すると、??に示される光子密度分布が表示される。横軸はデバイス内の位置を、縦軸は波長を表している。この図は、材料界面での反射に起因する干渉効果も含め、多層構造内を伝播する光子の空間分布を示している。
「吸収された光子分布」タブを選択すると、??に示す吸収光子分布図が表示される。デバイス内部の全光場強度を示す光子密度プロットとは異なり、この図では物質によって実際に吸収され電子・正孔対へと変換された光子のみが示されている。
局所的な吸収光子密度は\[ G_{\mathrm{opt}}(x,\lambda)=\alpha(x,\lambda)\,\phi(x,\lambda) \]によって与えられる。ここで\(\alpha\)は波長依存の光吸収係数、\(\phi\)は局所的な光子密度である。光吸収の大部分は、電荷生成が主に起こるPM6:Y6活性層内で生じる。
波長が短い場合、光吸収係数が比較的大きいため、光子は表面界面付近で吸収される。一方、波長が長くなると光吸収長が増大し、光子は活性層のより深い部分まで侵入した後に吸収される。また、薄膜キャビティ構造内において進行波と後退波が干渉することで、強い定在波パターンも観測される。
多くの場合、全スペクトル分布を直接見るよりも、波長にわたって積分した総生成分布を調べる方が有用である。これは??に示されている生成率タブにおいても確認でき、そこではデバイス内部の位置の関数として体積当たりの電子・正孔対生成率の総量がプロットされている。
生成率プロット上に重ねて表示されている着色領域は、デバイス各層の電子エネルギー準位を表している。これらのエネルギー準位は材料データベースから取得されたもので、主にデバイス構造を視覚的に把握するための参考情報として用いられる。これらは光学的な転送行列法シミュレーション自体には直接的な影響を与えない。ただしPM6:Y6活性層に関しては、電荷輸送、再結合、および電荷抽出のシミュレーションを行う際に電気的なドリフト拡散モデルによってこれらのエネルギー準位が利用される。
6. 反射スペクトルおよび透過スペクトル
反射した光パワーが示されている ??伝送された光パワーについては、そこに示されている。 ??反射による損失は、活性層に結合される光エネルギーの量を減少させる。一方、透過光とは、吸収されることなくデバイスから放出される光子に相当する。
デバイス内部における光パワーは、反射・吸収・透過の各過程に分配される。この分配比率は\[ R(\lambda)+A(\lambda)+T(\lambda)=1 \]に示されており、\(R\)、\(A\)、\(T\)はそれぞれ反射光パワー、吸収光パワー、透過光パワーの割合を表す。この関係式はエネルギー保存則から直接導かれるものであり、転送行列法によるシミュレーション結果の整合性を確認するための有用な指標となる。
長波長領域においては、銀製の裏面電極が強力な光反射体として機能するため透過率はほぼゼロに近づく。その結果、長波長の光子の大部分は多層構造内で吸収されるか、あるいは前面表面へと反射され、活性層内における有効な光路長が増大することになる。
実際の有機太陽電池設計においては、光干渉現象や層の厚さ、反射損失などにわずかな変化が生じただけでも短絡電流密度に顕著な変動が現れる。このため転送行列法によるシミュレーションが広く用いられ、キャビティ効果や光学スペーサー、電極厚さを最適化することで活性層内での光吸収を最大化する試みが行われている。
8. 最大可能光電流の算出方法
これらの値は、生成されたすべての電子・正孔対が確実に収集されるという仮定の下で、光学モデルによって予測される最大の光電流を示しています。再結合、トラップ、輸送損失、および電荷抽出の非効率性といった要因は、この純粋に光学的な計算には含まれていません。
最大光電流密度は、波長および位置にわたって光生成率を積分することで求められる。\[ J_{SC,\mathrm{max}} = q \int \int G(x,\lambda)\, d\lambda\, dx \] ここで\(q\)は素電荷、\(G(x,\lambda)\)はデバイス内部における波長依存の体積生成率である。
最適化が十分に行われたデバイスにおいては、実測される短絡電流密度 \(J_{SC}\) はこの光学的な上限値にかなり近い値となることが多い。したがって、光学的シミュレーションは実現可能なデバイス性能の初期見積もりとして極めて有用である。もし光学的モデルにより最大光電流値が低く予測される場合、たとえ電気輸送モデルをいかに最適化しても高効率な太陽電池を実現することは不可能である。なぜなら、そのデバイスは根本的に十分な光を吸収していないからだ。
9. 照射スペクトルの変更方法
??に示されている光学関連の画面に戻り、灯台のアイコンとして表示される「光源」ボタンをクリックする。これにより??に示す光源エディタが開く。
デフォルト設定では、シミュレーションには標準的なAM1.5G太陽光スペクトルが用いられる。ただし、長波長側ではスペクトルが切り捨てられている点に留意されたい。多くの有機太陽電池シミュレーションにおいてはこれで十分であり、なぜなら活性層は赤外領域での吸収が比較的弱く、ほとんどの光電流生成は可視光領域で起こるからである。波長範囲を限定することで計算コストが削減され、シミュレーション速度も向上する。
実際には、多くの実験測定においては実際のAM1.5G太陽光ではなく実験室用の太陽光シミュレータが用いられる。特にキセノンランプは、その放射スペクトルが太陽光の可視領域をほぼ再現するため広く利用されている。実験で用いられる照射スペクトルを正確に合わせることは極めて重要であり、スペクトル強度の変化によって光生成分布や測定されるデバイス性能が大きく変動する可能性があるからである。
スペクトルのドロップダウンメニューから照明光源をキセノンに変更し、Rebuildボタンを押してください。??に示される生成されたスペクトルは、依然として長波長側で若干切り取られています。これは、光学波長メッシュの範囲がまだキセノンの全発光スペクトルを捉えるには不十分だからです。
波長範囲を拡大するため、光学設定画面に戻り ?? に示される光学メッシュエディタを開いてください。波長の上限値を約1500 nmまで増加させた後、Rebuildボタンを押してください。これにより光源設定画面に戻ると、?? に示される完全なキセノン発光スペクトルが得られるようになります。
照明スペクトルを変更した後、転送行列法によるシミュレーションを再実行してください。その後Generation rateタブに戻ると、??に示すような分布が得られます。荷電キャリアの生成分布は元のAM1.5G条件時とは大きく異なっており、これはキセノンランプが異なるスペクトル強度分布を生み出し、多層構造内部における光干渉条件も変化したためです。
こうした効果は、特に活性層の厚みが光の波長と同程度である薄膜型有機太陽電池において重要となる。そのため、照射光のスペクトルやキャビティ内での干渉現象、吸収長にわずかな変化が生じただけでも、短絡電流密度やデバイス全体の効率に測定可能な変化が現れるのである。
9. 光生成現象を電気モデルへ結合させる方法
転送行列法に基づくシミュレータによって算出された吸収光子分布は、ドリフト拡散法シミュレータ内において電気的生成項へと変換される。
この変換処理は ?? に示されている光子効率パラメータによって制御される。
物理的には、このパラメータは吸収された光子が自由な電子・正孔対を生成し、電流に寄与する確率を表している。
有機太陽電池においては、吸収された光子は最初に励起子を形成するため、この値は1よりも大幅に小さくなる可能性がある。このような励起子の解離や同時再結合、さらには電荷移動に伴う損失により、最終的な光電流が著しく減少することがある。
10. 光学スナップショットの確認
optical_snapshotsフォルダが含まれる出力ディレクトリ。
??に示される出力タブ内で、optical_snapshotsディレクトリをダブルクリックします。これにより光学スナップショットビューアが起動し、多層構造を持つ有機薄膜太陽電池全体において波長ごとの各種光学量を確認できるようになります。
これらの光学スナップショットには、光学メッシュ内の各波長における光子密度、光吸収係数、前方伝播電界 \(E^+\)、後方伝播電界 \(E^-\)、吸収された光子密度、および生成率の空間分布が記録されている。
??に示されている光子密度分布は、光学的干渉がデバイス内部における光の空間分布にどのように影響を与えるかを示している。また??に示されている吸収分布は、多層構造における波長依存性のある吸収係数\(\alpha\)に対応している。
??および??に示されている電界分布は、転送行列法ソルバーによって算出された進行波および後退波の光波に対応する。これらの電界間の干渉作用により定在波パターンが形成され、それが薄膜有機太陽電池内部における極めて不均一な光生成分布の原因となる。
11. まとめ
このチュートリアルにおいて、転送行列法を用いてPM6:Y6系有機太陽電池内部における光吸収、干渉現象、およびキャリア生成の様子をモデル化しました。また、光子密度分布、吸収された光子の分布状況、生成率分布、反射・透過スペクトル、さらには光学的モデルから予測される最大光電流値についても検討しました。
また、照射光のスペクトルを変化させることで、多層構造内部における光場分布や光生成分布がどのように変化するかも調査しました。こうした効果は特に薄膜型有機太陽電池において重要であり、光干渉やキャビティ構造のわずかな変化でもデバイス性能に顕著な影響が現れるからです。
転送行列法は、多層構造を持つ太陽電池デバイスをモデル化するための計算効率に優れ、物理的にも妥当性の高い手法である。この手法では必要な材料パラメータが比較的少なく、かつ光吸収や光電流の限界値に関する詳細な情報も得られるため、より複雑な電荷輸送や再結合モデルを導入する前段階として、有機太陽電池デバイスの最適化において広く用いられている。
OghmaNanoにおいては、転送行列法ソルバーによって算出された光生成分布をドリフト拡散モデルに直接結合させることが可能であり、これにより光吸収、電荷生成、輸送、再結合、および電荷抽出の各プロセスを自己整合的にシミュレーションできる。