A discrete formulation for three-dimensional winding number

この論文は、3 次元多様体上の滑らかな写像に対する巻き数（winding number）を、偶然の縮退や対称性によって生じる縮退を含む系にも直接適用可能なロバストな離散化手法（$\theta$-ギャップに基づく）として定式化し、実用的な単純なフラックスと厳密な整数量子化を保証する修正フラックスの 2 種類を提案するものである。

要約

この論文は、物理学の難しい概念である「3 次元の巻き数（ウィンドイング・ナンバー）」を、コンピュータで計算しやすいように工夫した新しい方法を提案するものです。

専門用語を抜きにして、**「宇宙の地図を描く」**という物語に例えて説明しましょう。

1. 何の話？「ねじれた宇宙の地図」

まず、この論文が扱っているのは、**「3 次元の空間（X）」という巨大なキャンバスに描かれた、「ねじれた模様（g）」**です。
この模様は、ある場所から別の場所へ移動するたびに、少しずつ回転したりねじれたりしています。

• 巻き数（Winding Number）： この模様全体が、空間全体で「何回、完全にねじれたか」を表す整数です。
  • 例：1 回ねじれたら「1」、2 回ねじれたら「2」。ねじれていなければ「0」。
  • この数は、超伝導体や素粒子の理論など、物理の重要な分野で「この物質はどんな性質を持っているか」を判別する**「指紋」**のような役割を果たします。

2. 従来の方法の「落とし穴」

この「ねじれ具合」を計算するには、通常、空間を小さな点（格子点）の集まりとして近似し、それぞれの点で模様の状態を調べる必要があります。

• 昔の方法（バンド追跡）：
  以前の計算方法は、**「それぞれのねじれ（バンド）を、隣の点まで手をつなぎながら追いかける」**というものでした。
  • 問題点： もし、2 つ以上のねじれが**「偶然重なる（縮退）」場所があったり、「対称性によって必ず重なる」**場所があったりすると、誰が誰の手を握っているか分からなくなってしまいます。まるで、混雑した駅で友達と待ち合わせをするとき、誰が誰だか分からなくなってしまうようなものです。この場合、計算が複雑になりすぎて、正しく答えが出せなくなることがありました。

3. 新しい方法の「魔法の穴（θ-ギャップ）」

この論文の著者（塩崎健氏）は、**「個別のねじれを追いかける必要はない！」**という発想で、全く新しい方法を考え出しました。

• アイデア：「ねじれのない安全地帯」を作る
  著者は、**「θ-ギャップ（θ-gap）」という概念を使います。
  これは、「ねじれの模様の中に、特定の角度（θ）の場所が『空いている（存在しない）』こと」**を意味します。

  • アナロジー：
    Imagine 円盤（ねじれ模様）を回しているところを想像してください。
    従来の方法は、円盤上の「赤い点」や「青い点」を一つずつ追いかけていました。
    しかし、新しい方法は、**「円盤のどこかに『赤い点』が絶対に存在しない『安全な穴』を見つける」**ことに集中します。

    空間の小さな区画（立方体）ごとに、「この区画では、この角度のねじれは存在しない（穴が開いている）」と確認できれば、その区画全体を安全に扱えるようになります。

• なぜこれがすごいのか？
  ねじれが重なり合ったり（縮退）、複雑に絡み合ったりしても、「特定の角度に穴が開いている」という事実さえあれば、個別のねじれを追いかける必要がなくなります。
  「混雑した駅で、誰が誰だか分からない」状態でも、「改札口（穴）が空いている場所」さえ分かれば、人々はスムーズに通り抜けることができます。 これにより、計算が非常にシンプルで頑丈（ロバスト）になりました。

4. 2 つの計算ツール：「簡易版」と「厳密版」

この新しい方法では、2 つの計算ツール（フラックス）を提案しています。

1. 簡易版（Φp）：

   • 特徴： 非常にシンプルで、その場所の「穴」の情報だけで計算できます。
   • 性能： 格子（点の網）を細かくすればするほど、ほぼ確実に正しい整数の答えが出ます。
   • 使い道： 実用的な計算には、この「簡易版」で十分で、とても使いやすいです。

