[論文レビュー] Angle-Dependent {\it Ab initio} Low-Energy Hamiltonians for a Relaxed Twisted Bilayer Graphene Heterostructure

Shiang Fang, Stephen Carr+3 more|arXiv (Cornell University)|Jul 31, 2019

Graphene research and applications被引用数 13

ひとこと要約

本稿では、密度汎関数理論（DFT）とWannier関数変換、k·p展開を組み合わせることで、歪みのないねじれステッキングバイレイヤーグラフェン（tBLG）の角度依存のab initio低エネルギーハミルトニアンを提示する。この手法は、角度依存の擬似ゲージ場と運動量に依存する層間散乱を捉え、魔法の角度近くでのフラットバンドおよび粒子・空孔非対称性の正確なモデル化を可能にし、ねじれ2次元ヘテロ構造における相互作用電子理論の構築に不可欠である。

ABSTRACT

We present efficient angle-dependent low-energy Hamiltonians to describe the properties of the twisted bilayer graphene (tBLG) heterostructure, based on {\it ab initio} calculations of mechanical relxation and electronic structure. The angle-dependent relaxed atomic geometry is determined by continuum elasticity theory, which induces both in-plane and out-of-plane deformations in the stacked graphene layers. The electronic properties corresponding to the deformed geometry are derived from a Wannier transformation to local interactions obtained from Density Functional Theory calculations. With these {\it ab initio} tight-binding Hamiltonians of the relaxed heterostructure, the low-energy effective theories are derived from the projections near Dirac cones at K valleys. For twist angles ranging from 0.7$^\circ$ to 4$^\circ$, we extract both the intra-layer pseudo-gauge fields and the inter-layer coupling terms in the low-energy Hamiltonians, which extend the conventional low-energy continuum models. We further include the momentum dependent inter-layer scattering terms which give rise to the particle-hole asymmetric features of the electronic structure. Our model Hamiltonians can serve as a starting point for formulating physically meaningful, accurate interacting electron theories.

研究の動機と目的

• 機械的リラクゼーション効果を含む、角度依存の高精度な低エネルギー有効ハミルトニアンを、ねじれステッキングバイレイヤーグラフェン（tBLG）に対して開発すること。
• 従来の連続体モデルの限界を克服し、ab initioで得られた電子構造とリラクゼーションに起因する擬似ゲージ場を組み込むこと。
• 運動量に依存する層間散乱と粒子・空孔非対称性を捉えることで、物理的に意味のある相互作用電子理論の構築を可能にすること。
• 高精度かつ物理的に透明な多スケールフレームワークを提供し、tBLGおよび関連するファンデルワールスヘテロ構造を体系的にモデル化すること。
• 再現可能性を支援し、他の2次元ヘテロ構造への応用を可能にするために、導出されたモデルスクリプトを公開すること。

提案手法

• VASPを用いたab initio DFT計算を実施し、PBE汎関数に加え、ファンデルワールス補正にSCAN+rVV10を用いてtBLGのリラクゼーションされた原子配置を決定する。
• 9×9グリッドを用いた一般化されたスタッキングフォールトエネルギー（GSFE）計算により、面内および面外の変形を含むリラクゼーションされた原子位置を特定する。
• DFTバンドの最大局在化Wannier関数（MLWF）変換を用いてab initioタイトバインディングハミルトニアンを構築し、物理的に透明でパラメータフリーな記述を実現する。
• K谷の近傍におけるk·p展開を適用し、角度依存の擬似ゲージ場と層間結合を捉える低エネルギー有効ハミルトニアンを導出する。
• 実験的に観測された粒子・空孔非対称性を再現するために、運動量に依存する層間散乱項を導入する。
• 0.7°から4°の範囲のねじれ角度においてモデルを検証し、魔法の角度付近でのフラットバンドやバンドギャップといった既知の特徴と整合性を保つ。

実験結果

リサーチクエスチョン

• RQ1 ねじれステッキングバイレイヤーグラフェンにおける機械的リラクゼーションは、異なるねじれ角度において低エネルギー電子構造にどのように影響を与えるか？
• RQ2 リラクゼーションされたtBLGの低エネルギーハミルトニアンに現れる角度依存の擬似ゲージ場および層間結合項は何か？
• RQ3 運動量に依存する層間散乱項は、電子スペクトルにおける粒子・空孔非対称性にどのように寄与するか？
• RQ4 DFTおよびWannier変換から導かれたab initioタイトバインディングモデルは、実験で観測されたフラットバンド特徴およびバンドギャップを正確に再現できるか？
• RQ5 導出されたハミルトニアンは、ねじれ2次元ヘテロ構造における相互作用電子理論の構築にどの程度有効な基盤を提供できるか？

主な発見

• リラクゼーションされたtBLG構造は、格子エネルギーの最小化に起因し、顕著な面内および面外変形を示すが、その変形パターンはねじれ角度に強く依存する。
• 低エネルギーハミルトニアンに角度依存の擬似ゲージ場が出現し、再正規化されたディラックコーン分散およびバンドトポロジーの記述に不可欠である。
• 運動量に依存する層間散乱項が、電子構造における粒子・空孔非対称性の主な起因であると特定され、実験的観測と整合する。
• モデルはフラットバンドの電子および空孔側における単粒子ギャップを正確に再現しており、リラクゼーション効果に敏感であり、実験データと一致する。
• 従来の連続体モデルを拡張し、フィッティングなしに角度依存のab initioパラメータを含めるため、相関状態の予測精度が向上する。
• このフレームワークは一般化可能であり、オープンソーススクリプト（https://github.com/stcarr/kp_tblg で公開）として実装されており、トリレイヤーグラフェンや遷移金属ジ chalcogenideスーパーラティスなど、他の2次元ヘテロ構造への応用が可能である。