[论文解读] Economical Quasi-Newton Self Consistent Field Solver

Samuel A. Slattery, Kshitijkumar A. Surjuse+1 more|PubMed|Jul 2, 2023

Matrix Theory and Algorithms参考文献 6被引用 6

一句话总结

该论文提出QUOTR，一种新型拟牛顿轨道优化器，通过使用带信赖域约束的预处理有限内存BFGS算法，并利用海塞矩阵的低秩结构，高效求解信赖域子问题，从而最小化梯度评估次数。与标准Roothaan-Hall/DIIS求解器相比，该方法在收敛速度和鲁棒性方面表现更优，尤其适用于小HOMO-LUMO能隙或具有开壳层特征的复杂体系。

ABSTRACT

We present an efficient quasi-Newton orbital solver optimized to reduce the number of gradient evaluations and other computational steps of comparable cost. The solver optimizes orthogonal orbitals by sequences of unitary rotations generated by the (preconditioned) limited-memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) algorithm equipped with trust-region step restriction. The low-rank structure of the L-BFGS inverse Hessian is exploited when solving the trust-region problem. The efficiency of the proposed "Quasi-Newton Unitary Optimization with Trust-Region" (QUOTR) solver is compared to that of the standard Roothaan-Hall approach accelerated by the Direct Inversion of Iterative Subspace (DIIS), and other exact and approximate Newton solvers for mean-field (Hartree-Fock and Kohn-Sham) problems.

研究动机与目标

• 开发一种更稳健且高效的平均场电子结构方法轨道优化器，以减少对启发式收敛参数的依赖。
• 解决Roothaan-Hall/DIIS求解器的局限性，例如在小HOMO-LUMO能隙、开壳层构型或复杂电子结构体系中出现的收敛缓慢或行为异常问题。
• 在保持超线性收敛速度的前提下，最小化福克矩阵评估次数（NF），使其效率与标准DIIS加速求解器相当或更优。
• 通过采用带信赖域控制的拟牛顿优化方法，实现对非Aufbau解和大规模体系的可靠收敛。

提出的方法

• 采用预处理有限内存BFGS（L-BFGS）算法生成幺正轨道旋转。
• 对L-BFGS步长施加信赖域约束，以确保稳定性和全局收敛性。
• 利用L-BFGS逆海塞矩阵的低秩结构，高效求解信赖域子问题。
• 采用幺正优化框架，确保迭代更新过程中轨道的正交性。
• 集成信赖域步长限制，防止过大且不稳定的更新，从而提升鲁棒性。
• 通过梯度和海塞矩阵近似进行能量最小化，优化轨道参数，避免完整的对角化计算。

实验结果

研究问题

• RQ1 在具有挑战性的电子结构问题中，采用信赖域控制的拟牛顿方法是否能在收敛速度和鲁棒性方面优于标准Roothaan-Hall/DIIS求解器？
• RQ2 在不牺牲收敛速率或稳定性的情况下，轨道优化中福克矩阵评估次数可被减少到何种程度？
• RQ3 QUOTR求解器在HOMO-LUMO能隙趋近于零、开壳层构型以及d-和f-元素体系中的表现如何，这些体系中标准SCF求解器常会失效？
• RQ4 L-BFGS海塞矩阵的低秩结构是否可被有效利用，以加速轨道优化中信赖域子问题的求解？
• RQ5 所提出的方法在大规模体系和非LCAO表示下是否仍能保持高效性和鲁棒性？

主要发现

• QUOTR求解器在福克矩阵评估次数（NF）上优于标准Roothaan-Hall/DIIS求解器，尤其在具有挑战性的体系中表现更优。
• 对于HOMO-LUMO能隙趋近于零的体系，QUOTR能可靠收敛，而DIIS则可能失效或表现出不稳定行为。
• 该求解器在开壳层体系及过渡金属体系（d-和f-元素）中表现出稳健性能，包括标准SCF方法发生发散的情况。
• 信赖域机制可有效防止过大且不稳定的更新，显著提升收敛鲁棒性，且无需手动调节参数。
• 低秩海塞矩阵近似可高效求解信赖域子问题，降低每步迭代的计算开销。
• QUOTR实现了与牛顿法相当的超线性收敛速率，同时每步迭代成本低于基于精确海塞矩阵的求解器。