Abstract
Vibrational levels of polyatomic molecules are analysed with Van Vleck perturbation theory to connect experimental energy levels to computed molecular potential energy surfaces. Vibrational matrix elements are calculated from a quartic potential function via second-order Van Vleck perturbation theory, a procedure that treats both weak and strong interactions among vibrational states by approximately block-diagonalising the vibrational Hamiltonian. A clear and complete derivation of anharmonic and resonance constants as well as general expressions for both on- and off-diagonal matrix elements of the transformed Hamiltonian is presented. The equations are written in partial fraction form and as a constant multiplied by a harmonic oscillator matrix element to facilitate removing the effect of strongly interacting resonant states both in analytical formulae and in computer code. The derived equations are validated numerically, and results for the isotopomers of formaldehyde (H2CO, HDCO, D2CO) are included. The implications of the equations on zero-point energy calculations and experimental fits are discussed. The VPT2+K method is defined by these results for use in fitting and calculating vibrational energy levels.
Acknowledgements
We gratefully acknowledge John L. Davisson for his computation of molecular force fields and Dr Robert W. Field and Dr Peter R. Taylor for many insightful discussions. This work was supported by Howard Hughes Medical Institute and by National Science Foundation grants CHE-0624602 and CHE-1039925.