Full configuration interaction quantum Monte Carlo for accurate $\textit{ab initio}$ nuclear structure calculations: algorithms and calculation details

기술 요약: 정확한 Ab Initio 핵 구조 계산을 위한 전 구성 상호작용 양자 몬테카를로(FCIQMC)

문제 정의
Ab initio 핵 이론은 카이랄 유효 장론에서 유도된 미시적 해밀토니안으로부터 유한 핵을 기술하는 것을 목표로 한다. 다체 전개 방법(예: Coupled Cluster, IMSRG, MBPT)은 더 무거운 핵에 대한 계산을 확장하지만, 이는 정량화하기 어려운 잔여 다체 불확실성을 초래하는 절단(truncation)에 의존한다. 반면, No-Core Shell Model (NCSM)은 전 구성 상호작용(Full Configuration Interaction, FCI) 해에 직접 접근할 수 있는 경로를 제공하지만, 지수 함수적으로 증가하는 힐베르트 공간의 제약을 받는다. 변형(deformation), 클러스터링 또는 열린 껍질 구조로 인해 발생하는 강한 고차 상관관계가 존재하는 시스템에서, 명시적으로 다체 해밀토니안 행렬을 저장하지 않고도 큰 모델 공간 내에서 FCI 극한에 도달할 수 있는 확률론적 솔버(stochastic solver)에 대한 절박한 필요성이 존재한다.

방법론
본 연구는 원래 양자 화학을 위해 개발된 Full Configuration Interaction Quantum Monte Carlo (FCIQMC) 알고리즘을 핵 구조 계산에 적응시킨 것이다. 이 방법은 선택된 대칭 섹터(입자 수, 패리티 및 각운동량 투영) 내에서 슈뢰딩거 방정식을 풀기 위해 허수 시간(imaginary time)에서의 확률론적 투영 접근법을 채택한다.

설명된 주요 알고리즘 구성 요소는 다음과 같다:

• 허수 시간 전파 (Imaginary-Time Propagation): 파동 함수는 슬레이터 행렬식(Slater determinants)에 분포된 부호가 있는 워커(signed walkers)의 집단으로 표현된다. 전파 과정은 스폰(spawning), 사멸/복제(death/cloning), 그리고 애니힐레이션(annihilation) 단계를 포함하며, 이를 통해 페르미온의 부호 구조(fermionic sign structure)를 확률론적으로 확립한다.
• 들뜸 생성 (Excitation Generation): 균일한 생성보다 효율을 높이기 위해, 저자들은 해밀토니안 행렬 원소가 큰 들뜸을 우선적으로 샘플링하는 사전 계산된 히트 배스(heat-bath) 생성기를 활용한다.
• 이니시에이터 근사 (Initiator Approximation, i-FCIQMC): 거대한 힐베르트 공간에서의 부호 문제(sign problem)를 극복하기 위해, 약하게 점유된 행렬식으로부터의 스폰을 제한한다. 이는 유한한 워커 편향(finite-walker bias)을 도입하며, 이는 워커의 수가 증가함에 따라 체계적으로 감소한다.
• 적응형 시프트 (Adaptive Shift, AS) 보정: 비-이니시에이터 행렬식에 국소적 시프트를 할당하여, 워커 수가 무한대로 접근함에 따라 정확한 극한을 보존하면서 FCI 극한으로의 수렴을 가속화한다.
• 관측량 추정량 (Observable Estimators):
  • 투영/시행 추정량 (Projected/Trial Estimators): 참조 행렬식 또는 압축된 시행 파동 함수를 사용하여 에너지를 계산하는 데 사용된다.
  • 순수 추정량 (Pure Estimators): 반지름과 같은 비가환 관측량에 필수적이다. 이는 이차 추정량에 내재된 편향을 제거하기 위해 두 개의 독립적인 워커 복제본을 사용하는 축약 밀도 행렬(Reduced Density Matrix, RDM) 샘플링을 통해 계산된다.
• 들뜬 상태 (Excited States): 다수의 상태를 동일한 대칭 섹터 내에서 직교화된 다중 상태 전파를 통해 해결함으로써 들뜬 상태로 확장된다. 여기서 여러 워커 앙상블이 병렬로 전파되며 각 단계에서 직교화된다.

주요 기여

