技術要約：格子場理論における相転移の数値的研究

本論文は、格子場理論における臨界現象の解明と、シミュレーションの信頼性確保を目的とした、二つの補完的な数値的研究を提示するものである。本研究は、微細な臨界挙動を識別するための高精度なツールの需要と、モンテカルロ・サンプリングの熱力学的整合性を検証する必要性に正面から取り組んでいる。研究は、二次元一般化XY（gXY）モデルに対するテンソルネットワーク解析と、コンパクトU(1)格子ゲージ理論のための構成温度推定器の開発という、二つの異なる部分に分かれている。

第I部：一般化XYモデルのテンソル繰り込み研究

問題と背景
Korshunov、Lee、Grinsteinによって提案された二次元一般化XY（gXY）モデルは、変形パラメータ $\Delta$ を通じて異なる相の間を補間する。このモデルは、整数渦（強磁性相）と半整数渦（ネマティック相）の間の競争を特徴とし、これらは無秩序な常磁性相によって隔てられている。相構造には、整数および半整数渦の両方に関するベレジンスキー・コステリッツ・ザウス（BKT）転移と、その間に介在するイジング型転移が含まれる。モンテカルロ法や行列積状態（MPS）法を用いた先行研究では、これらの相の存在は確立されているものの、相境界の精密なマッピングや、転移線が交わる多重臨界領域の性質の特定には課題が残されていた。

手法
著者は、大きなボンド次元を扱うためにGPU加速を実装した高次テンソル繰り込み群（HOTRG）法を採用している。そのアプローチは以下の通りである：

1. 特性展開（Character Expansion）： gXYモデルのボルツマン重みを双対変数（特性展開）を用いて展開し、分配関数をテンソルネットワークの縮約へと変換する。
2. HOTRGによる粗視化： ネットワークは、高次特異値分解（HOSVD）を用いて反復的に粗視化される。これにより、複数の方向に対して同時に支配的なもつれ（エンタングルメント）の寄与を保持しながら、ボンド次元を切り詰めることが可能となり、標準的なTRGよりも高い精度が得られる。
3. 不純物テンソル技術： 磁化や感受性などの熱力学的観測量を抽出するために、外部対称性の破れを記述する不純物テンソル（外部場 $h$ および $h_1$）をネットワーク内に挿入する。これにより、数値的なノイズを導入する可能性のある有限差分近似に頼ることなく、自由度の微分を計算できる。
4. 外挿： 二次元においてはマーミン・ワグナーの定理により外部場なしでの自発的対称性の破れが禁じられているため、臨界温度は感受性のピークをゼロ磁場極限へと外挿することによって決定される。

主な結果

• 相図のマッピング： 本研究は、$(\Delta, T)$ 平面における相図の精密なマッピングを提供する。強磁性（整数渦）、ネマティック（半整数渦）、および無秩序相の三つの異なる相を特定した。
• 転移線： これらの境界は、二つのBKT転移（整数および半整数）と、一つのイジング転移によって特徴付けられる。
• 多重臨界点： 数値結果は、転移線の逐次的な合流シナリオを示唆している。半整数BKT線とイジング線はまず $\Delta \approx 0.34$ 付近で合流し、続いてBKT線との合流は $\Delta \approx 0.36(2)$ で起こる。多重臨界点は $\Delta = 0.36(2)$ および $T_c = 0.716(3)$ に位置する。
• 推定器の比較： BKT型の転移においては、比熱のピークは臨界温度を系統的に過大評価することが示された。対照的に、磁化率のゼロ磁場外挿は、より信頼性の高い推定器となる。
• パフォーマンス： GPU加速されたHOTRGの実装は、ボンド次元 $D=91$ およびフーリエカットオフ $m \in [-50, 50]$ を用いて結果を計算することに成功し、熱力学的極限において高い精度を達成した。

意義
本研究は、注意深く制御されたGPU加速HOTRGが、複雑な相構造を解明するための非常に正確かつ効率的な枠組みを提供し、従来のモンテカルロ法やMPS研究を改善できることを実証している。これは、gXYモデルにおける多重臨界領域の決定的な数値的特性付けを与えるものである。

---

第II部：格子ゲージ理論における構成温度推定器

問題と背景
標準的な正準モンテカルロ・シミュレーション、特にハイブリッドモンテカルロ（HMC）においては、アルゴリズムの正当性を検証するための直接的な保存則が不足している。分子動力学ではエネルギーのような保存量が用いられるが、正準アンサンブルでは、未知のシステムや参照データのない状態点において、観測量を既知の熱力学的特性と比較することに困難が伴う。さらに、格子ゲージ理論における離散化効果やアルゴリズムによるアーティファクトは、熱力学的観測量や臨界温度を歪ませる可能性がある。

