First-principles calculation of electron-phonon spectral functions for defects using phonon interpolation

技術要約：フォノン補間を用いた欠陥の電子-フォノン分光関数の第一原理計算

問題の定義
ダイヤモンド中の負電荷を持つ窒素空孔（NV）中心のようなワイドバンドギャップ半導体における点欠陥は、電子-フォノン結合の影響を強く受けた光学スペクトルを示します。これらのスペクトルを正確に予測するためには、フォノン分解された電子-フォノン結合密度（フーガン・ライスの因子）および結果として得られる光学線形形状関数を計算する必要があります。しかし、直接的なab initio（第一原理）計算は、アクセス可能なスーパーセルの有限なサイズによって厳しく制限されます。$N$個のイオンを含むスーパーセルは、$3N-3$個の振動モードしか生成せず、その結果、フォノンスペクトルのサンプリングが疎になります。この疎なサンプリングは、低エネルギー領域に不可欠な長波長の音響連続体や、光学フォノン領域におけるファン・ホーブ特異点付近の微細構造を捉えることができず、これらは正確な実験的線形形状を再現するために不可欠です。これまでの手法では、埋め込み戦略や機械学習ポテンシャルを利用してより大きな系へのアクセスを試みてきましたが、それらは複雑なモデル構築を伴うか、あるいは学習済みポテンシャルの転移性に依存することが多くあります。

手法
著者らは、局所的な欠陥の物理と、振動連続体のサンプリングをデカップル（分離）するフォノン補間法を提案しています。この手法の核心は、ギャップを持つホスト内における局在化した欠陥状態間の電荷保存的な遷移に対して、遷移誘起力が欠陥の周囲に強く局在するという観察に基づいています。

1. 遷移力の再構成（Transition-Force Reconstruction）: 有限なスーパーセル（周期的な欠陥配列を表す）の緩和された変位場を使用する代わりに、本手法では局在化した遷移力源（$F^{tr}$）を再構成します。これは、同一の原子配置における初期状態と最終状態のボルン・オッペンハイマー力（Born-Oppenheimer forces）の差から導かれます。この力は、有限のスーパーセル（例：$4\times4\times4$）内での密度汎関数理論（DFT）を用いて計算されます。
2. 局在性と収束性（Locality and Convergence）: 遷移力は欠陥からの距離に対して指数関数的に減衰することが示されています。スーパーセルが欠陥軌道の周期的な像の重なりを抑制できるほど十分に大きければ、遷移力は急速に収束します。これにより、管理可能なサイズのab initioスーパーセルから得られた力場が、孤立した欠陥の正確なソース項として機能することが可能になります。
3. ハイパーセル補間（Hypercell Interpolation）: 本手法は、欠陥を含むスーパーセルを$M_1 \times M_2 \times M_3$回繰り返すことで「ハイパーセル」を構築します。力定数行列（動力学行列）は、ab initioスーパーセルからこのより大きなハイパーセルへと補間されます。これらの補間された動力学行列を密な$q$点グリッド上で対角化することにより、本手法は密な振動モード（フォノン）のセットを生成します。
4. 分光関数の計算（Spectral Function Calculation）: 局在化した遷移力が、この密な振動連続体に結合されます。部分的なフーガン・ライス因子（$S_{q,\nu}$）は、遷移力を補間されたフォノン固有ベクトルに投影することで計算されます。これにより、滑らかで連続的なフーガン・ライス分光密度$S(\hbar\omega)$および光学線形形状関数（$L_{em}$および$L_{abs}$）を、極めて微細なスペクトル分解能で算出することが可能となり、実質的に数百万個の原子（例：$32\times32\times32$ハイパーセル $\approx$ 1700万個の原子）をシミュレートすることと同等の効果が得られます。

主な貢献

