A Jacobi-Wilson description coupled to a block-Davidson algorithm: An efficient scheme to calculate highly excited vibrational levels

We present a new approach based on the block-Davidson scheme which provides eigenvalues and eigenvectors of highly excited (ro) vibrational states of polyatomic molecules. The key ingredient is a prediagonalized-perturbative scheme applied to a subspace of a curvilinear normal-mode basis set. This approach is coupled to the Jacobi vector description recently developed by our group [C. Leforestier, A. Viel, F. Gatti, C. Munoz, and C. Iung, J. Chem. Phys. 114, 2099 (2001)], and applied to the HFCO and H2CO molecules, which represent the main difficulties of such calculations for any available method. The first one presents a significant state density because of its low symmetry and the presence of a fluorine atom, while strong resonances and intermode couplings occur in H2CO. This study establishes the robustness, the numerical efficiency, and the versatility of the method which is compared to the regular Lanczos and Davidson schemes. It is also shown that the eigenvectors can be obtained within a given accuracy easily set by the user. This point constitutes one of the main advantages of the method as very few potential-energy surfaces achieve an accuracy of the order of a wave number for highly excited states. Furthermore, this method allows one to restrict the calculations to selected energy levels based on their zero-order descriptions. (C) 2005 American Institute of Physics.