[논문 리뷰] Quantum Monte Carlo determinantal algorithm without Hubbard-Stratonovich transformation: a general consideration

A. N. Rubtsov|arXiv (Cornell University)|2003. 02. 12.

Theoretical and Computational Physics인용 수 4

한 줄 요약

이 논문은 상호작용하는 페르미온을 위한 연속시간 양자 몽테카를로 행렬식 알고리즘을 제안하며, 허버드-스트라토니쇼프 변환을 회피하기 위해 작용을 가우시안 부분과 섭동 상호작용 부분으로 분리한다. 이 방법은 상호작용의 거듭제곱에 대한 분할 함수의 급수 전개를 사용하며, 부호 문제와 계산 비용을 분석하여, 전통적인 방법에 비해 국소적이지 않은 상호작용과 시간에 국한되지 않은 상호작용에 대해 유리한 스케일링을 보인다.

ABSTRACT

Continuous-time determinantal algorithm is proposed for the quantum Monte Carlo simulation of the interacting fermions. The scheme does not invoke Hubbard-Stratonovich transformation. The fermionic action is divided into two parts. One of them contains the interaction and certain additional terms; another one is purely Gaussian. The first part is considered as a perturbation. Terms of the series expansion for the partition function are generated in a random walk process. The sign problem and the complexity of the algorithm are analyzed. We argue that the scheme should be useful particularly for the systems with non-local interaction.

연구 동기 및 목표

• 상호작용하는 페르미온을 위한 양자 몽테카를로 알고리즘을 개발하여 허버드-스트라토니쇼프 변환을 피한다.
• 특히 비국소적이고 장거리 상호작용에 대해 계산 비효율성과 부호 문제를 해결한다.
• 지연된 상호작용이나 소산이 있는 시스템과 같은 시간에 국한되지 않은 상호작용을 가진 시스템의 시뮬레이션을 가능하게 한다.
• 국소적 상호작용에 대해 표준 방법과 유사한 계산 스케일링을 유지하면서 비국소적 경우에서 성능을 향상시킨다.
• 부호 문제를 최소화하기 위해 최적화된 매개변수를 갖춘 연속시간 행렬식 QMC의 일반적인 프레임워크를 제공한다.

제안 방법

• 페르미온 작용이 가우시안 부분 $ S_0 $ 와 섭동 부분 $ W $ 로 분리되며, $ S_0 $ 는 운동 에너지와 매개변수화된 이동 $ \alpha $ 를 포함한다.
• 분할 함수는 상호작용의 거듭제곱에 대한 급수 전개로 전개되며, 상호작용 정점 위의 랜덤 워크를 통해 항이 생성된다.
• 행렬식 $ D $ 는 $ S_0 $ 로부터 유도된 전체 전파함수 $ G_0 $ 를 사용하여 계산된다.
• 매개변수 $ \alpha $ 는 평균 무게 $ \overline{|W|} $ 를 최소화하도록 최적화되어 부호 문제와 자기상관 시간을 감소시킨다.
• 이 방법은 연속된 허수 시간에서 작동하여 시간 이산화에 의한 체계적 오차를 피한다.
• 알고리즘은 연속시간 프레임워크에서 제안되었지만, 적분을 합으로 대체함으로써 이산 시간 형태로 유도할 수 있다.

실험 결과

연구 질문

• RQ1 행렬식 QMC 알고리즘을 허버드-스트라토니쇼프 변환 없이 구성할 수 있는가?
• RQ2 비국소적 상호작용에 대해 이러한 알고리즘의 계산 스케일링은 기존 방법과 어떻게 비교되는가?
• RQ3 최적의 $ \alpha $ 선택을 통해 부호 문제를 효과적으로 완화할 수 있는가?
• RQ4 시간에 국한되지 않은 상호작용을 가진 시스템에 대해 이 방법이 유리한 스케일링을 유지하는가?
• RQ5 이 알고리즘은 장거리 또는 지연된 상호작용을 가진 시스템에 적용 가능한가?

주요 결과

• 한 번의 몽테카를로 스텝당 계산 비용은 일반 허버드 유형 모델에 대해 $ (\beta U N)^2 $ 으로 스케일링되고, 해밀토니안 시스템에 대해서는 $ N^2 $ 으로, 기존의 허버드-스트라토니쇼프 기반 방법과 동일한 스케일링을 보인다.
• 비국소적 상호작용의 경우, 이 알고리즘은 국소적 상호작용과 거의 동일한 계산 비용을 유지한다. 반면 기존 방법은 $ \beta N^2 / \tau_0 $ 로 스케일링된다.
• 부호 문제의 수준은 허버드 모델의 히르슈 방법과 유사하지만, $ \alpha \approx \overline{n} $ 로 선택함으로써 완화될 수 있으며, 반발력 상호작용의 경우 최적값은 약 1 근처이다.
• 반만포화 상태에서 $ \alpha = 1/2 $ 이면 홀수 $ k $ 에 대해 $ \Omega_k = 0 $ 이 되어 대칭에 의해 유도된 부호 상쇄가 발생하며, 이는 랜덤 워크에서 특별한 대우가 필요함을 나타낸다.
• 보조 장이 없고 연속시간 형식이므로, 지연 상호작용과 같은 시간에 국한되지 않은 상호작용의 시뮬레이션을 가능하게 한다.
• 이 알고리즘은 일반적이며, 적분을 합으로 대체함으로써 이산 시간으로 변환할 수 있으며, 동일한 스케일링 행동을 유지한다.