Deep-learning Hamiltonian reveals twist-tunable flat bands and nonlinear photocurrents in SrTiO3 moire bilayers

技術要約：深層学習ハミルトニアンが明らかにする、SrTiO$_3$ モアレ二層におけるねじれ制御可能な平坦バンドと非線形光電流

問題提起
モアレ物理学を複雑な酸化物へと拡張することは、電子状態を操作するための新たな道を切り拓くものである。しかし、酸化物モアレ系に必要とされる巨大なスーパーセルは、第一原理計算による系統的な解析を計算量的に極めて困難にしている。グラフェンや遷移金属ディカルコゲナイドのようなファンデルワールス材料では、ねじれ角の制御によって平坦バンドや相関状態の生成に成功しているが、SrTiO$_3$ のようなペロブスカイト酸化物に同様のエンジニアリングを適用することは困難である。SrTiO$_3$ は、格子、軌道、電荷、およびスピンの自由度が強く結合しているため、単純な連続体モデルやタイトバインディングモデルへの簡略化が難しい。さらに、密度汎関数理論（DFT）の計算コストはシステムサイズに対して急激に増大するため、ねじれ角に依存する電子構造および光学応答を系統的に記述することを阻んでいる。

手法
計算上のボトルネックを克服するため、著者らは DFT と E(3) 等変深層学習ハミルトニアン・フレームワークである DeepH-E3 を組み合わせた。

• 学習データの生成: OpenMX パッケージを用い、様々な層間スライディング構成を持つ未ねじれ二層 SrTiO$_3$ スーパーセルに対して DFT 計算を行った（計623構成）。これらの構造は、モアレ超格子に見られる局所的なスタッキング環境をサンプリングしている。
• モデルの学習: E(3) 等変グラフニューラルネットワークである DeepH-E-3 モデルを用いて、原子構造からハミルトニアン行列要素を直接予測するように学習させた。この際、スピン軌道相互作用を組み込み、空間対称性を保持している。モデルは DFT 参照データに対して検証され、ハミルトニアン要素において 0.26 meV という平均絶対誤差（MAE）を達成し、バンド構造においても優れた一致を示した。
• 適用: 学習済みモデルを、$53.13^\circ$ から $8.80^\circ$（最大 2,380 個の原子を含むスーパーセルに対応）までの範囲の整合的なねじれ二層 SrTiO$_3$ 系に適用した。
• 応答計算: 予測されたハミルトニアンを HopTB パッケージと連携させ、誘電応答、第二高調波発生（SHG）、シフト電流、および固有スピンホール伝導率（SHC）を計算した。

主な貢献と結果

1. ねじれ制御可能な平坦バンド:
   本研究は、ねじれ角が減少するにつれて価電子帯最大値（VBM）が系統的に進化することを明らかにしている。

• 大角領域（$\theta > 20^\circ$）では、バンドは比較的分散的であり、弱いフォールディングを示す。
• 中間および小角領域（$\theta \lesssim 15^\circ$）へと角度が減少すると、VBM バンドは著しく平坦化する。
• 研究された最小角（例：$8.80^\circ$）では、ほぼ分散のない平坦バンドが出現し、報告された正孔質量は未ねじれ二層と比較して 1 桁以上増加している。
• グラフェン系において平坦バンドが層間ホッピングの干渉から生じるのに対し、SrTiO$_3$ における平坦化は、スタッキングに依存した静電的および軌道的なハイブリダイゼーションに起因する。バンド幅と有効質量は非単調な変化を示しており、これはねじれ角そのものよりも、離散的な整合スタッキング幾何学に対する敏感さを示唆している。

2. 光学的および輸送応答への差異的な影響:
   著者らは、ねじれが未ねじれ二層に対する様々な物理的応答にどのように影響するかを評価した。

• 線形誘電応答: 誘電関数（実部および虚部の両方）は、計算された角度範囲全体で概ね変化しない。光学ギャップと主要な吸収特徴の変化は極めて小さく、この領域における線形光学構造がモアレ変調に対して堅牢であることを示唆している。
• スピンホール伝導率 (SHC): 固有のスピン・ベリー曲率に支配される SHC は、ねじれ角に対してほぼ不変である。ピーク位置と振幅はわずかな非単調変化を示すのみであり、大角から中間角の範囲において、積分されたスピン・ベリー曲率がモアレ摂動に対して鈍感であることを示している。
• 非線形光学応答 (SHG およびシフト電流): これらの応答は、ねじれによって強く影響を受ける。
  • 第二高調波発生 (SHG): ねじれは界面の反転対称性を破るため、未ねじれ二層の弱い応答と比較して、SHG の強い増強をもたらす。
  • シフト電流: シフト電流（第二次非線形光電流）は、明確なねじれ角依存性を示す。未ねじれ二層は弱い応答を示す一方で、ねじれた構造は顕著な増強を示す。シフト電流のピーク振幅は、ねじれ角が $53.13^\circ$ から $22.62^\circ$ に減少するにつれて大幅に増加する。

意義と主張
本論文は、ねじれエンジニアリングが、酸化物モアレ系（具体的には二層 SrTiO$_3$）の電子および光電子特性を制御するための効果的なメカニズムとして機能することを主張している。主な知見は以下の通りである：

• 電子構造: ねじれ角は、低エネルギー電子構造を系統的にチューニングでき、小角極限においてほぼ平坦な価電子帯を誘起できる。
• 光電子制御: 線形応答およびスピン輸送は、研究された範囲内ではねじれに対して比較的鈍感であるが、非線形光学応答は高度にチューナブルである。特に、シフト電流はねじれ角に対して最も明確な依存性を示しており、酸化物系における非線形光電流を制御するための新たな手段を提供する。
• 手法の妥当性: 本研究は、従来の DFT では実用的な規模ではない、複雑で強く結合した酸化物材料において、深層学習ハミルトニアン（DeepH-E3）を用いた大規模な電子構造計算が可能であることを実証した。

著者らは、ねじれた酸化物におけるすべての結合された応答の完全な取り扱いは依然として未解決の方向性であるが、本研究は、複雑な酸化物におけるねじれエンジニアリング現象を探索するために機械学習を使用することの実現可能性を成功裏に確立したと結論付けている。