[論文レビュー] On the Nature of Self-Consistency in Density Functional Theory

Nick Woods|arXiv (Cornell University)|Mar 5, 2018

Advanced Chemical Physics Studies参考文献 63被引用数 5

ひとこと要約

本稿は、CASTEPソフトウェア内での数学的および物理的アプローチを用いて、Kohn-Sham密度汎関数理論における自己無撞着性加速を調査する。Marks & Luke (2011) や Banerjee et al. (2015) の高度な混合手法を実装したが、CASTEPのデフォルト手法との間に顕著な改善は認められなかった。一方で、実空間およびフーリエ空間における物理的に根拠のある誘電率予め条件付けフレームワークを提案し、KS系の線形応答モデルを用いることで、収束の系統的向上が可能であることを示した。

ABSTRACT

A thesis providing a pedagogical introduction to the problem of achieving self-consistency in density functional theory. Contained is an introduction to the framework of Kohn-Sham density functional theory, leading then to the considerations required to solve the equations of Kohn-Sham density functional theory. Specifically, a focus is placed on where current self-consistent field methodology is inefficient and/or fails to converge to a solution. As such, this review spans sub-disciplines such as numerical analysis of linear and non-linear systems, linear response theory, and general electronic structure theory. Toward the end of the thesis, certain contemporary methods for achieving self-consistency from literature are outlined, and a novel, computationally efficient preconditioning strategy is proposed. This work is implemented in the CASTEP software.

研究の動機と目的

• Kohn-Sham密度汎関数理論における自己無撞着解の収束速度を向上させること。
• Marks & Luke (2011) および Banerjee et al. (2015) が提唱した最近の数学的混合手法が、CASTEPソフトウェア環境内での有効性を評価すること。
• 実空間またはフーリエ空間における不均一なKS誘電率応答に基づく物理的に根拠のある予め条件付けフレームワークの開発および検証すること。
• 系の線形応答に関する物理的洞察が、純粋に数学的混合戦略を越えて収束を系統的に改善できるかを評価すること。
• 平面波DFTコードに適用可能な一般化可能で実装に依存しない自己無撞着性加速フレームワークを提供すること。

提案手法

• Kohn-Sham DFT問題を物理的構造を参照せずに、非線形系としてブラックボックスとして扱い、数学的混合技術を適用する。
• CASTEP内に Marks & Luke (2011) および Banerjee et al. (2015) の混合スキームを導入し、密度およびポテンシャルの混合に適用する。
• 自己無撞着サイクルの予め条件付けを実行するため、実空間またはフーリエ空間における不均一なKohn-Sham誘電率関数をモデル化するフレームワークを開発する。
• 線形応答理論を用いて誘電率モデルを情報提供し、密度／ポテンシャル更新に対する物理的に根拠のある予め条件付けを可能にする。
• 15種類以上の代表的物質を用いたテストスイートを用いて、さまざまな混合および予め条件付け戦略における収束性能をベンチマークする。
• 残差ノルムと反復回数を用いて収束行動を比較し、数学的および物理的アプローチの両方を評価する。

実験結果

リサーチクエスチョン

• RQ1 最近提唱された数学的混合手法（Marks & Luke (2011) および Banerjee et al. (2015)）は、15種類以上の物質を対象とした範囲で、CASTEPのデフォルト混合手法よりも収束を顕著に改善するか？
• RQ2 線形応答理論に基づく物理的に根拠のある誘電率予め条件付けフレームワークは、Kohn-Sham DFTにおける収束を向上させることができるか？
• RQ3 表現の選択（実空間対フーリエ空間）が、誘電率予め条件付けフレームワークの性能にどのように影響するか？
• RQ4 系の応答に関する物理的洞察は、純粋にアルゴリズム的な混合手法を越えて自己無撞着性をどの程度向上させられるか？
• RQ5 DFTにおける反復的混合と物理的モデリングを組み合わせた場合、収束に系統的かつ一般化可能な改善が得られるか？

主な発見

• Marks & Luke (2011) および Banerjee et al. (2015) の数学的混合手法は、15種類以上の物質を対象としたテストスイートにおいて、CASTEPのデフォルト混合手順と比較して顕著な系統的改善を示さなかった。
• 線形応答理論に基づく誘電率予め条件付けフレームワークの実装は成功しており、柔軟性を備えており、実空間およびフーリエ空間の両方でモデル化が可能であった。
• 予め条件付けフレームワーク内での単純な誘電率モデルは、収束向上の測定可能な可能性を示し、物理的洞察が反復スキームを強化できることを示唆した。
• 本研究は、純粋に数学的混合手法では限界があるものの、誘電率応答のような物理的原則を組み込むことで、より効果的な予め条件付けの道筋が開かれる可能性があることを確認した。
• フレームワークは一般化可能で実装に依存しないため、CASTEPを超えて平面波DFTコードにおける自己無撞着性加速に広く適用可能であると示唆された。
• 新しい手法の反復回数や残差ノルムに関する定量的ベンチマークは提供されていないが、定性的な傾向から、今後の最適化の可能性が示唆された。