A Least-Squares Commutator in the Iterative Subspace Method for Accelerating Self-Consistent Field Convergence

An improved DIIS method using a versatile residual matrix to accelerate SCF starting from a crude guess

The minimization of the commutator of the Fock and density matrices as the error matrix in the direct inversion of the iterative subspace (CDIIS) developed by Pulay is a powerful self-consistent field (SCF) acceleration technique for the construction of optimum Fock matrix, if initiated with a fair initial guess. In this work, we present an alternative minimized error matrix to the commutator in the CDIIS, namely the residual or the gradient of the energy-functional for a Slater determinant subject to the orthonormality constraints among orbitals, representing the search for a newly improved Fock matrix in the direction of the residual in the direct inversion of the iterative subspace (RDIIS). Implemented in the computational chemistry package GAMESS, the RDIIS is compared with the standard CDIIS and the second order SCF orbital optimization (SOSCF) for tested molecules started with a crude guess. As a result, the RDIIS stably and efficiently performs the SCF convergence acceleration. Furthermore, the RDIIS is considerably independent on the subspace size with the concentrated linear coefficients accounting proportionally for the Fock matrices close to the current iteration.

Galván I, Lindh R. A Story of Three Levels of Sophistication in SCF/KS-DFT Orbital Optimization Procedures

In this work, three versions of self-consistent field/Kohn-Sham density functional theory (SCF/KS-DFT) orbital optimization are described and benchmarked. The methods are a modified version of the geometry version of the direct inversion in the iterative subspace approach (which we call r-GDIIS), the modified restricted step rational function optimization method (RS-RFO), and the novel subspace gradient-enhanced Kriging method combined with restricted variance optimization (S-GEK/RVO). The modifications introduced are aimed at improving the robustness and computational scaling of the procedures. In particular, the subspace approach in S-GEK/RVO allows the application to SCF/KS-DFT optimization of a machine learning technique that has proven to be successful in geometry optimizations. The performance of the three methods is benchmarked for a large number of small- to medium-sized organic molecules, at equilibrium structures and close to a transition state, and a second set of molecules containing closed- and open-shell transition metals. The results indicate the importance of the resetting technique in boosting the performance of the r-GDIIS procedure. Moreover, it is demonstrated that already at the inception of the subspace version of GEK to optimize SCF wave functions, it displays superior and robust convergence properties as compared to those of the standard state-of-the-art SCF/KS-DFT optimization methods.

Iterative subspace algorithms for finite-temperature solution of Dyson equation

One-particle Green’s functions obtained from the self-consistent solution of the Dyson equation can be employed in the evaluation of spectroscopic and thermodynamic properties for both molecules and solids. However, typical acceleration techniques used in the traditional quantum chemistry self-consistent algorithms cannot be easily deployed for the Green’s function methods because of a non-convex grand potential functional and a non-idempotent density matrix. Moreover, the optimization problem can become more challenging due to the inclusion of correlation effects, changing chemical potential, and fluctuations of the number of particles. In this paper, we study acceleration techniques to target the self-consistent solution of the Dyson equation directly. We use the direct inversion in the iterative subspace (DIIS), the least-squared commutator in the iterative subspace (LCIIS), and the Krylov space accelerated inexact Newton method (KAIN). We observe that the definition of the residual has a significant impact on the convergence of the iterative procedure. Based on the Dyson equation, we generalize the concept of the commutator residual used in DIIS and LCIIS and compare it with the difference residual used in DIIS and KAIN. The commutator residuals outperform the difference residuals for all considered molecular and solid systems within both GW and GF2. For a number of bond-breaking problems, we found that an easily obtained high-temperature solution with effectively suppressed correlations is a very effective starting point for reaching convergence of the problematic low-temperature solutions through a sequential reduction of temperature during calculations.