[논문 리뷰] Scalable and predictive spectra of correlated molecules with moment truncated iterated perturbation theory

Oliver J. Backhouse, Alejandro Santana‐Bonilla외 1명|arXiv (Cornell University)|2021. 06. 07.

Spectroscopy and Quantum Chemical Studies참고 문헌 94인용 수 19

한 줄 요약

이 논문은 상호작용 분자의 전체 단체 스펙트럼 함수를 계산하기 위한 확장 가능하고 예측 가능하며 자기일관된 many-body 섭동 이론 방법— moment-truncated iterated perturbation theory (AGF2)—를 소개한다. 스펙트럼 모멘트 보존를 통해 효과적 자가에너지 동역학을 압축함으로써, 이 방법은 O[N⁵] 스케일링을 달성하고, 대규모 병렬 계산(약 1000 오비탈까지)을 가능하게 하며, 아르테미신과 같은 복잡한 시스템의 전자 여기 상태 스펙트럼을 정확히 해석함으로써 오랫동안 지속된 그들의 최저 에너지 여기 상태에 대한 할당의 모순을 해결한다.

ABSTRACT

A reliable and efficient computation of the entire single-particle spectrum of correlated molecules is an outstanding challenge in the field of quantum chemistry, with standard density functional theory approaches often giving an inadequate description of excitation energies and gaps. In this work, we expand upon a recently-introduced approach which relies on a fully self-consistent many-body perturbation theory, coupled to a non-perturbative truncation of the effective dynamics at each step. We show that this yields a low-scaling and accurate method across a diverse benchmark test set, capable of treating moderate levels of strong correlation effects, and detail an efficient implementation for applications up to $\sim1000$ orbitals on parallel resources. We then use this method to characterise the spectral properties of the artemisinin anti-malarial drug molecule, resolving discrepancies in previous works concerning the active sites of the lowest energy fundamental excitations of the system.

연구 동기 및 목표

• 강한 상호작용을 갖는 분자의 전체 단체 스펙트럼 함수를 계산하기 위한 확장 가능하고 정확한 방법을 개발하는 것.
• 표준 DFT 및 섭동 방법의 국한된 점을 해결함으로써 상호작용 시스템에서의 여기 에너지와 밴드 갭을 기술하는 것.
• 복잡한 분자 시스템에 대한 대규모 시뮬레이션(약 1000 오비탈까지)을 가능하게 하며, 효율적인 병렬 구현을 제공하는 것.
• 항말라리아 약물인 아르테미신 분자에서 최저 에너지 여기 상태의 공간적 국소화에 대한 문헌상의 모순을 해결하는 것.

제안 방법

• 이 방법은 보조 스펙트럼 함수 접근법(AGF2)을 완전히 자기일관되게 적용하여, 보조 자유도를 통해 두 번째 순서 자가에너지 문제를 정적 고유값 문제로 재구성한다.
• 첫 번째 스펙트럼 모멘트를 보존하는 모멘트 절단 근사법을 사용하여 정확성을 유지하면서 계산 비용을 감소시킨다.
• 효과적 자가에너지는 Σpq(ω) = Σk (vpk v†qk)/(ω − Ek) 형태의 극점 합으로 표현되어, 디슨 방정식의 행렬 고유값 문제로 재구성된다.
• 알고리즘은 오픈소스 PySCF 프레임워크에 구현되어 있어 효율적인 병렬 처리와 대규모 응용을 가능하게 한다.
• 입자 및 정공 여기 상태를 별도로 나누는 것을 피함으로써, 모든 에너지 스케일에서 전체 여기 스펙트럼에 동시에 접근할 수 있다.
• 이 방법은 자가일관된 ADC(2) 방법의 확장과 동치이며, 동적 자가에너지 성분의 체계적 절단 덕분에 더 높은 스펙트럼 정확도를 제공한다.

실험 결과

연구 질문

• RQ1 모멘트 절단 및 자기일관된 many-body 섭동 방법이 상호작용 분자의 스펙트럼에 대해 높은 정확도와 낮은 계산 스케일링을 동시에 달성할 수 있는가?
• RQ2 다양한 분자 시스템에서 모멘트 절단 AGF2 방법의 스펙트럼 정확도가 ADC(2), GW, EOM-CCSD와 같은 기존 방법과 비교하여 어떻게 되는가?
• RQ3 이 방법은 아르테미신과 같은 복잡한 분자에서 여기 상태 특성의 할당에 대한 애매함을 어느 정도 해결할 수 있는가?
• RQ4 병렬 아키텍처에서 이 방법을 사용하여 신뢰할 수 있는 스펙트럼 정확도로 처리할 수 있는 최대 시스템 크기(오비탈 수 기준)는 얼마인가?
• RQ5 모멘트 절단 전략은 상호작용 시스템에서 전체 동적 자가에너지 기술과 비교해 스펙트럼 해상도를 향상시키는가?

주요 결과

• AGF2 방법은 효율적인 병렬 처리를 통해 O[N⁵] 계산 스케일링을 달성하여 약 1000 오비탈까지의 시스템에서의 시뮬레이션을 가능하게 한다.
• GW100 벤치마크 세트에서 AGF2는 전하 여기 상태에 대해 ADC(2) 및 GW와 유사한 정확도를 보이며, EOM-CCSD보다 더 뛰어난 스케일링 성능을 확보한다.
• 이 방법은 아르테미신에서 최저 에너지 여기 상태에 대한 문헌상의 모순된 할당을 성공적으로 해결하였으며, 실험적 피크 전자 스펙트로스코피와 일치하는 바에 따라 활성 부위로 엔도퍼옥사이드 다리를 규명하였다.
• 모멘트 절단 자가에너지 근사는 엄격한 보존 근사 조건을 위반하더라도 핵심 스펙트럼 특징을 유지하며, 비절단 자기일관된 접근보다 정확도가 향상된다.
• PySCF에 구현된 이 방법은 재현 가능한 대규모 계산을 가능하게 하며, 에너지 의존 스펙트럼 함수에 직접적인 액세스를 제공한다.
• 완전한 자기일관성 덕분에 기저 상태 선택에 민감하지 않으며, 강한 상호작용 영역에서의 강건성을 향상시킨다.