О вычислении потенциала в многоатомных системах

Предлагается численный метод нахождения потенциала многоатомной системы в прямом пространстве. Отличительная особенность метода состоит в разделении электронной плотности $\rho$ и потенциала $\varphi$ на две части: $\rho=\rho_0+\hat\rho$, $\varphi=\varphi_0+\hat\varphi$, где $\rho_0$ — сумма сферических атомных плотностей, а потенциал $\varphi_0$ порождается плотностью $\rho_0$. Потенциал $\hat\varphi$ находится путем решения уравнения Пуассона. Граничные условия получены путем разложения обратной величины расстояния между двумя точками в ряд по полиномам Лежандра. Для обеспечения точности метода расчетная область разбивается на многогранники Вороного и применяются асимптотические оценки итераций при замене характеристической функции гладкими приближениями. Для численного решения уравнения Пуассона использованы двухсеточный метод и Фурье-преобразование. Получена оценка точности метода $O(h^{\gamma-1})$, где $h$ — шаг сетки, $\gamma$ — фиксированное число, большее $1$. Погрешность метода проанализирована на модельной двухатомной задаче. Библ. 19. Фиг. 2. Табл. 1.