[論文レビュー] Dual simulation of finite density lattice QED at large mass

Michael Kniely, Christof Gattringer|arXiv (Cornell University)|Feb 3, 2015

Quantum Chromodynamics and Particle Interactions参考文献 1被引用数 4

ひとこと要約

本稿では、大質量における2つのdegenerateなウィルソンフェルミオン流体を有する有限密度格子QEDの二重モンテカルロシミュレーションを、ヘッジング展開を要約することでフェルミオンループを導出し、ゲージ場を整数値のプラケット占有数に双対化することにより提示する。主な貢献は、正の総符号を持つ準平面的ループのサブクラスを同定することであり、これにより符号問題のないシミュレーションが可能となり、無限大のループ縮退のおかげで凝縮相における粒子密度が飽和しないことが明らかになった。

ABSTRACT

We discuss a mapping of lattice QED with two flavors and a chemical potential to dual variables, which are surfaces for the gauge fields and loops for the fermions. The gauge fields are completely dualized and the corresponding dual variables are integer valued plaquette occupation numbers with constraints that lead to a structure of surfaces that are either closed or bounded by fermion loops. The fermion loops are obtained from a resummed hopping expansion (large mass expansion) of the determinant of the Wilson-Dirac operator. The loops can come with both positive and negative signs. We identify a sub-class of loops, which we refer to as quasi-planar loops, where the total sign is positive. For this sub-class a dual Monte Carlo simulation is possible and we discuss its implementation and some results. In particular we address condensation phenomena at finite chemical potential.

研究の動機と目的

• ゲージおよびフェルミオン自由度を両方双対化することで、有限密度格子QEDにおける複素作用問題を克服すること。
• フェルミオン行列式の大きな質量（ヘッジング）展開から生じるフェルミオンループの構造を分析すること。
• パウリ排他的性とトレース構造からの総符号が正であるサブクラスとしての準平面的ループを同定し、符号問題のないモンテカルロシミュレーションを可能にすること。
• 双対表現を用いて有限化学ポテンシャルにおける凝縮現象を研究すること。
• ヘッジング展開の切り捨てが、凝縮相における粒子密度の飽和の出現に与える影響を調査すること。

提案手法

• プラケット占有数を用いてゲージ場を完全に双対化し、整数制約を満たす閉じた表面またはフェルミオンに接続された表面が得られる。
• ウィルソン＝ディラック作用素の行列式に要約されたヘッジング展開を適用し、後退しない・繰り返しのないループが生成され、関連する重みと符号が得られる。
• パウリ排他的性とトレース構造からの総符号が正である準平面的ループをサブクラスとして同定し、正定値モンテカルロサンプリングを可能にする。
• プラケット占有数の変更を用いて適格性を維持しながら、ループ移動、空間的移動、局所的変形の更新を含む双対メトロポリスアルゴリズムを実装する。
• 双対表現を用いて、化学ポテンシャルによる微分を用いて粒子数密度などの物理量を計算する。
• $N_S^3 \times N_T$ 格子上でシミュレーションを実施し、$N_S, N_T$ は4から16の範囲で変動させ、平衡化に $O(10^4)$ スイープ、測定に $O(10^4)$ 回の測定を行い、ジャックナイフ誤差推定を用いる。

実験結果

リサーチクエスチョン

• RQ1 ゲージおよびフェルミオン自由度を両方双対化した場合、有限密度格子QEDの完全な双対表現を構築できるか？
• RQ2 フェルミオンループの符号が双対形式において果たす役割は何か？ また、総符号が正であるサブセットのループを同定できるか？
• RQ3 双対ゲージ表面の制約とフェルミオンループ統計が符号問題に与える影響は何か？
• RQ4 ループクラスの選択（例：線形対比して準平面的）が、有限化学ポテンシャルにおける凝縮現象の出現に与える影響は何か？
• RQ5 ヘッジング展開の切り捨てが、凝縮相における粒子数密度の飽和にどの程度影響を及えるか？

主な発見

• 線形ループモデルでは、凝縮相において粒子数密度が $n = 2$ に飽和する。これは、各ループにおける最大二重占有と整合的である。
• 準平面的ループモデルでは、粒子数密度の飽和は観察されず、これはこのようなループの無限大の縮退が、高化学ポテンシャルでも完全な充填を妨げるためである。
• 凝縮転移は、より高い温度（$N_T$ が小さい）でなだらかに洗練され、これはクロスオーバー転移の予想に一致する。
• 粒子数密度は、$\beta = 0.9$ の閉じ込め相および $\beta = 1.1$ のクーロン相の両方で、化学ポテンシャルに単調に増加するが、転移が発生する $\mu$ 値は異なる。
• ヘッジング展開における微妙なキャンセルが重要である：準平面的モデルにおける飽和の欠如は、物理的上限 $2N_S^3$ があるにもかかわらず、このようなループの無限大の数に起因する。
• 双対シミュレーションフレームワークは、準平面的ループサブクラスにおいて符号問題なしに凝縮の物理を正しく捉えており、今後の有限密度QEDの研究に有効であることが妥当化された。