有機エレクトロニクス（OLED）の界面設計：正孔注入障壁の算出#
有機エレクトロニクス（OLED）デバイスの動作は、電極からの電荷注入、輸送、そして再結合による発光という一連のプロセスで成り立っています。デバイスの低消費電力化には電極界面の注入障壁の低減が、また所望の色を得るには発光層のエネルギーギャップ制御が不可欠です。本事例では、第一原理計算ソフトウェアAdvance/PHASEの固体表面解析機能と、量子化学計算ソフトウェアNWChemの高精度な分子励起状態計算機能を組み合わせ、デバイスの「入口（注入）」と「出口（発光）」の両面を評価するワークフローを紹介します。
Keywords: 有機EL (OLED), 第一原理計算, NWChem, 仕事関数, 正孔注入障壁, TD-DFT, アントラセン
理論背景：本解析における評価指標の物理的意味#
OLEDの評価において、Energy Level Alignmentが重要です [1]。 本解析では、以下の指標を用いて物理的に妥当なバンドアライメントを構築します。
- 仕事関数： 金属表面のフェルミ準位（）の絶対値。正孔注入の「入口」の高さを示します。解析には平面波・擬ポテンシャルコードAdvance/PHASEを使用します。
- 有機固体のイオン化ポテンシャル： 孤立分子のHOMO（軌道エネルギー）を用いるだけでは、気相の近似値に過ぎません。本手法では、PCM（連続誘電体モデル）を用いてアモルファス有機固体内の誘電環境（）を模し、中性状態とカチオン状態の全エネルギー差から「固体のイオン化ポテンシャル（）」を算出するSCF法を採用しています。解析にはオープンソースの量子化学ソフトウェアNWChem [2] を使用します。
- LUMO*： 単なる空軌道のエネルギーではなく、算出されたHOMOレベル（）に、TD-DFTで求めた「最低励起エネルギー（光学ギャップ）」を加算することで、励起子効果を取り込んだ電子側の基準レベルを算出します。厳密には、分子の構造緩和（ストークスシフト）が存在するため、このHOMOとLUMO*の差は実際の発光エネルギーよりもやや大きくなります。しかし、キャリア注入の障壁やデバイスの基本的なエネルギーアライメントを議論する上では、このLUMO*が電子側の基準レベルとして広く用いられます。解析にはNWChem [2] を使用します。
Advance/PHASEによる電極の絶対電位計算#
陽極材料である金電極に対して、その最安定面であるAu(111)面を対象として、7層のスラブモデルとESM (Effective Screening Medium)法 [3] を用いて仕事関数を算出しました（図1）。得られた仕事関数5.22 eV （実験値：5.26 eV [4]）は、真空準位（0 eV）から見たフェルミ準位の深さを示しており、有機層への正孔注入の起点となります。
図1. 左：Au(111)面のスラブモデル（視認性のため、4x4スーパーセルを示しています。半透明表示の中間層をバルク層として固定しています）。右：Au(111)面の面内平均ポテンシャルプロファイル。
NWChemによる分子特性と発光エネルギーの算出#
発光層のモデル分子としてアントラセンを用い、長距離補正汎関数（CAM-B3LYP）および基底関数 6-311+G** を用いたTD-DFT計算を行いました。得られた吸収スペクトル（図2）の第1励起状態は3.62 eV (342.3 nm) であり、これが本計算における理論上の光学ギャップ（発光エネルギーの上限）となります。
図2. NWChemによるシミュレーション吸収スペクトル。最低励起エネルギーを光学ギャップとして採用します。
物理的解釈の補足：発光波長とストークスシフト
孤立したアントラセン分子の垂直吸収エネルギー（342.3 nm）は近紫外領域に位置しますが、実際のデバイス（固体薄膜）では、分子の励起状態における構造緩和（ストークスシフト）や周囲の分極効果により、発光はより長波長側（青色領域）へシフトします。本解析では、材料固有の特性である「光学ギャップ」を基準に議論を進めます。
界面エネルギーアライメントの統合評価#
Advance/PHASEとNWChemの計算結果を統合し、界面のエネルギーバンド図を作成しました（図3）。
図3. Au(111) / アントラセン界面のエネルギー準位図。真空準位を0 eVとしています。
【解析結果の考察】
NWChemを用いたSCF法による計算の結果、アントラセンの気相でのイオン化ポテンシャル（7.33 eV）に対し、PCMを用いた固体状態（アモルファス）のイオン化ポテンシャル()は 6.15 eVと算出されました。この差（1.18 eV）が、薄膜内での正孔を安定化させる固体分極エネルギー（）に相当します。
図3のエネルギー準位図において、HOMOレベルは真空準位から見て （-6.15 eV）の位置に定義されます。これとAdvance/PHASEで得られたAu(111)の仕事関数（5.22 eV）を比較することで、実際のデバイス環境に即した正孔注入障壁は 0.93 eVと見積もられました。気相の軌道エネルギー（HOMO）をそのまま用いた場合の過大評価が、誘電体モデルの導入により物理的に正しく補正されています。
また、光学ギャップ（3.62 eV）を基にしたLUMO*準位は-2.53 eVとなり、発光エネルギーとキャリア輸送の両面を矛盾なく記述できる、実践的な統合界面モデルが構築されました。
補足：真空準位の接続（Schottky-Mott則）と界面双極子
図3のエネルギー図は、金属と有機層が相互作用しない理想的な極限（Schottky-Mott則）を仮定し、真空準位を平坦に接続しています。しかし、実際の有機/金属界面では、電荷移動や「push-back効果（金属表面の電子雲の押し返し）」によって界面双極子（Interface Dipole, ）が形成され、真空準位にシフトが生じることが知られています（参考文献[1]）。本事例の手法は、「一次スクリーニング・モデル」として有効ですが、さらに厳密な界面アライメントを評価する場合は、界面双極子の計算が必要となります。
まとめ#
本事例では、第一原理計算ソフトウェアAdvance/PHASEとNWChemを連携させ、電極の「仕事関数」と分子の「光学ギャップ」を組み合わせることで、OLED界面の注入・発光特性を包括的に評価しました。この手法は、新規材料のスクリーニングにおいて、熱力学的な注入条件と光学的な発光条件を同時に評価するための有効な手段となります。
本解析の詳細や、研究への適用可能性に関するご相談はこちら
お問い合わせ参考文献#
- H. Ishii, K. Sugiyama, E. Ito, and K. Seki, "Energy Level Alignment and Interfacial Electronic Structures at Organic/Metal and Organic/Organic Interfaces", Adv. Mater. 11, 605 (1999).
- E. Aprà et al., "NWChem: Past, present, and future", J. Chem. Phys. 152, 184102 (2020).
- M. Otani and O. Sugino, "First-principles calculations of charged surfaces and interfaces: A plane-wave nonrepeated slab approach", Phys. Rev. B 73, 115407 (2006).
- G. V. Hansson and S. A. Flodström, "Photoemission study of the bulk and surface electronic structure of single crystals of gold", Phys. Rev. B 18, 1572 (1978).
関連ページ#
- 第一原理計算ソフトウェア Advance/PHASE
- 解析分野：ナノ・バイオ
- 産業分野：材料・化学