SpectrumSDT: A program for parallel calculation of coupled rotational-vibrational energies and lifetimes of bound states and scattering resonances in triatomic systems

Abstract

We present SpectrumSDT – a program for calculations of energies and lifetimes of bound rotational-vibrational states below and scattering resonances above the dissociation threshold on a global potential energy surface of a triatomic system, which may include stable molecules, weekly-bound van-der-Waals complexes, and unbound atom + diatom scattering systems. Large-amplitude vibrational motion is treated explicitly using hyper-spherical coordinates. Three options for the rotational-vibrational interaction are supported: uncoupled (symmetric top rotor), partially coupled (to include interaction between several nearest states only) and full-coupled (vibrating asymmetric-top rotor). In addition to energies and lifetimes, SpectrumSDT is able to integrate ro-vibrational wave functions over the user-defined regions of potential energy surface, which helps to classify these states. In this release of the code, SpectrumSDT is limited to ABA-type molecules with wave functions that do not extend into the regions near Eckart singularities. Program summaryProgram title: SpectrumSDTCPC Library link to program files:https://doi.org/10.17632/9gftxjs4yk.1Developer's repository link:https://github.com/IgorGayday/SpectrumSDTLicensing provisions: GNU General Public License 3 (GPL)Programming language: FortranNature of problem: Calculations of energies and lifetimes of bound rotational-vibrational states below and scattering resonances above the dissociation threshold on a global potential energy surface of a triatomic system, which may include stable molecules, weekly-bound van-der-Waals complexes, and unbound atom + diatom scattering systems.Solution method: A Hamiltonian matrix is built in APH coordinates using an optimal 2D basis set, adjusted for a given problem. A complex absorbing potential (CAP) is added to define a boundary condition for calculation of scattering resonances above the dissociation threshold. The eigenstates of the Hamiltonain matrix are found using a state-of-the-art iterative eigensolver.Restrictions: The present version is restricted to ABA-molecules and wave functions that do not extend into linear and equilateral triangle configurations.Additional comments including restriction and unusual features: Probabilities and lifetimes of the wave functions can be calculated in user-defined regions on the PES, which allows to analyze their localization properties (i.e. automatic isotopomer assignment) and obtain channel-specific lifetimes for scattering resonances.