Unified description of exciton, phonon and plasmon dispersions in 2D materials from an optical conductivity approximation

技術要約：2D材料におけるエキシトン、フォノン、およびプラズモンの分散に関する統一的記述

問題提起
二次元（2D）材料において、長距離クーロン相互作用の変調された挙動（$v_{2D}(q) \sim 2\pi/q$）は、エキシトン、フォノン、プラズモンを含む縦方向の基本励起の分散を根本的に変化させる。LO/TO分裂が$\Gamma$点で発生するバルクの3D結晶とは異なり、2D材料では$q=0$においてこの分裂の崩壊が見られる。その結果、縦モードは非解析的な線形分散（$\propto q$）を示す一方で、横モードは放物線型の分散（$\propto q^2$）を維持する。

これらの分散を正確に記述するには、通常、高密度の運動量格子上で応答関数を計算する必要があるが、これはベテ・サルピーター方程式（BSE）やGW近似のような多体摂動論的手法においては計算コストが極めて高い。さらに、2Dのマクロなクーロン電場の影響下にある、電子的（エキシトン、プラズモン）およびイオン的（フォノン）な縦方向励起の異なる振る舞いを同時に考慮した統一的な理論枠組みが欠けていた。

手法：光学伝導度近似（OCA）
著者らは、有限運動量（$q$）の密度密度応答関数 $\chi(q, \omega)$ を、ゼロ運動量における光学伝導度 $\sigma(\omega) \equiv \sigma(q=0, \omega)$ に結びつける統一的な枠組みである「光学伝導度近似（OCA）」を導入する。

この導出は、連続の式と、誘起電流および全マクロ電場の間の関係に基づいている。主要な理論的洞察は、可約な密度密度応答 $\chi(q, \omega)$ はマクロなハートリーポテンシャルによる特異な$q$依存性を示すが、伝導度 $\sigma(q, \omega)$ は本質的にこのマクロなクーロン発散から自由であるという点にある。したがって、$\sigma(q, \omega)$ は$\Gamma$点近傍ではるかに緩やかな運動量依存性を示す。

OCAは、$\sigma(q, \omega) \approx \sigma(0, \omega)$ と仮定する。この近似を、連続の式から導かれた以下の厳密な関係式に代入することで、OCAの式が得られる：
$$\chi(q, \omega) = \frac{\sigma(q, \omega)|q|^2}{-i\omega - v_{2D}(q)\sigma(q, \omega)|q|^2}$$
$$\chi_{OCA}(q, \omega) = \frac{\sigma(\omega)|q|^2}{-i\omega - v_{2D}(q)\sigma(\omega)|q|^2}$$
ここで、$v_{2D}(q) = 2\pi/q$ である。この定式化により、有限$q$の損失関数（$-\text{Im}[\chi]$に比例）を、$q=0$の光学伝導度のみを用いて計算することが可能になる。伝導度は、ローレンツ型の共鳴の和としてモデル化された電子的な寄与（$\sigma_{el}$）と振動的な寄与（$\sigma_{ph}$）に分解される。

主な結果
著者らは、単層および二層の六方晶窒化ホウ素（hBN）とグラフェンに対してOCAを適用し、有限運動量格子を用いて行われたab-initio計算（エキシトンにはBSE、フォノンにはDFPT、プラズモンにはGW-BSE）と比較した。

1. hBNのフォノンとエキシトン： OCAは、単層hBNにおける縦光学（LO）フォノンおよび低エネルギーの輝くエキシトン（AおよびB枝）のエネルギーおよび強度の分散を正確に再現した。これらのモードの線形分散およびB枝とA枝の強度比の急速な増加は、有限$q$での$\sigma$を計算することなく捉えられている。
2. グラフェンのプラズモン： 本フレームワークは、グラフェンの$\pi$プラズモンの線形分散を正常に再現した。結果は、この分散がバンド構造の分散によるものではなく、2Dクーロン相互作用が$q=0$の光学伝導度に作用することによる$q$依存性によって駆動されていることを裏付けている。
3. モデル解析： $\sigma(\omega)$を単一のローレンツ型と非吸収性バックグラウンドで近似することにより、著者らは運動量依存のエネルギー $\Omega(q)$ および残差 $R(q)$ の解析的な表現を導出した。彼らは、線形分散係数が$q=0$における光学伝導度の残差に直接比例することを証明した。
4. 限界： OCAは$q=0$において厳密であり、$q$が小さい領域でも精度を保つ。しかし、暗い（横方向の）励起、例えばダークエキシトンや横光学（TO）フォノンについては記述できない。これらはマクロな縦方向電場と結合しないためである。さらに、本手法は、グラフェンの励起スペクトルの開始（onset）を捉えることはできない。この開始はディラックコーンの線形バンド構造によって決定される純粋なバンド構造の効果であるが、プラズモン分散は集団的な効果である。

意義および主張
本論文は、2Dの縦方向励起の統一的な記述を提供し、それらの鋭い線形分散が、基礎となるバンド構造の分散ではなく、主に2Dのマクロなクーロン相互作用の形式（$v_{2D} \sim 1/q$）の結果であることを示していると主張している。

本研究の意義は二重にある：

• 計算効率： OCAは、高価な計算を伴う高密度運動量格子を用いずに、正確な低運動量電子エネルギー損失分光（EELS）データを生成するための効率的なスキームを提供する。これは、単一の$q=0$計算を用いて、ab-initioの点間を効果的に補間するものである。
• 逆手法： このフレームワークは逆手法として機能し、研究者が低運動量EELSスペクトルから完全な複素光学伝導度を抽出することを可能にする。これにより、標準的な光学測定における高エネルギー到達範囲や試料純度の要件に関する制限を回避し、広範なエネルギー領域にわたる光学特性へのアクセスが可能になる。

著者らは、このスキームが半導体、半金属、および金属系に適用可能であり、減少次元におけるプラズモン、エキシトン、およびフォノンの分散に対して一貫した物理的描像を提供すると述べている。