10min:

KRISTY L. MARDIS AND EDWIN L. SIBERT, Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin - Madison, Madison, WI 53706.

Ab initio methods are now capable of calculating quartic force fields which accurately reproduce experimental observables such as band heads, _o and inertial constants, B. We present a method of refining these ab initio force field constants, f ^ab. The accuracy of the f ^ab permits the approximation of a figure of merit function, \beginequation \chi²_y = \sum_i\left(y_i ^expt - y_i\right)²/w_i²,\hspace3mm where~ y = _o~ or~ B, \endequation as a Taylor Series truncated at second order, \beginequation \chi²_y ~ A_y ^ab +\Upsilon_y^T( f-f ^ab) + \frac12( f-f ^ab)^T F_ y( f-f ^ab). \endequation Often, the computational effort of minimization is dominated by the derivative calculation. However, calculating _o's and B's by Van Vleck perturbation theory facilitates the fast calculation of the Hessian. The merit function is then minimized analytically by diagonalization of this Hessian. This method of improving the force constants has several advantages. First, there is great latitude in choosing a merit function. We choose \beginequation \chi² = (1 + \alpha)\chi²__o + \alpha\chi²_B. \endequation However, additional terms, such as different isotopic species, can be easily added to the merit function. Second, the calculation of the Hessian provides insight into the contours of the potential energy surface (PES). Third, small eigenvalues, which correspond to long, flat valleys in the PES, can be neglected. Finally, the near quadratic nature of \chi² permits the minimization to occur in a single step. Specific applications of this method to CO₂ and H₂CO are presented.