Kozin Igor N, Law Mark M, Tennyson Jonathan, Hutson Jeremy M
Department of Chemistry, University of Aberdeen, Meston Walk, Aberdeen AB24 3UE, United Kingdom.
J Chem Phys. 2005 Feb 8;122(6):064309. doi: 10.1063/1.1850471.
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 approximately 23 000 cm(-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.
一种用于精确计算四原子分子转动和振动激发态的通用、全维计算方法得到了进一步发展。所得的计算机程序可以串行和并行模式运行,特别适用于执行大幅度运动和异构化的分子。本文给出了该方法在乙炔/亚乙烯基异构化体系中的应用。使用基于正交卫星矢量的坐标系进行了六维的大规模计算,确定了乙炔和亚乙烯基态在势能极小值以上约23000cm⁻¹以下的振动项值和波函数。这使得能够对异构化势垒及以上的乙炔和亚乙烯基态进行表征。这些计算采用了更广泛的振动基组,因此考虑的态密度比迄今报道的该体系的任何变分计算都要高得多。将计算得到的态密度与通过实验确定的态密度进行比较表明,我们的计算是目前该体系最符合实际情况的。实际上,与该体系以往的任何研究相比,已经收敛了更多的态。基于HC-CH双原子-双原子坐标对乙炔低激发态的计算结果与基于正交卫星矢量的计算结果几乎相同。还与基于HH-CC双原子-双原子坐标的计算进行了比较。