[논문 리뷰] A New View on Density Corrected DFT: Can One Get a Better Answer for a Good Reason?

Devin J. Hernandez, Adam Rettig외 1명|arXiv (Cornell University)|2023. 06. 26.

Catalysis and Oxidation Reactions인용 수 10

한 줄 요약

본 논문은 kappa-OOMP2로부터 얻은 향상된 밀도를 사용한 density corrected 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.

연구 동기 및 목표

• DC-DFT의 HF 및 kappa-OOMP2 밀도를 사용하여 SIE에 민감한 문제들에서 효과를 평가한다.
• 총 오차를 밀도 주도 오차와 함수주도 오차로 분해하여 부정확성의 원인을 식별한다.
• 장벽 높이, 할로겐 결합, 결합 파괴 문제를 평가하여 더 나은 밀도가 결과를 향상시키는지 악화시키는지 이해한다.
• Coulson-Fischer 지점에서 DC-DFT 포텐셜 에너지 표면의 도함수 거동과 불연속성을 조사한다.

제안 방법

• 세 가지 밀도(자체수렴 원래 밀도, HF 밀도, κ-OOMP2 밀도)를 사용하여 단일점 에너지를 계산한다.
• 다양한 밀도함수(FSPW92, PBE, r2SCAN, PBE0, B3LYP, ωB97X-V, ωB97M-V)를 BH76, SIE4x4, Bauzá 할로겐 결합 집합, 및 SIE-민감 케이스 데이터세트에서 평가한다.
• κ-OOMP2를 “더 나은” 밀도의 대리로 사용하고, 오차를 밀도 주도 및 함수주도 기여로 분해한다.
• DC-DFT 변형(HF 기반 및 κ-OOMP2 기반)을 자기일치 DFT와 비교하여 오차 상쇄 및 성능을 평가한다.
• 다이폴 모멘트와 포텐셜 에너지 표면을 분석하여 DC-DFT에서의 비 헬만-파인만(Non-Hellmann-Feynman) 기여를 설명한다.

실험 결과

연구 질문

• RQ1 DC-DFT에서 κ-OOMP2 밀도 사용이 HF 밀도에 비해 SIE에 민감한 문제에서 정확성을 향상시키는가 아니면 악화시키는가?
• RQ2 다양한 함수들에 대해 HF-DFT에서 밀도 주도 오차와 함수 주도 오차가 어느 정도 상쇄되는가?
• RQ3 DC-DFT 변형은 BH76의 장벽 높이, Bauzá 세트의 할로겐 결합, SIE4x4의 결합 파괴 문제에서 어떻게 성능을 보이는가?
• RQ4 Coulson-Fischer 지점에서 DC-DFT 포텐셜 에너지 표면의 도함수 불연속성은 SCF 해와 비교해 어떤 차이가 있는가?
• 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 밀도는 제한된 개선만 제공하며, 극단적 SIE 사례에서 dc-DFT의 한계를 강조한다.