[논문 리뷰] Implementation of the HMC algorithm on the tempered Lefschetz thimble method

Masafumi Fukuma, Nobuyuki Matsumoto외 1명|arXiv (Cornell University)|2019. 12. 31.

Theoretical and Computational Physics참고 문헌 37인용 수 24

한 줄 요약

이 논문은 라티스 양자장 이론에서의 수치적 부호 문제를 해결하기 위해, 혼합 몬테카를로(HMC) 알고리즘을 온도 조절된 레프셰츠 끈(method) (TLTM)에 구현한다. 특히 페르미온계에서의 문제를 다루기 위해, 유동 표면 위에서 분자역학를 수행하기 위해 RATTLE 알고리즘을 적응시켰고, 페르미온 행렬식의 영점에 대한 강력한 처리 방법을 개발함으로써, 2×2 허버드 모델 테스트에서 메트로폴리스 알고리즘 대비 계산 비용을 70% 감소시키고 자가상관 시간을 크게 낮추는 데 성공하였다.

ABSTRACT

The tempered Lefschetz thimble method (TLTM) is a parallel-tempering algorithm towards solving the numerical sign problem, where the system is tempered by the antiholomorphic gradient flow to tame both the sign and ergodicity problems simultaneously. In this paper, we implement the hybrid Monte Carlo (HMC) algorithm for transitions on each flowed surface, expecting that this implementation on TLTM will give a useful framework for future computations of large-scale systems including fermions. Although the use of HMC in Lefschetz thimble methods has been proposed so far, our crucial achievement here is that HMC is implemented on TLTM so as to work within the parallel-tempering algorithm in TLTM, especially by developing an algorithm to handle zeros of fermion determinants in the course of the molecular-dynamics process. We confirm that the algorithm works correctly by applying it to the sign problem of the Hubbard model on a small lattice, for which the TLTM is known to work with the Metropolis algorithm. We show that the use of HMC significantly reduces the autocorrelation times with less computational times compared to the Metropolis algorithm.

연구 동기 및 목표

• 페르미온계에서 고자기상관과 느린 혼합 문제를 야기하는 메트로폴리스 알고리즘의 한계를 극복하기 위해.
• 대규모 페르미온계 시뮬레이션에서 효율성을 향상시키기 위해 HMC를 TLTM에 구현하기 위해.
• 분자역학를 수행하는 동안 발생하는 페르미온 행렬식의 영점 처리 방법을 개발하기 위해.
• 소규모 격자 크기를 가진 허버드 모델을 기준으로 하여 HMC의 TLTM 구현을 검증하기 위해.
• HMC를 TLTM에 적용함으로써, 부호 문제가 있는 시스템에서 메트로폴리스 알고리즘 대비 계산 비용과 자가상관 시간을 감소시킬 수 있는지 보여주기 위해.

제안 방법

• 유동 표면 Σₜ 위에서 분자역학를 수행하기 위해 RATTLE 알고리즘을 적응시켜, 속도-버레트 적분을 통한 제약 조건 유지 보장을 수행한다.
• Σₜ 위에서 제약 조건이 부여된 운동 방정식을 풀기 위해 적응형 단계 크기 제어와 재스케일링을 적용한 뉴턴-라프슨 반복 기법을 도입한다.
• Dₜ(행렬식 영점으로 향하는 영역)에 진입하는 해를 피하기 위해 Dₜ에의 근접도를 정지 조건으로 사용하는 기준을 도입한다.
• 뉴턴 반복이 실패하거나 Dₜ에 진입할 경우, 상세균형을 유지하기 위해 운동량 반전(Ψ)을 후속 조치로 사용한다.
• HMC 동역학을 TLTM 프레임워크에 통합하여, 다양한 유동 시간 t에 걸친 병렬 온도 조절을 가능하게 한다.
• Nₛ=2×2 및 Nₜ=5를 가진 허버드 모델에 알고리즘을 적용하여 정확한 해와 메트로폴리스 기반 TLTM의 결과를 비교한다.

실험 결과

연구 질문

• RQ1 HMC는 TLTM 프레임워크 내에서 페르미온계의 샘플링 효율성을 향상시키기 위해 성공적으로 구현될 수 있는가?
• RQ2 유동 끈 표면 위에서 분자역학를 수행할 때, 반전성을 해치지 않으면서 페르미온 행렬식의 영점을 어떻게 처리할 수 있는가?
• RQ3 부호 문제가 있는 시스템에서 HMC는 메트로폴리스 알고리즘 대비 자가상관 시간과 계산 비용을 줄일 수 있는가?
• RQ4 반화성 이격점에서 벗어난 상태에서 HMC 확장된 TLTM 방법은 안정적이고 정확한가?
• RQ5 페르미온 행렬식 영점 근처 영역을 순환할 때, 알고리즘은 상세균형과 에르고딕성을 유지할 수 있는가?

주요 결과

• 2×2 허버드 모델에서 TLTM에 HMC를 적용함으로써, 독립된 구성 상태당 계산 비용이 메트로폴리스 알고리즘의 약 30%로 감소하였다.
• HMC를 적용한 결과 자가상관 시간이 크게 감소하여 빠른 혼합과 향상된 샘플링 효율성을 나타냈다.
• 허버드 모델에서의 전자 밀도에 대해 정확한 결과를 정확히 재현하여 알고리즘의 정확성을 확인하였다.
• 재스케일링과 영역 탐지 기반의 반복적 뉴턴 해법을 통해 페르미온 행렬식의 영점 처리가 강력하고 안정적이었다.
• 해가 발산되거나 금지된 영역에 진입할 경우 운동량 반전을 사용함으로써 상세균형과 반전성을 유지하였다.
• 이 프레임워크는 확장 가능하며, 자유도가 더 높은 더 큰 시스템에서는 더욱 높은 효율성 향상을 기대할 수 있다.