[論文レビュー] The generalized fermion-bag approach

Shailesh Chandrasekharan, Anyi Li|arXiv (Cornell University)|Nov 22, 2011

Markov Chains and Monte Carlo Methods参考文献 4被引用数 4

ひとこと要約

本論文は、一般化されたフェルミオン・バッグ手法を導入し、四フェルミオン格子場理論における符号問題を解消する新しいモンテカルロ法を提案する。この手法は、符号問題を引き起こす未対応のフェルミオン自由度を「バッグ」に分離することで、質量ゼロのフェルミオンを正確に含む3次元質量ゼロのチラリング模型の、大規模シミュレーションを初めて可能にした。強-弱い結合双対性を活用することで、強い結合および弱い結合の両方で効率的であり、$40^3$ 格子上での正確な質量ゼロのフェルミオンを用いた大規模なシミュレーションが可能となり、臨界指数と量子臨界点 $U_c = 0.2608(2)$ を高精度で得た。

ABSTRACT

We present a new approach to some four-fermion lattice field theories which we call the generalized fermion bag approach. The basic idea is to identify unpaired fermionic degrees of freedom that cause sign problems and collect them in a bag. Paired fermions usually act like bosons and do not lead to sign problems. A resummation of all unpaired fermion degrees of freedom inside the bag is sufficient to solve the fermion sign problem in a variety of interesting cases. Using a concept of duality we then argue that the size of the fermion bags is small both at strong and weak couplings. This allows us to construct efficient algorithms in both these limits. Using the fermion bag approach, we study the quantum phase transition of the 3D massless lattice Thirrring model which is of interest in the context of Graphene. Using our method we are able to solve the model on lattices as large as $40^3$ with moderate computational resources. We obtain the precise location of the quantum critical point and the values of the critical exponents through this study.

研究の動機と目的

• 正確に質量ゼロのフェルミオンを含む四フェルミオン格子場理論をシミュレートできる新しいモンテカルロ法の開発。特に、ハイブリッド・モンテカルロ（HMC）などの従来手法が長年にわたり直面する符号問題を克服することを目的とする。
• 強相関フェルミオン系における量子臨界点のシミュレーションという課題に取り組む。特に、グラフェンなどの凝縮系物理学に関連するモデルを対象とする。
• フェルミオン・バッグのサイズが系の体積に比べて小さくなるように保証する双対性を利用することで、強い結合および弱い結合の両方で効率的なシミュレーションを実現すること。
• 大規模な格子上での有限サイズスケーリングを用いて、3次元質量ゼロのチラリング模型における臨界指数と量子臨界点の正確な位置 $U_c$ を算出すること。

提案手法

• 符号問題を引き起こす未対応のフェルミオン自由度を特定し、局所化された「フェルミオン・バッグ」に集約する。一方、ペアになったフェルミオンは有効的にボソン的であり、符号問題に寄与しない。
• 各バッグ内に存在するすべてのフェルミオン自由度の再結合（リサム）を実行することで、正の重みが得られ、効率的なマルコフ連鎖更新が可能になる。
• この手法は、ディラック行列 $D$、逆伝播関数行列 $G$、およびバッグに含まれないサイトを除いた縮小行列 $W$ を含む行列式の恒等式に基づく強-弱い結合双対性に依存している。
• フェルミオン・バッグのサイズは、臨界点でさえも通常 $\sim V/8$ の程度であり、系全体の体積に比べて非常に小さいことが示された。これは計算の効率性を保証する。
• 局所的更新は、バッグ内のフェルミオン自由度数 $N_b$ に対して $O(N_b^3)$ のスケーリングを示すが、全系の体積に依存しないため、高速なシミュレーションが可能になる。
• 臨界指数と臨界結合定数 $U_c$ を抽出するために、観測量 $\chi^{-1}L^{2-\eta}$、$\langle q_\chi^2 \rangle$、$R_f L^{2+\eta_\psi}$ に対して有限サイズスケーリング解析を適用する。

実験結果

リサーチクエスチョン

• RQ1 フェルミオン・バッグ手法は、正確に質量ゼロのフェルミオンを含む四フェルミオン格子場理論、特に3次元質量ゼロのチラリング模型において、符号問題を完全に解消できるか？
• RQ2 フェルミオン・バッグのサイズは結合強度にどのように依存するか？また、強い結合および弱い結合の両方で小さいままであるか？
• RQ3 量子臨界点においても、最小限の臨界遅れ（critical slowing down）を伴う効率的なサンプリングが可能か？
• RQ4 3次元質量ゼロのチラリング模型の量子臨界点における臨界指数 $\nu$、$\eta$、$\eta_\psi$ の正確な値は何か？
• RQ5 本手法を用いて正確に質量ゼロのフェルミオンを含む大規模シミュレーション（$40^3$）が可能か？また、質量のあるフェルミオンを用いた従来のHMCベースの研究と比較して、結果はどうなるか？

主な発見

• フェルミオン・バッグ手法は、3次元質量ゼロのチラリング模型における符号問題を効果的に解消し、中程度の計算リソースで $40^3$ 格子上での正確な質量ゼロのフェルミオンを含むシミュレーションを可能にした。
• 臨界結合定数は $U_c = 0.2608(2)$ と高精度に決定され、結合フィットにおける $\chi^2/\text{d.o.f.} = 1.3$ を達成した。
• 臨界指数は $\nu = 0.85(1)$、$\eta = 0.65(1)$、$\eta_\psi = 0.37(1)$ であり、普遍的スケーリング行動と整合的である。
• 動的臨界指数 $z \approx 0$ が観測され、スイープ単位で測定しても臨界遅れが見られないことが示された。これは $O(N_b^3)$ のスケーリングにもかかわらず、計算の高速性を示している。
• 量子臨界点における平均フェルミオン・バッグサイズは、約 $V/8$ であり、系全体の体積に比べて非常に小さいことが確認された。
• 従来のHMC手法は、符号問題と小さな固有値特異性のため、大規模格子上での正確な質量ゼロのフェルミオンのシミュレーションに失敗しているが、本手法はそれらを凌駒した。