Review of the tight-binding method applicable to the properties of moiré superlattices

本レビューは、様々なモアレ超格子の電子、輸送、および光学特性をモデリングするための原子論的タイトバインディング法と数値的手法に関する包括的な理論的・実践的ガイドを提供すると同時に、それらと有効低エネルギー連続体モデルとの関連性を明確にするものである。

要約

二枚の透明な、ハニカム構造（グラフェンのような）を持つプラスチックのシートを想像してみてください。これらを完全に重ね合わせると、一枚のシートのように見えます。しかし、片方のシートをわずかに回転させたり、少し引き伸ばしたりすると、パターンが一致しなくなります。その結果、個々のパターンよりもはるかに大きな、渦巻くような干渉パターン、すなわち**モアレ・スーパーラティス（モアレ超格子）**が形成されます。

これは、窓のスクリーンを二枚重ねて持ち、片方をひねる様子を想像すると分かりやすいでしょう。スクリーンの小さな穴よりもはるかに大きな、ゆっくりと動く巨大な波のパターンが現れます。原子の世界では、この「巨大な波」がある場所で、電気抵抗ゼロで電気が流れる超伝導や、物質が磁性を持つといった、非常に不思議で魔法のような物理現象が起こります。

しかし、この巨大な原子の波を研究することは、コンピュータにとって悪夢です。なぜなら、このパターンの一つの「単位」には数千の原子が含まれているからです。この巨大な群衆の中のすべての原子の挙動を計算しようとするのは、スタジアムにいる一人ひとりに質問して、全員の動きを予測しようとするようなものです。それには時間がかかりすぎ、膨大なメモリを必要とします。

この論文は、タイトバインディング（強結合）法と呼ばれる特定の「近道」に関するガイドブックです。以下に、簡単な比喩を用いて、この論文がどのように説明しているかを記します。

1. 問題点：多すぎる原子

論文では、小さな原子グループを研究するための強力なツール（密度汎関数理論、またはDFT）は存在するものの、これら巨大なモアレ・パターンに対しては遅すぎると指摘しています。一方で、単純な数学モデル（連続体モデル）は高速ですが、原子が物理的にどのように位置を変えたり緩和したりするかといった、微細な詳細を見落としてしまいます。

2. 解決策：タイトバインディングによる「近所」マップ

タイトバインディング法は、近所の地図のようなものです。スタジアム全体の物理学を一度に計算するのではなく、ある原子がそのすぐ隣にある近隣の原子（あなたのすぐ隣に座っている人）とどのように相互作用しているかだけを見ます。

• 仕組み: 原子は、主に隣人が誰であるか、そしてそれらがどれくらい離れているかによって、その挙動が決まると仮定します。
• 利点: これにより、個々の原子の詳細（原子が押しつぶされているか、あるいは引き伸ばされているか）を保持しつつ、数千もの原子を扱うのに十分な速さを維持できます。これは「ゴルディロックス（ちょうど良い）」ゾーンです。単純すぎず、かといって遅すぎもしません。

3. ツールキット：材料ごとの異なるマップ

論文では、これら三つの主要な材料に対して、どのように「近所の地図」を構築するかをレビューしています。

• グラフェン（炭素のハニカム構造）: マップは比較的単純で、炭素原子の間を電子がどのように「ホッピング（跳躍）」するかに焦点を当てています。論文では、マップ内の原子間の「距離」を調整することで、科学者がいつ「マジックアングル」の超伝導体になるかを正確に予測できることを示しています。
• TMDs（遷移金属ダイカルコゲニド）: これらは金属と他の元素が組み合わさった複雑なサンドイッチのようなものです。物理を正しく捉えるためには、より詳細なマップ（11種類の異なる「軌道」または電子の経路を使用するもの）が必要になります。
• hBN（六方晶窒化ホウ素）: これは他の材料の滑らかなベッド（基板）としてよく使われます。論文では、グラフェンの炭素原子と、このベッドのホウ素および窒生原子との間の相互作用をどのようにマッピングするかを説明しています。

4. 数学の処理：「ランダムウォーク」のトリック

モアレ・パターンが巨大になり（数百万の原子を含む場合）、近所のマップですら直接解くには大きすぎる場合があります。論文では、線形スケーリング法（カーネル多項式法など）と呼ばれる巧妙なトリックを紹介しています。

• 比喩: スタジアムにいる全員の平均身長を知りたいとします。全員を測る必要はありません。代わりに、数人のランダムな人々を選んで測定し、統計的な公式を使って群衆全体の平均を推測するのです。
• 結果: これにより、科学者は標準的なコンピュータ上で、光が材料とどのように相互作用するか、あるいは電気がどのように流れるかといったことを、数百万の原子を伴うシミュレーションとして計算できるようになります。

5. 「緩和」の魔法

