[논문 리뷰] Addressing Strong Correlation by Approximate Coupled-Pair Methods with Active-Space and Full Treatments of Three-Body Clusters

Ilias Magoulas, Jun Shen외 1명|arXiv (Cornell University)|2021. 11. 27.

Advanced Chemical Physics Studies참고 문헌 222인용 수 9

한 줄 요약

이 논문은 활성 공간과 세 개의 전자 집단(3체, T3) 집단 응축을 동시에 고려하는 새로운 근사 결합 쌍(ACP) 방법을 소개한다. 또한 점유 및 비점유 오비탈 수에 기반한 (T2)² 다이어그램의 적응적 스케일링을 포함한다. 이 방법은 기존의 CCSD 및 CCSDT가 수렴하지 못하는 강한 상관관계를 가지는 시스템—H6, H10 고리, H50 사슬—에서 안정적으로 수렴하며, 전체 CI 수준의 정확도를 유지하면서 계산 비용을 감소시킨다.

ABSTRACT

When the number of strongly correlated electrons becomes larger, the single-reference coupled-cluster (CC) CCSD, CCSDT, etc. hierarchy displays an erratic behavior, while traditional multi-reference approaches may no longer be applicable due to enormous dimensionalities of the underlying model spaces. These difficulties can be alleviated by the approximate coupled-pair (ACP) theories, in which selected $(T_2)^2$ diagrams in the CCSD amplitude equations are removed, but there is no generally accepted and robust way of incorporating connected triply excited ($T_3$) clusters within the ACP framework. It is also not clear if the specific combinations of $(T_2)^2$ diagrams that work well for strongly correlated minimum-basis-set model systems are optimum when larger basis sets are employed. This study explores these topics by considering a few novel ACP schemes with the active-space and full treatments of $T_3$ correlations and schemes that scale selected $(T_2)^2$ diagrams by factors depending on the numbers of occupied and unoccupied orbitals. The performance of the proposed ACP approaches is illustrated by examining the symmetric dissociations of the $ ext{H}_6$ and $ ext{H}_{10}$ rings using basis sets of the triple- and double-$ζ$ quality and the $ ext{H}_{50}$ linear chain treated with a minimum basis, for which the conventional CCSD and CCSDT methods fail.

연구 동기 및 목표

• 강한 상관관계를 가지며 많은 전자가 얽힌 전자들로 구성된 시스템에서 기존의 단일 참고 CCSD 및 CCSDT 방법이 실패하는 문제를 해결하기 위해.
• 근사 결합 쌍(ACP) 형식 내에서 세 개의 전자 집단(T3) 집단 응축을 통합하기 위한 견고한 프레임워크를 개발하기 위해.
• 강한 상관관계 영역에서 에너지 수렴성과 정확도에 미치는 활성 공간 대비 전체 T3 처리의 영향을 조사하기 위해.
• 점유 오비탈 수(no)와 비점유 오비탈 수(nu)에 기반한 (T2)² 다이어그램 스케일링이 큰 시스템에서 안정성과 정확도 향상에 얼마나 효과적인지 평가하기 위해.

제안 방법

• 특정 (T2)² 다이어그램을 CCSD 앰리티드 방정식에서 제거하면서도 T3 집단 응축을 유지하는 새로운 ACP 기법을 제안한다.
• 두 가지 변형을 도입한다: 하나는 T3 집단을 전체적으로 처리하고, 다른 하나는 활성 공간 기반 T3 집단 응축을 사용한다.
• 점유 오비탈 수(no)와 비점유 오비탈 수(nu)의 비율에 기반해 (T2)² 다이어그램에 적응적 스케일링 인자를 적용한다.
• CCSD, CCSDt, CCSDT 프레임워크를 기초로 하되, ACP 근사에 의해 수정한다.
• 스핀 적응형이고 수직 결합된 양자역학적 방법을 사용해 이중 전자 전이(T2) 집단 연산자를 적용하여 전자 상관관계를 처리한다.
• 결과를 전체 구성상태 상호작용(FCI) 및 대규모 밀도 행렬 재정렬 그룹(LDMRG) 데이터와 비교하여 검증한다.

실험 결과

연구 질문

• RQ1 CCSD 및 CCSDT가 수렴하지 못하는 강한 상관관계 시스템에서, T3 집단을 전체 또는 활성 공간 기반으로 처리하는 것이 ACP 방법의 안정성을 향상시킬 수 있는가?
• RQ2 no/(no+nu) 및 nu/(no+nu) 비율에 따라 (T2)² 다이어그램을 스케일링함으로써 ACP 방법의 정확도와 수렴성이 어떻게 향상되는가?
• RQ3 H50와 같은 큰 시스템에서, (T2)² 다이어그램 제거 전략의 성능은 기저함수 크기에 따라 달라지며, 특히 큰 기저함수에서의 영향은 어떠한가?
• RQ4 더 큰 기저함수에서 (T2)² 다이어그램 제거 패atters가 최적인지, 아니면 재가중 조정이 필요할까?
• RQ5 ACP 방법은 전체 CI의 지수적 스케일링을 피하면서도 H6, H10, H50 시스템의 대칭적 분리 상황에서 거의 FCI 정확도를 달성할 수 있는가?

주요 결과

• ACCSD(1, 3+4/2) 방법은 DCSD와 동일하며, H6(cc-pVTZ) 및 H10(DZ) 고리에서 거의 FCI 정확도를 달성하고 안정적인 수렴을 보였다.
• STO-6G 기저함수를 사용한 H50 사슬에서는 CCSD가 2.0 보هر를 초과하면 수렴하지 못하지만, ACCSD(1, 3+4/2)는 수렴하고, ACCSDT(1, 3+4/2)는 1.8 보هر를 초과해도 수렴하였다.
• ACP 방법에서 전체 T3 처리는 활성 공간 T3 처리보다 더 낮은 오차를 보였으며, 특히 분리 상태에서 뚜렷한 이점이 있었다.
• no/(no+nu) 및 nu/(no+nu) 비율에 따라 (T2)² 다이어그램을 스케일링함으로써 모든 시험된 시스템에서 안정성과 정확도가 크게 향상되었다.
• H6 및 H10 시스템에서는 평균 절대 오차(MUE)가 1 mEh 이하로 유지되었고, H50 시스템에서는 10 mEh 이내의 오차를 유지를 하였다.
• 제안된 ACP 기법은 강한 상관관계 시스템에서 에너지 정확도와 수렴성 측면에서 표준 CCSD 및 CCSDT를 모두 능가하였다.