Jump-Diffusion Stochastic Quantization for Euclidean Lattice Field Theories

技術要約：ユークリッド格子場理論におけるジャンプ・ディフュージョン・ストキャスティック量子化

問題提起
本論文は、ストキャスティック量子化を用いたユークリッド量子場理論（QFT）の数値シミュレーションにおける根本的な失敗モード、すなわち、作用の障壁によるエルゴード性の喪失について論じている。従来のストキャスティック量子化では、場 $\phi$ は確率偏微分方程式によって支配される連続的な拡散過程（ランジュバン動力学）に従って進化する。この手法は局所的な観測量には効率的であるが、作用 $S[\phi]$ が明確に異なるセクター（例：トポロジカル・セクターや相共存領域）を隔てる高い障壁を持つ場合、苦戦することになる。

このようなシナリオでは、連続的な経路が障壁を横切る確率は $\exp[-\Delta S_{\text{barrier}}]$ に比例してスケールする。システムが連続極限や臨界点に近づくにつれ、これらの障壁は現実的なシミュレーション時間内では乗り越えられないものとなる。これは「トポロジカル・フリーズ（topological freezing）」を引き起こす。サンプラーは単一のセクターに閉じ込められ、収束しているように見えても、正しいアンサンブル平均を得ることができなくなる。既存の解決策（テンパリング、マルチカノニカル法、または幾何学的戦略である開境界など）は、多くの場合、大幅なシステムの変更や複雑な再重み付けを必要とする。

手法
著者は、基礎となる確率過程を純粋な拡散から**ジャンプ・ディフュージョン（跳躍拡散）**へと拡張することにより、ストキャスティック量子化の自然な一般化を提案している。このフレームワークは、非連続的な経路、具体的にはレヴィ・フライト（Lévy flights）を取り入れるものであり、場 $\phi$ が微小な拡散ステップに加えて有限のジャンプを行うことを可能にする。

