Hessian operator

The last algorithm converges much more rapidly than the steepest descent method (5.29). The main difficulty is that it is a rather complicated problem to calculate the inverse quasi-Hessian operator. 

Note that the first term on the left-hand side of the last equation can be treated as a Hessian operator applied to Aa ... 

A few comments on the layout of the book. Definitions or common phrases are marked in italic, these can be found in the index. Underline is used for emphasizing important points. Operators, vectors and matrices are denoted in bold, scalars in normal text. Although I have tried to keep the notation as consistent as possible, different branches in computational chemistry often use different symbols for the same quantity. In order to comply with common usage, I have elected sometimes to switch notation between chapters. The second derivative of the energy, for example, is called the force constant k in force field theory, the corresponding matrix is denoted F when discussing vibrations, and called the Hessian H for optimization purposes. 

More variables are retained in this type of NLP problem formulation, but you can take advantage of sparse matrix routines that factor the linear (and linearized) equations efficiently. Figure 15.5 illustrates the sparsity of the Hessian matrix used in the QP subproblem that is part of the execution of an optimization of a plant involving five unit operations. 

Lastly, the gold is well washed, dried, and mixed with its own weight of bisulphatn of petassa, or a small quantity of borax and nitre, and fused in a Hessian crucible exposed ta a very strong beat. By this operation, the last portions of chloride of silver are removed, and a button of perfectly pure gold is obtained. 

Each such null vector may be considered an invariant or symmetry of the thermodynamic system, because it corresponds to an operation (change of extensive variables Xt) that produces no response in any intensive state variable and thus leaves the thermodynamic state unaltered (Sidebar 7.2). As described in Sidebar 10.3, these invariants also play a role somewhat analogous to overall rotations and translations ( null eigenmodes of the Hessian matrix) in the theory of molecular vibrations. 

In this section we shall compute the energy gradient and the Hessian matrix corresponding to the energy expression (3 25). We introduce the variation of the Cl coefficients by operating on the MCSCF state I0> with the unitary... 

The order of the operators in (4 5) is not arbitrary, since they do not commute. The reverse order, however, leads to more complicated expressions for the Hessian matrix, and since the final result is independent of the order, we make the more simple choice given in (4 5). The energy corresponding to the varied state (4 5) will be a function of the parameters in die unitary operators, and we can calculate the first and second derivatives of this function... 

The formulas above give the gradient and the Hessian in terms of matrix elements of the excitation operators. They can be evaluated in terms of one-and two-electron integrals, and first and second order reduced density matrices, by inserting the Hamiltonian (3 24) into equations (4 9), (4 11), and (4 13)-(4 15). Note that transition density matrices and are needed for the evaluation of the Cl coupling matrix (4 15). 

We shall in this section derive the explicit expressions for the elements of the gradient vector and the Hessian matrix. The derivation is a good exercise in handling the algebra of the excitation operators fey and the reader is suggested to carry out the detailed calculations, where they have been left out in the present exposition. 

Let us finally take a closer look at the orbital Hessian matrix, H(00). The calculation now involves the evaluation of commutators between the Hamiltonian and products of excitation operators according to equation (4 14). In spite of the rather tiring algebra, the result takes a surprisingly simple form ... 

The MCSCF gradient expression was first given by Pulay (1977). The MCSCF Hessian and first anharmonicity expressions were derived by Pulay (1983) using a Fock-operator approach, and by Jprgensen and Simons (1983) and Simons and Jorgensen (1983) using a response function approach. 

The Cl gradient expression was derived and implemented by Krishnan et al. (1980) and Brooks et al. (1980). The generalization to MRCI is due to Osamura et al. (1981, 1982a,b). The Hessian expression was derived by Jorgensen and Simons(1983) and implemented by Fox et al. (1983). Recently, a more efficient implementation has been reported by Lee et al. (1986). MRCI derivative expressions up to fourth order have been derived by Simons et al. (1984). The introduction of the Handy-Schaefer technique (Handy and Schaefer, 1984) greatly improved the efficiency of Cl derivative calculations. The calculation of Cl derivatives within the Fock-operator formalism has recently been reviewed by Osamura et al. (1987). 

---