Tensor-network formulation of QCD in the strong-coupling expansion

本論文は、任意の次元・色数・フレーバー数を持つ格子 QCD の強結合展開を、局所テンソル（数値部分とグラスマン部分から構成）のネットワークとして定式化し、$2\times2$ 格子における解析的計算結果と、より大規模な格子向けに開発された新しい手法「order-separated GHOTRG」の初期結果を提示するものである。

要約

この論文は、「素粒子の世界（クォークやグルーオン）がどう振る舞うか」を、新しい計算手法を使って解き明かそうとする研究です。

専門用語を避け、日常の比喩を使ってわかりやすく説明しましょう。

1. 問題の正体：「見えない迷路」と「複雑な計算」

まず、この研究が取り組んでいるのは**QCD（量子色力学）**という分野です。これは、原子核を構成する「クォーク」という小さな粒子たちが、どうやって集まって物質を作っているかを説明する理論です。

しかし、この理論をコンピュータでシミュレーションしようとすると、大きな壁にぶつかります。

• 壁： 「化学ポテンシャル（物質の密度）」が高い状態を計算しようとすると、計算式の中に**「マイナスとプラスが激しく入り乱れて、結果がゼロになってしまう」**という現象（符号問題）が起きます。
• 例え： 暗闇の迷路で、足元の石が「プラス（登れる）」と「マイナス（落ちる）」を交互に示していますが、その数が多すぎて、どちらが正しい道かわからなくなってしまうような状態です。これまでの計算方法（モンテカルロ法）では、この迷路の先（高密度の状態）に進むことができませんでした。

2. 新しい武器：「テンソル・ネットワーク」という巨大なパズル

そこで、この論文の著者たちは**「テンソル・ネットワーク」**という新しいアプローチを使いました。

• 比喩： 巨大なパズルを想像してください。
  • 従来の方法は、パズルのピースを一つずつランダムに置いて、全体像がどうなるか推測する「試行錯誤」でした。
  • 新しい方法は、パズルのピースを**「論理的なルールに従って、すべてつなぎ合わせて、完成図を直接描く」**方法です。

彼らは、QCD の複雑な計算式を、「小さなブロック（テンソル）」の集まりとして書き換えました。

• このブロックには、**「数字」と「グラスマン数（フェルミオンの性質を表す特殊な数）」**という 2 つの部品が入っています。
• これらをすべてつなぎ合わせると、最終的に「パズルの完成図（分配関数）」が得られ、そこから物理的な性質（質量やエネルギーなど）が読み取れるようになります。

3. 強結合展開：「弱い力」から「強い力」へ

これまでの研究では、粒子同士の相互作用が「無限に強い」極限の場合しか扱えていませんでした。しかし、現実の世界はもっと複雑です。

• 今回の突破： 彼らは、相互作用の強さを少しずつ変えていく**「強結合展開」**という手法を取り入れました。
• 例え： 音楽の楽譜を、最初は「ド」だけ（無限結合）で考えていたのが、次に「ド・レ・ミ・ファ…」と音階を少しずつ足していって、よりリアルなメロディ（現実の物理現象）を再現しようとする試みです。
• 彼らは、この新しい「パズル（テンソルネットワーク）」を使って、2×2 という小さな盤面（格子）で、このメロディを**4 番目の音まで（βの 4 乗まで）**正確に計算することに成功しました。

4. 重要な発見：「計算の仕方」で結果が変わる

面白い発見がありました。同じパズルから答えを出すとき、「どう計算するか」によって、答えの精度が全く違うということです。

• 方法 A（直感的な計算）： 「パズルの完成図そのもの」から答えを出す方法。
  • 結果：小さな盤面では合いますが、盤面が大きくなると、遠くまで見通せなくなります。
• 方法 B（対数をとった計算）： 「パズルの完成図の『対数』（複雑さを整理した形）」から答えを出す方法。
  • 結果：盤面が大きくなっても、モンテカルロ法（従来の試行錯誤法）のデータと非常に良く一致します。

結論： 大きな世界を正しく見るには、単にパズルを解くだけでなく、**「その構造を整理して（対数化して）から解く」**のが正解でした。

5. 今後の展望：「OS-GHOTRG」という新技術

