Multilevel Plaquette-Space Sampling for Lattice Gauge Theories with Local Constraint Solves

技術要約：格子ゲージ理論におけるマルチレベル・プラケット空間サンプリング

問題提起
格子ゲージ理論（LGT）は、粒子物理学および核物理学における第一原理計算に不可欠であり、ウィルソン作用によって定義される高次元の確率分布のサンプリングを必要とする。従来のマルコフ連鎖モンテカルロ（MCMC）法は、格子間隔が連続極限に近づく際、あるいは相関が長距離に及ぶ際（弱結合時）、クリティカル・スローイングダウン（臨界減速）やトポロジカルな凍結の問題に直面する。生成モデル（例：ノーマライジング・フロー）は、これらの問題を緩和するための有望な代替策として台頭しているが、既存のアプローチは通常、リンク空間においてサンプリングを行う。リンク空間では、作用は非局所的であり、弱結合に向かうにつれて相関が長距離化するため、生成モデルが変数間の必要な協調性を学習することが困難になる。

プラケット空間で直接サンプリングすることは、ウィルソン作用がプラケット変数において局所的であり、ゲージ冗長性が大部分除去されるため、より自然な表現を提供できる。しかし、一般的なプラケット構成は、必ずしも有効なリンク構成に対応するわけではない。有効な構成は、ビアンキ恒等式（幾何学的整合性制約）を満たす必要があり、これによりプラケット場は低次元の多様体に制限される。高次元において、これらの制約を満たすことは、計算コストの高い、かつ高度に結合した線形システムを解くことを意味する。

手法：マルチレベル・プラケット空間サンプラー（PSS）
著者らは、ビアンキ恒等式を厳密に満たしつつ、扱いやすいサンプリング密度を維持しながら、有効なプラケット構成を生成する**マルチレベル・プラケット空間サンプラー（PSS）を提案する。その核心となる革新は、グローバルに結合した制約問題を一連の局所的な解法へと変換する粗から密への分解（coarse-to-fine factorization）**である。

1. 変数の分解: この手法は、「自由な」プラケット変数と「決定された」プラケットを区別する。決定されたプラケットは、ビアンキ恒等式およびグローバルな閉鎖条件（ホロノミー）によって固定される。
2. マルチレベル構造:
   • 粗レベル（Coarse Level）: 少数のプラケット変数とホロノミーを生成し、一部のビアンキ恒等式を満たす。
   • 精細化レベル（Refinement Levels）: 次の各精細レベルにおいて、粗いプラケットはより細かいプラケットへと分解される。極めて重要な点は、粗いプラケットの値が、対応する精細レベルでのウィルソンループの値として保持されることである。
   • 局所制約解法: この保持戦略により、粗いレベルで満たされた制約は自動的に継承される。各精細化ステップにおいて、すべての決定されたプラケットは、高々4つの新たに生成された変数にのみ依存する。各ステップにおける制約問題のサイズは、全格子サイズに依存しない。
3. 制約補完写像 ($R_i$): 決定論的、単射的、かつ体積保存的な写像が、生成された自由変数と粗い場を、完全な有効プラケット構成へと変換する。この写像は、結果の密度が厳密であることを保証し、重要度再重み付けを可能にする。
4. 生成モデル: 各レベルの自由変数は、粗いレベルおよびホロノミーに条件付けられた条件付きノーマライジング・フロー（具体的には、U(1)に対する有理二次スプラインおよびSU(2)に対する等変フロー）を用いてサンプリングされる。

主な貢献

• プラケット空間サンプリングの一般化: 本研究は、プラケット空間サンプリングを、二次元U(1)理論から**四次元U(1)および二次元SU(2)**ゲージ理論へと拡張した。
• 局所制約を備えたマルチレベルサンプラー: ビアンキ恒等式を単純な局所的手順を用いて厳密に満たす階層的サンプラーを提案している。このアプローチは、グローバルな線形ソルバーの計算スケーリング問題を回避しながら、粗いレベルを通じて長距離相関を捉えることができる。
• 実証的検証: 本論文は、マルチレベルPSSが、長距離相関を特徴とする物理的に重要な領域において、シングルレベルおよびマルチレベルのリンク空間ベースラインを大幅に上回る性能を示すことを実証している。

結果
手法は以下の3つの設定で評価された：

1. 2D U(1)（格子サイズによるスケーリング）: 一定の物理量ラインにおいて、マルチレベルPSSは、リンク空間のベースラインよりも格子サイズ ($L$) に対して著しく優れたスケーリングを示した。$L=64$ において、リンク空間のフローは深刻な重み退化（有効サンプルサイズ、ESS $\approx 4 \times 10^{-5}$）に苦しんだが、マルチレベルPSSは 0.10 のESSを維持した。ESSの優位性は、$L=8$ での $2.0\times$ から、$L=32$ では $11.9\times$ へと拡大した。
2. 4D U(1)（弱結合）: リンク空間のフローが劣化する弱結合領域（質量のないクーロン相）において、マルチレベルPSSは増大する優位性を示した。PSSとリンク空間ベースラインの間のKLダイバージェンスの差は、大きな $\beta$ において2桁を超えた。
3. 2D SU(2)（相関長の増大）: 相関長 $\xi/a$ が増加する（格子サイズの約40%に達する最大6.6まで）につれ、リンク空間のフローは急速に効率を失った。シングルレベルのプラケットサンプラーはより堅牢であったが、やはり長距離相関において悪化した。これに対し、マルチレベルPSSは全領域を通じて効率を維持した。

すべてのケースにおいて、マルチレベルPSSは正確な観測量（平均プラケット、トポロジカル感受率）を算出し、高いESSを維持したが、ベースラインはしばicht分布を正しく捉えられなかったり、高い分散に苦しんだりすることが多かった。

意義と主張
本論文は、プラケット変数（作用を局所的にする）と、粗から密への階層構造（制約構造を局所的に保つ）の組み合わせが、困難な領域――2D U(1における連続極限への接近、4D U(1)における弱結合への移行、および2D SU(2)における相関長の増大――において、ますます有益になることを主張している。

著者らは、この戦略により、グローバルソルバーの計算スケーリングのペナルティを受けることなく、厳密な制約処理が可能になると強調している。彼らは、これらの結果が、より大きな四次元格子および非アーベルゲージ理論への手法の拡張を動機付けるものであると述べている。また、現在の実験は比較的小さな格子に限定されており、稀な欠陥構成（例：弱結合時のモノポール）については、依然としてアンダーサンプリングされている可能性があり、将来的な研究では、より表現力豊かなモデルや、離散的なトポロジカル自由度の明示的な取り扱いが必要であると、謙虚に限界を認めている。