[논문 리뷰] DDEC6: A Method for Computing Even-Tempered Net Atomic Charges in Periodic and Nonperiodic Materials

Thomas A. Manz, Nidia Gabaldon Limas|arXiv (Cornell University)|2015. 12. 24.

Advanced Chemical Physics Studies참고 문헌 5인용 수 13

한 줄 요약

DDEC6는 전체 전자 밀도를 제약 최적화 접근법을 사용해 분할하여, 주기적 및 비주기적 물질에서 화학적으로 일관되고 이 trasferable하며 기저함수 집합에 영향을 받지 않는 네트워크 원자 전하(NAC)를 계산하는 데 사용되는 새로운 방법이다. 선형 스케일링, 강건한 수렴성, 그리고 다양한 물질에서 전자 이동 경향, 코어 수준 이동, 전기쌍극모멘트를 정확히 재현함으로써 기존의 DDEC 및 히르슈펠드 방법보다 안정성과 정확도에서 뛰어나다.

ABSTRACT

Net atomic charges (NACs) are widely used in all chemical sciences to concisely summarize key information about the partitioning of electrons among atoms in materials. Although widely used, there is currently no atomic population analysis method suitable for being used as a default method in quantum chemistry programs. To address this challenge, we introduce a new atoms-in-materials method with the following nine properties: (1) exactly one electron distribution is assigned to each atom, (2) core electrons are assigned to the correct host atom, (3) NACs are formally independent of the basis set type because they are functionals of the total electron distribution, (4) the assigned atomic electron distributions give an efficiently converging polyatomic multipole expansion, (5) the assigned NACs usually follow Pauling scale electronegativity trends, (6) NACs for a particular element have good transferability among different conformations that are equivalently bonded, (7) the assigned NACs are chemically consistent with the assigned atomic spin moments, (8) the method has predictably rapid and robust convergence to a unique solution, and (9) the computational cost of charge partitioning scales linearly with increasing system size. Across a broad range of material types, the DDEC6 NACs reproduced electron transfer trends, core electron binding energy shift trends, and electrostatic potentials across multiple system conformations with excellent accuracy compared to other charge assignment methods. Due to non-nuclear attractors, Bader's quantum chemical topology could not assign NACs for some of these materials. The DDEC6 method alleviates the bifurcation or runaway charges problem exhibited by earlier DDEC variants and the Iterative Hirshfeld method. These characteristics make the DDEC6 method ideally suited for use as a default charge assignment method in quantum chemistry programs.

연구 동기 및 목표

• 양자화학 소프트웨어 프로그램에서 신뢰성 있고 이 trasferable하며 물리적으로 의미 있는 기본 원자 인구 분석 방법을 개발하는 것.
• 이전의 DDEC 및 반복 히르슈펠드 방법에서 나타나는 분기 문제와 과도한 전하 문제를 해결하는 것.
• 기저함수 집합의 선택과 무관하게 네트워크 원자 전하를 확보하기 위해 총 전자 밀도에서 직접 유도하는 것.
• 네트워크 원자 전하가 폴링 전기음성도 순서를 따르며 원자 스핀 모멘트와 화학적으로 일관되게 유지되도록 하는 것.
• 대규모 시스템, 특히 주기적 물질을 포함한 시스템에서 선형 계산 스케일링과 빠르고 유일한 수렴성을 달성하는 것.

제안 방법

• DDEC6는 물리적 일관성을 강제하는 제약 최적화 프레임워크를 사용해 총 전자 밀도를 원자 중심 분포로 분할한다.
• 이 방법은 각 원자당 정확히 한 개의 전자 분포를 할당하며, 코어 전자는 해당 원자에 정확히 할당된다.
• 네트워크 원자 전하는 총 전자 분포의 함수로 계산되어 기저함수 집합에 독립적이게 된다.
• 매우 매끄럽고 잘 국소화된 원자 전하 분포 덕분에 다원자 다극모멘트 전개가 효율적으로 수렴된다.
• 알고리즘은 네트워크 원자 전하와 스핀 모멘트 사이의 화학적 일관성을 강제하여 물리적으로 비합리적인 전하-스핀 분리 문제를 방지한다.
• 볼록 최적화 공식화 덕분에 수렴이 보장되고 빠르며, 계산 비용은 시스템 크기와 선형으로 증가한다.

실험 결과

연구 질문

• RQ1 다양한 물질에서 수렴성, 유일성, 기저함수 집합 선택에 영향을 받지 않는 전하 분할 방법을 개발할 수 있는가?
• RQ2 DDEC 프레임워크는 이전 버전에서 관찰된 분기 문제와 과도한 전하 문제를 어떻게 제거할 수 있는가?
• RQ3 DDEC6로 계산된 네트워크 원자 전하는 전자 이동 및 코어 수준 이동의 실험적 경향을 어느 정도 정확히 재현하는가?
• RQ4 DDEC6는 다양한 분자 구형과 주기적 시스템에서 이 trasferable하고 화학적으로 일관된 전하를 생성할 수 있는가?
• RQ5 비핵적 최소점이 존재하는 시스템을 포함한 대규모 또는 주기적 시스템에서도 이 방법은 선형 스케일링과 빠른 수렴성을 유지하는가?

주요 결과

• DDEC6는 네트워크 원자 전하가 기저함수 유형과 무관하게 공식적으로 독립적이며, 총 전자 밀도에서 직접 유도된다는 점에서 유의미하다.
• 이전의 DDEC 및 히르슈펠드 방법에서 나타나는 수렴 불안정성과는 달리, DDEC6는 빠르고 강건한 유일한 해로 수렴한다.
• DDEC6 전하는 동일한 결합 환경을 가진 다양한 분자 구형에서 뛰어난 이 trasferability를 보인다.
• 계산된 네트워크 원자 전하는 폴링 전기음성도 순서를 따르며 원자 스핀 모멘트와 화학적으로 일관된다.
• DDEC6는 다양한 물질에서 전자 이동 및 코어 수준 결합 에너지 이동의 실험적 경향을 정확히 재현한다.
• Bader의 양자화학적 위상 구조가 실패하는 비핵적 최소점이 존재하는 시스템에서도 DDEC6는 네트워크 원자 전하를 성공적으로 할당한다.