Hessian в Gamess

вопросы строения молекул и квантовой химии
Ответить
Аватара пользователя
AlexKsen
Сообщения: 602
Зарегистрирован: Пн мар 07, 2011 8:43 am

Hessian в Gamess

Сообщение AlexKsen » Чт ноя 01, 2018 8:34 pm

Добрый вечер.
Помогите разобраться, почему не стали рассчитываться частоты :dontknow: .
hessian.txt
У вас нет необходимых прав для просмотра вложений в этом сообщении.

Аватара пользователя
Гесс
Сообщения: 13055
Зарегистрирован: Ср фев 15, 2012 11:19 pm

Re: Hessian в Gamess

Сообщение Гесс » Чт ноя 01, 2018 10:18 pm

Дык это, два экрана от конца аутпута -

строка 389389 ENERGY DID NOT CONVERGE...ABORTING HESSIAN
отматываем на окончание scf - так как аутпут забит орбиталями то это почти самое начало
строка 3748 SCF IS UNCONVERGED, TOO MANY ITERATIONS
проблему можно пытаться решить технически разнообразно, понижением критериев сходимости, судя скажем по энергии оно очень хорошо сходилось первые шагов 15 и достигло весьма неплохих цифр. Можно поменять конвержер, скажем я там в аутпуте вижу DIIS=T DEM=F SOSCF=F, возможно с SOSCF оно бы свелось. Ну и в большинстве кодов есть бронебойный Ньютон-Рафсон, который работает почти всегда.
Впрочем тут есть и другая проблема, сама система.
Во первых - вы хотите гессиан - а это оно типа так уже оптимизирована? Я так понимаю это HF/6-31? Спорный выбор на мой вкус, но главное что мне ненравится - заряд минус 5. Такие полианионы на мой взгляд 1) изрядно афизичны и 2) будут иметь проблемы со сходимостью (поэтому я удивлен что вам удалось его оптимизировать). На мой взгляд сульфогруппы надо забить протонами и тогда возможно проблемность снизится.

Аватара пользователя
AlexKsen
Сообщения: 602
Зарегистрирован: Пн мар 07, 2011 8:43 am

Re: Hessian в Gamess

Сообщение AlexKsen » Чт ноя 01, 2018 10:40 pm

Гесс, спасибо за ответ.
Геометрия системы оптимизирована. На первом этапе попробую SOSCF=.Т.

$BASIS GBASIS=N31 NGAUSS=6 NDFUNC=1 NPFUNC=1 DIFFSP=.TRUE. DIFFS=.TRUE. $END
$SCF DIRSCF=.T. DIIS=.T. DAMP=.T. SHIFT=.T. SOSCF=.T. $END
$CONTRL SCFTYP=RHF RUNTYP=HESSIAN DFTTYP=B3LYP
ICHARG=-5 MULT=1 MAXIT=100 $END
$FORCE METHOD=ANALYTIC VIBANL=.TRUE. $END
$SYSTEM TIMLIM=1200 MWORDS=300 $END

Это было 6-31++G(d,p), попробую b3lyp/6-31++G(d,p).

Аватара пользователя
Гесс
Сообщения: 13055
Зарегистрирован: Ср фев 15, 2012 11:19 pm

Re: Hessian в Gamess

Сообщение Гесс » Чт ноя 01, 2018 10:50 pm

В этом инпуте у вас DFTTYP=B3LYP, в прошлом такого небыло, оптимизировано оно у вас в чем?
Если оно как то оптимизировалось то странно что у него теперь возникают проблемы со сведением scf в том же(?) методе. Тогда можно было бы прочитать волновую функцию хотя я ужене помню где ее гамесс хранит.
Впрочем независимо от того сведется оно или нет - заряд минус 5 мне не нравится от слова совсем.

Аватара пользователя
AlexKsen
Сообщения: 602
Зарегистрирован: Пн мар 07, 2011 8:43 am

Re: Hessian в Gamess

Сообщение AlexKsen » Чт ноя 01, 2018 10:53 pm

Оптимизировано в b3lyp/6-31++G(d,p).
В первый раз проглядел.

Аватара пользователя
AlexKsen
Сообщения: 602
Зарегистрирован: Пн мар 07, 2011 8:43 am

Re: Hessian в Gamess

Сообщение AlexKsen » Пт ноя 02, 2018 9:48 am

SOSCF=.Т. - не помогло решить проблему.

Подскажите, что еще можно попробовать сделать?
У вас нет необходимых прав для просмотра вложений в этом сообщении.

Аватара пользователя
Гесс
Сообщения: 13055
Зарегистрирован: Ср фев 15, 2012 11:19 pm

Re: Hessian в Gamess

Сообщение Гесс » Пт ноя 02, 2018 10:52 am

:shock:
То есть вас не смутило что этот аутпут в 100 с лишним раз меньше чем предыдущий и у вас не возникло желания посмотреть "а че там отличного от предыдущего"?
Я вам могу сказать - он такой маленький потому что там нет орбиталей. А орбиталей там нет потому что scf не только не сошелся, но даже и не начинался. А не начался он потому что строка 3496 - -SOSCF- AND -DIIS- ARE MUTUALLY EXCLUSIVE CONVERGERS вы попробовали одновременно включить два взаимоисключающих конвержера. Если включаете SOSCF - то выключайте DIIS.
С другой стороны непонятно зачем эти пляски. Да, у вас не сошелся HF/6-31++G(d,p), но ведь b3lyp/6-31++G(d,p) - как то сходился - соответственно у него есть шанс сейчас свестись опять без игр с конвержерами.
Мое замечание про заряд я так понимаю вы игнорите, ну дело ваше. :dontknow:

Аватара пользователя
AlexKsen
Сообщения: 602
Зарегистрирован: Пн мар 07, 2011 8:43 am

Re: Hessian в Gamess

Сообщение AlexKsen » Пт ноя 02, 2018 10:59 am

Гесс, спасибо за помощь.
Ваше замечание по поводу заряда, я обязательно учту. Только для начала попробую с зарядом -5, ведь оптимизация то прошла успешно.

Аватара пользователя
Гесс
Сообщения: 13055
Зарегистрирован: Ср фев 15, 2012 11:19 pm

Re: Hessian в Gamess

Сообщение Гесс » Пт ноя 02, 2018 11:02 am

да оптимизировать наверное можно и с зарядом минус 10, смысл только сомнителен.

Ответить

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

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

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