分子集合体と固体の光吸収スペクトルと光電子特性を計算するためのエキシトニック・ハミルトン人

過去20年間で、凝集した有機分子から作られた固体およびフィルムは、光電子デバイスで複数の用途を発見しました。このような材料に基づくデバイスは、軽量化、柔軟性、低消費電力、インクジェット印刷を使用した安価な生産の可能性など、多くの魅力的な特性を持っています。有機発光ダイオード(OLD)に基づくディスプレイは、液晶ディスプレイを携帯電話、ラップトップ、テレビ、およびその他の電子機器^{1、2、3、42,}の最先端として置き¹換えています。^3,4照明アプリケーションの OLD の重要性は、今後⁴年で増加すると予想されます。有機太陽光発電デバイスの性能は着実に改善しており、単接合型有機太陽電池⁵に対して最近報告された電力変換効率は16%を超えています。有機材料はまた、光ファイバー通信などの他の技術を破壊する可能性があり、その使用により、15 THz以上の^66,7以上の非常に高い帯域幅を有する電気光学変調器の開発が可能^になる。

オプトエレクトロニクスの用途に向けて固体分子材料を最適化する上での大きな課題は、その特性が一般的に材料のナノスケール構造に大きく依存していることである。製造プロセスでは、化学気相成長、光学活性分子の⁸個のテンプレート化などの制御された成長技術(すなわち、ポリマーマトリックス^99、10)、熱アニーリング^{11、12、,12}などによって、ある程度材料のナノ構造を定義することができます。¹⁰しかし、ナノスケールの疾患は、ほとんどの分子材料に固有であり、通常は完全に排除することはできません。したがって、有機光電子材料の合理的な設計には、物質の特性にどのような影響を与えるかを理解し、最適な性能を得るためにそれを設計する方法を見つけることが不可欠です。

分子材料の障害の程度は、通常、バンド理論で記述できる電子構造を有する周期的な結晶構造の摂動として扱うには大きすぎる。一方、バルク材料やフィルムの特性を再現するためにシミュレーションに含まれなければならない分子の数は、密度関数論^(DFT)13、14および時間依存密度機能理論(TD-DFT)15,16のような量子化学的方法を第一原理に使用するには大きすぎる。^{,14 15,16}オプトエレクトロニクスに応用されている有機分子は、典型的には、比較的大きなπ共役系を有する。多くはドナーとアクソクターグループを持っています。このような分子における正しい電荷移動挙動を捕捉することは、それらの光電子特性を計算するために不可欠であるが、TD-DFT^{17、18、19、2018,19,20}の¹⁷長距離補正ハイブリッド機能を使用してのみ達成することができる。このような機能を用いる計算は、システムの大きさに応じて超直線的にスケールし、現在のところ、〜104未満の原子基底関数を使用して記述できる個々の有機分子または小分子集合体の光電子特性を⁴モデル化する場合にのみ実用的である。多数のクロモフォアで構成される無秩序な材料を記述できるシミュレーション方法は、これらのシステムをモデル化するのに非常に有用であろう。

分子材料における分子間相互作用の大きさは、多くの場合、材料を構成する個々の分子間のエネルギーパラメータ(固有エネルギーや励起エネルギーなど)の変動の順序に匹敵するか、または小さい。このような場合、マルチスケールモデリングは、大きな障害分子系^21,22,23,22の光電子特性を理解し^最適化23するための最も有望なアプローチである。このアプローチでは、第一原理量子化学的方法(通常DFTおよびTD-DFT)を用いて、物質を構成する個々の分子の特性を正確に計算します。バルク分子材料を表すのに十分な大きさの材料サンプルのHamiltonian(おそらく、周期的な境界条件を採用することによって)は、個々の分子について計算されたパラメータを使用して構築される。このHamiltonianは、大きな分子集合体、薄膜、またはバルク分子材料の光電子パラメータを計算するために使用することができます。

励起子モデルは、分子材料の励起状態が励起基に表現される多尺度モデルのサブクラスである: クーロンの引力^24,2525によって結合される電子正²⁴孔対。多くの励起状態のプロセスをモデル化する場合、電子と正孔が同じ分子に局在するフレンケル励起子²⁶のみを含むだけで十分である。電荷移動励起子は、電子と正孔が異なる分子に局在する場合、場合によっては含める必要がある(例えば、ドナー・アクセプター系における電荷分離をモデル化する場合^)27,28。,28励起子モデルは、個々の分子に対する第一原理計算のみを使用してパラメータ化できるマルチスケールモデルですが、分子間相互作用を引き続き説明しています。彼らが説明できる2つの主要な相互作用タイプは、(a)分子間のエキシトン結合または分子間の偏光または分子間の移動能力を特徴付ける分子間の励起的な結合であり、(b)周囲の分子の電荷分布によって分子上の電子密度の静電分極化である。我々は、これらの要因の両方が、光吸収スペクトル²⁹および第1超極極性³⁰のような分子集合体の光学的および電気光学的特性をモデル化するために重要であることを示した。

本論文では、大きな分子集合体やバルク分子材料の光学スペクトルやその他の光電子特性の計算に使用できる、励起子モデルをパラメータ化するためのプロトコルを紹介する。興奮性のハミルトンは、密結合のハミルトン^{24、25であると25}仮定され、

ここでε_iは、物質中のi^番目の分子の励起エネルギーであり、b_ijは、i^番目とj^番目の分子の間の励起的結合であり、â_i^†と â _iはそれぞれ、材料中のi^番目の分子の励起状態に対する生成および消滅演算子である。エキシトン・ハミルトンのパラメータは、材料を構成する個々の分子に対して行われるTD-DFT計算を使用して見つかります。これらのTD-DFT計算では、材料内の他のすべての分子の電荷分布は、分子の電子密度の静電分極化を考慮して原子点電荷の静電的埋め込みによって表される。個々の分子に対する励起エネルギー ε_iは、TD-DFT計算出力から直接取り出される。分子間の励起結合は_、bij、b遷移密度キューブ法³¹を用いて計算され、ガウシアン³²におけるTD-DFT計算の出力から得られた相互作用分子の地盤から励起状態への移行密度を有し、Multiwfn多機能波関数アナライザ³³を用いて後処理される。バルク分子固形物の特性をシミュレートするために、周期的境界条件を、Hamiltonianに適用してもよい。

現在のプロトコルでは、ユーザがガウス³²およびMultiwfn³³プログラムにアクセスできる必要があります。このプロトコルはGaussian 16、リビジョンB1、Multiwfnバージョン 3.3.8 を使用してテストされていますが、これらのプログラムの他の最近のバージョンでも動作するはずです。さらに、このプロトコルは、カスタム C++ ユーティリティと、カスタム Python 2.7 および Bash スクリプトを使用 https://github.com/kocherzhenko/ExcitonicHamiltonianします。この計算は、Unix/Linux ファミリのオペレーティング システムを実行しているマシンで実行することを目的としています。