[논문 리뷰] An efficient algorithm for numerical computations of continuous densities of states

Kurt Langfeld, Biagio Lucini외 2명|arXiv (Cornell University)|2015. 09. 28.

Theoretical and Computational Physics참고 문헌 21인용 수 7

한 줄 요약

이 논문은 연속 자유도를 가진 시스템, 예를 들어 격자 양성거시학 이론과 같은 시스템에서 연속 상태 밀도를 효율적으로 계산하기 위한 개선된 로그 선형 근사화(Logarithmic Linear Relaxation, LLR) 알고리즘을 제시한다. 히스토그램을 피하고 지수적 오차 억제를 가능하게 함으로써, 전통적인 중요도 샘플링이 실패하는 강한 준정적 안정성 영역에서도 $20^4$에 이르는 격자 크기까지 임계 결합 상수와 비열 피크 값과 같은 열역학적 관측량을 고정밀도로 추정할 수 있다.

ABSTRACT

In Wang-Landau type algorithms, Monte-Carlo updates are performed with respect to the density of states, which is iteratively refined during simulations. The partition function and thermodynamic observables are then obtained by standard integration. In this work, our recently introduced method in this class (the LLR approach) is analysed and further developed. Our approach is a histogram free method particularly suited for systems with continuous degrees of freedom giving rise to a continuum density of states, as it is commonly found in Lattice Gauge Theories and in some Statistical Mechanics systems. We show that the method possesses an exponential error suppression that allows us to estimate the density of states over several orders of magnitude with nearly-constant {\it relative} precision. We explain how ergodicity issues can be avoided and how expectation values of arbitrary observables can be obtained within this framework. We then demonstrate the method using Compact U(1) Lattice Gauge Theory. A thorough study of the algorithm parameter dependence of the results is performed and compared with the analytically expected behaviour. We obtain high precision values for the critical coupling for the phase transition and for the peak value of the specific heat for lattice sizes ranging from $8^4$ to $20^4$. Our results perfectly agree with the reference values reported in the literature, which covers lattice sizes up to $18^4$. Robust results for the $20^4$ volume are obtained for the first time. This latter investigation, which, due to strong metastabilities developed at the pseudo-critical coupling, so far has been out of reach even on supercomputers with importance sampling approaches, has been performed to high accuracy with modest computational resources. Other situations where the method is expected to be superior to importance sampling techniques are pointed out.

연구 동기 및 목표

• 연속 자유도를 가진 시스템, 예를 들어 격자 양성거시학 이론에서 연속 상태 밀도를 계산하기 위한 히스토그램을 사용하지 않고도 효율적인 방법을 개발하는 것.
• 강한 준정적 안정성, 일阶 상전이 또는 거친 자유 에너지 표면을 가진 시스템에서 중요도 샘플링의 한계를 해결하는 것.
• LLR 프레임워크 내에서 관측량 기대값의 수렴성을 엄밀히 증명하여 이 방법의 이론적 기반에 빈틈을 메우는 것.
• 특히 표준 방법이 실패하는 큰 체적 영역에서, Compact U(1) 격자 양성거시학 이론에서 이 방법의 강건성과 정확성을 입증하는 것.
• $8^4$에서 $20^4$에 이르는 격자 크기 범위에서 임계 결합 상수와 비열 피크 값에 대한 고정밀도 결과를 제공하며, $20^4$에서의 첫 성공적 계산을 포함한다.

제안 방법

• LLR 알고리즘은 에너지 히스토그램 없이 로그 선형 완화 기반으로 상태 밀도 $\rho(E)$를 반복적으로 개선한다.
• 현재 $\rho(E)$ 추정치의 역수 비례 가중치를 가진 확률적 업데이트 절차를 사용하여 에너지 표면의 균일한 탐색을 보장한다.
• 기준 척도와 체적 스케일링 분석을 통해 역온도 $a_k = \partial \ln \rho / \partial E$가 체적에 영향을 받지 않음을 확인함으로써 열역학적 일관성을 검증한다.
• 모멘트 전개 기법을 적용하여 $\rho(E)$의 체적 스케일링을 유도하며, 모멘트가 $v = V / \xi^4$에 대해 선형으로 스케일링됨을 보여, 이 방법이 열역학적 기대와 일치함을 확인한다.
• 관측량 기대값은 $\rho(E)$에 대한 수치적 적분을 통해 계산되며, 반복적 개선 과정에서 지수 수렴을 통해 오차 억제가 달성된다.
• 특히 일阶 상전이 근처에서 준정적 상태에 갇히는 것을 방지하기 위해 샘플링 가중치를 동적으로 조정함으로써 에르고딕성을 보장한다.

실험 결과

연구 질문

• RQ1 LLR 알고리즘이 연속 시스템에서 몇 개의 주요 스케일에 걸쳐 거의 일정한 상대 오차로 상태 밀도를 고정밀도로 추정할 수 있는가?
• RQ2 LLR 방법은 Lattice Gauge Theories의 일阶 상전이와 같은 강한 준정적 안정성을 가진 시스템에서 어떻게 성능을 발휘하는가?
• RQ3 LLR 방법은 지수적 오차 억제를 보이며, 큰 격자에서 계산 효율성 면에서 중요도 샘플링을 능가할 수 있는가?
• RQ4 상태 밀도와 그 도함수의 체적 스케일링 행동은 어떻게 되며, 열역학적 기대와 일치하는가?
• RQ5 중요도 샘플링이 준정성으로 인해 실패하는 $20^4$ 격자 크기에서 LLR 방법이 임계 결합 상수와 비열 피크 값을 신뢰성 있게 계산할 수 있는가?

주요 결과

• LLR 알고리즘은 지수적 오차 억제를 실현하여 상태 밀도의 몇 개의 주요 스케일에 걸쳐 거의 일정한 상대 정밀도를 달성한다.
• 이번에 처음으로 $20^4$ 격자 크기에서 임계 결합 상수와 비열 피크 값에 대한 강력한 결과를 확보하였으며, 이는 이전까지 중요도 샘플링이 실패한 강한 준정적 안정성 영역이었다.
• $8^4$에서 $20^4$에 이르는 격자 크기에서 계산된 임계 결합 상수와 비열 피크 값이 $18^4$까지 문헌 값과 완벽하게 일치하여 방법의 정확성을 검증한다.
• LLR 방법에서 상전이 근처에서 상관 시간은 체적에 대해 최대 제곱수준으로 증가하는 반면, 중요도 샘플링에서는 기대되는 지수적 증가와 대비하여 우수한 확장성(스케일러빌리티)을 보여준다.
• 역온도 $a_k = \partial \ln \rho / \partial E$가 체적에 영향을 받지 않음을 확인하여 이 방법의 열역학적 일관성을 뒷받침한다.
• 논문에서 엄밀히 증명한 바와 같이, 임의의 관측량 기대값은 정확한 수렴을 보이며, 이는 LLR 접근법의 이론적 기반에 있던 빈틈을 메우는 결과이다.