Higher Berry curvature, second Chern numbers and magnetoelectric coupling in crystalline insulators

技術要約：結晶性絶縁体における高次ベリー曲率、第2チェルン数、および磁気電気結合

問題提起
トポロジカル不変量は、量子多体系の基底状態を分類するために不可欠である。ベリー曲率（2形式）およびそれに付随するチェルン数は、ゼロ次元の基底状態（例：ブリルアンゾーンにおけるブロッホ状態）の族を特徴付けることができるが、近年の理論的発展は、より空間的に広がった系への一般化を示唆している。この枠組みにおいて、$d$ 次元の空間における短距離量子もつれ状態の族は、より高次の微分形式（例：1次元系における3形式）によって特徴付けられる。具体的には、3次元空間でパラメータ化された1次元系の場合、関連する不変量は、高次ベリー曲率から導出されるディクミエ・ドゥーディ・カプスティン・スポディニコ（DDKS）数 $\nu$ である。

汎用的な相互作用系または自由フェルミオン系において、これらの高次不変量を明示的に量子化された形で計算するための手法には、依然として大きな空白が存在する。特定の自由フェルミオンモデル（例：4次元ディラックモデル）については解析的な表現が存在するものの、4次元における第2チェルン数 ($ch_2$) を計算するための標準的な数値的手法は、格子ゲージ理論や実空間の磁場応答に依存することが多い。これらの標準的な手法は、「明示的な量子化」を欠いていることが多く、すなわち、構成によって整数値を返すのではなく、連続極限に到達するか、あるいは広範なメッシュの精緻化を必要とする。さらに、3次元空間におけるDDKS数の物理的解釈、特に軌道磁気電気結合およびチェルン・サイモンズ・アキシオン角 ($\theta_{CS}$) との関係は、未解決の課題である。

手法
著者らは、4次元（4D）チェルン絶縁体の格子モデルを、3次元のブリルアンゾーン ($T^3$) によってパラメータ化された、並進対称性を持つ無限の1次元（1D）鎖の族として書き換えることで、これらの問題に対処している。

1. モデルの再定式化： 4Dディラック・ハミルトニアンから出発し、著者らは一つの運動量成分 ($k_4$) に関して部分フーリエ変換を行い、パラメータ $\lambda = (\lambda_1, \lambda_2, \lambda_3)$ に依存する1Dハミルトニアンの族 $\{H_\lambda\}$ を得た。これらのハミルトニアンは、単位セルあたり2つの軌道を持つ1D鎖を記述する。
2. 数値フレームワーク (iMPS)： 1D鎖の基底状態は、無限行列積状態（iMPS）を用いて近似される。著者らは、高次ベリー曲率を計算するために、Fukui-Hatsugai-Suzukiの手法を一般化した、文献[13]で提案されたアルゴリズムを利用している。
3. 離散高次ベリー曲率：
   • パラメータ空間 $T^3$ は、グリッドに分割され、四面体に三角分割される。
   • 各頂点において、iMPS基底状態がiDMRGにより得られる。
   • 隣接するiMPSテンソル間の遷移関数は、混合転送行列の主要な固有ベクトルを用いて計算される。
   • 三角形の面に対して、ボンド次元の変化を扱うためにシュミット・スペクトルで重み付けされた、ゲージ不変な位相 $\phi(\Delta_2)$ が定義される。
   • 四面体に対する高次ベリー曲率 $F(\Delta_3)$ は、その境界上のこれらの位相の向き付き和として定義される。
   • DDKS数 $\nu$ は、$T^3$ 全体の三角分割にわたる $F(\Delta_3)$ の総和として得られ、明示的に量子化された整数となる。
4. 比較分析： 著者らは、標準的な格子ゲージ理論の手法（文献[34]）を用いて同じモデルの第2チェルン数 $ch_2$ を計算し、その相図を比較している。
5. 磁気電気結合： 本研究では、DDKS数と、3Dにおける軌道磁気電気結合を支配するチェルン・サイモンズ・アキシオン角 $\theta_{CS}$ との関係を調査している。著者らは、ポンピングサイクルに沿った $\theta_{CS}$ の変化（解析的に導出）と、「高次ベリー位相」$\theta_{HBP}$（iMPS構成から導出）の変化を比較している。

