[论文解读] DC-DFT for Open Shells: How to Deal with Spin Contamination

Hayoung Yu, Suhwan Song+3 more|arXiv (Cornell University)|Jul 20, 2023

Advanced Chemical Physics StudiesPhysics and Astronomy被引用 3

一句话总结

本文研究了在开放壳体系的密度修正密度泛函理论（DC-DFT）中，无限制哈特ree-fock（UHF）密度的自旋污染问题。结果表明，使用限制性开壳体系HF（ROHF）密度替代UHF密度可显著降低平均绝对误差，尤其在自旋污染情况下，该方法在GGA、杂化和meta-GGA泛函中均提升了精度，DC-DFT泛函的误差最高可降低至二分之一。

ABSTRACT

Density functional theory (DFT) is widely used to predict chemical properties, but its accuracy is limited by functional approximations and their approximate self-consistent densities. Density-corrected DFT (DC-DFT) is the study of the errors due to densities and Hartree-Fock DFT (HF-DFT) uses HF densities to improve energetics. With increasing use of HF-DFT, the question of how to address strong spin contamination in the HF calculation becomes increasingly important. We compare two different open-shell HF densities across 13 different DFT functionals and two DC-DFT methods. For significant spin contamination, ROHF densities outperform UHF densities by as much as a factor of 3, depending on the energy functional, and ROHF-DFT improves over self-consistent DFT for most of the tested functionals. We refine the DC(HF)-DFT algorithm, recommending ROHF-DFT in cases of severe spin contamination.

研究动机与目标

• 解决在开放壳体系中应用密度修正密度泛函理论（DC-DFT）时，无限制哈特ree-fock（UHF）密度的自旋污染问题。
• 评估在存在自旋污染的情况下，限制性开壳体系HF（ROHF）密度是否比UHF密度提供更精确的能级计算结果。
• 评估不同密度泛函近似（DFA）在自旋污染和自旋未污染情况下的表现，当与ROHF和UHF密度配对时。
• 通过识别ROHF密度在何种条件下表现更优，为开放壳体系中DC-DFT的实际应用提供指导。
• 将HF-DFT和DC-DFT的适用范围扩展至具有显著自旋污染的复杂开放壳体系，如自由基和过渡金属配合物。

提出的方法

• 本研究基于GMTKN55数据库评估自洽（sc）、UHF和ROHF密度，将反应分类为自旋污染（Δ⟨Ŝ²⟩ > 0.1）或自旋未污染。
• 在DC-DFT框架下，将总误差分解为泛函误差（ΔEF）和密度驱动误差（ΔED），以隔离密度质量的影响。
• 对每种泛函，计算开放壳体系子集的加权总平均绝对偏差（WTMAD-2），分别针对sc-、UHF-和ROHF-DFT结果。
• 采用2 kcal/mol的密度敏感性阈值（基于Sim等人）来区分密度敏感与不敏感的情况。
• 比较13种泛函的表现：4种GGA、4种meta-GGA、5种杂化泛函，以及2种全HF-DFT泛函（HF-r²SCAN-DC4和BL1p）。
• 所有计算均使用Ahlrichs def2-QZVPPD基组，并基于Simulations of Chemistry Framework，使用自定义Python代码实现CUHF计算。

实验结果

研究问题

• RQ1 在存在显著自旋污染的开放壳体系中，将DC-DFT中的UHF密度替换为ROHF密度是否能降低误差？
• RQ2 在自旋污染情况下，UHF与ROHF密度在不同泛函类型（GGA、meta-GGA、杂化、DC-DFT）中的表现差异如何？
• RQ3 当使用UHF密度时，自旋污染在多大程度上降低了DC-DFT的精度？ROHF密度能否有效缓解这一问题？
• RQ4 在自旋污染体系中，使用ROHF密度替代UHF密度时，是否存在与泛函类型相关的误差降低趋势？
• RQ5 即使存在自旋污染，是否可以可靠地推荐在DC-DFT中使用ROHF密度于开放壳体系？

主要发现

• 在自旋污染情况下，使用ROHF密度使WTMAD-2误差相比UHF密度最高降低至二分之一，尤其在HF-r²SCAN-DC4等DC-DFT泛函中表现显著。
• 对于UBL1p，自旋污染情况下的WTMAD-2误差达25.42 kcal/mol，而使用ROBL1p后降至5.40 kcal/mol，显示出显著改善。
• 对于杂化泛函PBE0，UHF-DFT在自旋污染情况下的WTMAD-2为13.07 kcal/mol，而ROHF-DFT将其降低至5.24 kcal/mol。
• 在自旋未污染情况下，ROHF-DFT的误差仅略高于UHF-DFT（例如，PBE泛函下分别为6.13与5.20 kcal/mol），表明在非污染体系中使用ROHF的代价极小。
• 双杂化泛函BL1p在UHF中受自旋污染影响最严重，但ROHF-DFT将其自旋污染误差从25.42 kcal/mol降至5.40 kcal/mol。
• 本研究证实，UHF密度中的自旋污染会导致较大的密度驱动误差（ΔED），而用ROHF密度替代UHF密度可有效抑制开放壳体系中的此类误差。