[论文解读] A New View on Density Corrected DFT: Can One Get a Better Answer for a Good Reason?

Devin J. Hernandez, Adam Rettig+1 more|arXiv (Cornell University)|Jun 26, 2023

Catalysis and Oxidation Reactions被引用 10

一句话总结

本文分析使用 kappa-OOMP2 提供的改进密度的密度校正 DFT（DC-DFT），结果表明误差抵消常常推动 HF-DFT 的成功，而更高质量的密度反而会削弱低阶泛函，对 SIE 和导数不连续性有影响。

ABSTRACT

Despite its widespread use, density functional theory (DFT) has several notable areas of failure; perhaps the most well-studied of these failures is self-interaction error (SIE). Density corrected DFT (DC-DFT) was proposed as a potential solution to systems where SIE causes traditional DFT to fail. The Hartree-Fock (HF) density is then used for cases where the DFT energy is suitable but the self-consistent density is erroneous. In this study, we investigate the utility of the higher quality orbital optimized MP2 densities in DC-DFT for barrier heights and halogen bonded complexes. For functionals such as PBE and r$^2$SCAN, find that these densities yield worse results than the HF density due to favorable cancellation between the density-driven and functional-driven errors, confirming a recent study. Error decomposition reveals functional driven error, not density driven error, to be the primary cause of inaccuracy in DFT calculations where SIE is prominent. We therefore advise caution when using HF-DFT, because the only rigorous way to remove large functional-driven errors in lower rungs of Jacob's ladder is by climbing to higher rungs that include exact exchange. We recommend that better functionals be improved by using a better density in SIE-sensitive cases. Examples support the value of this variant of DC-DFT. We also emphasize that DC-DFT potential energy surfaces have first derivative discontinuities at Coulson-Fischer points, in contrast to the second derivative discontinuities in SCF solutions.

研究动机与目标

• 评估在 SIE 敏感问题中使用 HF 和 kappa-OOMP2 密度的 DC-DFT 的有效性。
• 将总误差分解为密度驱动和泛函驱动分量，以识别不准确的来源。
• 评估势垒高度、卤键结合和断键问题，以了解更高质量的密度何时改善或降低结果。
• 研究 DC-DFT 势能面在 Coulson-Fischer 点处的导数行为和不连续性。

提出的方法

• 使用三种密度进行单点能量计算：自洽本征密度、Hartree-Fock (HF) 密度，以及 κ-OOMP2（κ-OOMP2）密度。
• 在 BH76、SIE4x4、Bauzá 卤键集合以及 SIE 敏感案例等数据集上，使用一系列泛函（SPW92、PBE、r2SCAN、PBE0、B3LYP、ωB97X-V、ωB97M-V）评估能量。
• 将 κ-OOMP2 作为“更好”密度的代理，将误差分解为密度驱动和泛函驱动贡献。
• 比较 DC-DFT 变体（基于 HF 和基于 κ-OOMP2）的结果与自洽 DFT，以评估误差抵消和性能。
• 分析偶极矩和势能面，以说明 DC-DFT 中非亥姆霍兹-费米贡献。

实验结果

研究问题

• RQ1 在相对于 HF 密度，使用 κ-OOMP2 密度的 DC-DFT 在 SIE 敏感问题上是改善还是恶化准确性？
• RQ2 在不同泛函下，HF-DFT 的密度驱动与泛函驱动误差在多大程度上互相抵消？
• RQ3 DC-DFT 变体在 BH76 势垒、高级泛函的 Bauzá 卤键集以及 SIE4x4 断键问题上的表现如何？
• RQ4 与 SCF 解相比，Coulson-Fischer 点处的 DC-DFT 势能面导数不连续性有哪些？
• RQ5 κ-OOMP2 密度能否提供一种可靠、经济的近似真密度用于分析误差成分？

主要发现

• HF-DFT 常因误差抵消带来对势垒高度的优于自洽 DFT 的表现，尤其在低阶泛函中更为显著。
• κ-OOMP2 密度在 NSP 情况下显著改进偶极矩，相较于 HF 和自洽 DFT；但对低阶泛函 HF-DFT 的性能造成退化，因为失去了误差抵消。
• 对于 BH76 势垒，HF-DFT 在 PBE 和 r2SCAN 上显著降低 MAE，而 κ-OOMP2-DFT 可能恶化或仅略有改善，表明密度驱动误差并非唯一的不准确来源。
• 分解结果显示在 SIE 敏感问题中，泛函驱动误差占主导地位；HF 的密度驱动误差可能很大，高阶泛函可以缓解但无法消除。
• 卤键结合（Bauzá 集）显示 HF-DFT 在较高阶泛函上通常与自洽 DFT 相当，κ-OOMP2-DFT 仅提供温和提升，表明在这些情况下对误差抵消的依赖较少。
• SIE4x4 的断键问题显示出较大的泛函误差，HF-DFT 通过密度驱动纠正在一定程度上缓解，而 κ-OOMP2 密度提供的改进有限，凸显 dc-DFT 在极端 SIE 情况下的局限性。