2. 厳密版（˜Φp）：

   • 特徴： 簡易版の計算結果に、少しだけ「補正」を加えたものです。
   • 性能： 数学的に**「絶対に整数になる」**ことを保証します。
   • 使い道： 理論的に完璧な証明が必要な場合や、簡易版で少しズレが生じる可能性がある場合に使います。

5. 結論：なぜこれが重要なのか？

この論文は、**「ねじれた物理現象を、コンピュータで簡単に、かつ間違いなく計算できる新しいルール」**を提案しました。

• 従来の課題： ねじれが重なり合うと計算が破綻していた。
• 今回の解決： 「穴（ギャップ）」を見つけるという視点を変えれば、重なり合っても計算できる。
• 実用性： 複雑な物質（超伝導体など）の性質を、コンピュータシミュレーションで正確に予測できるようになります。

まとめると：
「ねじれた世界を数える」ために、以前は「一人一人のねじれを名前を呼んで追いかける」大変な作業が必要でしたが、この論文では**「ねじれが通れない『穴』を見つけ、その穴を基準に全体を数える」**という、賢くて簡単な方法を見つけたのです。これにより、どんなに複雑な物質でも、その「指紋（トポロジカルな性質）」を正しく読み取れるようになりました。

技術概要

以下は、Ken Shiozaki 氏による論文「A discrete formulation for three-dimensional winding number (3 次元巻き数の離散定式化)」の技術的な要約です。

1. 問題の背景と目的

3 次元の向き付けられた閉多様体 $X$ からユニタリ群 $U(N)$ への滑らかな写像 $g: X \to U(N)$ に対して定義される3 次元巻き数（winding number） $W_3$ は、時間反転対称性を持つ 3 次元超伝導体のトポロジカル不変量や、非可換格子ゲージ理論におけるインスタントン数の計算など、物理学の重要な分野で中心的な役割を果たしています。

従来の離散化手法（格子近似）では、数値計算のために写像 $g$ が格子点上で定義されることが一般的です。しかし、既存の手法（Höckendorf らによるもの）は、隣接する格子点間で固有状態（固有ベクトルと固有値）を追跡・一致させる（matching）プロセスに依存しています。このアプローチは、偶然の縮退や対称性によって強制された縮退（degeneracies）が存在する系において、バンドの追跡が複雑化し、ロバスト性に欠けるという課題がありました。

本研究の目的は、縮退が存在する場合でも直接的に適用可能で、ロバストな 3 次元巻き数の離散計算手法を開発することです。

2. 手法と定式化

本研究は、**$\theta$-ギャップ（$\theta$-gap）**の概念に基づいた新しい離散定式化を提案しています。

• $\theta$-ギャップの定義:
  行列 $g \in U(N)$ が、ある実数 $\theta \in [0, 2\pi)$ に対して固有値 $e^{i\theta}$ を持たないとき、その行列は $\theta$-ギャップを持つと定義されます。
• 局所的な 2 形式 $B_\theta$ の構成:
  3 形式 $H = \frac{1}{12\pi} \text{Tr}[(g^{-1}dg)^3]$ は局所的に完全微分 $H=dB_\theta$ と書けます。ここで $B_\theta$ は対角化行列 $\gamma$ と対角行列 $\Lambda$ を用いて構成され、$\theta$ の選択に依存する項 $R_\theta$ と独立な項 $Q$ からなります。
• 格子への適用:
  多様体 $X$ を立方格子 $L$ で近似します。各セル（キューブ）$c$ に対して、その内部のすべての点で $\theta_c$-ギャップが存在するように $\theta_c$ を選択します（固有値の分布に基づき、ギャップの最も広い部分の中央などを $\theta_c$ とします）。