理論的基礎は、クルーグ（Courrège）の定理とレヴィ・イトー分解（Lévy-Itô decomposition）に基づいている。これらは、最も一般的な連続時間マルコフ過程は、決定論的なドリフト、ガウス型拡散、およびジャンプの3つの成分で構成されることを示している。確率分布 $P(\phi, t)$ の進化は、一般化されたチャップマン・コルモゴロフ（CK）方程式、あるいは一般化されたフォッカー・プランク方程式によって記述される：
$$\partial_t P = \underbrace{-\partial_\phi [b P] + \partial_\phi^2 P}_{\text{拡散}} + \underbrace{\int d\phi' [W(\phi|\phi')P(\phi') - W(\phi'|\phi)P(\phi)]}_{\text{ジャンプ}}$$
ここで $W(\phi'|\phi)$ はジャンプの遷移レート密度である。ターゲットとなる分布 $e^{-S[\phi]}$ への収束を保証するために、ジャンプ項は詳細釣合い条件を満たさなければならない。

論文では、これらのジャンプ更新を構築するための3つの具体的な戦略を概説している：

1. 加法的プロポーザル（クラスA）： 固定された対称なプロポーザル密度 $q_0(\eta)$ を用いたメトロポリス受容ステップを組み合わせる。これは、標準的なメトロポリス・ヘイスティングス法を、ストキャスティック量子化の枠組みの中でジャンプとして適用したものに相当する。
2. 可逆写像（クラスB）： 測度を保存する（ヤコビアン $J_T = 1$）決定論的かつ可逆な変換 $T$（例：グローバルなシフトや回転）を用いる。これらは、作用の変化 $\Delta S$ に基づくメトロポリス因子によって受理される。
3. 状態依存型の写像ファミリー（クラスC）： 離散系の文献から着想を得た戦略であり、パラメータ $m$（例：位置やサイズ）によってパラメータ化された可逆写像のファミリー $\{T_m\}$ を用いる。一様な選択ではなく、パラメータ $m$ が現在の状態に依存して確率 $p(m|\phi)$ で選択される。論文では、「局所的にバランスの取れた」重みが導出されており、そこでは選択確率が $e^{-\Delta S_m/2}$ に比例する。このアプローチは、作用のコストを順方向と逆方向の移動の間に効果的に分配し、困難な遷移に対する受理率を大幅に向上させる。

これらのジャンプの効率性は、受理推定値の導出によって分析される。測度保存型の対称ジャンプの場合、受理率は $erfc(\sigma/ \sqrt{8})$ となる（ここで $\sigma$ は作用コスト $\Delta S$ の標準偏差である）。状態依存型のジャンプの場合、関連するコスト変数は $\Delta S_{\text{eff}} = \ln(Z'/Z)$ となる（ここで $Z$ は写像のファミリーにわたる正規化因子である）。

主な結果
本手法は、2つのベンチマーク・システムでテストされた：

1. 傾斜二重井戸（Tilted Double Well）：

   • セットアップ： 2つの局所解を障壁で隔てるポテンシャル $S(x) = \gamma(x^2-1)^2 + cx$ を持つ1次元系。
   • 結果： 片方の井戸に初期化された標準的なランジュバン動力学は、トラップされたままとなり、正しいグローバル平均のサンプリングに失敗した（$\langle x \rangle_{\text{diff}} \approx 0.985$ に対し、真の値は $\approx -0.765$）。「反転（flip）」写像にガウスノイズを加えたジャンプ・ディフュージョン・アプローチは、障壁を頻繁に横切り、正しい期待値を回復することに成功した（$\langle x \rangle_{\text{jumps}} \approx -0.755$）。

2. 2D $U(1)$ ゲージ理論：

   • セットアップ： トラス上のコンパクト $U(1)$ ゲージ理論。これは連続極限（大きな $\beta$）付近でのトポロジカル・フリーズで知られている。トポロジカル電荷 $Q$ は、コスト $\sim 2\beta$ の障壁を介して変化する。
   • グローバル更新： 格子全体にフラックスを均一に広げることで、トポロジカル電荷を $Q=\pm 1$ シフトさせる「定数フラックス」ジャンプが構築された。この更新は、標準的なランジュバン動力学が完全に凍結（遷移ゼロ）していた $\beta=12.5$ においても、トポロジカル・セクター間の遷移を維持することに成功した。ジャンプ・ディフュージョン・サンプラーは、トポロジカル感受率 $\chi_t$ およびプラケット平均の厳密な解析結果を再現した。
   • 局所（ワインディング）更新： マター場が存在する現実的な理論においてグローバル更新がコスト高になることを考慮し、「ワインディング（winding）」更新（$2\pi$ の捻れを伴う局所ゲージ変換）がテストされた。
     • 非情報型（Uninformed）： ランダムに配置されたワインディング更新は、ワインディングのサイズ $L_w$ が減少するにつれて受理率が低下し、$L_w$ が小さい場合にフリーズが発生した。
     • 状態依存型（State-Informed）： クラスCの戦略（局所的な作用コストに基づいてワインディングの位置を選択する）を用いることで、受理率が劇的に改善された。驚くべきことに、状態依存型選択を用いた最小のワインディング（$L_w=2$）は、非情報型のグローバル更新よりも高い効率を達成し、正しいトポロジカル・サンプリングを回復した。

意義と主張
本論文は、ジャンプを組み込むことでエルゴード性崩壊を克服できる、ストキャスティック量子化の完全なマルコフ的一般化を提供すると主張している。

• エルゴード性の回復： 主要な貢献は、ジャンプ・ディフュージョン過程が、従来の拡散が失敗する領域（例：トポロジカル・フリーズ）において、作用の風景にある障壁を効率的に探索し、エルゴード性を回復できることを示したことである。
• 理論的枠組み： 一般化されたフォッカー・プランク方程式を通じて、詳細釣合いを満たすジャンプ更新を設計するための厳密な理論的基礎を確立した。
• 状態依存型選択による効率化： ジャンプのパラメータを状態依存の選択問題として扱うこと（クラスC）が、個々のジャンプ写像の高いコストを軽減できるという重要な洞察を与えた。これにより、マター場を持つ理論において極めて重要となる、局所的な更新であっても効率的なサンプリングが可能になる。
• 今後の展望： 著者は、本研究をより複雑な理論への足掛かりとして位置づけている。このフレームワークは、$CP(N-1)$ モデル（ゲージ場とマター場を含む）への適用、そして最終的にはフルQCDへの適用を意図している。また、複雑な平面における特異点を回避するための複素ランジュバン・シミュレーションや、拡散が定義されていない離散系への応用可能性についても言及している。

本研究は、すべての幾何学において直ちに臨界減速（critical slowing down）を解決すると主張するものではなく、ジャンプとカーネルの自由度（例：フーリエ加速）を組み合わせることが有望な将来の方向性であることを示唆している。結果は、ジャンプ・ディフュージョンが格子場理論におけるエルゴード性回復において有効であることの概念実証として提示されている。