Efficient local classification of parity-based material topology

技術要約：パリティに基づく物質トポロジーの効率的な局所分類

問題提起
結晶性物質におけるトポロジカル相の分類は、 Chern数や $Z_2$ 指数といった、並進対称性に依存する運動量空間の手法を通じて確立されている。しかし、これらの従来の手法は、並進対称性が存在しない非周期系（例：準結晶、アモルファス材料、無秩序系）には適用できない。$Z_2$ トポロジカル絶縁体（クラス AII）に対する既存の実空間代替手法（例：Bott指数、スピンBott指数などの投影位置演算子）は、大きな制限に直面している。具体的には、これらはハミルトニアンにスペクトルギャップが存在することを前提としており、場合によってはスピン投影にもギャップを必要とするため、現実的な無秩序モデルで一般的なギャップレス系には適用できない。また、射影演算子ベースの手法は密行列（dense matrices）の使用を強いるため、低エネルギー近似に限定され、そのような近似を事前に特定する必要がある。したがって、スペクトルギャップ、スピン分解、または低エネルギー近似に依存せずに、パリティに基づくトポロジーを分類できる、数値的に効率的な実空間フレームワークが切実に求められている。

手法
著者らは、**スペクトル・ロカライザー（spectral localizer）**と、新規でスケーラブルな疎行列数値アルゴリズムを組み合わせたフレームワークを提示している。

1. スペクトル・ロカライザー・フレームワーク:
   この手法は、システムのハミルトニアン $H$ と位置演算子 $X, Y$ から構成されるスペクトル・ロカライザー演算子 $L_{(x,E)}$ を利用する。2次元のクラス AII 系（フェルミオンの時間反転対称性 $T^2 = -I$ を持つもの）において、このフレームワークは、特定の基底において $L_{(x,E)}$ を純虚数の歪対称行列へと変換する性質を利用する。局所的なトポロジカル不変量は、この歪対称行列のパフィアン（Pfaffian）の符号によって決定される：
   $$\xi^L_{(x,E)} = \text{sign}[\text{Pf}(Q^\dagger L_{(x,E)} Q)] \in Z_2$$
   極めて重要な点は、この手法はシステムがバルク・スペクトルギャップやギャップのあるスピンセクターを持つことを必要としないことである。これは、エネルギー分解された局所不変量とともに、トポロジカルな保護を定量化し、$\mu \to 0$ のときに局在境界状態の存在を示す局所ギャップ尺度 $\mu_{(x,E)}$ を提供する。

2. 疎行列パフィアン符号アルゴリズム:
   大規模な物理系に対してこのフレームワークを計算可能にするため、著者らは、完全な行列式やパフィアンを計算することなく、大きな疎な歪対称行列のパフィアンの符号を計算するスケーラブルなアルゴリズムを開発した。

   • 因子分解: この手法は、疎な歪対称行列 $S$ を、単位対角成分を持つ下三角行列 $L$ と、歪対称なブロック対角行列 $T$ を用いた分解 $S = L T L^\top$ に依拠している。
   • 符号計算: パフィアンの符号は、$L$ の対角成分の符号と $T$ の非対角ブロックの符号から導出される：$\text{sign}[\text{Pf}(S)] = \text{sign}[\det(L)] \cdot \text{sign}[\text{Pf}(T)]$。
   • 実装: アルゴリズムは、Trilinos フレームワークの一部である Tacho ソフトウェアパッケージを利用しており、2段階の置換戦略を採用している。第一に、最大カーディナリティまたは重み付きマッチングアルゴリズム（$P_s$）を用いて、安定な $2\times2$ ピボットを確保するために大きな数値的成分を副対角線上に移動させる。第二に、フィルイン（fill-in）を最小限に抑えるために、入れ子状の解体アルゴリズム（METIS による $P_n$）を用いて行列を再順序付けする。このアプローチにより、標準的なワークステーション上で数百万の自由度を持つ行列を扱うことが可能となり、実行時間のスケーリングにおいて従来の密行列アルゴリズムを大幅に凌駕する。

主な結果
本フレームワークは、2つの異なる非周期系に適用された：

1. 準結晶量子スピンホール（QSH）系:
   著者らは、反対のChern数を持つ2つの結合された準結晶層（Penroseタイリングに基づく）のヘテロ構造をモデル化した。

   • ギャップ系: 手法は、ギャップ領域における非自明な $Z_2$ トポロジーと、それに伴うエッジ状態を正確に特定した。
   • ギャップレス系: 枠組みは、一方がギャップを持ち、もう一方がギャップレスであるギャップレス・ヘテロ構造に適用された。グローバルなスペクトルギャップが存在しないにもかかわらず、局所指数 $\xi^L$ と局所ギャップ $\mu$ は、トポロジカル相を正しく診断し、局所ギャップが閉じた箇所での「リーキー（leaky）」な境界局在状態を予測した。これは、射影ベースの手法が失敗する場面でも本手法が機能することを示している。

2. 脆弱なトポロジーを持つフォトニック準結晶:
   手法は、$C_2T$ 対称性（ボゾンの時間反転対称性）を持つ2次元フォトニック準結晶（Ammann-Beenkerタイリング）に適用された。

   • システムは有限差分法を用いて離散化され、約 $3.2 \times 10^6$ 次元の疎行列となった。
   • アルゴリズムは、標準的なワークステーション上で1〜2時間で、広範な周波数範囲にわたって脆弱なトポロジカル不変量を計算することに成功した。
   • 特定されたトポロジカル相転移は、局所ギャップの消失と一致しており、疎行列因子分解戦略の正確性を裏付けている。

意義と主張
本論文は、スペクトル・ロカライザーのフレームワークと効率的な疎行列数値手法の組み合わせが、パリティに基づくトポロジカル相を診断するための統一的かつ堅牢なツールを提供すると主張している。

• 限界の克服: 従来のアプローチとは異なり、本手法は並進対称性、バルク・スペクトルギャップ、または補助的な演算子（スピン投影など）のギャップに依存しない。ハミルトニアン全体に直接作用するため、低エネルギー近似の必要性を排除している。
• スケーラビリティ: 疎行列パフィアン符号アルゴリズムの開発により、従来の密行列線形代数手法では到達不可能であった、数百万の自由度を持つ現実的な材料モデルの解析が可能となった。
• 物理的洞察: 本フレームワークは、単なる二値分類だけでなく、トポロジカルな保護の固有の尺度（局所ギャップ）と、ギャップが閉じる場所での状態の局在化を通じたバルク・境界対応への直接的なつながりを提供する。
• 幅広い適用性: 著者らは、これらの手法が、従来の運動量空間による診断が適用できない、ねじれたファンデルワールス系や大規模なメタマテリアルといった複雑なプラットフォームにおけるトポロジカル相の系統的な発見において価値を持つと期待している。

本研究は、実空間におけるトポロジカル分類と、大規模な疎演算子のホモトピー類との間の直接的な数学的接続を確立しており、周期性のないトポロジカル物質を解析するための実用的な経路を提示している。