как построить молекулярную орбиталь?

[alexph]

Здравствуйте!

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

нормировка: dx1x2=C*x1*x2*Exp[-a*r^2] - C=(2048*a^7/pi^3)^(1/4),
чтобы при интегрировании (dx1x2)^2 давало 1

gaussianHELP - keyword Gen -
по поводу преобразований декартовыских и чистых угловых функций Гауссиан предлагает посмотреть
H. B. Schlegel and M. J. Frisch, Int. J. Quant. Chem. 54, 83 (1995).

[Osmiy]

Спасибо, alexph.
Я уже решил работать с 5D без преобразования в 6D. У меня в программе просто была реализована функция подсчета интегралов перекрывания АО, в которых угловые части были в виде произведения степеней координат, а тут такой сюрприз с 5D-функциями.
Теперь у меня рекурсивно происходит разложение, типа таково:
S((xx-yy)*(2zz-xx-yy)*R1*R2)=S(xx*(2zz-xx-yy)*R1*R2) - S(yy*(2zz-xx-yy)*R1*R2)=
=2*S(xx*zz*R1*R2) - S(xx*xx*R1*R2) - S(xx*yy*R1*R2) - 2*S(yy*zz*R1*R2) + S(yy*xx*R1*R2) + S(yy*yy*R1*R2).
Оставшиеся интегралы уже подсчитываются аналитически.