主な貢献と結果

• $ch_2$ の明示的に量子化された計算： 著者らは、高次ベリー曲率とiMPSを用いたアルゴリズムを通じて、DDKS数が4Dチェルン絶縁体のトポロジカルな相図を、第2チェルン数 $ch_2$ の相図と正確に一致して生成することを実証した。この一致は、標準的な格子ゲージ理論による $ch_2$ の計算が顕著な有限グリッド誤差に苦しむ領域を含む、異なる質量領域 ($m$) においても成立する。iMPS法は、粗いグリッド（例：$4^3$）上でも正確な整数を与えるのに対し、格子ゲージ法は収束するために非常に細かいグリッド（例：$10^4$ 点）を必要とする。
• Künneth公式による理論的接続： 本論文は、Künneth公式を用いて、4D第2チェルン数と3D DDKS数の間の厳密な結びつきを確立している。4Dトーラス $T^4$ を $S^1 \times T^3$ として分解することにより、著者らは $T^4$ 上の第2チェルン類が、$H^1(S^1, \mathbb{Z})$ の生成元と $\nu$ によって測定される $H^3(T^3, \mathbb{Z})$ の類との積に対応することを示した。したがって、この特定の分解において、$ch_2 = \nu$ となる。
• 4Dディラックモデルにおける磁気電気結合： 4Dディラックモデルにおいて、著者らはポンピングサイクルに沿ったチェルン・サイモンズ・アキシオン角の変化 $\Delta\theta_{CS}$ が、高次ベリー位相の変化 $\Delta\theta_{HBP}$ によって正確に再現されることを見出した。これは、BZスライスの第1チェルン数が消失する系において、高次ベリー位相と磁気電気応答の間に直接的な対応関係があることを示唆している。
• ホップ絶縁体における反例： 著者らは、非自明なホップ不変量を持つ2バンドモデルであるホップ絶縁体を用いて、この対応関係をテストしている。この場合、チェルン・サイモンズ角はホップ不変量によって量子化されているが、計算された高次ベリー位相 $\theta_{HBP}$ は量子化されておらず、$\theta_{CS}$ とも一致しない。著者らは、この不一致を、ホップ不変量の「デリケートな」性質、すなわち、それは運動量空間における2バンドモデルに対してのみ厳密に定義されるという条件が、熱力学的極限における実空間のiMPS記述では保持されない可能性があることに帰している。

意義と主張
本論文は、高次ベリー曲率とiMPSを用いて、4D格子モデルの第2チェルン数を計算する、初の明示的に量子化された数値的計算を提供したと主張している。これは、高次のベリー曲率が高次の広がりを持つ状態の族のトポロジカル不変量のプローブとして機能し得るという理論的枠組みを検証するものである。

著者らは、結果を量子幾何学およびトポロジカルな性質の文脈の中に位置づけており、高次ベリー位相と磁気電気結合の対応関係は4Dディラックモデルでは成立するものの、ホップ絶縁体の例が示すように普遍的ではないことを指摘している。彼らは、本手法が、単一粒子ブロッホ記述が破綻する相互作用系において特に適していることを強調しているが、現在は自由フェルミオンのベンチマークに焦点を当てている。本研究は、合成次元や準周期的な駆動が、これらの高次元トポロジカル現象を実現し、探査するための実験プラットフォームとして有望であることを浮き彫りにしている。

論文は、3次元における高次ベリー曲率と物理的応答の関係、特に $\theta_{CS}$ への結合の普遍性に関する問題は依然として未解決であると控えめに結論付けている。これは、高次ベリー位相が標準的なチェルン・サイモンズ応答とは異なる幾何学的情報を捉えている可能性を示唆しており、非線形輸送や基底状態の物理における役割についてさらなる調査を求めている。