[논문 리뷰] Accelerating Hartree--Fock exchange calculation using the TURBOMOLE program system: different techniques for different purposes

Arnim Hellweg, Uwe Huniar|arXiv (Cornell University)|2016. 10. 25.

Molecular spectroscopy and chirality인용 수 4

한 줄 요약

이 논문은 TURBOMOLE 프로그램 시스템에서 허트리-폭크 교환을 가속화하기 위해 RI-J(쿠론에 대한 해상도-항정), RI-K(교환에 대한 해상도-항정), 및 반수치적 K(seminumerical K, senex)라는 세 가지 계산 기법을 평가한다. rik(교환에 대한 RI)는 큰 기본 집합을 가진 소형에서 중형 분자에 가장 효율적이며, rij(쿠론에 대한 RI)는 대부분의 시스템에서 직접 SCF를 능가하지만, 매우 큰 시스템에서는 메모리 제약으로 인해 RI가 열세를 보인다.

ABSTRACT

Hartree--Fock theory is one of the most ancient methods of computational chemistry, but up to the present day quantum chemical calculations on Hartree--Fock level or with hybrid density functional theory can be excessively time consuming. We compare three currently available techniques to reduce the computational demands of such calculations in terms of timing and accuracy.

연구 동기 및 목표

• TURBOMOLE 프로그램 시스템에서 허트리-폭크 교환 계산을 가속화하기 위한 세 가지 현대적 기법—rij, rik, senex의 효율성과 정확도를 비교하기 위해.
• 다양한 분자 시스템과 기본 집합 크기에 따라 최적의 방법을 선택하는 데 혼란을 겪는 사용자들을 해결하기 위해.
• 시간, 정확도, 시스템 크기를 기반으로 한 방법 선택에 실용적인 지침을 제공하기 위해.
• 유기 사슬, 포도당, 안토니움 클러스터 등 다양한 시스템에서 이러한 방법의 성능을 벤치마킹하기 위해.
• 반직접적 및 잘라내기 RI 변형이 계산 비용과 수치 오차에 미치는 영향을 평가하기 위해.

제안 방법

• TURBOMOLE 프로그램 시스템을 활용하여 세 가지 기법—rij(쿠론 부분에 대한 RI), rik(교환 부분에 대한 RI), senex(교환 적분의 반수치적 평가)—를 구현하고 성능을 벤치마킹한다.
• 적분 스크리닝과 해상도-항정(RI) 근사법을 적용하여 이차전자 적분의 형식적 스케일링을 N⁴에서 N² 이하로 감소시킨다.
• 다중극 가속 RI-J(MARI-J)를 사용하여 쿠론 포텐셜 평가를 효율적으로 수행하고 스케일링을 N으로 감소시킨다.
• 교환 부분은 정확한 스크리닝(rik) 또는 공간 격자 위에서의 반수치적 적분(sexex)을 통해 평가하며, 후자의 경우 표준 DFT 격자를 사용한다.
• 성능와 정확도를 비교하기 위해 직접 SCF(dscf), 반직접적 SCF(적분을 디스크에 저장), 잘라내기 RI(감소된 보조 기본 집합) 변형을 비교한다.
• def-SV(P)에서 QZVPP에 이르기까지 다양한 기본 집합을 사용하여 아밀로스 사슬에서 Asₙ 클러스터(n=4,8,12)에 이르기까지 다양한 시스템에서 성능 벤치마크 계산을 수행한다.

실험 결과

연구 질문

• RQ1 다양한 시스템 크기와 기본 집합 유형에서 허트리-폭크 교환 계산에 대해 rij, rik, senex 중 어느 기법이 가장 효율적인가?
• RQ2 큰 산란 기본 집합을 가진 분자에서 rik의 성능은 직접 SCF(dscf) 및 rij와 비교해 볼 때 어떻게 되는가?
• RQ3 반직접적 및 잘라내기 RI 변형이 계산 속도와 수치 정확도에 어떤 영향을 미치는가?
• RQ4 senex는 큰 기본 집합을 가진 하이브리드-DFT 또는 후-HF 계산에서 속도와 정확도 사이의 타협 가능성을 제공할 수 있는가?
• RQ5 RI 방법의 메모리 오버헤드가 성능 이득을 상쇄하기 시작하는 시스템 크기는 어느 정도인가?

주요 결과

• rij는 테스트된 모든 시스템에서 직접 SCF(dscf)를 항상 뛰어넘는 성능을 보이며, 중형 분자에서 가장 큰 속도 향상을 보였다.
• rik는 소형에서 중형 분자에 큰 기본 집합(예: def2-TZVPPD)을 사용할 경우 가장 효율적인 방법으로, 포도당의 MP2 계산에서 HF 계산 시간을 최대 98%까지 감소시켰다.
• 큰 기본 집합을 가진 안토니움 클러스터(예: QZVPP)에서는 rik가 여전히 dscf보다 빠르지만, 더 작은 시스템에 비해 성능 이점이 감소한다.
• senex는 일반적인 DFT에서 작은 기본 집합을 사용할 경우 너무 느리고, 고정밀도 파동함수 이론에 적합하지 않아 정확도가 떨어지지만, 이전에는 불가능했던 두 성분 상대론적 하이브리드-DFT 계산을 가능하게 한다.
• 반직접 dscf는 소형 시스템에서만 직접 dscf보다 빠르며, 더 큰 시스템에서는 더 느리고 추가적인 I/O 오버헤드를 유발한다.
• 잘라내기 rik는 전체 rik보다 빠르지만 성능 향상은 미미하고 오차가 증가하므로 고정밀도 파동함수 이론에는 부적합하다.