Алгоритмы решения краевых задач для атомных тримеров в коллинеарной конфигурации методом Канторовича

Модель атомных тримеров с парными молекулярными взаимодействиями в коллинеарной конфигурации формулируется в виде двумерной краевой задачи в якобиевских полярных координатах. Последняя сводится методом Канторовича к одномерной краевой задаче для системы ОДУ второго порядка, используя разложение искомого решения по угловым базисным функциям, зависящим от гиперрадиуса, как от параметра. Представлены алгоритмы решения параметрической краевой задачи методом конечных элементов и вычисления асимптотических разложений параметрических угловых базисных функций и эффективных потенциалов системы ОДУ при больших значениях параметра. Эффективность алгоритмов подтверждается сравнением асимптотических решений параметрической задачи на собственные значения и эффективных потенциалов с их численными значениями, полученных методом конечных элементов при больших значениях параметра. Применимость алгоритмов демонстрируются на примере расчетов асимптотических разложений решений параметрической краевой задачи и эффективных потенциалов, и собственных значений энергии связи тримера бериллия в коллинеарной конфигурации.

Algorithms for solving the boundary-value problems for atomic trimers in collinear configuration using the Kantorovich method

The model of atomic trimers with molecular pair interactions for collinear configuration is formulated as a 2D boundary-value problem (BVP) in the Jacobi and polar coordinates. The latter is reduced to a 1D BVP for a system of second-order ordinary differential equations (ODEs) by means of the Kantorovich method using the expansion of the desired solutions over a set of angular basis functions, parametrically dependent on the (hyper)radial variable. The algorithms for solving the 1D parametric BVP by means of the finite element method (FEM) and calculating the asymptotes of the parametric angular functions and effective potentials of the system of ODEs at large values of the parameter are presented. The efficiency of the algorithms is confirmed by comparing the calculated asymptotic solutions and effective potentials with those of the parametric eigenvalue problem obtained by applying the FEM at large values of the parameter. The applicability of the algorithms is demonstrated by calculating the asymptotic expansions of the parametric BVP solution, effective potentials and sets of binding energies for the beryllium trimer in the collinear configuration.