4.5 Converging SCF Calculations

4.5.1 Introduction

As for any numerical optimization procedure, the rate of convergence of the SCF procedure is dependent on the initial guess and on the algorithm used to step towards the stationary point. Q-Chem features a number of SCF optimization algorithms which can be selected via the $rem variable SCF_ALGORITHM, including:

Methods that are based on extrapolation or interpolation:

  • The highly successful DIIS procedures. These are the default (except for restricted open-shell SCF calculations) and are available for all orbital types (see Section 4.5.3). The damping349 and level-shifting technique842, 383, 330 can also be invoked together with DIIS (R, U only).

  • ADIIS (the augmented DIIS algorithm developed by Hu and Yang,404 available for R and U only).

Methods that make use of orbital gradient:

  • Direct Minimization (DM), which has been re-implemented as simple steepest descent with line search, and is available for all orbital types. DM can be invoked after a few DIIS iterations.

  • Geometric Direct Minimization (GDM) which is an improved and highly robust version of DM and is the recommended fall-back when DIIS fails. Like DM, It can also be invoked after a few iterations with DIIS to improve the initial guess. GDM is the default algorithm for restricted open-shell SCF calculations and is available for all orbital types (see Section 4.5.7).

  • GDM_LS: It is essentially a preconditioned (using orbital energy differences as the preconditioner) L-BFGS algorithm with line search, available for R, U, RO and OS_RO.

  • GDM_QLS: GDM with quadratic line search, available for R and U.

  • Square Gradient Minimization (SGM) family: SGM is a GDM-inspired method for excited state orbital optimization. Currently, three variants of this approach are supported: SGM (for RO and and OS_RO), SGM_LS (for R, U, RO and OS_RO) and SGM_QLS (for R and U), which are based upon the corresponding GDM approaches. For further details, see Section 4.5.12.

Methods that require orbital Hessian:

  • NEWTON_CG/NEWTON_MINRES (solve 𝐇𝐝=-𝐠 for the update direction with CG/MINRES solvers).

  • SF_NEWTON_CG (the “saddle-free" version of NEWTON_CG).

The analytical orbital Hessian is available for R/U/RO/G/CR unless special density functionals (e.g., those containing VV10) are used, while the use of finite-difference Hessian is available for all orbital types by setting FD_MAT_VEC_PROD = TRUE.

In addition to these algorithms, there is also the maximum overlap method (MOM) which ensures that DIIS always occupies a continuous set of orbitals and does not oscillate between different occupancies. MOM can also be used to obtain higher-energy solutions of the SCF equations (see Section 7.5). The relaxed constraint algorithm (RCA), which guarantees that the energy goes down at every step, is also available via the old SCF code (set GEN_SCFMAN = FALSE). Nevertheless, the performance of the ADIIS algorithm should be similar to it.

Since the code in GEN_SCFMAN is highly modular, the availability of different SCF algorithms to different SCF (orbital) types is largely extended in general. For example, the old ROSCF implementation requires the use of the GWH guess and the GDM algorithm exclusively. Such a limitation has been eliminated in GEN_SCFMAN based RO calculations.