[논문 리뷰] Efficiency of a Multi-Reference Coupled Cluster method

Emmanuel Giner, Grégoire David외 2명|arXiv (Cornell University)|2015. 09. 10.

Advanced Chemical Physics Studies참고 문헌 58인용 수 7

한 줄 요약

이 논문은 결합된 다전자 상태를 고려한 비수축 다전자 결합 클러스터(MRCC) 방법을 제시하며, 기존의 CASSDCI보다 정확도를 향상시키기 위해 ${\hat{T}}^2$ 연산자를 통해 삼중 및 사중 전이를 추정함으로써 MRSDCI 행렬을 반복적으로 보정한다. 이 방법은 결합 깨짐 및 다전자 상태 시스템에서 화학적 정확도(오차 ~1–2 mH)를 달성하며, 전체적으로 병렬 처리가 가능하고 기준 상태 및 전이 선택에 대한 유연성을 제공한다.

ABSTRACT

The multi-reference Coupled Cluster method first proposed by Meller et al (J. Chem. Phys. 1996) has been implemented and tested. Guess values of the amplitudes of the single and double excitations (the ${\hat T}$ operator) on the top of the references are extracted from the knowledge of the coefficients of the Multi Reference Singles and Doubles Configuration Interaction (MRSDCI) matrix. The multiple parentage problem is solved by scaling these amplitudes on the interaction between the references and the Singles and Doubles. Then one proceeds to a dressing of the MRSDCI matrix under the effect of the Triples and Quadruples, the coefficients of which are estimated from the action of ${\hat T}^2$. This dressing follows the logics of the intermediate effective Hamiltonian formalism. The dressed MRSDCI matrix is diagonalized and the process is iterated to convergence. The method is tested on a series of benchmark systems from Complete Active Spaces (CAS) involving 2 or 4 active electrons up to bond breakings. The comparison with Full Configuration Interaction (FCI) results shows that the errors are of the order of a few milli-hartree, five times smaller than those of the CASSDCI. The method is totally uncontracted, parallelizable, and extremely flexible since it may be applied to selected MR and/or selected SDCI. Some potential generalizations are briefly discussed.

연구 동기 및 목표

• 강한 상관관계를 가진 시스템을 위한 크기 일致성, 크기 집중성, 유연성을 갖춘 다전자 결합 클러스터 방법을 개발하기 위해.
• 단일 전이 및 이중 전이와의 기준 상태 상호작용에 기반해 보정된 진폭을 사용하여 MRCC의 다중 기반 문제를 해결하기 위해.
• 반복적인 행렬 보정을 통해 $\hat{T}^2$ 연산자로 삼중 및 사중 전이 기여도를 추정함으로써 CASSDCI보다 정확도를 향상시키기 위해.
• 상태별 형식을 통해 자극 상태 및 선택된 기준 상태 공간에의 적용을 가능하게 하기 위해.
• 계산 비용 효율성, 비수축성, 병렬 처리가 가능한 고정밀 양자화학 계산을 위한 프레임워크를 제공하기 위해.

제안 방법

• 단일 및 이중 전이에 대한 추정 진폭($\hat{T}$)을 MRSDCI 계수에서 추출함으로써 방법을 시작한다.
• 기준 상태와 단일/이중 전이 간의 상호작용 강도에 따라 이러한 진폭을 스케일링함으로써 다중 기반 문제를 해결한다.
• MRSDCI 행렬은 $\hat{T}^2$의 작용을 통해 삼중 및 사중 전이의 기여도를 추정함으로써 반복적으로 보정된다.
• 보정은 중간 효과 해밀토니안 형식을 따르며, 크기 일치성 유지 및 상관 에너지 추정 향상에 기여한다.
• 각 반복 단계에서 보정된 행렬을 대각화하고, 진폭이 안정화될 때까지 반복한다.
• 이 방법은 완전히 비수축적이며 기준 상태 및 전이 선택에 대해 매우 민감하며, 가장 비용이 많이 드는 단계(삼중/사중 전이 계산)에서 본질적으로 병렬 처리가 가능하다.

실험 결과

연구 질문

• RQ1 상태별 비수축 MRCC 방법이 다전자 상태 시스템에서 CASSDCI보다 더 높은 정확도를 달성할 수 있는가?
• RQ2 기준 상태와 단일/이중 전이 간의 상호작용에 기반한 물리적으로 타당한 진폭 스케일링을 통해 MRCC의 다중 기반 문제를 어떻게 해결할 수 있는가?
• RQ3 $\hat{T}^2$로 추정한 삼중 및 사중 전이를 행렬에 반복적으로 보정함으로써 상관 에너지 추정은 어느 정도 향상될 수 있는가?
• RQ4 이 형식은 크기 일치성과 정확도를 유지하면서 자극 상태로 확장될 수 있는가?
• RQ5 선택적 기준 상태 및 전이 공간 절단을 적용할 경우 계산적 트레이드오프와 성능 향상은 어떠한가?

주요 결과

• 이 방법은 풀 CI와 비교해 약 1–2 밀리하르트리의 오차를 기록하며, 이는 CASSDCI의 오차보다 약 5배 작다.
• 결합 깨짐 과정 및 2 또는 4개의 활성 전자를 가진 시스템을 포함한 다양한 벤치마크 시스템에서 정확도가 유지된다.
• 이 방법은 완전히 비수축적이며 병렬 처리가 가능하며, 가장 비용이 많이 드는 단계(삼중/사중 전이)는 분산 계산에 적합하다.
• $\hat{T}^2$ 기반 반복적 보정 절차는 명시적인 삼중 또는 사중 진폭을 사용하지 않고도 고차 상관 효과를 효과적으로 포착한다.
• 이 방법은 선택된 기준 상태 공간 및/또는 선택된 SDCI 구성에 적용 가능하여 계산 비용을 감소시킬 수 있다.
• 상태별 또는 상태 평균 형식을 통해 자극 상태로의 확장이 가능하며, 계산 후 고유벡터의 수직화 처리가 가능하다.