[논문 리뷰] Renormalized Second-order Perturbation Theory for The Electron Correlation Energy: Concept, Implementation, and Benchmarks

Xinguo Ren, Patrick Rinke외 2명|arXiv (Cornell University)|2012. 12. 15.

Advanced Chemical Physics Studies참고 문헌 1인용 수 4

한 줄 요약

이 논문은 Kohn-Sham 프레임워크 내에서 RPA, SOSEX, 그리고 엄밀하게 다루어진 재규합 단일(excited) 상태(rSE)를 조합한 개선된 재규합 제2차 섭동이론(rPT2)을 소개한다. 이 방법은 공유결합, 비공유결합, 그리고 약한 결합 시스템 전반에서 균형 잡힌 정확도를 달성하며, 특히 희귀기체 이량체와 S22/S66 기준치에서 분산력에 의해 지배되는 상호작용에 대해 이전의 근사적 rSE 처리 방식보다 향상된 성능을 보인다.

ABSTRACT

We present a renormalized second-order perturbation theory (rPT2), based on a Kohn-Sham (KS) reference state, for the electron correlation energy that includes the random-phase approximation (RPA), second-order screened exchange (SOSEX), and renormalized single excitations (rSE). These three terms all involve a summation of certain types of diagrams to infinite order, and can be viewed as "renormalization" of the 2nd-order direct, exchange, and single excitation (SE) terms of Rayleigh-Schrödinger perturbation theory based on an KS reference. In this work we establish the concept of rPT2 and present the numerical details of our SOSEX and rSE implementations. A preliminary version of rPT2, in which the renormalized SE (rSE) contribution was treated approximately, has already been benchmarked for molecular atomization energies and chemical reaction barrier heights and shows a well balanced performance [Paier et al, New J. Phys. 14, 043002 (2012)]. In this work, we present a refined version of rPT2, in which we evaluate the rSE series of diagrams rigorously. We then extend the benchmark studies to non-covalent interactions, including the rare-gas dimers, and the S22 and S66 test sets. Despite some remaining shortcomings, we conclude that rPT2 gives an overall satisfactory performance across different chemical environments, and is a promising step towards a generally applicable electronic structure approach.

연구 동기 및 목표

• 약한 결합 시스템에서 표준 RPA가 기인하는 결합력 부족 문제를 수정할 수 있는 체계적으로 개선된 전자 상관 방법을 개발하기 위해.
• 이전 rPT2 버전에서 근사적으로 처리된 바 있었던 재규합 단일 상태(rSE) 기여를 엄밀하게 구현하기 위해.
• 공유결합 에너지, 반응율 장벽, 비공유결합 상호작용(예: 분산력)을 포함한 다양한 화학적 환경에서 rPT2를 기준치로 삼기 위해.
• 균형 잡힌 성능과 수렴 행동을 입증함으로써 rPT2가 전망 가능하고 일반적으로 적용 가능한 전자구조 방법이 될 수 있음을 확립하기 위해.

제안 방법

• 이 방법은 세 가지 무한도형 합산을 조합한다: RPA(Random-phase approximation), SOSEX(제2차 스크리닝된 교환), rSE(재규합 단일 상태), 모두 Kohn-Sham 기반의 제2차 Rayleigh-Schrödinger 섭동이론에서 유도된다.
• rSE 기여는 Fock 행렬의 준표준형 변환을 통해 계산되며, 고립된 상태와 가상 오비탈 블록을 별도로 대각화하여 에너지 표현식에서 역행렬을 효율적으로 평가할 수 있도록 한다.
• 최종 rSE 에너지는 변환된 오비탈에 대한 합으로 표현된다: $ E_{\text{c}}^{\text{rSE}} = \sum_{ia} \frac{\tilde{f}_{ai}\tilde{f}_{ia}}{\tilde{\epsilon}_{i}-\tilde{\epsilon}_{a}} $, 여기서 $ \tilde{f}_{ia} $ 및 $ \tilde{\epsilon} $ 는 재정의된 전이 진폭과 고유값을 의미한다.
• 이 방법은 교환-상관 기능(예: PBE)을 사용한 Kohn-Sham 기반을 사용하며, 상관 에너지를 SCF 이후 보정으로 평가한다.
• 기저함수 수렴성은 Dunning 스타일의 cc-pVXZ/aug-cc-pVXZ 및 FHI-aims NAO 타이어 기저함수를 사용하여 평가하며, BSSE 보정을 적용한다.
• 구현은 원자화 에너지, 반응율 장벽, 비공유결합 상호작용 세트(S22, S66, 희귀기체 이량체)를 대상으로 한 기준치를 통해 검증된다.

실험 결과

연구 질문

• RQ1 엄밀한 재규합 단일 상태(rSE) 기여 처리가 희귀기체 이량체와 같은 약한 결합 시스템의 결합 에너지 곡선에서 비물리적인 거동을 제거할 수 있는가?
• RQ2 통합된 rPT2 프레임워크 내에서 RPA, SOSEX, rSE를 포함함으로써 공유결합, 전이 상태, 분산력에 의해 지배되는 시스템을 포함한 다양한 화학적 환경에서 균형 잡힌 정확도를 달성할 수 있는가?
• RQ3 S22 및 S66 테스트 세트에 대해 rPT2 방법은 비공유결합 상호작용 에너지에서 다른 RPA 기반 방법과 비교해 어떻게 성능을 보이는가?
• RQ4 특히 산란 함수를 포함한 기저함수의 선택, 특히 aug-cc-pV5Z 또는 하이브리드 타이어 3/4 + 산란 함수 기저함수(t3/4+a5Z-d)는 rPT2의 결합 에너지 예측 수렴성과 정확도에 어떤 영향을 미치는가?

주요 결과

• 엄밀한 rSE 처리를 통한 개선된 rPT2 방법은 이전의 근사적 rSE 버전에서 발생했던 희귀기체 이량체의 비물리적 거동을 완전히 제거한다.
• rPT2는 공유결합 에너지(원자화 에너지), 반응율 장벽, 비공유결합 상호작용 전반에서 균형 잡힌 성능을 보이며, 일관성 면에서 표준 RPA 및 다른 RPA 기반 방법보다 뛰어난 성능을 보인다.
• S22 및 S66 테스트 세트에서 rPT2는 RPA보다 유의미하게 향상된 정확도를 보이며, 평균 절대편차가 다른 고급 후-RPA 방법과 비교해 유사하거나 더 낫다.
• 산란 함수의 포함, 특히 aug-cc-pV5Z 또는 하이브리드 타이어 3/4 + 산란 함수 기저함수(t3/4+a5Z-d)는 메탄 이량체와 같은 분산력에 의해 지배되는 시스템의 정확한 기술에 필수적이다.
• 기저함수 수렴성은 표준 cc-pVXZ보다 FHI-aims NAO 타이어 기저함수를 사용할 경우 더 빠르며, 동일한 기저 크기에서 aug-cc-pVXZ는 더 뛰어난 정확도를 제공한다.
• rPT2 방법은 다양한 시스템에서 뛰어난 강건성과 신뢰성을 보이며, 일반적으로 적용 가능한 전자구조 접근법으로서의 잠재력을 확인한다.