Calculating energy levels of isomerizing tetra-atomic molecules. II. The vibrational states of acetylene and vinylidene

A general, full-dimensional computational method for the accurate calculation of rotationally and vibrationally excited states of tetra-atomic molecules is further developed. The resulting computer program may be run in serial and parallel modes and is particularly appropriate for molecules executing wide-amplitude motions and isomerizations. An application to the isomerizing acetylene∕vinylidene system is presented. Large-scale calculations using a coordinate system based on orthogonal satellite vectors have been performed in six dimensions and vibrational term values and wave functions for acetylene and vinylidene states up to ≈23000cm−1 above the potential minimum have been determined. This has permitted the characterization of acetylene and vinylidene states at and above the isomerization barrier. These calculations employ more extensive vibrational basis sets and hence consider a much higher density of states than in any variational calculations reported hitherto for this system. Comparison of the calculated density of states with that determined empirically suggests that our calculations are the most realistic achieved for this system to date. Indeed more states have been converged than in any previous study of this system. Calculations on lower lying excited states of acetylene based on HC–CH diatom-diatom coordinates give nearly identical results to those based on orthogonal satellite vectors. Comparisons are also made with calculations based on HH–CC diatom-diatom coordinates