Пытаюсь руками посчитать значение орбитальных функций из базиса 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) + ... )
На выходе получаю совпадающие с пакетом 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) )
Разбирался ли кто-то с точным представлением функций в таком базисе?