[論文レビュー] Remarks on the construction of worm algorithms for lattice field theories in worldline representation

Mario Giuliani, Christof Gattringer|arXiv (Cornell University)|Feb 15, 2017

Scientific Research and Discoveries参考文献 13被引用数 3

ひとこと要約

本稿では、世界線表現における格子場理論に対して一般化されたウジ虫アルゴリズムを提示する。このアルゴリズムでは、各サイトにおける重み係数がそのサイトのすべてのフラックスに依存する。アーマチュアパラメータ A を導入し、ウジ虫の開始および終了確率を制御することで、ウジ虫の長さとシミュレーション効率の最適化を可能にする。この手法は詳細つり合いを満たし、4次元相対論的ボーズガスにおいて数値的に妥当性が確認された。A は平均ウジ虫長と1ステップあたりのフラックス変化に顕著な影響を与えることが示された。

ABSTRACT

We introduce a generalized worldline model where the partition function is a sum over configurations of a conserved flux on a d-dimensional lattice. The weights for the configurations of the corresponding worldlines have factors living on the links of the lattice, as well as terms which live on the sites x and depend on all fluxes attached to x. The model represents a general class of worldline systems, among them the dual representation of the relativistic Bose gas at finite density. We construct a suitable worm algorithm and show how to correctly distribute the site weights in the various Metropolis probabilities that determine the worm. We analyze the algorithm in detail and give a proof of detailed balance. Our algorithm admits the introduction of an amplitude parameter A that can be chosen freely. Using a numerical simulation of the relativistic Bose gas we demonstrate that A allows one to influence the starting and terminating probabilities and thus the average length and the efficiency of the worm.

研究の動機と目的

• サイト依存重み係数を有する格子場理論におけるウジ虫アルゴリズムの理論的基盤の欠如に応えること。
• すべてのフラックスに依存するサイト重みを、メトロポリス受容確率に正しく組み込むウジ虫アルゴリズムを構築すること。
• ウジ虫の開始および終了確率を制御するアーマチュアパラメータ A を導入し、ウジ虫長とシミュレーション効率の最適化を可能にすること。
• リンク重みおよびサイト重みを有する一般世界線モデルにおいて、提案されたアルゴリズムの詳細つり合いを証明すること。
• 4次元相対論的ボーズガスの数値シミュレーションを通じて、アルゴリズムの性能と調整可能性を示すこと。

提案手法

• モデルは、d次元格子上での保存フラックスを有する世界線系を一般化し、リンク重みおよびすべてのフラックスに依存するサイト重みを含む。
• ウジ虫が格子を移動し、リンク上のフラックスを変化させ、メトロポリスステップを用いて重みを更新するウジ虫アルゴリズムを構築する。
• サイト重みはメトロポリス確率に配分され、ウジ虫の開始時には分母に、終了時には分子に現れることで、正しい統計的重みが保証される。
• アーマチュアパラメータ A を開始および終了確率に導入し、小さなまたは大きなサイト重みに対応するとともに、ウジ虫長の調整が可能になる。
• 遷移確率を個々のウジ虫ステップに分解し、各ステップとその逆ステップを一致させることで、詳細つり合いを証明する。これにより、正しい平衡分布が保証される。
• ウジ虫長、フラックス変化、ステップ効率に及ぼす A の影響を調査するため、4次元相対論的ボーズガスの数値シミュレーションを実施する。

実験結果

リサーチクエスチョン

• RQ1 すべてのフラックスに依存するサイト依存重み係数を、ウジ虫アルゴリズムのメトロポリス確率に正しく組み込む方法は何か？
• RQ2 アーマチュアパラメータ A はウジ虫の開始および終了確率をどのように制御し、ウジ虫長にどのような影響を与えるか？
• RQ3 複雑なサイト重みを有する状況下でも、提案されたウジ虫アルゴリズムが詳細つり合いを満たすことを証明できるか？
• RQ4 アーマチュアパラメータ A は1ウジ虫ステップあたりの平均フラックス変化および全体のシミュレーション効率にどのように影響を与えるか？
• RQ5 アーマチュアパラメータ A は、格子場理論の異なる結合定数領域において、性能最適化にどの程度利用可能か？

主な発見

• アーマチュアパラメータ A により、ウジ虫の開始および終了確率を独立に制御でき、平均ウジ虫長の調整が可能になる。
• A が増加するにつれて、1ステップあたりの正規化フラックス変化量（Δflux / N_steps）は減少する。これは、長いウジ虫が必ずしも1ステップあたりのフラックス変化を大きくするわけではないことを示している。
• 1リンクあたりの正規化フラックス変化量（Δflux / 4V）は A に従って増加する。これは、A が増加するにつれてより長いウジ虫が生成されることを確認している。
• 受理ステップ数および総フラックス変化量における有限サイズ効果は、ウジ虫が周期的境界をぐるりと回るためであると帰属づけられる。
• 提案されたアルゴリズムは、リンク重みおよびサイト重みを有する一般世界線モデルにおいて詳細つり合いを維持する。これは、その使用に理論的基盤を提供する。
• この手法は他の系およびウジ虫アルゴリズムの変種へも一般化可能であり、A を用いた性能最適化の体系的手段を提供する。