Resurrecting the coherent state variational algorithm for large $N$ gauge theories

本論文は、フェルミオンの有無にかかわらずSU($N$)格子ゲージ理論に適用可能な新しい実装を導入することにより、大$N$ゲージ理論に対するコヒーレント状態変分法の実現可能性を再評価し、無限二次元空間格子上のハミルトニアン・ヤン・ミルズ理論に関する初期の結果を提示するものである。

要約

あなたは、非常に巨大で、信じられないほど複雑なパズルを解こうとしているところだと想像してください。このパズルは、宇宙を繋ぎ止めている基本的な力（具体的には、陽子や中性子の中のクォークを束ねる強い力）を表しています。このパズルはあまりに巨大で、ピースの数は無限にあり、コンピュータを使って一つひとつのピースを解こうとするのは、スプーンで海を飲もうとするようなものです。

何十年もの間、物理学者たちはこれを解くために「モンテカルロ・シミュレーション」と呼ばれる手法を用いてきました。これは、目隠しをしたハイカーが山の中をあてもなく歩き回り、いつか最も低い谷（理論の基底状態）にたどり着くことを期待してランダムに足を進めるようなものです。この方法は機能しますが、非常に時間がかかり、遠くから山を見渡そうとすると（複雑さが無限大になる「大きなN」極限において）、非常に厄介なものになります。

新しいアプローチ：「コヒーレント状態」による地図

ローレンス・G・ヤフェによるこの論文は、このパズルを解くための異なる方法を提案しています。著者は、ランダムにさまようのではなく、「コヒーレント状態」という数学的概念に基づいた「地図」を使うことを提案しています。

このパズルを、混沌とした混乱状態としてではなく、滑らかな風景として捉えてください。「大きなN」極限（パズルの複雑さが無限大になる場合）では、量子的な奇妙さが消え去り、その風景は「古典的」になります。これは、霧に包まれた混沌とした夜（量子）と、晴れ渡った明るい昼（古典）の違いのようなものです。

著者の手法は、この滑らかな風景における絶対的な最低点（最小値）を見つけることです。一度谷の底を見つければ、その周囲にある丘の形を簡単に把握することができます。これにより、物理学者は粒子（グルーボール）の質量や、それらがどのように衝突するかといった、従来の「さまよい歩き」の手法では非常に困難だった事柄を計算できるようになります。

道具：「ゴルディオン（Gordion）」

これを行うために、著者は「ゴルディオン」と名付けられた新しいコンピュータプログラムを構築しました。この名前は、誰も解くことができなかった絡まった結び目（ゴルディオンの結び目）に直面したアレクサンダー大王の伝説にちなんだ巧みなものです。アレクサンダーは、糸を一本ずつ解こうとする代わりに、剣でその結び目を一気に切り裂きました。

同様に、「ゴルディオン」プログラムは、無限のパズルの糸を一本ずつ解こうとはしません。代わりに、「ループ・リスト」戦略を用います。それは最も重要なループ（粒子が辿る経路）に焦点を当て、それ以外を無視することで、実質的に複雑さを「切り抜ける」のです。

何を発見したのか？

著者は、この新手法をいくつかのシナリオでテストしました。

1. 単純なテストケース： 彼らは、非常に小さく単純なパズル（一つの「プラケット」または正方形のループ）から始めました。プログラムは完璧に動作し、既知の正確な回答と一致しました。これは、「剣」が鋭利であり、「地図」が正確であることを証明しました。
2. 2次元格子（平坦な世界）： 彼らは、二次元の格子にこれを適用しました。数学的な簡略化をあまり行わなくても、プログラムは正しい答えに非常に近い値を出しました。これは、通常は非常に難しいとされる領域（弱結合）においても同様です。
3. 3次元格子（現実世界のシミュレーション）： 彼らは、2+1次元の格子（二つの空間次元と一つの時間次元）を試しました。これはより困難な課題です。プログラムは強い相互作用に対してはうまく機能しましたが、相互作用が弱くなるにつれて苦戦し始めました。

限界：「切り捨て（Truncation）」の問題

主な課題は、通常のデスクトップコンピュータで実行するために、プログラムがパズルのピースの一部を無視しなければならないことです。これは「切り捨て」と呼ばれます。

• 比喩： 複雑な絵画を、大きな筆致の色だけをリストアップして説明しようとしていると考えてみてください。最初はこれで上手くいきます。しかし、ズームインしていく（あるいは、物理現象がより微細になっていく）と、細部を見落としてしまいます。
• 結果： プログラムは、「塗料」が厚く大胆なとき（強結合）には見事に機能します。しかし、塗料が薄く詳細になっていく（弱結合）につれて、近似がずれ始めます。プログラムは時として、物理的に不可能な結果（確率が100%を超えるなど）を生成することがあり、これは、利用できるピースが足りなくなったことを示しています。

