[論文レビュー] Adaptive Finite Element Approximations for Kohn-Sham Models

Huajie Chen, Xiaoying Dai+3 more|arXiv (Cornell University)|Feb 27, 2013

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

ひとこと要約

本稿では、密度汎関数理論におけるKohn-Sham方程式を解くための適応型有限要素法を提示している。残差に基づく後験的誤差推定器とD"orflerのマーク戦略を用いることで、収束性と準最適な複雑度を保証する。この手法により、原子核や化学結合付近で電子密度の変動が急激な系に対しても、効率的かつ高精度な電子構造計算が可能になり、堅牢な適応性を有する。

ABSTRACT

The Kohn-Sham equation is a powerful, widely used approach for computation of ground state electronic energies and densities in chemistry, materials science, biology, and nanosciences. In this paper, we study the adaptive finite element approximations for the Kohn-Sham model. Based on the residual type a posteriori error estimators proposed in this paper, we introduce an adaptive finite element algorithm with a quite general marking strategy and prove the convergence of the adaptive finite element approximations. Using D{\" o}rfler's marking strategy, we then get the convergence rate and quasi-optimal complexity. We also carry out several typical numerical experiments that not only support our theory,but also show the robustness and efficiency of the adaptive finite element computations in electronic structure calculations.

研究の動機と目的

• 密度汎関数理論におけるKohn-Shamモデルのための適応型有限要素アルゴリズムを開発し、電子構造計算の精度と効率を向上させること。
• Kohn-Shamフレームワークに生じる非線形固有値問題に特化した後験的誤差推定器を確立すること。
• D"orflerのマーク戦略を用いて、適応型有限要素近似の収束性と準最適な複雑度を証明すること。
• 典型的な分子系を用いた数値実験を通じて、理論的知見を検証すること。

提案手法

• Kohn-Sham方程式の有限要素離散化における要素残差とジャンプ項に基づく残差型後験的誤差推定器を提案する。
• メッシュの細分化のための一般化されたマーク戦略を導入し、収束性と最適な複雑度を達成するためにD"orflerの戦略を用いる。
• 適合メッシュの細分化による適応型有限要素離散化を採用し、二分法による更新により形状の正規性と入れ子の有限要素空間を保証する。
• 変分式とエネルギーノルムの誤差推定を用いて、解の誤差と局所誤差指標を結びつける。
• 誤差推定器と振動項を含む後験的境界を導出し、誤差推定の信頼性と有効性を保証する。
• 反復的に離散問題を解き、誤差指標を計算し、細分化のための要素をマークし、メッシュを更新する適応アルゴリズムを設計する。

実験結果

リサーチクエスチョン

• RQ1 適応型有限要素法は、密度汎関数理論における非線形Kohn-Sham固有値問題に対して収束するか？
• RQ2 Kohn-Shamモデルのための適応型有限要素近似の収束速度と計算複雑度は何か？
• RQ3 局所的電子密度変動を持つ系において、残差に基づく後験的誤差推定器はメッシュ細分化を効果的に誘導できるか？
• RQ4 多様な分子系、特に多様な電子構造を有する系においても、適応型手法は耐障害性と効率性を維持できるか？

主な発見

• 誤差と誤差推定器を含む収縮性の性質により、エネルギーノルムにおける適応型有限要素法の収束性が証明された。
• D"orflerのマーク戦略のもとで、準最適な複雑度が達成された。これは、自由度の数が近似誤差に対して最適な割合で増加することを保証する。
• 数値実験により、さまざまな分子系における基底状態エネルギーと電子密度の計算において、適応型手法の耐障害性と効率性が確認された。
• 後験的誤差推定器は信頼性があり、有効性を有しており、誤差と振動項を含む境界により、正確な誤差の指標が得られた。
• エネルギーノルムにおける収束速度は線形であり、収縮係数は1未満に厳密に保たれ、安定的かつ単調な誤差低減が示された。
• 局所的なメッシュ細分化により、原子核や化学結合付近での電子密度の急激な変化が効果的に解像された。