Явная формула для гауссовских атомных орбиталей в def2-svp базисе

[Kuzma]

Добрый день!

Пытаюсь руками посчитать значение орбитальных функций из базиса def2-svp исходя из предположения, что функции считаются как

Код: Выделить всё

\xi(\bar (r-Ra)) = N(m, c1, c2, c3, ...) Ylm(\phi, \theta)*( c1*exp(-a1*r^2) + c2*exp(-a2*r^2) + c3*exp(-a3*r^2) + ... )

Где N - нормализирующая константа.

На выходе получаю совпадающие с пакетом pyscf числа для s-орбиталей, но несовпадающие для p и d, причем несовпадение похоже экспоненциальное по |r|.
В некоторых статьях нашел, что может использоваться

Код: Выделить всё

\xi(\bar (r-Ra)) = N(m, c1, c2, c3) F(r^2) * Ylm(\phi, \theta)*( c1*exp(-a1*r^2) + c2*exp(-a2*r^2) + c3*exp(-a3*r^2) )

где F может быть как константой, так и из большого класса функций.
Разбирался ли кто-то с точным представлением функций в таком базисе?

[Kuzma]

Кажется, сам ответил на свой вопрос и F_m(r) = |r|^m