「因子分解」の試み

著者は、欠けているピースを補うための巧妙なトリックを試みました。もし大きなループが二つの小さなループから成っているならば、その大きなループの値は、二つの小さなものの積になるはずだ、という仮定です。彼らはこれを「因子分解」と呼びました。

しかし、その結果は期待外れでした。この仮定が役に立つこともありましたが、多くの場合、状況を悪化させるか、あるいは何も変えないものでした。これは、複雑なスープの味を、単に二つの材料の味を掛け合わせることで推測しようとするようなもので、必ずしも全体の味わいを捉えることはできません。

結論

この論文は、「コヒーレント状態」のアプローチが、これらの無限のパズルを見るための強力な新しい方法であることを結論づけています。これにより、物理学者はランダムなシミュレーションによる統計的なノイズを回避し、理論の「無限」バージョンに直接取り組むことができます。

現在のバージョン（標準的なデスクトップコンピュータで動作するもの）は、まだ最も困難な3次元のパズルを解き明かしてはいませんが、コンセプトが機能していることは証明されました。著者は、より優れたコンピュータ（スーパーコンピュータ）と、欠けているピースを扱うより賢明な方法があれば、この手法が、現在の方法よりもはるかに直接的に、粒子の散乱や崩壊を正確に計算するという、現在では不可能な問題を解決できる可能性があると示唆しています。

要するに、著者は新しい剣（ゴルディオン）を研ぎ澄ませ、それが単純な結び目を完璧に切り裂くことを示しました。それはより大きな結び目をも切り裂き始めていますが、仕事をやり遂げるためには、より大きな手（計算能力）と、より鋭い刃（より優れた近似）を必要としています。

技術概要

技術要約：大規模Nゲージ理論におけるコヒーレント状態変分アルゴリズムの復活

問題提起
本論文は、大規模 $N$ 極限における非可換ゲージ理論（具体的には $SU(N)$および$U(N)$）を数値的に解く際の課題に取り組んでいる。大規模$N$極限では、量子場理論のダイナミクスはコヒーレント状態の位相空間上での古典的な最小化問題へと簡約されるが、非自明な格子ゲージ理論に対してこの問題を解くことは依然として計算量的に困難である。モンテカルロ法などの既存の手法は、実時間特性（崩壊幅や散乱振幅など）の抽出に制限があり、確率的サンプリングを必要とする。また、標準的なハミルトニアン・截断法は、2+1次元および3+1次元におけるヒルベルト空間の指数関数的な増大に苦慮する。著者は、有限体積効果や統計的サンプリングの分散なしに、基底状態の性質、励起スペクトル、および散乱振幅を計算することを目的として、$N=\infty$ 極限に直接アクセスするためのコヒーレント状態変分法の実現可能性を再評価している。

手法
本論文は、1980年代に定式化されたコヒーレント状態変分アルゴリズムの新しい実装として、"Gordion" と名付けられた統一されたC++プログラムについて記述している。この手法は、以下のコアコンポーネントに基づいている。

1. 状態表現： アルゴリズムは「ループ・リスト」截断スキームを使用する。状態は、ウィルソンループ（およびフェルミオン双線形が存在する場合はそれら）の期待値の有限集合によって表現される。これらの観測量は、「強結合次数」に基づいて選択される。これは、ある観測量を省略した場合に最初に誤差が生じる強結合展開の次数によって定義される。
2. 変分パラメータ： ウィルソンループの期待値を直接変分パラメータとして使用する代わりに（これはユニタリ制約を破り、複雑な最小化境界を招く）、リーマン正規座標を用いる。これらは、無限次元のコヒーレンス群 $G$（ループと電場の反エルミートな組み合わせの指数関数）の生成子の係数である。状態の変形は、これらの群要素を作用させることで生成され、これにより状態が物理的な領域内に留まることが保証される。
3. アルゴリズムのステップ：
   • 記号的生成： プログラムは、保持された観測量とコヒーレンス群の生成子の間の交換子を記号的に評価し、大規模 $N$ 位相空間上の測地線方程式を定義する。
   • エネルギー／自由エネルギー最小化： リーマン正規座標に関する大規模 $N$ ハミルトニアン（またはユークリッド自由エネルギー）の勾配と曲率（ヘッセ行列）を計算する。
   • ニュートン反復： 最小値はニュートン反復によって求められる。これは、最小値の位置を予測するために線形方程式を解き、期待値を更新するために測地線方程式を積分する。
   • スペクトル抽出： ハミルトニアン理論の場合、曲率行列とラグランジュ・ブラケット（逆ポアソン括弧）を含む一般化固有値問題を解くことにより、微小振動の周波数が決定され、グルーボールまたはメソンの質量スペクトルが得られる。