• 定式化: 本論文は、光学線形の調和フーガン・ライス理論に直接つながるフォノン補間手順を定式化しています。この手法は、変位場の代わりに局在化した遷移力源を用いることで、補間された動力学行列を用いて振動連続体をサンプリングすることを可能にします。
• 効率性: 本手法は、大規模なスーパーセルの直接的なab initio計算に伴う膨大なコストを回避しながら、局所的な欠陥物理に対するab initioの記述精度を維持します。また、大規模な埋め込みモデルの構築や複雑なカットオフ力定数行列の管理を必要とする埋め込み法よりも、計算の導入が容易です。
• 実装: 著者らは、標準的なab initioコード（Phonopy経由）と連携して、力定数の後処理および振動特性の補外を行う公開Pythonソフトウェアパッケージ（pysideband）を提供しています。

結果
本手法は、ダイヤモンド中のNV$^-$中心、具体的には$^3A_2$基底状態と$^3E$励起状態の間の光学遷移に対してベンチマークが行われました。

• 力の局在化: $2\times2\times2$から$4\times4\times4$のスーパーセルにおける遷移力の解析により、強い局在性が確認されました。収束指標（参加数、RMS局在半径）は$4\times4\times4$スーパーセルで安定しており、これを孤立欠陥のソースとして用いることの妥当性が検証されました。
• 分光密度: $32\times3\times32$ハイパーセルへの補間により、滑らかなフーガン・ライス分光密度が得られました。結果は以下を示しました：
  • 63 meV付近に中心を持つ支配的な結合バンド。
  • 低エネルギーの音響領域における正しい線形スケーリング（$S(\hbar\omega) \propto \hbar\omega$）を実現し、非物理的な発散を回避。
  • 光学フォノンバンド端付近のファン・ホーブ特異点に関連する微細構造の解像。これらは疎なスーパーセル計算では通常見落とされるものです。
• 線形形状:
  • 放出（Emission）: 基底状態のフォノンを用いた計算による放出スペクトルは、低温での実験データと極めて良好に一致し、ゼロフォノン線、63 meVのフォノンレプリカ、およびサイドバンドのエンベロープを再現しました。
  • 吸収（Absorption）: 基底状態のフォノンを用いた計算では、特に最初のフォノンピークにおいて実験的な吸収プロファイルとの不一致が見られました。しかし、励起状態のフォノンを用いて吸収を再計算したところ、一致が大幅に改善され、実験的に観察される主要サイドバンドの平坦で弱く分裂した構造が回収されました。これは、吸収が$^3E$状態のヤーン・テラー活性に対して敏感であることを浮き彫りにしています。
• 収束性: 全フーガン・ライス因子（$S_{tot}$）は、ハイパーセルのサイズに従って系統的に収束することが冪乗則（power-law）アプローチによって示され、無限系の極限への補外が可能であることが示されました。

意義と主張
本論文は、フォノン補間法が、局在化した欠陥に対して収束した電子-フォノン分光関数を得るための、効率的かつ系統的に改善可能な経路を提供すると主張しています。これは、標準的なab initio欠陥計算における振動サンプリングの疎さという特有の制限に対処するものです。

著者らは、本アプローチが以下の特徴を持つことを強調しています：

• 局所的な欠陥緩和に対するab initioの記述を保持しつつ、振動サンプリングの制限を克服している。
• 機械学習ポテンシャルのアプローチとは異り、かつそれを補完するものであり、その精度はポテンシャルの学習誤差ではなく、基礎となる第一原理の力源によって制御されている。
• 大規模な埋め込みモデルの構築や複雑なカットオフ行列の管理を必要としないため、埋め込み法よりも展開が容易である。

本研究は限界についても認めており、現在の実装は並行モードおよび等周波数近似、ならびに調和近似に依存しています。したがって、放出の記述には成功しているものの、吸収スペクトルの記述には励起状態の表面およびヤーン・テラー活性を慎重に考慮する必要があり、最終状態のポテンシャルエネルギー面が支配的な役割を果たす系に対しては、完全なドゥシュンスキー（Duschinsky）モード混合および明示的なビブロン結合の取り扱いが必要であることを示唆しています。