[论文解读] A systematic benchmark of the ab initio Bethe-Salpeter equation approach for low-lying optical excitations of small organic molecules

Fabien Bruneval, Samia M. Hamed+1 more|arXiv (Cornell University)|May 18, 2015

Spectroscopy and Quantum Chemical Studies被引用 13

一句话总结

本研究使用Thiel的基准分子集，对28种小有机分子的低激发态光学激发进行从头算Bethe-Salpeter方程（BSE）基准测试。结果表明，BSE单重激发能对DFT平均场起始点高度敏感，其中含高比例精确交换的杂化泛函——尤其是此类泛函——能与高阶波函数方法实现定量一致；而三重激发态则因BSE核和GW准粒子能级中的系统性偏差而导致准确性较低。

ABSTRACT

The predictive power of the ab initio Bethe-Salpeter equation (BSE) approach, rigorously based on many-body Green's function theory but incorporating information from density functional theory, has already been demonstrated for the optical gaps and spectra of solid-state systems. Interest in photoactive hybrid organic/inorganic systems has recently increased, and so has the use of the BSE for computing neutral excitations of organic molecules. However, no systematic benchmarks of the BSE for neutral electronic excitations of organic molecules exist. Here, we study the performance of the BSE for the 28 small molecules in Thiel's widely-used time-dependent density functional theory benchmark set [M. Schreiber et al. J. Chem. Phys. 128, 134110 (2008)]. We observe that the BSE produces results that depend critically on the mean-field starting point employed in the perturbative approach. We find that this starting point dependence is mainly introduced through the quasiparticle energies obtained at the intermediate GW step, and that with a judicious choice of starting mean-field, singlet excitation energies obtained from BSE are in excellent quantitative agreement with higher-level wavefunction methods. The quality of the triplet excitations is slightly less satisfactory.

研究动机与目标

• 系统评估从头算Bethe-Salpeter方程（BSE）对小有机分子中性激发态的预测准确性。
• 通过分析对DFT平均场起始点的依赖性，识别BSE激发能中系统性误差的来源。
• 将BSE性能与更高阶波函数方法（如CCSD(T)）进行比较，评估其在单重和三重激发态中的可靠性。
• 确定BSE方法是否能在有限分子体系中实现光学带隙和激发能的定量准确性。
• 研究GW准粒子步骤和屏蔽库仑相互作用对小分子BSE结果的影响。

提出的方法

• 采用微扰GW信息框架实施BSE方法，其中准粒子能级在GW级别计算，电子-空穴相互作用通过屏蔽库仑相互作用W处理。
• 研究采用多种DFT起始点，包括局域泛函（PBE）和杂化泛函（B3LYP、BHLYP、tCAM-B3LYP），以评估其对BSE结果的影响。
• 通过Bethe-Salpeter方程形式化求解粒子-空穴通道中的两粒子格林函数，计算激发能。
• BSE核由GW自能和屏蔽库仑相互作用W构建，W采用随机相位近似（RPA）计算。
• 结果以Thiel基准集中28种小有机分子的高精度ΔSCF-CCSD(T)参考数据为基准进行验证。
• 通过将偏差分解为来自准粒子能级（GW）和BSE核的贡献，执行误差分析，以分离主要误差来源。

实验结果

研究问题

• RQ1 与高阶参考数据相比，BSE计算的小有机分子单重和三重激发能的准确性如何？
• RQ2 DFT平均场起始点的选择在多大程度上影响BSE对分子激发态的结果？
• RQ3 BSE激发能中的主要误差来源是什么——GW准粒子能级还是BSE核（W和电子-空穴相互作用）？
• RQ4 当使用最优起始泛函时，BSE能否实现与波函数方法相当的单重激发能定量准确性？
• RQ5 BSE中的系统性偏差如何影响单重-三重激发能级差的预测？与TD-DFT相比如何？

主要发现

• BSE方法计算的单重激发能对DFT平均场起始点高度敏感，其中局域泛函（如PBE）导致结果较差。
• 杂化泛函，特别是高精确交换含量的BHLYP和tCAM-B3LYP，能与ΔSCF-CCSD(T)参考数据实现极佳的定量一致。
• BSE结果对起始泛函的依赖性主要源于GW准粒子能级的变化，而非BSE核或屏蔽库仑相互作用W。
• 即使采用最优起始点，BSE仍存在残余误差：单重激发能被高估约0.2 eV，三重激发能被低估约0.4 eV。
• 由于单重激发能被高估且三重激发能被低估，BSE系统性地高估单重-三重能级差约0.6 eV，尽管该误差量级与最佳TD-DFT方案相当。
• 尽管单重激发能结果良好，但BSE准确性的物理基础并不完全可靠，其准确性依赖于误差抵消，与TD-DFT类似，而非源于HOMO-LUMO能隙的根本正确性。