Smolyak representations with absorbing boundary conditions for reaction path Hamiltonian model of reactive scattering

Abstract

In this work, we present the efficient combination of Smolyak representations with time independent quantum mechanical approach using absorbing boundary conditions for the cumulative reaction probability calculations of a multidimensional reactive scattering problem. Our approach uses both kinds of Smolyak representations (finite basis and grid) which drastically reduces the size of the basis representation for the cumulative reaction operator. The cumulative reaction probability is thus obtained by solving the eigenvalue problem within the context of reaction path Hamiltonian using the compact Smolyak basis combined with an iterative Lanczos algorithm. Benchmark calculations are presented for reactive scattering models with a linear reaction coordinate and applied to a 25D model highlighting the efficiency of the present approach for multidimensional reactive processes.