[논문 리뷰] Efficient and scalable Path Integral Monte Carlo Simulations with worm-type updates for Bose-Hubbard and XXZ models

Nicolas Sadoune, Lode Pollet|arXiv (Cornell University)|2022. 04. 26.

Physics of Superconductivity and Magnetism인용 수 4

한 줄 요약

이 논문은 보즈-후버드 및 XXZ 모델에 대한 경로적분 몬테카를로 시뮬레이션을 위한 효율적이고 확장 가능한, C++14 호환 오픈소스 웜 알고리즘 구현을 제시한다. 이 알고리즘은 시스템 크기와 역온도에 대해 거의 선형 스케일링을 달성하며, 초유체 상에서 자가상관 시간이 거의 1에 가까워져 강한 상관관계를 가지는 양자 시스템의 임계 현상에 대한 정확한 시뮬레이션을 가능하게 한다.

ABSTRACT

We present a novel and open-source implementation of the worm algorithm, which is an algorithm to simulate Bose-Hubbard and sign-positive spin models using a path integral representation of the partition function. The code can deal with arbitrary lattice structures and assumes spin-exchange terms, or bosonic hopping amplitudes, between nearest-neighbor sites, and local or nearest-neighbor interactions of the density-density type. We explicitly demonstrate the near-linear scaling of the algorithm with respect to the system volume and the inverse temperature and analyze the autocorrelation times in the vicinity of a U(1) second order phase transition. The code is written in such a way that extensions to other lattice models as well as closely-related sign-positive models can be done straightforwardly on top of the provided framework.

연구 동기 및 목표

• 서명이 양수인 양자 모델을 위한 확장 가능하고 오픈소스의 웜 알고리즘 구현을 개발하는 것.
• 근처 이웃 터널링과 밀도-밀도 상호작용을 포함한 임의의 격자 구조에서의 시뮬레이션을 가능하게 하는 것.
• 보스-후버드 및 XXZ 모델에 대해 시스템 크기와 역온도에 대해 거의 선형 스케일링을 달성하여 임계 영역 연구를 효율적으로 수행하는 것.
• 미래의 다른 격자 모델과 관련된 서명이 양수인 시스템으로의 확장에 적합한 모듈식이고 확장 가능한 프레임워크를 제공하는 것.

제안 방법

• 알고리즘은 포크 상태를 사용한 경로적분 표현과 운동에너지에 대한 강한결합 전개를 사용한다.
• 단일 입자 그린 함수를 샘플링하기 위해 INSERTWORM-GLUEWORM 및 MOVEWORM 업데이트를 포함한 웜 유형의 업데이트 방식을 적용한다.
• 시간 순서를 유지하기 위해 다체 양자장 이론에서 유도된 순서대로 터널링 사건을 저장하기 위해 C++의 벡터-리스트 데이터 구조를 사용한다.
• ALPSCore 라이브러리에 의존하며, 정적 격자 정의와 XML 기반 격자 입력을 모두 지원하여 유연성을 확보한다.
• 입자 수, 에너지, 밀도-밀도 상관관계, 도약 수, 그리고 허수시간 그린 함수와 같은 주요 물리량은 전용 추정기로 계산된다.
• 코드는 높은 성능을 위해 설계되었으며, 메모리 사용량은 터널링 사건 수에 비례하고, 단일 노드 및 슈퍼컴퓨터에서 효율적으로 스케일링된다.

실험 결과

연구 질문

• RQ1 보즈-후버드 및 XXZ 모델에 대해 웜 알고리즘을 어떻게 구현할 수 있을까? 이때 시스템 크기와 역온도에 대해 거의 선형 스케일링이 보장되어야 한다.
• RQ2 U(1) 두 번째 계급 전이점 근처에서 자가상관 시간의 행동은 어떠한가? 이는 시스템 크기와 역온도에 따라 어떻게 변화하는가?
• RQ3 특히 임계 영역에서, 작은 시스템에서 큰 시스템으로 확장할 때 성능과 메모리 효율성이 유지되는가?
• RQ4 유지보수성, 확장 가능성, 성능 측면에서 이 구현은 이전의 웜 알고리즘 코드와 비교해 어떻게 다른가?

주요 결과

• 알고리즘이 시스템 부피와 역온도에 대해 거의 선형 스케일링을 보이며, 대규모 시뮬레이션을 가능하게 한다.
• 초유체 상에서 자가상관 시간이 거의 1에 가까워지며, 웜이 장거리 상관관계를 효과적으로 탐색할 수 있음을 시사한다.
• 메모리 사용량은 터널링 사건 수에 비례하며, 중간 β 값에서 96×96 시스템에 대해 약 70MB로 추정되며, 변동성에 따라 두 배로 증가한다.
• 코드는 높은 성능을 유지하며, 3.1 GHz i5 프로세서를 탑재한 단일 iMac 노드에서 2×10^7 건 이상의 업데이트를 초당 수행한다.
• 이 구현은 매우 확장 가능하며, 정적 또는 XML 기반 정의를 통해 임의의 격자를 지원하고, 새로운 모델에 쉽게 적용할 수 있다.
• 보즈-후버드 및 XXZ 모델에 대한 자동화된 테스트를 통해 정확도가 정확 대각화 결과와 비교하여 확인되었다.