1. 형식론의 적응: 본 논문은 핵 해밀토니안을 위한 FCIQMC의 구체적인 구현, 즉 3체 핵력을 다루는 방식과 핵 상호작용에 요구되는 특정 들뜸 생성 체계를 상세히 설명한다.
2. 결정론적 FCI와의 검증: 작은 모델 공간($e_{\text{max}}=2$)에서 $^4$He, $^8$Be, $^{12}$C, $^{16}$O에 대해, 확률론적 FCIQMC 결과(에너지 및 제곱 점양성자 반지름 $R_p^2$)가 통계적 불확실성 내에서 정확한 결정론적 FCI 값을 재현함을 보여준다.
3. 체계적 오차 제어: 이니시에터 근사에서의 잔여 유한-워커 편향이 어떻게 제어되고 정량화될 수 있는지 입증한다. 저자들은 큰 모델 공간($e_{\text{max}}=10$)에서 거듭제곱 법칙 외삽($O(N_w) = O_\infty + \lambda N_w^{-\alpha}$)을 사용하는 전략을 도입한다. 여기서 지수 $\alpha$는 빠르게 수렴하는 투영 에너지와 더 느리게 수렴하는 순수 에너지 사이의 일관성을 요구함으로써 제약된다.
4. 들뜬 상태로의 확장: 본 논문은 명시적인 대각화 없이 동일한 대칭 섹터 내에서 다수의 상태를 해결함으로써 낮은 들뜸 상태의 스펙트럼을 계산할 수 있는 능력을 보여준다.

결과

• 벤치마킹: $e_{\text{max}}=2$ 공간에서, FCIQMC 에너지와 반지름은 FCI 벤치마크와 일치한다. 적응형 시프트 보정은 유한-워커 편향을 크게 줄이고 이니시에터 임계값($n_{\text{init}}$)에 대한 의존성을 약화시킨다.
• 큰 모델 공간 ($e_{\text{max}}=10$):
  • $^4$He의 경우, 계산은 워커 수($N_w$)가 증가함에 따라 NCSM 결과와 일치하며 직접 수렴한다.
  • $^8$Be, $^{12}$C, $^{16}$O의 경우, 투영 에너지는 빠르게 수렴하는 반면, RDM 기반 순수 추정량(에너지 및 반지름)은 더 느린 수렴성을 보인다. 에너지 제약 거듭제곱 법칙 외삽을 적용함으로써, 저자들은 절단된 다체 전개 방법(MBPT, CC, IMSRG)을 넘어서는 상관관계를 포착하는 결과를 얻는다.
  • 외삽 전략은 결과가 적응형 시프트 오프셋 매개변수($\tilde{\Delta}$)에 독립적이며, 서로 다른 이니시에터 임계값 간에 일관됨을 보여줌으로써 검증된다.
• 들뜬 상태: $^6$Li에 대해, 이 방법은 다섯 개의 가장 낮은 상태($1^+, 3^+, 0^+, 2^+_1, 2^+_2$)를 성공적으로 분해한다. 들뜸 스펙트럼은 낮은 워커 수에서도 올바르게 정렬되며, 다체 절단($N_{\text{max}}$)에 대해 NCSM 결과보다 빠르게 수렴한다.

의의
본 논문은 FCIQMC를 힐베르트 공간의 크기로 인해 결정론적 방법이 실패하는 핵 시스템에서 FCI 극한에 접근할 수 있는, 체계적으로 개선 가능한 확률론적 솔버로 위치시킨다. 이는 두 가지 주요 역할을 수행한다:

1. 정확한 솔버: 큰 모델 공간에서 접합된 상관관계를 포착하여 지상 상태 에너지 및 반지름에 대한 고정밀 벤치마크를 제공한다.
2. 벤치마크 도구: 절단된 다체 전개 방법들의 잔여 다체 불확실성을 평가하기 위한 참조를 제공함으로써, 정확한 해에 근접한 기준을 제시한다.

이 연구는 FCIQMC가 닫힌 껍질 및 열린 껍질 핵을 동일한 프레임워크 내에서 다룰 수 있고, 들뜬 상태 분광학으로 확장될 수 있음을 확인하며, 보고된 물리 결과가 담긴 동반 서신 [80]의 알고리즘적 및 수치적 토대를 구축한다. 저자들은 향후 연구로서 명시적인 3체 상호작용, 더 큰 단일 입자 기저, 그리고 집단 및 전기약 상호작용 관측량에 대한 응용을 언급하였다.