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.
Andreana M. Rosnik and William F. Polik, "VPT2+K spectroscopic constants and matrix elements of the transformed vibrational Hamiltonian of a polyatomic molecule with resonances using Van Vleck perturbation theory", Molecular Physics, Vol. 214, No. 2, pp. 261-300, DOI:10.1080/00268976.2013.808386 (2014).