Files in this item



application/pdf5277.pdf (23kB)
(no description provided)PDF


Title:Variational Rovibrational Calculation For Linear Tetraatomic Molecules: I. The C8v4 Approach
Author(s):Schröder, Benjamin
Subject(s):Theory and Computation
Abstract:Variational rovibrational calculations for small molecules with up to 3 atoms are nowadays a routine task.\footnote{\ Tennyson, \textit{J. Chem. Phys.} \textbf{145}, 120901 (2016).} By employing accurate potential energy surfaces it is possible to achieve so called “spectroscopic accuracy” of 1 cm$^{-1}$ for small linear molecules.\footnote{\ Schröder and Sebald, \textit{J. Chem. Phys.} \textbf{144}, 044307 (2016); Makhnev~\textit{et al.}, \textit{J. Phys. Chem. A} \textbf{122}, 1326 (2018).} However, rovibrational calculations for linear molecules with more than 3 atoms are still rather scarce.\footnote{\ Mladenovi\'{c}, \textit{J. Chem. Phys.} \textbf{141}, 224304 (2014); Chubb~\textit{et al.}, \textit{J. Chem. Phys.} \textbf{149}, 14101 (2018).} A new variational method (C8v4) will be presented that is able to obtain rovibrational term energies and wave functions of tetraatomic linear molecules based on Watsons isomorphic Hamiltonian for linear molecules.\footnote{\ Watson, \textit{Mol. Phys.} \textbf{14}, 465 (1970).} The rovibrational wavefunction is expanded in (symmetrized) products of harmonic oscillator (1D and 2D) and rigid-rotor functions. The intricacies related to the vibrational angular momentum in linear molecules require a careful study of symmetry properties. Kinetic energy matrix elements in the chosen basis set can be evaluated in a fast mixed numerical/analytical fashion. The main computational bottlenecks are the integration of the potential energy matrix and the diagonalisation of the rovibrational Hamiltonian. An efficient hybrid MPI/OMP parallelization has been implemented, exploiting the block structure of the Hamiltonian, that yields a significant reduction of the computational time. The diagonalization employs a $K$-contraction scheme that reduces the number of basis functions in the final rovibrational Hamilonian by at least an order of magnitude. Example benchmark calculations for acetylene (HCCH) are presented. The converged ($N_\text{vib} = 2785$) results for HCCH up to $J_\text{max} = 2$ accurately reproduce earlier variational calculations\footnote{\ Bramley and Handy, \textit{J. Chem. Phys.} \textbf{98}, 1378 (1993).} in both $g$ and $u$ symmetry.
Issue Date:2021-06-23
Publisher:International Symposium on Molecular Spectroscopy
Genre:Conference Paper / Presentation
Date Available in IDEALS:2021-09-24

This item appears in the following Collection(s)

Item Statistics