小さな盤面（2×2）では手計算でもできましたが、現実の大きな盤面（8×8 やそれ以上）になると、パズルのピースが多すぎて、従来の方法では「どの音階（展開係数）がどこに含まれているか」を区別できなくなってしまいます。

そこで、彼らは**「OS-GHOTRG（順序分離型 GHOTRG）」**という新しい技術を導入しました。

• 比喩： 大勢の合唱団の中で、特定の「ソプラノの声」だけを聞き分けるための、超高性能なマイクのようなものです。
• これを使えば、大きな盤面でも、必要な「音階（展開係数）」だけを正確に抽出して計算できるようになります。

まとめ

この論文は、**「素粒子の高密度状態という、これまで解けなかった難問を、新しいパズル解法（テンソルネットワーク）で解き明かすための道筋を示した」**という画期的な成果です。

• 何をした？ 複雑な計算式を「ブロックパズル」に変えて、新しいルールで解けるようにした。
• 何がわかった？ 大きな世界を正しく計算するには、「整理してから解く（対数展開）」のが重要だとわかった。
• これから？ 「順序分離型」という新しいツールを使って、より大きな世界（現実的な格子サイズ）の計算が可能になる。

これは、宇宙の始まりや中性子星の内部など、極限状態の物質を理解するための、非常に重要な第一歩です。

技術概要

この論文は、非ゼロの化学ポテンシャルを持つ格子 QCD（量子色力学）の強結合展開を、テンソルネットワーク法を用いて定式化したものである。著者らは、ゲージ自由度とクォーク自由度を積分消去し、分配関数を局所的なテンソル（数値部分とグラスマン部分からなる）の完全な縮約として書き換えることに成功している。

以下に、論文の技術的要点を問題設定、手法、主要な貢献、結果、意義の観点から詳細に要約する。

1. 問題設定と背景

• 符号問題（Sign Problem）: 格子 QCD におけるモンテカルロシミュレーションは、化学ポテンシャル $\mu$ が存在する場合、ディラック演算子の行列式が複素数になるため「符号問題」に直面し、計算が困難になる。
• 既存手法の限界: 再重み付け、$\mu$ に関するテイラー展開、虚数化学ポテンシャルからの解析接続、複素ランジュバン法などの既存手法は、$\mu/T > 1$ の領域（高密度領域）には適用できない。
• 強結合展開と双対変数: 無限結合極限（$g \to \infty$）における双対変数（占有数）の定式化は有望視されているが、有限結合（強結合展開）への拡張は未解決であった。
• 目的: 任意の次元、カラー数、フレーバー数、および非ゼロ化学ポテンシャルに対して、強結合展開（$\beta \sim 1/g^2$ の逆数展開）をテンソルネットワーク形式で定式化し、大規模格子での計算を可能にすること。

2. 手法と定式化

論文では、以下のステップを経て分配関数をテンソルネットワークに変換している。

A. テイラー展開と双対変数の導入

• フェルミオン作用とゲージ作用（ウィルソン・プランケット作用）を含むボルツマン因子を、リンク上のホッピング項とプランケット項に対して個別にテイラー展開する。
• これにより、リンク上のフェルミオンホッピング数（$k, \bar{k}$）とプランケット上の占有数（$n$）という「双対変数（占有数）」が導入される。

B. ゲージリンクの積分とカラー自由度の縮約

• 各リンク上の SU($N_c$) ゲージ場 $U$ に関する積分を実行する。
• 従来の双対変数定式化では現れる非局所的な符号因子を回避するため、著者らは Gagliardi と Unger、Borisenko らが提案した効率的な積分公式（Weingarten 関数を用いたもの）を採用する。
• 重要な簡略化: 両端にグラスマン変数が存在することを利用し、カラー添字の和（置換群 $S_p$ 上の和）を大幅に削減する。これにより、カラー自由度を完全に消去し、局所的なテンソル構造を得る。

C. グラスマン変数の積分と補助変数の導入

• 元のグラスマン変数（$\psi, \bar{\psi}$）を積分消去するために、各リンクに新しい色・フレーバーを持たない補助グラスマン変数（$\chi$）を導入する。
• これにより、元の非局所的なホッピング項が、局所的なグラスマン変数と数値テンソルの積に変換される。
• 質量項もテイラー展開され、サイトごとの質量占有数 $w$ として扱われる。

