[論文レビュー] A novel density of state method for complex action systems

Biagio Lucini, Kurt Langfeld|arXiv (Cornell University)|Nov 1, 2014

Quantum many-body systems参考文献 13被引用数 4

ひとこと要約

本論文は、複素作用を有する系、特に有限化学ポテンシャルにおける $ζ(3)$ スピン模型での符号問題に対処するため、LLRアルゴリズムに基づく新しい状態密度法を提案する。一般化された状態密度を計算し、半解析的統合を用いることで、極めて深刻な符号問題を正確に解消する。$μ=2$ においても、キャンセルが $10^{-16}$ に達するが、正確な双対アルゴリズム結果と一致する。

ABSTRACT

Recently, a new and efficient algorithm (the LLR method) has been proposed for computing densities of states in statistical systems and gauge theories. In this talk, we explore whether this novel density of states method can be applied to numerical computations of observables in systems for which the action is complex. To this purpose, we introduce a generalised density of states, in terms of which integrals of oscillating observables can be determined semi-analytically, and we define a strategy to compute it with the LLR method. As a case study, we apply these ideas to the Z(3) spin model at finite density, finding a remarkable agreement of our results for the phase twist with those obtained with the worm algorithm for all explored chemical potentials, including values for which there are cancellations over sixteen orders of magnitude. These findings open new perspectives for dealing with the sign problem on physically more relevant systems.

研究の動機と目的

• 複素作用を有する統計的およびゲージ理論における物理量の計算のための、信頼性の高い手法を開発すること。標準的なモンテカルロ法は符号問題のため失敗する。
• 状態密度計算のためのLLRアルゴリズムを、振動的で複素数値の経路積分を有する系へと拡張すること。
• 双対形式による正確な結果が得られる物理的に意味のあるモデル、すなわち有限密度における $ζ(3)$ スピン模型に対して、本手法の正確性と精度を検証すること。
• 一般化された状態密度が十分な精度で計算された場合、振動的物理量の半解析的統合が可能であることを示すこと。

提案手法

• 複素経路積分が1次元の振動的積分に簡略化されるように、一般化された状態密度を定義する。
• LLRアルゴリズムを用いて、多数のオーダーにわたる高精度な状態密度 $\rho(n)$ を数値的に計算する。
• 状態密度の対数を多項式補間によりモデル化する：$\ln\rho(n) = \sum_{k=0}^{p} c_k n^{2k}$。これにより $n \to -n$ の下での対称性が保持される。
• 位相ねじれ $p(\mu)$ などの振動的物理量は、$\rho(n)$ を用いて半解析的に計算され、複素重みの直接的なモンテカルロサンプルを回避する。
• 正規化された積分を用いる：$p(\mu) = \frac{1}{V} \frac{I_1}{I_2} = \frac{1}{V} \frac{\sum_n \rho(n) n \sin(\kappa\sqrt{3}\sinh(\mu)n)}{\sum_n \rho(n) \cos(\kappa\sqrt{3}\sinh(\mu)n)}$。これにより高精度な評価が可能になる。
• 結果は、符号問題を有しない正確な双対アルゴリズムシミュレーションと照合される。

実験結果

リサーチクエスチョン

• RQ1 標準的な重要度サンプリングが失敗する複素作用系に、LLRに基づく状態密度法を拡張できるか？
• RQ2 数値的に決定された状態密度に基づく半解析的アプローチにより、振動的物理量はどの程度正確に計算できるか？
• RQ3 符号問題における極めて深刻なキャンセル、例えば $10^{-16}$ に達する場合、この手法はどの程度解消できるか？
• RQ4 強い符号問題の領域においても、$\ln\rho(n)$ の多項式アンザッツが安定的かつ正確に保たれるか？
• RQ5 位相ねじれなどの物理量に関して、本手法は正確なアルゴリズムと比較して、精度と信頼性において優れているか？

主な発見

• 本手法は、$\mathbb{Z}(3)$ スピン模型における位相ねじれ $p(\mu)$ を、$μ=2$ においても正確な双対アルゴリズムの統計誤差内に一致する精度で計算できた。このとき符号問題は極めて深刻である。
• $μ=2$ において、重なり因子 $O(\mu) \simeq 10^{-16}$ であり、極めて深刻なキャンセルが生じるが、本手法は高精度でこれを解消できた。
• $2p=8$ 項までの多項式補間において、$c_0$ と $c_2$ のみが顕著に非ゼロであったため、アンザッツの妥当性が確認された。
• 位相ねじれの分子 $I_1$ と分母 $I_2$ は、$μ=2$ で $10^{-16} - 10^{-19}$ にまで低下するが、その比は安定的かつ物理的に意味のある値を保った。
• 直接 $\rho(n)$ から $O(\mu)$ を再構築する手法は、強い符号問題では十分な精度を欠くが、本手法はそれを上回る性能を示した。
• 半解析的手法と、符号問題のない双対アルゴリズムを用いた直接モンテカルロシミュレーションとの間で良好な一致が観察され、本手法の有効性が全範囲の $\mu$ で裏付けられた。