[논문 리뷰] DC-DFT for Open Shells: How to Deal with Spin Contamination

Hayoung Yu, Suhwan Song외 3명|arXiv (Cornell University)|2023. 07. 20.

Advanced Chemical Physics StudiesPhysics and Astronomy인용 수 3

한 줄 요약

이 논문은 개방계 시스템에서 밀도보정된 DFT(DC-DFT)에 응용된 비제한형 하트리-폭크(UHF) 밀도에서 스핀 오염을 조사한다. UHF 밀도 대신 제한형 개방계 HF(ROHF) 밀도를 사용할 경우 평균 절대 오차가 크게 감소함을 보여주며, 특히 스핀 오염이 있는 경우에 효과적이다. 이는 GGA, 메타-GGA, 하이브리드 함수 모두에서 정확도를 향상시키며, DC-DFT 함수의 오차는 최대 2배까지 감소한다.

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.

연구 동기 및 목표

• 개방계 시스템에서 밀도보정된 DFT(DC-DFT)에 적용된 비제한형 하트리-폭크(UHF) 밀도에서 스핀 오염 문제를 해결하기 위해.
• 스핀 오염이 존재할 경우, 제한형 개방계 HF(ROHF) 밀도가 UHF 밀도보다 더 정확한 에너지 계산을 제공하는지 평가하기 위해.
• 스핀 오염이 있는 경우와 없는 경우에서 다양한 밀도 기능형 근사(DFA)가 ROHF 밀도와 UHF 밀도와 조합되었을 때의 성능을 평가하기 위해.
• ROHF 밀도가 더 우수한 결과를 얻는 조건을 규명하여, 개방계 시스템에서 DC-DFT를 적용하는 데 실용적인 지침을 제공하기 위해.
• 자기모멘트가 크게 오염된 라디칼 및 전이 금속 복합체와 같은 도전적인 개방계 시스템에 대해 HF-DFT 및 DC-DFT의 적용 범위를 확장하기 위해.

제안 방법

• 연구는 GMTKN55 데이터베이스를 사용하여 자성수렴(scf), UHF, ROHF 밀도를 평가하였으며, 반응을 스핀 오염(⟨Ŝ²⟩ > 0.1) 또는 스핀 무오염으로 분류하였다.
• DC-DFT 프레임워크를 적용하여 총 오차를 기능 오차(ΔEF)와 밀도 기반 오차(ΔED)로 분해함으로써 밀도 품질의 영향을 고립시켰다.
• 각 기능형에 대해 개방계 서브셋에 대한 sc-, UHF-, ROHF-DFT 결과의 가중 총 평균 절대편차(WTMAD-2)를 계산하였다.
• 밀도 민감도 임계값 2 kcal/mol(Shim 등에 기반)을 사용하여 밀도 민감성 있는 경우와 없는 경우를 구분하였다.
• 13종의 기능형을 비교하였으며, 이는 4종의 GGA, 4종의 메타-GGA, 5종의 하이브리드, 2종의 순수 HF-DFT 기능형(HF-r²SCAN-DC4 및 BL1p)을 포함한다.
• 모든 계산은 Ahlrichs의 def2-QZVPPD basis set를 사용하였으며, CUHF를 위한 커스텀 파이썬 코드를 포함한 Chemistry Simulation Framework를 활용하였다.

실험 결과

연구 질문

• RQ1 DC-DFT에서 UHF 밀도를 ROHF 밀도로 대체할 경우, 심한 스핀 오염이 있는 개방계 시스템에서 오차가 감소하는가?
• RQ2 스핀 오염이 있는 경우, 다양한 유형의 기능형(GGA, 메타-GGA, 하이브리드, DC-DFT)에서 UHF 밀도와 ROHF 밀도의 성능은 어떻게 다를까?
• RQ3 UHF 밀도를 사용할 경우 스핀 오염이 DC-DFT의 정확도를 얼마나 심각하게 떨어뜨리며, ROHF 밀도로 대체하면 이를 완화할 수 있는가?
• RQ4 스핀 오염이 있는 시스템에서 ROHF 밀도를 UHF 밀도로 대체할 경우, 기능형에 따라 오차 감소 추세가 존재하는가?
• RQ5 스핀 오염이 존재하더라도, DC-DFT에서 ROHF 밀도 사용이 개방계 시스템에 대해 신뢰성 있게 권장될 수 있는가?

주요 결과

• 스핀 오염이 있는 경우, ROHF 밀도를 사용하면 UHF 밀도 대비 WTMAD-2 오차가 최대 2배까지 감소하였으며, 특히 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 vs. 5.20 kcal/mol), 이는 비오염 시스템에서 ROHF를 사용하는 데 거의 손해가 없음을 시사한다.
• 이중 하이브리드 기능형 BL1p의 경우 UHF에서 스핀 오염으로 인한 오차가 가장 심각했으나, ROHF-DFT로 전환하면 오차가 25.42에서 5.40 kcal/mol으로 크게 감소하였다.
• 연구는 UHF 밀도에서의 스핀 오염이 큰 밀도 기반 오차(ΔED)를 유발하며, 이를 ROHF 밀도로 대체함으로써 개방계 시스템에서 이러한 오차를 효과적으로 억제할 수 있음을 확인하였다.