[论文解读] 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方程的电子结构计算数值实验，验证理论结果。

提出的方法

• 将非线性特征值问题表述为在 $ H^1_0(\bar{\Omega}) $ 中的变分问题，并对解的 $ L^2 $-范数施加约束。
• 在非结构化网格上使用 $ C^0 $-连续的分片线性元进行有限元Galerkin离散化。
• 基于离散系统的残量和离散解的对偶一致性，推导后验误差估计量。
• 采用基于残量泛函 $ 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所示。
• 通过局部网格加密，该方法成功捕捉了原子核和化学键附近的局域电子密度，提升了精度，同时避免了自由度的过度增加。
• 理论框架支持广泛的物理模型，包括玻色-爱instein凝聚体、TFW模型以及具有负核的Schrödinger-Newton系统。