手法
本論文では、Rughの統計力学の幾何学的定式化に基づく構成温度推定器を開発し、ベンチマークを行った。

1. 理論的基礎： この推定器は、運動量変数とは独立して、ユークリッド作用（$S$）の勾配とヘッセ行列から温度を導出する。逆構成温度 $\beta_{config}$ は以下のように定義される：
   $$\beta_{config} = \left\langle \nabla_\phi \cdot \left( \frac{\nabla_\phi S}{|\nabla_\phi S|^2} \right) \right\rangle$$
   ここで $\nabla_\phi$ はすべての格子場変数に関する勾配である。
2. ゲージ不変性： 共変微分を用いることで、ゲージ不変となるよう定式化されており、コンパクトゲージ理論への適用が可能である。
3. 診断への応用： 確率分布が $P[\phi] \propto e^{-S[\phi]}$ となる正しくサンプリングされたアンサンブルにおいて、推定器は $\beta_{config} = 1$ （自然単位系）を与えるはずである。この値からの逸脱は、サンプリングの不整合、熱化の遅れ、または数値的エラーを示す。
4. ベンチマーク： 推定器は、一、二、および四次元のユークリッド空間におけるコンパクトU(1)格子ゲージ理論を用いてテストされた。
   • 1Dおよび2D： これらは既知の解析解を持つ制御されたテストベッドとして機能する（1D/2D U(1)には相転移が存在しない）。
   • 4D： これは $\beta \approx 1.0$ 付近で弱い一次の閉じ込めークーロン相転移を伴う、非自明なベンチマークを提供する。
   • 近似： 大体積（特に4D）の場合、非対角のヘッセ行列項は無視できることが示されており、計算コストを削減するために対角項のみを用いた簡略化された推定器を使用できる。

主な結果

• 検証： 1Dおよび2DのU(1)理論において、測定された構成温度 $\beta_M$ は入力結合 $\beta$ （または $\beta_{config} \approx 1$）と一致しており、実装の妥当性と推定器の安定性が検証された。
• エラー検出： 推定器はアルゴリズムのエラーを検出することに成功した。メトロポリスの受理・拒否ステップに意図的な欠陥（誤った乱数範囲の使用）を設けたテストにおいて、標準的な観測量は問題を即座に明らかにできなかったが、$\beta_M$ は入力値から大きく逸脱した。
• 熱化のモニタリング： 推定器は熱化挙動を追跡し、熱平衡化の遅れやメタスタビリティの早期指標を提供する。
• 相転移に対する非感受性： 秩序パラメータ（例：プラケット感受率）とは異なり、構成温度は相転移において特異な挙動を示さず、臨界領域を通って滑らかに推移する。これは、これが秩序パラメータではなく、一貫性チェックとしての役割を担っていることを裏付けている。

意義
構成温度推定器は、格子場理論における熱的サンプリングを検証するための、ゲージ不変で運動量に依存しない診断ツールを提供する。これは、大規模なプロダクション・ランの前に、サンプリングの病理、数値的エラー、および遅い熱化を検出するための独立した手段を提供する。本研究は、大きな格子体積においては、精度を大きく損なうことなく対角ヘッセ行列項を無視することで、効率的に推定器を実装できることを示している。

---

全体の結論

本論文は、これら二つのアプローチが、臨界現象の数値的研究における補完的な課題に対処していると結論付けている：

1. テンソルネットワーク（HOTRG）： 競合する秩序を持つ系の複雑な相図と臨界挙動をマッピングするための、決定論的で高精度な経路を提供し、熱力学的極限を効果的に扱う。
2. 構成温度： モンテカルロ・シミュレーションにおける確率的サンプリングの信頼性を検証するための、低オーバーヘッドで独立した診断を提供する。

今後の方向性として、より高い分数化（$q \ge 3$）を持つ一般化XYモデルや、結合されたゲージ自由度へのHOTRGの適用、および非アーベルゲージ理論や有限温度QCDにおける大規模なハイブリッドモンテカルロ・シミュレーションへの構成温度モニターの統合が挙げられる。著者は、これらの手法は格子場理論に限定されるものではなく、精密な数値制御と信頼できる診断が不可欠な凝縮系物理学においても価値を持つと述べている。