[論文レビュー] Correcting hybrid density functionals to model Y6 and other non-fullerene acceptors

Tom Ward, Isabel Creed+2 more|arXiv (Cornell University)|Mar 5, 2026

Organic Electronics and Photovoltaics被引用数 0

ひとこと要約

この論文はY6および関連NFAsを正確にモデリングするためにスクリーン付きレンジ分離ハイブリッドを調整し、レンジ分離パラメータの減少がオフ・ザ・シェルフの機能を上回る励起状態予測を改善することを示している。

ABSTRACT

Recently developed fused-ring organic electron-acceptors such as Y6 have strong oscillator strength, good charge-carrier transport and low bandgaps. They therefore have enormous current technical application to optoelectronic devices, such as solar cells. Due to the large number of atoms involved in representative aggregates of these materials, we need an efficient electronic structure method to model them. Standard density functional theory poorly describe charge-transfer states, and were developed for vacuum calculations of individual molecules. In this work we tune a range-separated hybrid functional for Y6. We characterise representative dimers of the solid-state and show that Y6 dimers show the extensive solvatochromic effects are due, in part, to oscillator strength borrowing. We provide an explanation for the short optimally-tuned range-separation parameter, based in the Penn model for the frequency dependent dielectric of a semiconductor. We caution that standard range-separated hybrids are less accurate than global hybrids for these, and similar, materials. We show how reducing the range-separation length improves the accuracy of standard functionals, without an involved tuning process.

研究の動機と目的

• 低帯域ギャップ・高振動強度NFAs（例：Y6）の正確な量子化学モデリングを動機づける。
• Y6二量体のための最適に調整されたスクリーン付きレンジ分離ハイブリッド（OT-SRSH）を開発・検証する。
• OT-SRSHを標準的なグローバルハイブリッドおよび他のSRSH機能と比較する。
• 標準のレンジ分離ハイブリッドが不十分である理由を説明し、ペンモデルを介した誘電体スクリーニングへの調整と結びつける。
• モノマー、二量体、および小さな集合体間で調整パラメータの移行性を実証する。

提案手法

• Koopmansベースの損失関数 J(ω) によって HOMO/LUMO エネルギーと IP/EA を一致させるように LC-ωhPBE パラメータ（ω, α, β）を最適化する。
• α+β = 1/εr を制約として εr ≈ 6 の固体状態スクリーニングを反映させる。
• PCMを用いて誘電環境で調整し、モノマーと二量体間の移植性を評価する。
• TheoDOREで励起状態の特性を分析し、CT 対 FE の含有量（ωCT）と局在化（ωPR）に焦点を当てる。
• GW/BSE/MM結果や他のTDDFT機能（B3LYP, PBE0, CAM-B3LYP, wB97XD）と比較してベンチマークを取る。
• Y6モノマーと結晶構造から抽出した接触ペア二量体に対してGaussian16（6-31G(d,p)）でTDDFT計算を実施する。

実験結果

リサーチクエスチョン

• RQ1 最適に調整されたスクリーン付きレンジ分離ハイブリッドはY6二量体のCTおよびFE励起状態を正確に記述できるか。
• RQ2 OT-SRSHとオフザシェルのSRSHおよびグローバルハイブリッドは、Y6接触ペア二量体の状態の順序とエネルギー予測においてどう異なるか。
• RQ3 調整された機能はモノマーと二量体の間で一貫して移行し、集合体のスケーラブルなモデリングを可能にするか。
• RQ4 誘電モデルに導かれたレンジ分離パラメータの縮小は、フルOT-SRSHで予測される励起を再現できるか。
• RQ5 CT-FE混成と単結・三重項の順序付けは、固体状態の集合体におけるY6の光物理にどのような影響を及ぼすか。

主な発見

Table 1: Dimer	ω (a0^-1)	α	β
Default	N/A	0.40	1
D1 (Contact)	0.120	0.230	-0.063
D2 (Contact)	0.110	0.250	-0.083
D3 (Contact)	0.110	0.240	-0.073
D4 (Contact)	0.120	0.250	-0.083
D5 (Contact)	0.120	0.240	-0.073
D6 (Contact)	0.120	0.240	-0.073
D7 (Contact)	0.130	0.250	-0.083
D8 (Contact)	0.130	0.250	-0.083
D9 (Non-Contact)	0.130	0.250	-0.083
D10 (Non-Contact)	0.130	0.250	-0.083
Monomer	N/A	0.140	0.260	-0.093

• OT-SRSHはすべての接触ペア二量体（D1–D6）における最初の6つの励起状態の平均二乗誤差を最小化する。
• 非調整のレンジ分離ハイブリッドはCT含有励起に対してOT-SRSHおよびグローバルハイブリッドと比較して最大の誤差を生む。
• ωを低下させたCAM-B3LYPの調整はOT-SRSHの結果を再現する。範囲分離長さを短くすることで二量体光物理の本質を捉えられるが、完全調整は不要。
• 同じ調整パラメータはモノマーと二量体計算の間で移行可能で、励起状態特性は一貫して得られる。
• 三つの二量体クラスが出現：D1/D2はFE状態が近くCT様の最小対称性を持つCT様の最低単量体を示す；D4は混成CT-FEの単量体を示し、三重項はFE支配で局在化が変動する。
• 単結状態におけるCT-FE混成はCT状態の安定化をもたらす一方、三重項CT状態はこの混成によって不安定化し、観察された singlet-triplet の順序付けを説明する。
• CT-FE混成は固体状態でのsinglet状態のシフトに強く影響を与え、Y6接触ペア二量体で予測されるsinglet-triplet反転を支持する。
• Pennモデル由来のω ≈ 0.6 q_s および ω ≈ 0.6 sqrt(2 μ E_g) の関係は、低帯域ギャップNFAs（例：Y6、εr ~ 6、E_g ~ 1.3 eV）に対して小さな ω を選択する物理的根拠を提供する。
• ωを0.16 a0^-1へと調整したCAM-B3LYPはOT-SRSHとexcellent agreementを示し、これらの系ではωを下げるだけで正確な記述を得られることを示している。