[論文レビュー] Pawpyseed: Perturbation-extrapolation band shifting corrections for point defect calculations

Kyle Bystrom, Danny Broberg+3 more|arXiv (Cornell University)|Apr 25, 2019

Structural Analysis and Optimization被引用数 7

ひとこと要約

本稿では、DFTにおける点欠陥計算のための自動的・並列化された摂動外挿法バンドシフト補正を可能にする高性能なPythonパッケージpawpyseedを紹介する。一次摂動理論と波動関数重なり解析を活用することで、シリコン欠陥における遷移準位の予測を向上させ、ハイブリッド汎関数計算と同等の精度を達成しながらも、スクリーニングのための計算的実行可能性を維持する。

ABSTRACT

Significant progress has been made recently in the automation and standardization of ab initio point defect calculations. However, the task of developing, implementing, and benchmarking charge corrections for density functional theory (DFT) point defect calculations is still an open challenge. Here we present a high-performance Python package called pawpyseed, which can read PAW DFT wave functions and calculate the overlap between wavefunctions from different structures. Using pawpyseed, we implement a new band shifting correction derived from first order perturbation theory. We benchmark this method by calculating the transition levels of several point defects in silicon and comparing to experimental and hybrid functional results. The new band shifting method can shift single-particle energies to improve transition level predictions and can be automated and parallelized using pawpyseed, suggesting it could be a useful method for high-throughput point defect calculations.

研究の動機と目的

• 半局所的DFTにおける欠陥遷移準位の不正確さという長年の課題に取り組む。
• 高価なハイブリッド汎関数に依存せずに、欠陥のバンド端シフトを計算的に効率的かつ自動的に補正する手法を開発する。
• 既存のDFTワークフローと互換性のあるスケーラブルで並列化された補正フレームワークを提供することで、点欠陥の高スループットスクリーニングを可能にする。
• 欠陥状態のバンド端に対する相対的位置を補正することで、欠陥形成エネルギーおよび遷移準位の精度を向上させる。
• キーディフェクトにおける実験データおよびハイブリッド汎関数結果と比較して、新しいバンドシフト法のベンチマークを実施する。

提案手法

• 欠陥とバルクスーパセル状態間の波動関数重なりを用いた摂動理論に起因するバンドシフト補正を実装する。
• pawpyseedを用いてPAW DFT波動関数を読み込み、異なるスーパセル構成におけるKohn-Sham状態間の重なり積分を計算する。
• GGAとハイブリッド汎関数のバンド端の差分に基づいて補正を適用する：ΔE_VBM = E_VBM,hybrid − E_VBM,GGA および ΔE_CBM = E_CBM,hybrid − E_CBM,GGA。
• 状態nの密度が欠陥およびその最近接原子にどれほど局在しているかを示すχ_n（局在率）を用いて、補正の割合を調整する局在化重み付き補正を導入する。
• Pythonベースのpawpyseedフレームワーク内にCで最適化されたルーチンを統合し、計算的に高負荷な重なりおよび射影計算を並列化する。
• 欠陥状態の特定、バンド端抽出、遷移準位の再キャリブレーションを含む、補正パイプライン全体を自動化する。

実験結果

リサーチクエスチョン

• RQ1 摂動外挿法バンドシフト補正法は、点欠陥における半局所的DFTの遷移準位予測精度を向上させることができるか？
• RQ2 本手法のバンドシフト補正は、シリコン点欠陥において標準的なフェルミ準位再キャリブレーションとハイブリッド汎関数結果と比べてどのように異なるか？
• RQ3 χ_nによる局在化の考慮（局在化率）は、欠陥準位補正の精度をどの程度向上させるか？
• RQ4 この補正手法は、高スループット欠陥スクリーニングに効率的に並列化および自動化できるか？
• RQ5 本手法は、バナディティや電子的性質が異なる置換型欠陥など、多様な欠陥タイプにわたり、安定した性能を示すか？

主な発見

• 摂動外挿法バンドシフト補正は、シリコン欠陥における遷移準位予測の精度を顕著に向上させ、ハイブリッド汎関数および実験的ベンチマークとの誤差を低減した。
• シリコンバナディティおよび置換型欠陥（B, P, Cu, S）において、本手法はハイブリッド汎関数結果と約0.1–0.2 eVの誤差で遷移準位を予測し、高い精度を達成した。
• 2原子素単位格子を用いたPBE-GGAとHSE06ハイブリッド関数との比較で、バンドギャップ補正はΔE_VBM = −0.36 eVおよびΔE_CBM = 0.22 eVを示した。
• 局在化状態の射影によるシフト補正法（局在化率χ_nで補正を重みづけ）は、混合性を持つ状態の精度を向上させ、局在化した準位の過剰補正を低減した。
• pawpyseedパッケージにより、補正ワークフローの完全な自動化と並列化が可能となり、計算クラスタ全体で効率的なスケーリングが実現した。
• 本手法は、深レベル状態や混合局在化/原子的性質を示す状態を含む多様な欠陥タイプにわたり、標準的なフェルミ準位再キャリブレーションを上回る性能を示した。