X

Search Results

Searching....

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

9.1.3 Eigenvector-Following (EF) Algorithm

(July 14, 2022)

The work of Cerjan and Miller, 188 Cerjan C. J., Miller W. H.
J. Chem. Phys.
(1981), 75, pp. 2800.
Link
and later Simons and co-workers, 1098 Simons J. et al.
J. Phys. Chem.
(1983), 87, pp. 2745.
Link
, 61 Banerjee A. et al.
J. Phys. Chem.
(1985), 89, pp. 52.
Link
showed that there was a better step than simply directly following one of the Hessian eigenvectors. A simple modification to the Newton-Raphson step is capable of guiding the search away from the current region towards a stationary point with the required characteristics. This is

𝐡=-i(Fibi-λ)𝐮i (9.6)

in which λ can be regarded as a shift parameter on the Hessian eigenvalue bi. Scaling the Newton-Raphson step in this manner effectively directs the step to lie primarily, but not exclusively (unlike Poppinger’s algorithm 976 Poppinger D.
Chem. Phys. Lett.
(1975), 35, pp. 550.
Link
), along one of the local eigenmodes, depending on the value chosen for λ. References  188 Cerjan C. J., Miller W. H.
J. Chem. Phys.
(1981), 75, pp. 2800.
Link
, 1098 Simons J. et al.
J. Phys. Chem.
(1983), 87, pp. 2745.
Link
, 61 Banerjee A. et al.
J. Phys. Chem.
(1985), 89, pp. 52.
Link
all use the same basic approach of Eq. (9.6) but differ in the means of determining the value of λ.

The EF algorithm 54 Baker J.
J. Comput. Chem.
(1986), 7, pp. 385.
Link
uses the rational function approach presented in Refs.  61 Banerjee A. et al.
J. Phys. Chem.
(1985), 89, pp. 52.
Link
, yielding an eigenvalue equation of the form

(𝐇𝐠𝐠t0)(𝐡1)=λ(𝐡1) (9.7)

from which a suitable λ can be obtained. Expanding Eq. (9.7) yields

(𝐇-λ)𝐡+𝐠=0 (9.8)

and

𝐠t𝐡=λ (9.9)

In terms of a diagonal Hessian representation, Eq. (9.8) rearranges to Eq. (9.6), and substitution of Eq. (9.6) into the diagonal form of Eq. (9.9) gives

λ=-i(-Fi2bi-λ) (9.10)

which can be used to evaluate λ iteratively.

The eigenvalues, λ, of the RFO equation Eq. (9.7) have the following important properties: 61 Banerjee A. et al.
J. Phys. Chem.
(1985), 89, pp. 52.
Link

  • The (n+1) values of λ bracket the n eigenvalues of the Hessian matrix λi<bi<λi+1.

  • At a stationary point, one of the eigenvalues, λ, of Eq. (9.7) is zero and the other n eigenvalues are those of the Hessian at the stationary point.

  • For a saddle point of order m, the zero eigenvalue separates the m negative and the (n-m) positive Hessian eigenvalues.

This last property, the separability of the positive and negative Hessian eigenvalues, enables two shift parameters to be used, one for modes along which the energy is to be maximized and the other for which it is minimized. For a transition state (a first-order saddle point), in terms of the Hessian eigenmodes, we have the two matrix equations

(b1F1F10)(h11)=λp(h11) (9.11)
(b2F2𝟎𝟎bnFnF2Fn0)(h2hn1)=λn(h2hn1) (9.12)

where it is assumed that we are maximizing along the lowest Hessian mode 𝐮1. Note that λp is the highest eigenvalue of Eq. (9.11), which is always positive and approaches zero at convergence, and λn is the lowest eigenvalue of Eq. (9.12), which it is always negative and again approaches zero at convergence.

Choosing these values of λ gives a step that attempts to maximize along the lowest Hessian mode, while at the same time minimizing along all the other modes. It does this regardless of the Hessian eigenvalue structure (unlike the Newton-Raphson step). The two shift parameters are then used in Eq. (9.6) to give the final step

𝐡=-(F1b1-λp)𝐮1+i=2n(Fibi-λn)𝐮i (9.13)

If this step is greater than the maximum allowed, it is scaled down. For minimization only one shift parameter, λn, is used which acts on all modes.

In Eq. (9.11)) and Eq. (9.12) it was assumed that the step would maximize along the lowest Hessian mode, b1, and minimize along all the higher modes. However, it is possible to maximize along modes other than the lowest, and in this way potentially locate transition states for alternative rearrangements/dissociations from the same initial starting point. For maximization along the kth mode (instead of the lowest mode), Eq. (9.11) is replaced by

(bkFkFk0)(hk1)=λp(hk1) (9.14)

and Eq. (9.12) now excludes the kth mode, but includes the lowest mode. Since what was originally the kth mode is the mode along which the negative eigenvalue is required, then this mode will eventually become the lowest mode at some stage of the optimization. To ensure that the original mode is being followed smoothly from one cycle to the next, the mode that is actually followed is the one with the greatest overlap with the mode followed on the previous cycle. This procedure is known as mode following. For more details and some examples, see Ref.  54 Baker J.
J. Comput. Chem.
(1986), 7, pp. 385.
Link
.