この論文の重要なポイントの一つは、原子は静止した彫像ではなく、動き、心地よい位置に落ち着く（緩和する）ということです。

• 比喩: グリッド状に並んで立っている人々の群衆を想像してください。グリッドをひねると、真ん中にいる人々はスペースを節約するために身を寄せ合い、端にいる人々は広がったりするかもしれません。
• 知見: タイトバインディング法が特別なのは、この「身寄せ（緩和）」を考慮できる点です。論文では、もしこの緩和を無視すれば、誤った物理現象を得ることになると示しています。緩和を含めることで、「フラットバンド（電子が動きにくくなり、強く相互作用し始めるエネルギー準位）」を正確に予測でき、それが超伝導のようなエキゾチックな現象につながるのです。

6. 論文における実例

著者らは、この手法を二つの具体的なストーリーで実証しています。

• 12角形の結晶: 彼らは、12角形（十二角形）のパターンを形成する、ねじれたグラフェン構造を研究しました。このパターンは単純な繰り返しではないため、標準的な数学では通用しません。タイトバインディング法は、「ランダムウォーク」のトリックを用いることで、このユニークな形状において光や電気がどのように振る舞うかを予測することに成功しました。
• 閉じ込められたエキシトン: 彼らは、WSe2の層がねじれたグラフェンの上に乗っているシステムを調査しました。彼らは、グラフェンの原子の「身寄せ（緩和）」が、リドバーグ・エキシトン（一種の励起粒子）を捕まえ、保持する小さなトラップ（罠）を作り出す様子を示し、実験で見られる特定の信号を説明しました。

まとめ

要するに、この論文は、巨大でねじれた原子パターンを理解するための、特定の種類のコンピュータモデルを構築し、使用するためのマニュアルです。著者は、タイトバインディング法が、個々の原子の動きや緩和を見るのに十分な詳細さを持ちつつ、巨大なモアレ・スーパーラティスの膨大なサイズを扱えるほど十分に高速であるため、この仕事に最適なツールであると主張しています。それは、単純で高速な理論と、低速で超精密なシミュレーションとの間の架け橋となるものです。

技術概要

モアレ超格子特性に適用可能なタイトバインディング法のレビュー

問題提起
二次元（2D）材料を相対的な回転や格子不整合を伴って積層することで形成されるモアレ超格子は、非従来型超伝導、相関絶縁状態、トポロジカル効果を含む、エキゾチックな量子現象を探索するための多才なプラットフォームとして台頭している。しかし、これらの系の理論的モデリングは大きな課題に直面している。オングストロームスケールの材料とは異なり、モアレ系は長波長の干渉パターンにより、非常に多くの原子（しばしば数千から数百万個）を含む。さらに、その電子構造は格子緩和や原子再構成に決定的に依存する。密度汎関数理論（DFT）は高い精度を提供するが、その計算コストにより、大規模なモアレ超格子のシミュレーションには非効率である。逆に、連続体モデルは有効な低エネルギー記述を提供するが、局所的な無秩序、歪み、または特定の化学的効果を捉えるために必要な原子解像度が欠けている。したがって、原子論的シミュレーションの精度と連続体モデルの効率性の間の溝を埋める理論的枠組みが必要とされている。

手法
本レビューでは、モアレ材料をモデリングするための中心的なツールとしてのタイトバインディング（TB）法に焦点を当てる。本論文は、グラフェン系、遷移金属ジカルコゲニド（TMD）系、および六方晶窒化ホウ素（hBN）系の主要な材料ファミリーに対する原子論的TBハミルトニアンの構築を体系的に分類している。

1. ハミルトニアンの構築:

   • グラフェン: 標準的なアプローチでは、単一粒子 $p_z$ 軌道モデルを利用する。ホッピングパラメータ（$t_{ij}$）は、原子間距離と減衰因子に依存するスレーター・ケスター（SK）関係を用いて決定される。このモデルは、緩和された原子位置に基づいてホッピング項を修正することにより、格子緩和を組み込んでいる。電子相互作用は、ハートリー・ポテンシャル（長距離クーロン相互作用のため）、ハートリー・フォック近似（交換相互作用を含む）、およびハバード $U$ 項（局所相関のため）を含む平均場近似を通じて扱われる。大規模なスーパーセルにおける自己整合計算の計算コストを削減するために、リスケーリング戦略が導入されている。
   • TMDs: TMDの複雑さのため、これらのモデルは通常、単層あたり11個の軌道基底（金属由来の5つの $d$ 軌道とカルコゲニド由来の3つの $p$ 軌道）を用いる。スピン軌道相互作用（SOC）はオンサイト項を介して組み込まれる。層間相互作用は、特定の軌道間（例：ホモバイレイヤーにおける $p$-$p$、ヘテロバイレイヤーにおける $p$-$d_{z^2}$）のホッピング項を加えることでモデル化され、そのパラメータはしばしばDFTデータからのワニエ変換によって導出される。
   • hBN: 回転積層hBNおよびグラフェン/hBNヘテロ構造のモデルは、ホウ素（B）および窒素（N）原子上の $p_z$ 軌道を利用する。これらのモデルは、BとNの原子間のオンサイトエネルギー差を考慮し、DFTの結果に適合させた様々なSKパラメトリゼーションを用いて層間ホッピングを扱う。

