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

вопросы строения молекул и квантовой химии
Ответить
Vlas
Сообщения: 13
Зарегистрирован: Пн мар 29, 2010 9:15 pm

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

Сообщение Vlas » Пн мар 29, 2010 9:28 pm

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

VTur
Сообщения: 7360
Зарегистрирован: Пт авг 31, 2007 1:36 pm

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

Сообщение VTur » Пн мар 29, 2010 10:43 pm

Перезапустите еще раз расчет с *.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 и наслаждайтесь его анализом
Последний раз редактировалось VTur Ср мар 31, 2010 11:51 pm, всего редактировалось 1 раз.
После отстоя требуйте долива

Vlas
Сообщения: 13
Зарегистрирован: Пн мар 29, 2010 9:15 pm

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

Сообщение Vlas » Пн мар 29, 2010 10:47 pm

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

Vlas
Сообщения: 13
Зарегистрирован: Пн мар 29, 2010 9:15 pm

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

Сообщение Vlas » Вт мар 30, 2010 10:20 am

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

VTur
Сообщения: 7360
Зарегистрирован: Пт авг 31, 2007 1:36 pm

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

Сообщение VTur » Вт мар 30, 2010 3:18 pm

В wfn декартовы волновые функции, где вместо сферических гармоник вводятся множители вида xkypzm. Они и определяют угловую зависимость.
После отстоя требуйте долива

Vlas
Сообщения: 13
Зарегистрирован: Пн мар 29, 2010 9:15 pm

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

Сообщение Vlas » Ср мар 31, 2010 12:00 am

Я поставил 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
Сообщения: 7360
Зарегистрирован: Пт авг 31, 2007 1:36 pm

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

Сообщение VTur » Ср мар 31, 2010 11:50 pm

Я в голове никакие формулы не держу. Только общий смысл. Касаемо выражения угловой зависимости xkypzm. Сумма 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). Уж не обессудьте.
После отстоя требуйте долива

Ответить

Вернуться в «квантовая химия и моделирование»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 2 гостя