[論文レビュー] Convergence of Adaptive Finite Element Approximations for Nonlinear Eigenvalue Problems

Hong Chen, Xuemin Gong+2 more|arXiv (Cornell University)|Jan 14, 2010

Advanced Numerical Methods in Computational Mathematics参考文献 31被引用数 5

ひとこと要約

本稿は、量子化学および材料科学でよく見られる非凸なエネルギー汎関数を伴う非線形固有値問題に対する適応型有限要素法の収束を確立する。後向き誤差推定と適応型メッシュ・リファインメント戦略を開発することで、著者らはこの手法が正確な基底状態解に収束することを証明した。微小構造計算における数値実験によってその有効性が検証された。

ABSTRACT

In this paper, we study an adaptive finite element method for a class of a nonlinear eigenvalue problems that may be of nonconvex energy functional and consider its applications to quantum chemistry. We prove the convergence of adaptive finite element approximations and present several numerical examples of micro-structure of matter calculations that support our theory.

研究の動機と目的

• 非凸なエネルギー汎関数を伴う非線形固有値問題に適応型有限要素法を適用した場合の収束理論の欠如に取り組む。
• 非局所項と変動するポテンシャルを含む非線形固有値問題に適した後向き誤差推定フレームワークを構築する。
• 特に原子核や化学結合付近の高い電子密度領域を的確に解像するように、局所的にメッシュをリファインする適応型有限要素アルゴリズムを設計する。
• 適切な正則性および安定性条件の下で、適応型アルゴリズムが正確な基底状態解に収束することを理論的に確立する。
• TFWモデルおよびGross-Pitaevskii方程式を用いた電子構造計算における理論的予測を、数値実験によって検証する。

提案手法

• 解の$L^2$ノルムに制約を課えた$ H^1_0(\bar{\Omega}) $における変分問題として非線形固有値問題を定式化する。
• 非構造化メッシュ上での$ C^0 $連続な線形有限要素を用いて、ガラーキン離散化を適用する。
• 離散系の残差と離散解の双対整合性に基づいて、後向き誤差推定器を導出する。
• エネルギー誤差と同等であることが示された、残差関数性$ F_h(λ_h, u_h) $の双対ノルムを含む残差に基づく誤差推定器を用いる。
• Dörfler型のマーク戦略を用いて収束を保証する、解→推定→マーク→リファインの適応型ループを実装する。
• 非局所項（例：畳み込み型非線形性）と非凸なエネルギー汎関数を扱うために、非線形関数解析と特別な技術を用いる。

実験結果

リサーチクエスチョン

• RQ1 非凸なエネルギー汎関数を伴う非線形固有値問題に対して、適応型有限要素法は正確な基底状態解に収束するか？
• RQ2 非局所的相互作用を含む非線形固有値問題に対して、後向き誤差推定器をどのように構築し、検証するか？
• RQ3 エネルギー汎関数の非凸性が存在する中で、適応型有限要素法が収束を保証する条件は何か？
• RQ4 適応型手法は、電子密度の急峻な勾配（原子核や化学結合付近）を効果的に解像できるか？
• RQ5 非線形固有値問題の文脈において、後向き誤差推定器と実際のエネルギー誤差の関係は何か？

主な発見

• 非凸なエネルギー汎関数を伴う非線形固有値問題に対し、適応型有限要素法は正確な基底状態解に収束することが示された。
• 双対ノルムにおける残差に基づく後向き誤差推定器がエネルギー誤差と同等であることが証明され、信頼性の高い適応型リファインが可能となった。
• 解が正則である、すなわち固有関数の直交補空間上でエネルギー汎関数の第二変分が強収縮性を満たすという仮定の下で収束が確立された。
• Thomas-Fermi-von WeizsäckerモデルおよびGross-Pitaevskii方程式に対する数値実験では、理論的予測と整合する収束率が得られ、図5.1、5.3、および5.5に示された。
• 局所的メッシュリファインによって原子核や化学結合付近の局在的電子密度を的確に捉えることができ、過剰な自由度を増やすことなく精度が向上した。
• 理論的枠組みは、ボーズ＝アインシュタイン凝縮、TFWモデル、負のカーネルを有するシュレーディンガー＝ニュートン系を含む広範な物理的モデルをカバーする。