Block diagonal method

Assuming that the diabatic space can be truncated to the same size as the adiabatic space, Eqs. (64) and (65) clearly define the relationship between the two representations, and methods have been developed to obtain the tians-formation matrices directly. These include the line integral method of Baer [53,54] and the block diagonalization method of Pacher et al. [179]. Failure of the truncation assumption, however, leads to possibly important nonremovable derivative couplings remaining in the diabatic basis [55,182]. 

The Newton-Raphson block diagonal method is a second order optimizer. It calculates both the first and second derivatives of potential energy with respect to Cartesian coordinates. These derivatives provide information about both the slope and curvature of the potential energy surface. Unlike a full Newton-Raph son method, the block diagonal algorithm calculates the second derivative matrix for one atom at a time, avoiding the second derivatives with respect to two atoms. 

Unconstrained optimization methods [W. II. Press, et. ah, Numerical Recipes The An of Scieniific Compulime.. Cambridge University Press, 1 9H6. Chapter 101 can use values of only the objective function, or of first derivatives of the objective function. second derivatives of the objective function, etc. llyperChem uses first derivative information and, in the Block Diagonal Newton-Raphson case, second derivatives for one atom at a time. TlyperChem does not use optimizers that compute the full set of second derivatives (th e Hessian ) because it is im practical to store the Hessian for mac-romoleciiles with thousands of atoms. A future release may make explicit-Hessian meth oils available for smaller molecules but at this release only methods that store the first derivative information, or the second derivatives of a single atom, are used. 

The method of vibrational analysis presented here can work for any polyatomic molecule. One knows the mass-weighted Hessian and then computes the non-zero eigenvalues which then provide the squares of the normal mode vibrational frequencies. Point group symmetry can be used to block diagonalize this Hessian and to label the vibrational modes according to symmetry. 

To complete this method of gradual diagonalization, we write the matrix Kh in a block-diagonal form ... 

The positional and thermal parameters, and a scale factor, were refined by block-diagonal least-squares methods. After six cycles with isotropic thermal parameters R was 0-136. Five further cycles with anisotropic parameters for the gold atom and isotropic parameters for the fluorines reduced R to 0-099, and t ee cycles with all atoms anisotropic completed the refinement, the final R being 0-090. Measured and calculated structure factors are listed in Table 1. 

---