[논문 리뷰] Spin-adapted selected configuration interaction in a determinant basis

Chilkuri, Vijay Gopal, Thomas Applencourt외 3명|arXiv (Cornell University)|2018. 12. 17.

Advanced Chemical Physics Studies인용 수 5

한 줄 요약

이 논문은 선택형 구성상태상호작용(SCI)에서 스핀 적응 파동함수를 계산적으로 효율적으로 생성하기 위한 알고리즘을 제시한다. 이 알고리즘은 임의의 스핀 상태를 갖는 디터미낸트 공간을 보완하여 $\hat{S}^2$ 연산자의 고유상태가 되도록 보장한다. 이 방법을 통해 대각화에 구성상태함수(CSF)를 사용할 수 있게 되어, 디터미낸트 기반 선택 방식을 변경하지 않은 채로 최대 4배의 메모리 절감이 가능해지며, 스핀 순수성과 함께 근접한 전자기준결합에너지(CI) 정확도와 개선된 외삽 안정성을 확보한다.

ABSTRACT

Selected configuration interaction (SCI) methods, when complemented with a second-order perturbative correction, provide near full configuration interaction (FCI) quality energies with only a small fraction of the Slater determinants of the FCI space. However, a selection criterion based on determinants alone does not ensure a spin-pure wave function. In other words, such SCI wave functions are not eigenfunctions of the $S^2$ operator. In many situations (bond breaking, magnetic system, excited state, etc), having a spin-adapted wave function is essential for a quantitatively correct description of the system. Here, we propose an efficient algorithm which, given an arbitrary determinant space, generates all the missing Slater determinants allowing one to obtain spin-adapted wave functions while avoiding manipulations involving configuration state functions. For example, generating all the possible determinants with 6 spin-up and 6 spin-down electrons in 12 open shells takes 21 CPU cycles per generated Slater determinant. The selection is still done with individual determinants, and one can take advantage of the basis of configuration state functions in the diagonalization of the Hamiltonian to reduce significantly the memory footprint.

연구 동기 및 목표

• 결합 단절, 자화 상태 및 자기적 시스템에서 정확도를 떨어뜨리는 디터미낸트 기반 선택형 구성상태상호작용(Sci) 방법의 스핀 오염 문제를 해결하기 위해.
• SCI 파동함수가 $\hat{S}^2$ 연산자의 고유상태가 되도록 하기 위해 필요한 모든 슬래터 디터미낸트를 생성함으로써, 디터미낸트 선택 과정을 변경하지 않는 방법을 개발하기 위해.
• 해밀토니안 대각화를 위한 구성상태함수(CSF) 기반으로의 효율적 전환을 가능하게 하여, 저장 메모리 사용량을 크게 줄이기 위해.
• 기존의 SCI 프레임워크인 CIPSI, SHCI, FCIQMC와의 호환성을 유지하여 광범위한 적용 가능성을 확보하기 위해.
• 기본 파동함수의 스핀 오염을 제거함으로써, $E_{\text{PT2}}$ 보정 및 에너지 외삽의 정확도를 향상시키기 위해.

제안 방법

• 알고리즘은 주어진 디터미낸트 집합에 대해 스핀 적응 부분공간을 닫기 위해 필요한 모든 누락된 슬래터 디터미낸트를 식별하고 추가하여, 결과 파동함수가 $\hat{S}^2$ 연산자의 고유상태가 되도록 보장한다.
• 디터미낸트 선택 이후의 후처리 단계로 작동하므로, 원래의 SCI 선택 기준과 계산 워크플로우를 유지한다.
• 디터미낸트 공간의 구조를 활용하여 필요한 추가 항목을 효율적으로 계산하며, 12껍질, 6α/6β 시스템에서 디터미낸트당 평균 21 CPU 사이클만 소요된다.
• 스핀 적응 후 파동함수는 CSF 기반으로 변환되어 대각화가 가능해지며, 이로 인해 메모리 사용량이 극적으로 감소한다.
• 선택 과정 중에 구성상태함수(CSF)를 직접 다루지 않아 효율성과 모듈성을 유지한다.
• 표준 SCI 알고리즘과 호환되며, 기저 파동함수의 스핀 오염을 제거함으로써 정확한 $E_{\text{PT2}}$ 보정이 가능해진다.

실험 결과

연구 질문

• RQ1 디터미낸트 기반 SCI 방법을 선택 알고리즘을 변경하거나 복잡한 CSF 기반 선택을 도입하지 않고도 스핀 적응으로 만들 수 있는가?
• RQ2 스핀 순수성을 확보하기 위해 필요한 최소한의 추가 디터미낸트를 생성하는 데 드는 계산 비용은 얼마인가?
• RQ3 스핀 적응 디터미낸트 선택 이후 CSF 기반으로 전환했을 때, 메모리 사용량은 어느 정도 감소할 수 있는가?
• RQ4 기본 파동함수의 스핀 오염이 결합 단절 과정에서 $E_{\text{PT2}}$ 외삽의 정확도에 어떤 영향을 미치는가?
• RQ5 제안된 방법이 CIPSI, SHCI, FCIQMC와 같은 기존 SCI 방법에 일반적으로 적용 가능한가?

주요 결과

• 알고리즘은 12껍질, 6α/6β 시스템에서 디터미낸트당 평균 21 CPU 사이클만으로 스핀 적응에 필요한 모든 누락된 디터미낸트를 생성하여, 선택 후 후처리 단계에서의 효율적 스핀 적응을 가능하게 한다.
• 스핀 적응 파동함수로 인해 스핀 오염이 감소하여 $E_{\text{PT2}}$ 외삽이 안정적이고 정확해지며, 질소 분자의 이완거리가 2 Å를 초과할 경우 비스핀 적응 계산에서 1 밀리하르트 이상의 진동이 관찰되었다.
• 디바이저드 대각화에서 CSF 기반 저장 방식을 사용할 경우, 디터미낸트 기반 저장 방식 대비 최대 4배의 메모리 절감이 가능하며, 특히 많은 수의 디터미낸트가 선택된 경우에 두드러진다.
• 선택된 디터미낸트 수가 증가할수록 $N_{\text{det}}/N_{\text{CSF}}$ 비율이 증가하여, 더 큰 CI 공간에서의 메모리 절감 효과가 커진다.
• 이 방법은 스핀 오염이 심각하게 영향을 미치는 도전적인 시스템, 예를 들어 분해 중인 N₂에서 정확한 에너지 외삽을 가능하게 한다.
• 이 방법은 CIPSI, SHCI, FCIQMC를 포함한 모든 SCI 방법에 일반적으로 적용 가능하며, 코드 수정이 최소화되고 계산 오버헤드가 거의 없다.