Как извлечь волновые функции из результатов расчета в GAMESS

вопросы строения молекул и квантовой химии
Osmiy
Сообщения: 28
Зарегистрирован: Сб май 31, 2008 9:04 pm

Re: Как извлечь волновые функции из результатов расчета в GAMESS

Сообщение Osmiy » Сб июл 19, 2008 6:53 pm

Уважаемые коллеги, не окажете ли еще чуть-чуть помощи.
С нормами примитивных функций FS,FP,... я разобрался, это
(2a/pi)^3/4 для s
(128a^5/pi^3)^1/4 для p
и т.д.
Нормы базисных функций Norm1S,Norm2S,... я не нашел, вычислить самому трудно- нужно отдельно для каждого типа оболочки и базиса.
Может кто покажет готовые выражения или хотя бы покажите как заставить GAMESS их показать.

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

Re: Как извлечь волновые функции из результатов расчета в GAMESS

Сообщение VTur » Ср июл 23, 2008 3:54 pm

Создайте *.wfn файл, по идее, в нём все, что для Вас нужно, есть.
После отстоя требуйте долива

Osmiy
Сообщения: 28
Зарегистрирован: Сб май 31, 2008 9:04 pm

Re: Как извлечь волновые функции из результатов расчета в GAMESS

Сообщение Osmiy » Ср июл 23, 2008 8:10 pm

В Docs конкретно про wfn-файлы ничего нет. Может вы имели Punch?
Если да, то там нет норм базисных функций, ну или я не знаю где они там находятся.

N_A_B
Сообщения: 133
Зарегистрирован: Вс ноя 04, 2007 10:43 am

Re: Как извлечь волновые функции из результатов расчета в GAMESS

Сообщение N_A_B » Чт июл 24, 2008 1:55 pm

Как я понимаю в данном случае - wfn-файл
это входной файл для AIM-анализа
содержит в себе информацию о волновых функциях.
Его можно создать, введя комманту aimpac=.true. ( могу ошибаться в написании, пишу на память)
Смотрите input.doc

После рассчета- смотрите pun-файл.
В нем будет строка begin бла-бла-бла
и до самого конца( там будет end бла-блабла)
Вот то что между двумя этими строками- и есть содержимое wfn-файла.
Буду дома- смогу написать нормально и точно.
жить будем.

Osmiy
Сообщения: 28
Зарегистрирован: Сб май 31, 2008 9:04 pm

Re: Как извлечь волновые функции из результатов расчета в GAMESS

Сообщение Osmiy » Чт июл 24, 2008 2:55 pm

N_A_B, да-да-да, так и есть.
Вобщем, я разобрался.
Я долго не мог получить одинаковое значениет ЭП, посчитанное мною и GAMESS.
Думал, что нужно считать нормы базисных функций. В итоге, они оказались равными 1.
Это я, блин, не те единицы измерения координат использовал!!! :oops:
Наконец-то, я теперь могу занятся любимиым делом- программированием.

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

Re: Как извлечь волновые функции из результатов расчета в GAMESS

Сообщение VTur » Пт июл 25, 2008 3:42 am

Osmiy писал(а):Электронную плотность буду считать сам. Сначала хочу написать программу визуализации ЭП, за одно можно правильность расчетов визуально проверить. Потом, возможно, буду писать программу по поиску критических точек и их параметров.
Если Ваша прога - визуализатор, то все-таки лучше создавать *.cub файл. В нем можно вывести массив (сетку, объем) точек в декартовых координатах, для каждой из которых будет посчитано (на выбор): значение электронной плотности, градиента эл. пл., лапласиан эл. пл., спиновая плотность, разность, сума и любые комбинации эл. плотностей. С таким массивом удобнее работать визуализатору. А также вся орбитальная муть. См. рукосуйство к Гауссиану.
После отстоя требуйте долива

Ответить

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

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

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