A method is introduced for the calculation of normal-mode vibrational frequencies of polyatomic molecules based on numerical differencing of analytical gradients in symmetry coordinates. This procedure requires a number of gradient evaluations equal to the largest number of symmetry coordinates belonging to any single irreducible representation of the molecular point group (plus a single gradient evaluation at the equilibrium configuration), which is fewer than the 3N-6 ( N atoms) gradient evaluations needed for schemes based on Cartesian or internal coordinates. While the proposed method will not generally be as efficient as procedures which involve the direct calculation of energy second derivatives analytically (as are now available for single-determinant wavefunctions) it appears to be equally accurate, and it should be the method of choice for frequency calculations involving multideterminant wavefunctions for which analytical second-derivative algorithms have yet to be developed. The method is illustrated by the calculation of equilibrium secondary deuterium-isotope effects on a number of reactions involving simple carbocations.
* National Science Foundation postdoctoral fellow.