Волновые функции из Гауссиана

[Vlas]

Подскажите, пожалуйста, как можно из Гауссиана взять волновые функции (электронные) для молекулы? Вотт, например, из прикрепленного файла

[VTur]

Перезапустите еще раз расчет с *.chk файлом.
1. или укажите GFINPUT и IOP(6/7=3), то в *.out файле
- получите свой базис в "Basis set in the form of general basis"
- получите орбитальные коэффициенты в "Molecular Orbital Coefficients " всех атомных орбиталей
- матрицу плотности "DENSITY MATRIX"
2. или укажите OUTPUT=WFN (не забудте 6d 10f)
- получите *.wfn файл со всеми экспонентами и коэффициентами
3. либо превратите *.chk в *.fck и наслаждайтесь его анализом

[Vlas]

Спасибо большое!

[Vlas]

В wfn - радиальные волновые функции? А как учесть угловую зависимость?

[VTur]

В wfn декартовы волновые функции, где вместо сферических гармоник вводятся множители вида x^ky^pz^m. Они и определяют угловую зависимость.

[Vlas]

Я поставил GFINPUT. Мне интересно что означает следущее:
-------------------
AO basis set:
Atom N1 Shell 1 S 6 bf 1 - 1 0.000000000000 0.000000000000 1.043182588633
0.6293480000D+04 0.1969788147D-02
0.9490440000D+03 0.1496128592D-01
0.2187760000D+03 0.7350053084D-01
0.6369160000D+02 0.2489367658D+00
0.1882820000D+02 0.6024594331D+00
0.2720230000D+01 0.2562017589D+00
...
----------------------------------
Расшифровка:
Atom N1-(фтом азота) , Shell 1 S 6 bf (6 примитивов эс-типа для оболочки 1S ), а все остальное я не понимаю:" 1 - 1 0.000000000000 0.000000000000 1.043182588633". Подскажите, что это? И еще: чтобы получить угловую зависимость используя декартовы функции, нужно x,y,z выразить через r,theta,phi. То есть:
x=r sin(theta) cos (phi)
y=r sin(theta) sin(phi)
z=r cos(theta)?
Или по-другому как-то?

[VTur]

Я в голове никакие формулы не держу. Только общий смысл. Касаемо выражения угловой зависимости x^ky^pz^m. Сумма k+p+m=l и только. l - орбитальное квантовое число. Т.е. для s орбиталей l=0 -> k=p=m=0.
Для p орбиталей l=1, их три варианта k=1 и p=m=0 или k=0 p=1 m=0, или k=0 p=0 m=1. Меняете экспоненты, получаете разные наборы.
Для d орбиталей l=2, т.е. k+p+m=2 и перебираются все варианты: 3p 6d 10f 15g и т.д.

В приведенных цифрах где-то указаны коэффициенты группировки гауссианов.
Да, вспомнил, что с GFINPUT желательно писать оверлей IOP(6/7=3). Уж не обессудьте.