X

Search Results

Searching....

9.1 Equilibrium Geometries and Transition-State Structures with Q-Chem

9.1.7 GDIIS

(July 14, 2022)

Direct inversion in the iterative subspace (DIIS) was originally developed by Pulay for accelerating SCF convergence. 996 Pulay P.
J. Comput. Chem.
(1982), 3, pp. 556.
Link
Subsequently, Csaszar and Pulay used a similar scheme for geometry optimization, which they termed GDIIS. 245 Csaszar P., Pulay P.
J. Mol. Struct. (Theochem)
(1984), 114, pp. 31.
Link
The method is somewhat different from the usual quasi-Newton type approach and is included in Optimize as an alternative to the EF algorithm. Tests indicate that its performance is similar to EF, at least for small systems, however, there is rarely an advantage in using GDIIS in preference to EF.

In GDIIS geometries generated in previous optimization cycles, 𝐱i, are linearly combined to find the “best” geometry for the current cycle

𝐱n=i=1mci𝐱i (9.33)

where the problem is to find the best values for the coefficients ci.

If we express each geometry by its deviation from the sought-after final geometry, 𝐱f, i.e., 𝐱f=𝐱i+𝐞i, where 𝐞i is an error vector, then it is obvious that if the conditions

𝐫=ci𝐞i (9.34)

and

ci=1 (9.35)

are satisfied, then the relation

ci𝐱i=𝐱f (9.36)

also holds.

The true error vectors 𝐞i are, of course, unknown. However, in the case of a nearly quadratic energy function they can be approximated by

𝐞i=-𝐇-1𝐠i (9.37)

where 𝐠i is the gradient vector corresponding to the geometry 𝐱i and 𝐇 is an approximation to the Hessian matrix. Minimization of the norm of the residuum vector 𝐫, Eq. (9.34), together with the constraint equation, Eq. (9.35), leads to a system of m+l linear equations

(B11B1m1Bm1Bmm1110)(c1cm-λ)=(001) (9.38)

where Bij=𝐞i|𝐞j is the scalar product of the error vectors 𝐞i and 𝐞j, and λ is a Lagrange multiplier.

The coefficients ci determined from Eq. (9.38) are used to calculate an intermediate interpolated geometry

𝐱m+1=ci𝐱i (9.39)

and its corresponding interpolated gradient

𝐠m+1=ci𝐠i (9.40)

A new, independent geometry is generated from the interpolated geometry and gradient according to

𝐱m+1=𝐱m+1-𝐇-1𝐠m+1. (9.41)

Note:  Convergence is theoretically guaranteed regardless of the quality of the Hessian matrix, as long as it is positive definite. The original GDIIS algorithm used a static Hessian (i.e. the initial guess Hessian) which was often a simple identity matrix. Updating the Hessian at each cycle generally results in more rapid convergence, and this is the default in Optimize.

Other improvements to the original method include limiting the number of previous geometries used in Eq. (9.33) by neglecting earlier geometries, and eliminating any geometries more than a certain distance from the current geometry (default = 0.3 a.u.).