Hybrid variation-perturbation method for calculating rovibrational energy levels of polyatomic molecules

Abstract

A procedure for calculation of rotation-vibration states of medium sized molecules is presented. It combines the advantages of variational calculations and perturbation theory. The vibrational problem is solved by diagonalizing a Hamiltonian matrix, which is partitioned into two sub-blocks. The first, smaller sub-block includes matrix elements with the largest contribution to the energy levels targeted in the calculations. The second, larger sub-block comprises those basis states which have little effect on these energy levels. Numerical perturbation theory, implemented as a Jacobi rotation, is used to compute the contributions from the matrix elements of the second sub-block. Only the first sub-block needs to be stored in memory and diagonalized. Calculations of the vibrational-rotational energy levels also employ a partitioning of the Hamiltonian matrix into sub-blocks, each of which corresponds either to a single vibrational state or a set of resonating vibrational states, with all associated rotational levels. Physically, this partitioning is efficient when the Coriolis coupling between different vibrational states is small. Numerical perturbation theory is used to include the cross-contributions from different vibrational states. Separate individual sub-blocks are then diagonalized, replacing the diagonalization of a large Hamiltonian matrix with a number of small matrix diagonalizations. Numerical examples show that the proposed hybrid variational-perturbation method greatly speeds up the variational procedure without significant loss of precision for both vibrational-rotational energy levels and transition intensities. The hybrid scheme can be used for accurate nuclear motion calculations on molecules with up to 15 atoms on currently available computers.

Molecular Physics (Handy Special Issue), in press

Cited by in corpus (8)

• The 2020 release of the ExoMol database: molecular line lists for exoplanet and other hot atmospheres
• ExoMol line lists XVIII. The high temperature spectrum of VO
• The ExoMol project: Software for computing large molecular line lists
• ExoMol line lists XXIV: A new hot line list for silicon monohydride, SiH
• ExoMol line lists XXXI: Spectroscopy of lowest eights electronic states of C${}_2$
• Ab initio effective rotational and rovibrational Hamiltonians for non-rigid systems via curvilinear second order vibrational Møller-Plesset perturbation theory
• Predicted Landé $g$-factors for open shell diatomic molecules
• A hybrid variational-perturbation calculation of the ro-vibrational spectrum of nitric acid