• フラックスの定義:
  巻き数は、各セルの境界（面、plaquette）上の積分の和として書き換えられます。
  $$W_3[g] = \frac{1}{2\pi} \sum_p \int_p (B_{\theta^-_p} - B_{\theta^+_p})$$
  ここで、$\theta^\pm_p$ は面 $p$ に隣接する 2 つのキューブのギャップパラメータです。この積分は、面の周回（ループ）上の線積分に変換され、さらに離散的なベリー位相として近似されます。

3. 2 つの離散フラックスの提案

著者は、実用的な「単純なフラックス」と、厳密な整数量子化を保証する「修正されたフラックス」の 2 つを提案しています。

1. 単純なフラックス ($\Phi_p$):

   • 各面の 4 つの頂点における対角化行列 $\gamma$ のみを用いて計算されます。
   • 式 (19) に示されるように、隣接する頂点間の固有ベクトルの内積（行列式）の積の偏角（Arg）として定義されます。
   • 特徴: 実用的で計算が容易ですが、厳密には整数値にならない可能性があります（$2\pi$ の曖昧さや、離散化誤差によるもの）。しかし、格子が十分に細かければ、ほぼ常に整数に量子化されることが示されています。

2. 修正されたフラックス ($\tilde{\Phi}_p$):

   • 厳密な整数量子化を保証するために導入されました。
   • 面 $p$ に隣接する 4 つのキューブの $\theta$-ギャップ情報を考慮し、エッジ（辺）ごとの寄与を再定義します。
   • 特定の配置（例えば、異なる $\theta$ 範囲の固有状態が混在する場合）において、非対角成分による位相の曖昧さを排除するため、$\theta$ 値をソートし直したブロック対角構造を持つユニタリ接続 $\tilde{A}$ を定義します。
   • 特徴: 構成上、$e^{2\pi i \tilde{W}^{\text{dis}}_3} = 1$ となり、計算結果 $\tilde{W}^{\text{dis}}_3$ は厳密に整数となります。

4. 数値検証と結果

提案された手法は、以下のモデルで検証されました。

• モデル計算:
  • 3 次元トーラス上のモデル $g_0$ において、解析的な巻き数と一致することを確認しました（$m$ の値に応じて $W_3 = -2, 1, 0$ となる領域を正しく再現）。
  • 縮退を含むモデル（$g_1 = g_0 \oplus e^{i\pi/2}g_0 \oplus u + v$）において、単純なフラックス $\Phi_p$ の和は整数にならず、修正フラックス $\tilde{\Phi}_p$ のみが厳密な整数値を与えることを確認しました。
• ランダムモデル:
  • 有限範囲のランダムホッピングを持つモデル $g(k)$ に対して、対称性制約 $Jg(k)^*J = g(k)$ を課すことで、特異点（vortex loop defects）を回避し、逆行列が存在するユニタリ行列を構成しました。
  • 格子サイズ $L$ を変化させた際、修正フラックス $\tilde{W}^{\text{dis}}_3$ はすべての格子サイズで厳密に整数（真の値 $-1$）を返し、単純フラックス$W^{\text{dis}}_3$ は格子が粗い段階では非整数値を示すが、格子を細かくするにつれて整数に収束することが確認されました。
  • 最大面フラックス $\max_p \Phi_p$ が格子サイズとともに 0 に収束し、離散化近似の精度向上を示しました。

5. 意義と結論

• ロバスト性の向上: 既存のバンド追跡法とは異なり、固有状態の一致付けを必要としないため、縮退（偶然の対称性によるものも含む）が存在する系においても直接適用可能です。
• 実用性と厳密性の両立: 実用的な計算には、格子が十分細かければほぼ整数になる「単純なフラックス」で十分ですが、厳密なトポロジカル不変量が必要な場合や粗い格子での計算には、「修正されたフラックス」が有効です。
• 将来の展望: この定式化は、インスタントン数（第 2 チェルン数）や、より一般的な対称空間への写像の次数など、他の高次元トポロジカル不変量の離散計算への拡張が期待されます。

本研究は、トポロジカル物質の数値解析や格子ゲージ理論において、3 次元巻き数を効率的かつ正確に評価するための強力なツールを提供するものです。