[論文レビュー] Lattice Deformation, Low Energy Models and Flat Bands in Twisted Graphene Bilayers

Niels R. Walet, F. Guinea|arXiv (Cornell University)|Mar 1, 2019

Graphene research and applications参考文献 26被引用数 7

ひとこと要約

本論文は、θ = 1.05°におけるねじれバイレイヤー石墨烯における原子スケールの格子緩和および環境依存の層間 hopping が、低エネルギー電子バンドおよびフラットバンド形成に顕著に与える影響を調査する。実際の力場ポテンシャルを用いたタイトバインディングモデルと洗練された連続体近似を用いることで、特に hopping パラメータが局所的な原子環境に依存する場合に、運動量に依存する hopping および層間結合における3次以上の調和成分を含めることが、正確なバンド構造記述に不可欠であることが示された。これにより、有効なミラクル角度のシフトおよび新しいディラック点の出現が生じる。

ABSTRACT

Twisted graphene bilayers show a complex electronic structure, further modified by interaction effects. The main features can be obtained from effective models, which make use a few phenomenological parameters. We analyze the influence of effects at the atomic scale, such as interlayer hopping and lattice relaxation, on the electronic bands. We assume that the twist angle and the size of the Moiré pattern is fixed, as it is usually the case in experiments. We obtain a strong dependence of the electronic structure on details of the models at the atomic scale. We discuss how to incorporate this dependence on effective models

研究の動機と目的

• 原子スケールの格子緩和および層間 hopping の詳細が、ねじれバイレイヤー石墨烯の電子バンド構造に与える影響を理解すること。
• 原子スケール効果を含めた場合に、標準的な連続体モデルが低エネルギー物理学をどれだけ正しく捉えられるかを評価すること。
• 環境依存の hopping パラメータおよび格子変形が、フラットバンドの出現および有効なミラクル角度に与える影響を特定すること。
• 運動量に依存する行列要素および高次調和成分を組み込むことで、タイトバインディング結果を正確に再現する洗練された連続体モデルを構築すること。
• 実験的に測定されたねじれ角度と理論的に定義されたミラクル角度の乖離を解明すること。

提案手法

• θ = 1.05°（共通のねじれ角）におけるねじれバイレイヤー石墨烯の面内格子緩和を、古典的力場モデル（LCBOP-I+KC および AIREBO-M+ILP）を用いて計算した。
• 原子スケールの相互作用に基づく層間 hopping パラメータ（距離および環境依存形）を用いたタイトバインディングモデルにより電子構造を計算した。
• 超格子調和成分を用いた層間 hopping の連続体展開を実施し、標準的な3調和モデルを越えて最大48調和まで含めた。
• 層間 hopping の行列要素の運動量依存性を組み込み、平均運動量および逆格子ベクトルの差（Gₘ - Gₙ）に依存する項を明示的に取り入れた。
• タイトバインディング結果と連続体近似を比較し、収束性および正確性を評価し、標準モデルがどこで失敗するかを特定した。
• 三角対称歪みおよび格子変形が各層のディラックコーンに与える影響を分析し、ディラック方程式への補正を加えた。

実験結果

リサーチクエスチョン

• RQ1 θ = 1.05°における格子緩和は、ねじれバイレイヤー石墨烯における AB、BA、および AA スタッキング領域の空間的分布にどのように影響するか？
• RQ2 環境依存の層間 hopping パラメータ（例：γ₃ ≠ γ₄）が、距離のみに依存するモデルと比較して、低エネルギーバンド構造にどの程度の影響を及ぼすか？
• RQ3 原子スケールの詳細（格子緩和や運動量に依存する hopping）を含めた場合、標準的な3調和連続体モデルが電子バンドを正確に記述できるか？
• RQ4 環境依存の hopping を有する場合に、タイトバインディング結果を定量的に再現するためには、連続体モデルでどの程度の調和数が必要か？
• RQ5 観測されたバンド特徴（新しいディラック点やバンドクロスイング）は、力場モデルおよび hopping モデルの選択にどのように依存するか？

主な発見

• LCBOP-I+KC および AIREBO-M+ILP 力場モデルによる予測によると、θ = 1.05°における格子緩和は、広範囲にわたる大きな AB および BA スタッキング領域を生じさせ、AA スタッキング領域は顕著に減少する。
• γ₃ と γ₄ を区別する環境依存の hopping パラメータは、距離のみに依存するモデルよりも、整列したバイレイヤー石墨烯およびグラファイトの電子的性質をより正確に再現する。
• 標準的な3調和連続体モデルは、特に hopping が局所的環境に依存する場合に、タイトバインディング計算で観測される粒子-hole 対称性の破れおよびバンドクロスイングを捉えられていない。
• 行列要素の運動量依存性（k および Gₘ - Gₙ に依存）と、3調和を超える高次調和成分（最大48調和）を連続体モデルに組み込むことで、タイトバインディング結果とほぼ完全に一致する。
• hopping パラメータが局所的環境に依存する場合、有効なミラクル角度がシフトし、新しいディラック点が出現する。これは、低エネルギー物理学が顕著に変化していることを示唆する。
• これらの結果は、ねじれバイレイヤー石墨烯における超伝導および絶縁相の出現が、現実的な原子スケールパラメータの範囲で安定しており、微調整されたモデル仮定に依存しない一般的メカニズムである可能性を示唆する。