ADIABATIC PSEUDOSPECTRAL METHODS FOR MULTIDIMENSIONAL VIBRATIONAL POTENTIALS

We describe a new algorithm for computing eigenvalues, spectral intensities, and selected eigenvectors of multidimensional vibrational potential surfaces. The method involves a synthesis of pseudospectral and sequential adiabatic reduction methods and merges the storage and computational advantages of the former with the improved basis set generated by the latter. The recursive residue generation method, which utilizes a Lanczos-based diagonalization procedure, is employed to calculate the observables. As a test case, we apply the method to computation of the infrared and stimulated emission pumping spectra for the HCN molecule and demonstrate a very large (one to three orders of magnitude) reduction in computational effort (for comparable accuracy) as compared to discrete variable representation (DVR)/adiabatic reduction or standard collocation approaches. We expect that this advantage will be increased considerably for larger (e.g., tetra-atomic) systems and will permit accurate basis set calculations on such systems to be carried out in a straightforward fashion.