2. 数値的手法:
   大規模なハミルトニアン行列（しばしば $N > 10^5$ 原子）を扱うため、本レビューでは標準的な対角化法を超えた手法について論じている。

   • 線形スケーリング法: カーネル多項式法（KPM）およびタイトバインディング伝搬法（TBPM）が強調されている。これらの確率論的アプローチは $O(N)$ でスケールするため、数百万個の原子を持つ系における状態密度（DOS）、光導電率、および輸送特性の計算を可能にする。
   • 機械学習: 学習データからab initio品質のTBハミルトニアンを構築するための、ディープラーニング（例：DeepH, HamGNN, DeepTB）を用いた最近の開発がレビューされている。これはパラメトリゼーションのプロセスを加速させる。

3. 連続体モデルとの接続:
   本論文は、原子論的TBハミルトニアンから低エネルギー有効連続体モデルを導出する詳細を述べている。モアレポテンシャルが層間トンネリングのフーリエ成分からどのように生じるかを説明している。本レビューは、局所的（運動量に依存しない）ポテンシャルがフラットバンドの出現を捉える一方で、緩和されたTBモデルで見られる粒子ホール非対称性などの特徴を再現するには、非局所的（運動量依存）なトンネリング項が必要であることを強調している。

主な貢献と結果

• 包括的な概要: 本論文は、グラフェン、hBN、およびTMDモアレ系に使用される特定のTBハミルトニアン、SKパラメータ、および緩和ポテンシャル（例：LAMMPSポテンシャル）に関する統一的なリファレンスを提供している。
• TBモデルの検証: 本レビューは、適切にパラメトリゼーションされ、緩和を含むTBモデルが、実験およびDFTの結果を再現することを示している。例えば、ねじれ二層グラフェン（TBG）のTB計算は、フェルミ速度が消失しフラットバンドが出現する「マジック角度」（$\theta \approx 1.05^\circ - 1.1^\circ$）を正常に予測している。
• 格子緩和の役割: 著者らは、原子緩和が単なる微小な補正ではなく、極めて重要な要因であることを強調している。TBGにおいて、緩和はAAドメインおよびABドメインの形成をもたらし、バンド構造を大きく変化させ、粒子ホール非対称性を誘起する。TMDやhBNにおいても、同様に緩和がバンドの平坦さとバンドギャップの大きさを決定する。
• 相互作用の取り扱い: 本レビューは、ハートリー、ハートリー・フォック、ハバードを含む平均場TBアプローチが、フラットバンド領域におけるモット絶縁状態や強磁性などの相関相の研究にどのように利用できるかを示している。
• ケーススタディ: TB法の有用性を実証するために、2つの具体的な例が提供されている：
  1. 十二角形グラフェン準結晶: 1000万個の原子を持つ系に対してTBPMを用いることで、層間相互作用に起因するDOSおよび光導電率の明確なピークを明らかにし、磁場中での新しいランダウ準位を予測している。
  2. WSe$_2$/TBGにおけるリドバーグ・モアレ励起子: TBを分子動力学法と組み合わせることで、緩和されたTBGのAA領域における電荷蓄積によるWSe$_2$内のリドバーグ励起子の空間閉じ込めを説明し、実験の反射率およびフォトルミネッセンススペクトルと一致させている。

意義と範囲
本論文は、タイトバインディング法をモアレ超立体の理論的研究のための中心的かつ実践的なツールとして位置づけている。その意義は、格子緩和、歪み、および化学的特異性をモデル化するために必要な原子レベルの解像度を保持しつつ、数千から数百万の原子を含む現実的なスーパーセルをシミュレートできるほど計算効率が高いという、独自のバランスにある。

著者らは、本レビューがモアレ系にTB法を適用する研究者のための「理論的および実践的なガイド」として機能すると述べている。これは、ab initio DFT計算と有効連続体モデルの間の溝を埋め、多体相互作用や外部場を組み込むための体系的な方法を提供している。本論文は、TB法は強力ではあるものの、非六方晶格子への手法の拡張、非直交基底に対する線形スケーリングアルゴリズムの改善、および新しいモアレ材料の自動発見のための機械学習の活用といった今後の課題があることを述べて締めくくっている。本レビューは、すべての相関問題を解決すると主張するものではなく、TBがフラットバンド関連現象の起源を理解するための、必要かつ正確な出発点を提供することを主張している。