[論文レビュー] An efficient algorithm for numerical computations of continuous densities of states

Kurt Langfeld, Biagio Lucini+2 more|arXiv (Cornell University)|Sep 28, 2015

Theoretical and Computational Physics参考文献 21被引用数 7

ひとこと要約

本稿では、連続的な自由度を持つ系、たとえば格子ゲージ理論のような系における連続的な状態密度の効率的計算のため、改善された対数線形緩和（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フレームワーク内での観測量期待値の収束を形式的に証明し、この手法の理論的基盤における空白を埋める。
• 特に、標準的手法が失敗する大規模な体積領域においても、強靭性と正確性を示す、コンパクトU(1)格子ゲージ理論への適用の検証。
• $8^4$ から $20^4$ の格子サイズにわたり、臨界結合定数および比熱ピーク値の高精度な結果を提供すること、特に$20^4$ での成功した計算を初実行すること。

提案手法

• LLRアルゴリズムは、エネルギーヒストグラムを必要とせず、対数線形緩和スキームを用いて状態密度 $\rho(E)$ を反復的に改善する。
• エネルギー状態を、現在の $\rho(E)$ の推定値の逆数に比例する重みでサンプリングする確率的更新手順を採用し、エネルギー空間の均一な探索を保証する。
• 基準スケールと体積スケーリング解析を用いて、逆温度 $a_k = \partial \ln \rho / \partial E$ が体積に依存しないことを確認し、熱力学的整合性を裏付ける。
• モーメント展開技術を適用して、$\rho(E)$ の系の体積 $V$ によるスケーリングを導出し、モーメントが $v = V / \xi^4$ に線形にスケーリングすることを示し、この手法が熱力学的期待と整合していることを確認する。
• 観測量の期待値は、$\rho(E)$ における数値積分によって計算され、反復的改善プロセスにおける指数的収束によって誤差抑制が達成される。
• 特に一次相転移付近で顕著なかく乱状態に陥りやすい状況においても、サンプリング重みを動的に調整することでエルゴード性を確保し、効果的にメタ安定状態からの脱出を実現する。

実験結果

リサーチクエスチョン

• RQ1 LLRアルゴリズムは、連続系において、数個のオーダーにわたり、ほぼ一定の相対誤差を維持しながら、状態密度の高精度推定を達成できるか？
• RQ2 LLR手法は、格子ゲージ理論における一次相転移のような強いかく乱状態を示す系において、どのように性能を発揮するか？
• 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アプローチの理論的基盤における空白を解消した。