[论文解读] Non-Koopmans Corrections in Density-functional Theory: Self-interaction Revisited

Ismaïla Dabo, Matteo Cococcioni+1 more|ArXiv.org|Jan 18, 2009

Advanced Chemical Physics Studies参考文献 1被引用 5

一句话总结

本文提出了一种无参数的非-Koopmans自相互作用校正（NK-SIC），相较于Perdew-Zunger自相互作用校正（PZ-SIC），在密度泛函理论中实现了更优的自相互作用误差抵消。该方法在预测电子性质方面表现出更优的准确性，特别是在氢链的电极化率预测上，优于标准DFT和PZ-SIC，其与CCSD(T)的平均绝对偏差降低至2.1 a.u.，而PZ-SIC为13.9 a.u。

ABSTRACT

In effective single-electron theories, self-interaction manifests itself through the unphysical dependence of the energy of an electronic state as a function of its occupation, which results in important deviations from the ideal Koopmans trend and strongly affects the accuracy of electronic-structure predictions. Here, we study the non-Koopmans behavior of local and semilocal density-functional theory (DFT) total energy methods as a means to quantify and to correct self-interaction errors. We introduce a non-Koopmans self-interaction correction that generalizes the Perdew-Zunger scheme, and demonstrate its considerably improved performance in correcting the deficiencies of DFT approximations for self-interaction problems of fundamental and practical relevance.

研究动机与目标

• 为解决局域和半局域DFT泛函中长期存在的自相互作用误差问题，该误差会扭曲电子能级和电离势。
• 量化并校正Kohn-Sham DFT中的非-Koopmans行为，其中本征值因非物理的自相互作用而偏离物理电离能。
• 通过改进自相互作用贡献的抵消方式，开发出优于Perdew-Zunger方案的更精确自相互作用校正方法。
• 在基准体系（包括H₂⁺、H₂和扩展的氢链）上验证新方法，以评估其在复杂电子响应性质预测中的能力。

提出的方法

• 提出一种非-Koopmans自相互作用校正（NK-SIC），通过改进Kohn-Sham势中自相互作用项的处理方式，对Perdew-Zunger方案进行推广。
• 将非-Koopmans贡献定义为 $ \Pi_{i\alpha} = \epsilon_{i\alpha} + \Delta E^{(0)}_{i\alpha} $，其中 $ \epsilon_{i\alpha} $ 为Kohn-Sham本征值，$ \Delta E^{(0)}_{i\alpha} $ 为非弛豫电离能，以量化自相互作用误差。
• 在Quantum-Espresso程序包中，采用预处理共轭梯度逐带算法，实现NK-SIC的自洽计算。
• 通过外加电场（0.005 a.u.）的有限差分响应方法评估纵向电极化率，以偶极矩变化作为电子离域化的探测手段。
• 在单电子和双电子体系（H⁺、H₂⁺、H₂）上验证该方法，并将其扩展至交替键长为2和3 bohr的周期性氢链。
• 采用CCSD(T)作为极化率计算的参考标准，对比HF、MP4、LSD、PZ-SIC和NK-SIC的结果。

实验结果

研究问题

• RQ1 与原始的Perdew-Zunger SIC相比，NK-SIC在Kohn-Sham DFT中如何改善自相互作用误差的抵消？
• RQ2 NK-SIC在分数电子占据体系中，对Kohn-Sham本征值的非-Koopmans行为校正程度如何？
• RQ3 NK-SIC能否准确描述扩展氢链的电极化率？该体系是电子结构理论中自相互作用误差的基准测试体系。
• RQ4 在预测氢二聚体电极化率方面，NK-SIC与后-Hartree-Fock方法（如CCSD(T)和MP4）相比表现如何？
• RQ5 NK-SIC是否在无需对Hartree自相互作用能进行经验缩放的情况下，实现了优于HF和PZ-SIC的准确性？

主要发现

• NK-SIC将氢链纵向电极化率的平均绝对偏差从LSD的39.9 a.u.降低至相对于CCSD(T)的2.1 a.u.，显著优于PZ-SIC的13.9 a.u.。
• 对于含有N个二聚体的氢链，NK-SIC预测的电极化率 $ \alpha_N $ 呈线性趋势，且在N = 2至6范围内始终与CCSD(T)保持在2.1 a.u.以内，表明其具有极高的预测准确性。
• NK-SIC的准确性优于HF，平均绝对偏差仅为2.1 a.u.，而HF为8.0 a.u.，尽管HF本身是无自相互作用的。
• 该方法在单电子密度的二阶展开中正确捕捉了Koopmans条件，表明其有效抑制了自相互作用误差。
• 在H₂⁺和H₂的离解过程中，NK-SIC与精确解和无自相互作用结果高度一致，证实其在基本电子结构问题中的鲁棒性。
• 与以往的SIC改进方法不同，NK-SIC为无参数方法，无需对Hartree自相互作用能进行经验下调，提供了更严格且系统化的校正。