D. テンソルネットワークの構築

• 最終的に、分配関数は局所的な数値テンソルとグラスマンテンソルの完全な縮約（Trace）として表現される。
• プランケット占有数 $n$ は、4 つのテンソルに共有されるため、直接のテンソル縮約には適さない。これを解決するため、プランケットを 4 つの「エッジ変数」に分解し、それらが等しくなるように制約（Kronecker デルタ）を課すことで、標準的なテンソルネットワーク構造（リンク変数のみで構成される）へ変換する。
• 初期テンソルのサイズを有限にするため、$\beta$ の展開次数を $n_{\text{max}}$ で切断する（切断近似）。

3. 主要な貢献と新規性

1. 強結合展開の一般化: 無限結合極限を超え、有限の $\beta$ に対する QCD の強結合展開を、任意の次元・カラー・フレーバーに対してテンソルネットワーク形式で定式化した初の試み。
2. OS-GHOTRG（Order-Separated GHOTRG）の提案:
   • 標準的な GHOTRG（Grassmann Higher-Order Tensor Renormalization Group）法では、数値的に分配関数 $Z$ を得るが、$\beta$ 展開係数 $Z_n$ を直接得ることは困難である。
   • 物理量（カイラル凝縮など）を高精度に計算するには、$Z$ ではなく $\ln Z$ または物理量自体を $\beta$ 展開して係数ごとに計算する必要がある（後述の結果参照）。
   • このため、展開係数 $Z_n$ を明示的に計算するための新しいアルゴリズム「順序分離型 GHOTRG（OS-GHOTRG）」を開発した（詳細は次号の論文に譲るが、本論文でその有効性を示した）。
3. 展開手法の比較と理論的解析:
   • 物理量の計算において、「$Z$ を展開してから対数をとる（展開 A）」と「$\ln Z$ を直接展開する（展開 B）」の 2 方法を比較。
   • 理論的解析と数値実験により、体積 $V$ が大きい場合や $\beta$ が大きい場合、展開 B（$\ln Z$ の直接展開）の方がモンテカルロ結果とよく一致し、安定であることを示した。

4. 結果

• 2x2 格子での解析的計算:
  • $N_c=3, N_f=1$ の 2 次元 QCD について、$n_{\text{max}}=4$ まで分配関数の展開係数 $Z_n$ を解析的に計算した。
  • カイラル凝縮 $\Sigma$ について、モンテカルロ（MC）データと比較。展開 B（$\ln Z$ の展開）が広い $\beta$ 範囲で MC データと一致し、展開 A は $\beta$ が大きくなると大きく乖離することが確認された。
• 8x8 格子での数値計算:
  • 開発した OS-GHOTRG 法を用いて 8x8 格子での計算を行った。
  • 標準的な GHOTRG や展開 A に比べ、展開 B を用いた OS-GHOTRG の結果がモンテカルロデータと非常に良く一致することを示した。
  • 標準的な GHOTRG 法では、不完全な高次項の影響により、大きな $\beta V$ 領域で物理量が一定値（プラトー）に収束する誤った挙動を示すことが確認された。

5. 意義と展望

• 符号問題の回避: 化学ポテンシャルがある場合でも、テンソルネットワーク法を用いることで数値的な不安定性（符号問題）を回避でき、強結合領域での QCD 相図の探索が可能になる。
• 高精度計算への道筋: 物理量を直接展開する手法（展開 B）と、展開係数を直接計算する OS-GHOTRG 法の組み合わせは、大規模格子での高精度な熱力学的観測量の計算を可能にする。
• 将来の展望: 本論文はシリーズの第 1 報であり、OS-GHOTRG 法の詳細と、より大規模な格子（4 次元など）での QCD 相図の計算が次の論文で行われる予定である。

結論:
この研究は、強結合 QCD をテンソルネットワークで扱うための堅固な定式化を提供し、特に物理量の計算における展開手法の重要性と、展開係数を直接扱う新しいアルゴリズム（OS-GHOTRG）の必要性を明らかにした。これは、化学ポテンシャル領域における QCD の非摂動的な研究に対する重要な進展である。