[论文解读] A Parallel Orbital-Updating Approach for Electronic Structure Calculations

Xiaoying Dai, Xingao Gong+2 more|arXiv (Cornell University)|May 1, 2014

Advanced Chemical Physics Studies参考文献 26被引用 8

一句话总结

本文提出了一种基于有限元离散化和Kohn-Sham方程迭代求解的并行轨道更新算法。通过将解分解为独立求解的线性系统（并行求解）和一个小规模特征值问题，该方法在分子体系的全势能计算中实现了高效率，数值实验验证了其有效性。

ABSTRACT

In this paper, we propose an orbital iteration based parallel approach for electronic structure calculations. This approach is based on our understanding of the single-particle equations of independent particles that move in an effective potential. With this new approach, the solution of the single-particle equation is reduced to some solutions of independent linear algebraic systems and a small scale algebraic problem. It is demonstrated by our numerical experiments that this new approach is quite efficient for full-potential calculations for a class of molecular systems.

研究动机与目标

• 开发一种可扩展的并行算法，用于求解电子结构计算中的Kohn-Sham方程。
• 通过将轨道更新分解为独立的线性求解和一个小特征子空间校正，降低非线性特征值问题的计算复杂度。
• 利用自适应有限元离散化和并行计算，实现高效的全势能计算。
• 提供一种基于变分原理和Galerkin投影的数值稳定且收敛的迭代格式。

提出的方法

• 该方法采用迭代轨道更新格式，每个轨道通过并行求解一个线性源问题独立更新。
• 在每次迭代中，基于当前电子密度更新有效势能，并通过求解 $ a(U^{(n)}; u^{(n+1/2)}_i, v) = \lambda^{(n)}_i (u^{(n)}_i, v) \quad \forall v \in V_{n+1} $ 计算新轨道。
• 在更新的试函数轨道张成的空间中求解一个小规模特征值问题，以将结果投影回正确的特征子空间：$ a(U^{(n+1/2)}; u^{(n+1)}, v) = \lambda^{(n+1)} (u^{(n+1)}, v) \quad \forall v \in \tilde{V}_{n+1} $。
• 在每一步中基于当前解应用自适应网格加密，以处理奇点和振荡行为。
• Galerkin投影 $ P_h $ 通过将解投影到有限元空间中，确保稳定性和收敛性。
• 通过适当选择初始猜测，该算法可兼容全势能和赝势设置。

实验结果

研究问题

• RQ1 并行轨道更新方法能否显著提升全势能电子结构计算的效率？
• RQ2 将轨道更新分解为独立线性求解和小特征子空间校正，对收敛性和精度有何影响？
• RQ3 该迭代格式的收敛行为如何受初始猜测质量和网格自适应性的影响？
• RQ4 在大规模分子体系中，该方法能否在降低计算成本的同时保持精度？
• RQ5 与传统自洽场（SCF）方法相比，该算法在可扩展性和鲁棒性方面表现如何？

主要发现

• 所提算法将Kohn-Sham特征值问题的求解转化为并行求解 $ N $ 个独立线性系统和一个小型特征值问题，实现了高并行效率。
• 数值实验表明，该方法在一类分子体系的全势能计算中具有高效性，显示出在超级计算机上的强可扩展性。
• 理论分析表明，计算特征子空间的误差受初始猜测误差、特征值偏差以及有限元空间中逼近误差的有界控制：$ d_1(u_i^{(1)}, V) \lesssim \sum_{k=1}^N \left( |\lambda_k - \lambda_k^{(0)}| + \inf_{v \in V_1} \|u_k - v\|_{1} + \|u_k - u_k^{(0)}\|_0 \right) $。
• 通过Galerkin投影和自适应网格加密，该方法保持了稳定性和收敛性，确保了最优逼近性质。
• 该方法对各种初始猜测（包括高斯型、Slater型及基于有限元的近似）均表现出鲁棒性。
• 对于具有强电子关联或奇异势能的体系，该方法特别有效，因为传统方法因病态性而难以处理。