4. 観測量の近似： 截断誤差を減らすための初期の試みとして、非保持の観測量を因子化スキーム（自己交差ループに対して $W_{\Gamma} \approx W_{\Gamma_1} W_{\Gamma_2}$）を用いて近似する手法が用いられたが、現在の実装では限定的な有用性に留まるとされている。

主な貢献

• "Gordion" の実装： 著者は、立方格子上の $U(N)$ ゲージ理論（基本フェルミオンの有無を問わない）を、ハミルトニアン形式およびユークリッド形式の両方で扱うことができる、現代的で効率的なC++実装を提示している。
• ベンチマーク： コードは、ユークリッドおよびハミルトニアンの両方の定式化において、厳密に解ける一プラケットモデルを用いて検証されており、変分パラメータ（生成子）の数が増えるにつれて厳密解に収束することが示されている。
• 2次元ユークリッド・ヤン・ミルズ理論： 非局所的な独立プラケット変数への簡約を用いずに、無限格子上の2次元ユークリッド・ヤン・ミルズ理論の結果が提示されている。本手法は、$\lambda \approx 0.7$ の弱結合領域まで、自由エネルギーおよび様々なループ期待値を厳密な結果と一致させることに成功している。
• 2+1次元ハミルトニアン・ヤン・ミルズ理論： 本手法を、無限格子上の2+1次元ハミルトニアン・ヤン・ミルズ理論に初めて適用した結果を提示している。基底状態エネルギー、および様々な対称性チャネル（$A_1, A_2, B_1, B_2, E$）における最低次グルーボール質量の計算を行っている。

結果

• 収束性： 一プラケットモデルにおいて、変分結果は厳密解に急速に収束する。3次から5つの生成子があれば、第三次相転移付近であっても、厳密な自由エネルギーとループ期待値を視覚的に区別できないレベルで再現できる。
• 2次元ユークリッド理論： 最大3つのプラケット演算子を含む次数6までの生成子を使用し、次数28までの観測量で截断した場合、$\lambda > 0.7$ において厳密な自由エネルギーおよびループ期待値を高い精度で再現した。$\lambda$ がさらに減少すると、截断誤差が顕著になる。
• 2+1次元ハミルトニアン理論：
  • 本手法は、基底状態エネルギーおよび小さなループの期待値を正常に計算している。
  • グルーボール質量は複数の対称性チャネルで抽出されている。結果は、期待される強結合挙動（質量が $\lambda$ に比例してスケールする挙動）を示している。
  • しかし、弱結合領域における収束は、2次元ユークリッドの場合よりも遅い。次数6の生成子による結果は、連続極限への外挿に必要な、原点への線形的な接近をまだ明確に示していない。
  • 限界： 計算は、観測量の数の急速な増大と測地線方程式のサイズによって制限されている。2+1次元理論において、次数6の生成子を用いた計算は、非物理的なループ期待値（$|W_\Gamma| \le 1$ を破るもの）の出現や、曲率固有値の複素化により、$\lambda \approx 1.4$ 以下では収束に失敗した。
• 観測量の近似： 因子化に基づく近似のテストが行われたが、結論は出ていない。これは計算上のストレージ要件を大幅に増加させたが、精度を一貫して向上させたり、高次の截断結果を模倣したりすることには至らなかった。

意義と主張
本論文は、コヒーレント状態変分法が、大規模 $N$ ゲージ理論（$N=\infty$）を直接研究するための実行可能な決定論的手法であることを主張している。その主な意義は、以下の点にある。

1. 標準的な格子シミュレーションに特有の有限サイズ効果を回避し、無限体積で直接扱えること。
2. 励起スペクトル（グルーボール/メソンの質量）に直接アクセスできること。また、原理的には、ユークリッド時間からの解析接続を介さずに、高次の微分を通じて崩壊幅や散測振幅を得られること。
3. ディラック行列式の複雑さを伴わずに、動力学的なフェルミオンを扱うこと。

著者は、2+1次元ハミルトニアン理論に対する現在の結果は「有望な初期の試み」であるが、截断誤差と計算リソース（デスクトップコンピュータで実施された）によって制限されていると控えめに結論付けている。本論文は、3+1次元QCDを解決した、あるいは2+1次元ヤン・ミルズの連続極限に到達したと主張するものではない。代わりに、本手法の実現可能性を確立し、物理的に重要な弱結合領域へと拡張するために対処すべき具体的な技術的課題――例えば、より優れた観測量選択基準、截断された観測量の改良された近似スキーム、および有限次元マスターフィールド表現の潜在的な有用性――を特定している。著者は、この古典的な最小化問題は困難ではあるものの、有限 $N$ に対する直接的なハミルトニアン截断や、量子コンピュータ上での実時間発展よりも、おそらく扱いやすいものであることを強調している。