Model of diatomic homonuclear molecule scattering by atom or barriers

Abstract

The mathematical model of quantum tunnelling of diatomic homonuclear molecules through repulsive barriers or scattering by an atom is formulated in the s-wave approximation. The 2D boundary-value problem (BVP) in polar coordinates is reduced to a 1D BVP for a set of second-order ODEs by means of Kantorovich expansion over the set of parametric basis functions. The algorithm for calculating the asymptotic form of the parametric basis functions and effective potentials of the ODEs at large values of the parameter (hyperradial variable) is presented. The solution is sought by matching the numerical solution in one of the subintervals with the analytical solution in the adjacent one. The efficiency of the algorithm is confirmed by comparing the calculated solutions with those of the parametric eigenvalue problem obtained by applying the finite element method in the entire domain of definition at large values of the parameter. The applicability of algorithms and software are demonstrated by the example of benchmark calculations of discrete energy spectrum of the trimer Be3 in collinear configuration.