Plasmons in twisted bilayer graphene across dispersive and flat bands

技術要約：分散型バンドおよびフラットバンド領域におけるツイスト二層グラフェンのプラズモン

問題提起
ツイスト二層グラフェン（TBG）の動的な誘電応答は、電荷ダイナミクス、集団励起、およびマジック角系における潜在的なペアリング機構を理解する上で極めて重要である。理論モデル（連続体近似またはタイトバインディング）は様々なプラズモンモードを予測してきたが、大角度およびマジック角近傍の角度を持つTBGにおける運動量依存の誘電応答の完全な原子論的初手計算（first-principles calculations）は、計算コストの面で極めて困難であった。モアレ超格子における原子数は、ツイスト角の減少に伴い急速に増加するため（第1のマジック角 $\theta \approx 1.1^\circ$ 付近では $\sim 10^4$ 個の原子に達する）、標準的な平面波（PW）密度汎関数理論（DFT）および時間依存密度汎関数理論（TDDFT）による計算は非常に過酷なものとなる。本研究は、電子バンドが分散型（大角度）から準フラット（小角度）のレジームへと遷移する際の、プラズモンスペクトルの原子論的な記述におけるこのギャップに対処するものである。

手法
著者らは、ランダム位相近似（RPA）内での時間依存密度汎関数理論（TDDFT）を用いて、マクロな誘電関数および電子エネルギー損失関数（ELOSS）を算出している。

• 構造モデル: 3つの整合的なTBG構成を分析している：
  1. 分散型バンド・レジーム（大角度）: 基準層間距離 $d_0 = 3.349$ Å における $\theta_1 = 21.7868^\circ$ および $\theta_2 = 13.1736^\circ$。
  2. フラットバンド・レジーム: 減少した「マジック距離」$\tilde{d}_7 = 2.587$ Å で処理された、より小さな角度 $\theta_7 = 4.4085^\circ$。この「マジック距離構成」は、実際の $\theta \approx 1.1^\circ$ の構造に必要な膨大な超格子を必要とせずに、マジック角TBGの準フラットバンド多様体を再現するために設計された制約付きモデルである。
• 計算フレームワーク: 計算には、PBE汎関数およびノルム保存型擬ポテンシャルを用いたQuantum ESPRESSOスイートを利用している。誘電応答はカスタムのインハウス・フレームワークを用いて計算される。$\theta_7$ 構成の巨大な超セルを扱うため、著者らは2次元応答を孤立させ、メモリ制約を管理するために、応答計算中の波動関数の運動エネルギーカットオフを制限するスラブ切断型クーロンカーネルを採用している。
• 解析: 研究は損失関数 $E_{LOSS}(q, \omega) = -\text{Im}[\epsilon_M(q, \omega)^{-1}]$ に焦点を当てている。プラズモン特性は、損失共鳴における誘電関数の実部と虚部の検討、特に $\text{Re}[\epsilon_M]$ のゼロ交差とそれに伴う小さな $\text{Im}[\epsilon_M]$ の観察を通じて評価される。

主な結果

1. 分散型バンド・レジーム（大角度）:

   • 固有応答: 電荷中性において、$\theta_1$ および $\theta_2$ 構造はグラフェンの半金属的なディラック特性を保持している。損失スペクトルは、グラファイトに対してレッドシフトした、強く分散する $\pi$ プラズモン（長波長で $\sim 4.4$ eV）によって支配されており、これは $\pi \to \pi^*$ 遷移に由来する。より弱い特徴は、モアレ折り畳みによって誘起されたvan Hove 特異点（VHS）に対応する。
   • ドープ応答: 電子ドープにより、$\sqrt{q}$ 分散を持つ二次元ディラックプラズモン（2DP）が活性化される。本研究は、$\theta_1$ と $\theta_2$ の構成間の近似的な幾何学的スケーリング関係 $\omega_{2DP}^{\theta_1}(q_1) \approx s_{1/2} \omega_{2DP}^{\theta_2}(q_2)$ （ここで $s_{1/2} \approx \sqrt{19/7}$ はモアレ格子定数の比）を確立している。このスケーリングは線形ディラックコーン領域内ではよく成立するが、非線形および異方的なバンド効果により、より高いドープレベルでは崩れる。
   • アコースティック様ブランチ: 単一粒子励起連続体の中に埋め込まれたアコースティック様励起（AE）ブランチが観察される。これらは激しくランダウ減衰しており、かつ運動量方向に非常に敏感であるため、明確に定義された層反対称アコースティックプラズモンとしての同定を妨げている。

2. フラットバンド・レジーム（マジック距離構成）:

   • 電子構造: 層間距離を $\tilde{d}_7$ に減少させることで、フェルミレベル付近に狭く準フラットな多様体（four bands）が形成され、顕著な状態密度（DOS）ピークが生じる。
   • 低エネルギー励起: 固有充填において、低エネルギーの特徴（フラットバンド励起、FBE）が観察される。しかし、誘電解析によれば、損失最大値において $\text{Re}[\epsilon_M]$ の明確なゼロ交差は見られず、これはこれが明確に定義されたプラズモンではなく、準フラット多様体および高DOS領域に関わるバンド内およびバンド間遷移によって駆動される集団的性質を持つ励起であることを示唆している。
   • バンド間プラズモン: 異なる有限エネルギーのピーク（光学的プラズモン、OP）が $\omega_{OP} \approx 0.40$ eV で特定されている。このモードは明確な誘電ゼロ交差を示し、準フラット多様体と隣接する高DOS領域との間のバンド間遷移に起因するとされている。
   • 充填依存性: FBEは準フラットバンドが完全に満たされると消失し、部分的な占有への依存性を裏付けている。一方、OPは、関連するバンド間遷移チャネルが利用可能であるため、すべての充填条件（固有、電子ドープ、正孔ドープ）において持続するが、そのエネルギーとスペクトル重みは変化する。

意義および主張
本論文は、分散型およびフラットバンドの両レジームにおけるTBGの誘電応答に関する、初の完全な原子論的TDDFT-RPAによる特性評価を提供すると主張している。

• 手法的な貢献: 本研究は、「マジック距離」構成が、実際のマジック角における膨大な $\sim 10^4$ 個の原子を必要とする超セルをシミュレートすることなく、マジック角TBGのプラズモンエネルギー尺度を推定するための計算可能な経路を提供することを実証している。
• 物理的洞察: 結果は、電荷ダイナミクスの階層性を浮き彫りにしている（グラフェン由来の応答のような堅牢な $\pi$ プラズモンは、大角度TBGにおいて持続する）。一方で、フラットバンド・レジームは、準フラット多様体の特定のトポロジーと高DOS領域との結合によって支配されている。
• 超伝導への影響: 著者らは、特定された集団モード、特にフラットバンド・レジームにおけるバンド間プラズモンが、マジック角TBGにおけるクーパー対形成へのプラズモン媒介の寄与を評価するための出発点を提供すると示唆している。彼らは、これらの周波数依存の遮蔽相互作用を超伝導DFTに組み込むことで、プラズモン媒介とフォノン媒介のペアリングメカニズムの相互作用を解明できる可能性があると提案している。

本論文は、マジック距離モデルに必要な構造圧縮が実験的に検証される必要があること、および現在の制約付きジオメトリには構造緩和効果が含まれていないことを述べ、実験的な実現可能性については控えめな姿勢を保っている。