Economical Quasi-Newton Self Consistent Field Solver
arXiv:2307.00560 · doi:10.1039/D3CP05557D
Abstract
We present an efficient quasi-Newton orbital solver optimized to reduce the number of gradient (Fock matrix) evaluations. The solver optimizes orthogonal orbitals by sequences of unitary rotations generated by the (preconditioned) limited-memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) algorithm incorporating trust-region step restriction. Low-rank structure of the inverse (approximate) Hessian is exploited not only in L-BFGS but also when solving the trust-region problem. The efficiency of the proposed ``Quasi-Newton Unitary Optimization with Trust-Region'' (QUOTR) method 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.
References in corpus (12)
- Recent developments in the PySCF program package
- Global rates of convergence for nonconvex optimization on manifolds
- An efficient implementation of two-component relativistic exact-decoupling methods for large molecules
- An assessment of initial guesses for self-consistent field calculations. Superposition of Atomic Potentials: simple yet efficient
- A trust-region augmented Hessian implementation for restricted and unrestricted Hartree-Fock and Kohn-Sham methods
- Orbital Optimization in Selected Configuration Interaction Methods
- Clustered Low-Rank Tensor Format: Introduction and Application to Fast Construction of Hartree-Fock Exchange
- Direct Energy Minimization Based on Exponential Transformation in Density Functional Calculations of Finite and Extended Systems
- Convergence analysis of adaptive DIIS algorithms with application to electronic ground state calculations
- Efficient Evaluation of Exact Exchange for Periodic Systems via Concentric Atomic Density Fitting
- A black-box, general purpose quadratic self-consistent field code with and without Cholesky Decomposition of the two-electron integrals
- Co-iterative augmented Hessian method for orbital optimization