고차 다체 섭동 이론을 위한 확률적 양자 몬테카를로

초록

본 논문은 PTQMC(perturbation theory quantum Monte Carlo)라는 새로운 확률적 알고리즘을 제안한다. 무작위 워커를 이용해 고차 MBPT 계수를 직접 샘플링함으로써 전통적인 지수적 비용을 회피하고, 리처드슨 짝짓기 모델에서 16차까지 정확히 재현한다. 또한 Padé 재합성 기법과 결합해 발산하는 섭동 급수를 안정적인 에너지 추정치로 변환한다. 마지막으로 파동함수 복잡도를 $e^{S}$라는 전역 지표로 정의해, 에너지 수렴만으로는 판단하기 어려운 섭동 전개의 유효성을 평가한다.

상세 분석

PTQMC는 전통적인 다체 섭동 이론이 직면한 두 가지 핵심 문제—구성 공간의 급격한 팽창과 고차 항의 수렴성 판단—를 확률적 샘플링으로 해결한다. 기본 아이디어는 $n$차 MBPT 에너지 항을 “경로 합” 형태로 표현하고, 각 경로를 무작위 워커가 순차적으로 탐색하도록 하는 것이다. 워커는 현재 구성 $|Φ_J⟩$에서 연결된 구성 $|Φ_I⟩$로 전이할 확률을 $|\langle Φ_I|H_1|Φ_J\rangle|/|\Delta_{I0}|$에 비례하게 부여한다. 부호는 행렬 원소와 현재 워커 가중치의 부호에 의해 결정된다. 이렇게 생성된 워커 집합은 $w^{(n)}_I$라는 부호가 있는 인구분포를 형성하고, 이는 정확한 MBPT 재귀식
$$w^{(n+1)}I=\frac{1}{\Delta{I0}}\sum_J\langle Φ_I|H_1|Φ_J\rangle w^{(n)}J-\sum{m=1}^n E^{(m)} w^{(n-m+1)}_I$$
을 무편향 추정한다. 워커가 동일 구성에 겹치면 부호 소멸(annihilation) 과정을 거쳐 통계적 분산을 감소시키며, 이는 FCIQMC와 구조적으로 유사하지만 목표가 전체 파동함수가 아니라 각 섭동 차수별 계수라는 점에서 차별화된다.

알고리즘의 복잡도는 $O(n N_w)$이며, 여기서 $n$은 섭동 차수, $N_w$는 워커 수이다. 이는 전통적인 연산량 $O(N^{2n})$에 비해 크게 개선된 형태다. 논문에서는 이 방법을 리처드슨 짝짓기 모델에 적용해 $g$(상호작용 강도) 파라미터를 다양하게 변조하면서 8차, 16차까지의 계수를 정확히 재현했다. 특히 $g=-1.2\sim-0.76$ 구간처럼 급격히 발산하는 경우에도 통계적 오차가 $N_w^{-1/2}$ 스케일을 따르며 안정적인 결과를 얻었다.

발산하는 원시 섭동 급수를 직접 사용하면 물리적 예측이 불가능하지만, PTQMC가 제공하는 고차 계수를 Padé 근사(