Hessian в Gamess
Hessian в Gamess
Добрый вечер.
Помогите разобраться, почему не стали рассчитываться частоты .
Помогите разобраться, почему не стали рассчитываться частоты .
У вас нет необходимых прав для просмотра вложений в этом сообщении.
Re: Hessian в Gamess
Дык это, два экрана от конца аутпута -
строка 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) будут иметь проблемы со сходимостью (поэтому я удивлен что вам удалось его оптимизировать). На мой взгляд сульфогруппы надо забить протонами и тогда возможно проблемность снизится.
строка 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) будут иметь проблемы со сходимостью (поэтому я удивлен что вам удалось его оптимизировать). На мой взгляд сульфогруппы надо забить протонами и тогда возможно проблемность снизится.
Re: Hessian в Gamess
Гесс, спасибо за ответ.
Геометрия системы оптимизирована. На первом этапе попробую 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).
Геометрия системы оптимизирована. На первом этапе попробую 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).
Re: Hessian в Gamess
В этом инпуте у вас DFTTYP=B3LYP, в прошлом такого небыло, оптимизировано оно у вас в чем?
Если оно как то оптимизировалось то странно что у него теперь возникают проблемы со сведением scf в том же(?) методе. Тогда можно было бы прочитать волновую функцию хотя я ужене помню где ее гамесс хранит.
Впрочем независимо от того сведется оно или нет - заряд минус 5 мне не нравится от слова совсем.
Если оно как то оптимизировалось то странно что у него теперь возникают проблемы со сведением scf в том же(?) методе. Тогда можно было бы прочитать волновую функцию хотя я ужене помню где ее гамесс хранит.
Впрочем независимо от того сведется оно или нет - заряд минус 5 мне не нравится от слова совсем.
Re: Hessian в Gamess
Оптимизировано в b3lyp/6-31++G(d,p).
В первый раз проглядел.
В первый раз проглядел.
Re: Hessian в Gamess
SOSCF=.Т. - не помогло решить проблему.
Подскажите, что еще можно попробовать сделать?
Подскажите, что еще можно попробовать сделать?
У вас нет необходимых прав для просмотра вложений в этом сообщении.
Re: Hessian в Gamess
То есть вас не смутило что этот аутпут в 100 с лишним раз меньше чем предыдущий и у вас не возникло желания посмотреть "а че там отличного от предыдущего"?
Я вам могу сказать - он такой маленький потому что там нет орбиталей. А орбиталей там нет потому что scf не только не сошелся, но даже и не начинался. А не начался он потому что строка 3496 - -SOSCF- AND -DIIS- ARE MUTUALLY EXCLUSIVE CONVERGERS вы попробовали одновременно включить два взаимоисключающих конвержера. Если включаете SOSCF - то выключайте DIIS.
С другой стороны непонятно зачем эти пляски. Да, у вас не сошелся HF/6-31++G(d,p), но ведь b3lyp/6-31++G(d,p) - как то сходился - соответственно у него есть шанс сейчас свестись опять без игр с конвержерами.
Мое замечание про заряд я так понимаю вы игнорите, ну дело ваше.
Re: Hessian в Gamess
Гесс, спасибо за помощь.
Ваше замечание по поводу заряда, я обязательно учту. Только для начала попробую с зарядом -5, ведь оптимизация то прошла успешно.
Ваше замечание по поводу заряда, я обязательно учту. Только для начала попробую с зарядом -5, ведь оптимизация то прошла успешно.
Re: Hessian в Gamess
да оптимизировать наверное можно и с зарядом минус 10, смысл только сомнителен.
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 18 гостей