Сходимость геометрии S1 состояния CITYP=CIS PC GAMESS и GAU.

[lmihael]

и Gaussian
Попробовал оптимизировать геометрию в переходном состояние для переноса протона в гидроксибензтиазоле.
Gamess вначале хорошо оптимизирует геометрию, но дойдя до точки близкой к минимому начинает плавать вокруг неё:
-1024.330
...
...
-1024.352
-1024.353
-1024.352
-1024.351
-1024.353
-1024.352
-1024.352
-1024.353
такого типа результат (H)
При попытке начать новый расчет используя лучшую геометрию - начинается такое же плавание выше исходного минимума, уменьшение DXMAX ничего не дает. Уменьшение CPTOL до 6.0E-05 привподит к очень плохой сходимости градиента, а плавание остается.
Использовал METHOD=QA и RFO. Вот мой input
$CONTRL INTTYP=HONDO EXETYP=RUN SCFTYP=RHF
mplevl=0 RUNTYP=optimize COORD=ZMT UNITS=ANGS MAXIT=60
CITYP=CIS NZVAR=69 $END
$statpt METHOD=qa NSTEP=200000 DXMAX=0.01 IFREEZ(1)=55 $end
$SYSTEM timlim=6000 MEMORY=199990000 $END
$CIS NSTATE=1 ISTATE=1 MULT=1 DAVCVG=8.0E-05 CPTOL=9.0E-06
RDCISV=.t. $END
$SCF DIRSCF=.t. NCONV=5 $END
$GUESS GUESS=moread norb=159 $END
$BASIS GBASIS=n311 NGAUSS=6 NDFUNC=1 NPFUNC=3 NFFUNC=0
DIFFSP=.f. DIFFS=.f. $END
$DATA
HBT same
C1
C
C 1 1.4942235
C 2 1.4065965 1 118.8767241
C 3 1.3789253 2 120.6498045 1 -0.0008245 0
C 4 1.4079475 3 120.4125144 2 -0.0007124 0
C 5 1.3848260 4 121.7650311 3 -0.0002337 0
O 2 1.2657370 1 119.4856602 5 179.9998858 0
C 1 1.4000767 2 116.6062359 3 -179.9998829 0
N 8 1.3433190 1 119.9859619 6 179.9998464 0
C 9 1.3440229 8 114.1727319 1 179.9994964 0
C 10 1.4083884 9 113.7757138 8 -0.0001631 0
S 11 1.7521329 10 110.0407232 9 -0.0001398 0
C 10 1.4052603 9 126.7363287 8 179.9999893 0
C 13 1.3723157 10 118.6024004 11 -0.0031797 0
C 14 1.4012406 13 121.5427405 10 -0.0047273 0
C 15 1.3828281 14 120.5424999 13 -0.0028969 0
H 3 1.0726443 2 117.1561683 7 0.0001888 0
H 4 1.0728354 3 120.1396246 2 179.9999055 0
H 5 1.0733369 4 119.3170237 3 -179.9999272 0
H 6 1.0732729 5 119.6081888 4 -179.9998363 0
H 7 1.2000000 2 108.0043145 3 -179.9996910 0
H 13 1.0729185 10 119.7830046 11 -179.9999497 0
H 14 1.0738758 13 119.3362869 10 -179.9998656 0
H 15 1.0730873 14 119.7792989 13 179.9996614 0
H 16 1.0733694 15 120.6276406 14 179.9998028 0
$END
в Gaussian'e в дальнейшем минимизируется нормально (хотя каждая итерация длится намного дольше) и оказывается что эта геометрия очень близка к оптимальной.

Подскажите какие использовать параметры для нормальной сходимости в GAMESS?

[VTur]

Попробуйте считать вторые производные ($STATPT HESS=CALC). Расчет каждого шага возрастёт в разы. Если очень долго добавьте IHREP=2

[lmihael]

Попробуйте считать вторые производные

Спасибо большое. Попробую.

[Alex Granovsky]

Ñòðàííûå Âû ëþäè. Âåäü åñëè âèäíî, ÷òî òî÷íîñòè ïî óìîë÷àíèþ íå õâàòàåò, çà÷åì åå äîïîëíèòåëüíî ïîíèæàòü.

Ïðèëàãàþ "ïðàâèëüíûé âõîäíîé ôàéë"


$CONTRL INTTYP=HONDO RUNTYP=optimize COORD=ZMT MAXIT=60 CITYP=CIS NZVAR=1
ICUT=10 $END
$SYSTEM timlim=6000 MEMORY=199990000 $END
$p2p p2p=.t. dlb=.t. $end
$BASIS GBASIS=n311 NGAUSS=6 NDFUNC=1 NPFUNC=3 $END
$SCF DIRSCF=.t. NCONV=5 $END
$CIS NSTATE=1 ISTATE=1 MULT=1 $END
$statpt METHOD=qa NSTEP=200000 $end
$zmat dlc=1 auto=1 nonvdw(1)=7,21,9,21 ifzmat(1)=1,7,21 $end
$DATA
HBT same
C1
C
C 1 1.4942235
C 2 1.4065965 1 118.8767241
C 3 1.3789253 2 120.6498045 1 -0.0008245 0
C 4 1.4079475 3 120.4125144 2 -0.0007124 0
C 5 1.3848260 4 121.7650311 3 -0.0002337 0
O 2 1.2657370 1 119.4856602 5 179.9998858 0
C 1 1.4000767 2 116.6062359 3 -179.9998829 0
N 8 1.3433190 1 119.9859619 6 179.9998464 0
C 9 1.3440229 8 114.1727319 1 179.9994964 0
C 10 1.4083884 9 113.7757138 8 -0.0001631 0
S 11 1.7521329 10 110.0407232 9 -0.0001398 0
C 10 1.4052603 9 126.7363287 8 179.9999893 0
C 13 1.3723157 10 118.6024004 11 -0.0031797 0
C 14 1.4012406 13 121.5427405 10 -0.0047273 0
C 15 1.3828281 14 120.5424999 13 -0.0028969 0
H 3 1.0726443 2 117.1561683 7 0.0001888 0
H 4 1.0728354 3 120.1396246 2 179.9999055 0
H 5 1.0733369 4 119.3170237 3 -179.9999272 0
H 6 1.0732729 5 119.6081888 4 -179.9998363 0
H 7 1.2000000 2 108.0043145 3 -179.9996910 0
H 13 1.0729185 10 119.7830046 11 -179.9999497 0
H 14 1.0738758 13 119.3362869 10 -179.9998656 0
H 15 1.0730873 14 119.7792989 13 179.9996614 0
H 16 1.0733694 15 120.6276406 14 179.9998028 0
$END

è Gaussian
Ïîïðîáîâàë îïòèìèçèðîâàòü ãåîìåòðèþ â ïåðåõîäíîì ñîñòîÿíèå äëÿ ïåðåíîñà ïðîòîíà â ãèäðîêñèáåíçòèàçîëå.
Gamess âíà÷àëå õîðîøî îïòèìèçèðóåò ãåîìåòðèþ, íî äîéäÿ äî òî÷êè áëèçêîé ê ìèíèìîìó íà÷èíàåò ïëàâàòü âîêðóã íå¸:
-1024.330
...
...
-1024.352
-1024.353
-1024.352
-1024.351
-1024.353
-1024.352
-1024.352
-1024.353
òàêîãî òèïà ðåçóëüòàò (H)
Ïðè ïîïûòêå íà÷àòü íîâûé ðàñ÷åò èñïîëüçóÿ ëó÷øóþ ãåîìåòðèþ - íà÷èíàåòñÿ òàêîå æå ïëàâàíèå âûøå èñõîäíîãî ìèíèìóìà, óìåíüøåíèå DXMAX íè÷åãî íå äàåò. Óìåíüøåíèå CPTOL äî 6.0E-05 ïðèâïîäèò ê î÷åíü ïëîõîé ñõîäèìîñòè ãðàäèåíòà, à ïëàâàíèå îñòàåòñÿ.
Èñïîëüçîâàë METHOD=QA è RFO. Âîò ìîé input
$CONTRL INTTYP=HONDO EXETYP=RUN SCFTYP=RHF
mplevl=0 RUNTYP=optimize COORD=ZMT UNITS=ANGS MAXIT=60
CITYP=CIS NZVAR=69 $END
$statpt METHOD=qa NSTEP=200000 DXMAX=0.01 IFREEZ(1)=55 $end
$SYSTEM timlim=6000 MEMORY=199990000 $END
$CIS NSTATE=1 ISTATE=1 MULT=1 DAVCVG=8.0E-05 CPTOL=9.0E-06
RDCISV=.t. $END
$SCF DIRSCF=.t. NCONV=5 $END
$GUESS GUESS=moread norb=159 $END
$BASIS GBASIS=n311 NGAUSS=6 NDFUNC=1 NPFUNC=3 NFFUNC=0
DIFFSP=.f. DIFFS=.f. $END
$DATA
HBT same
C1
C
C 1 1.4942235
C 2 1.4065965 1 118.8767241
C 3 1.3789253 2 120.6498045 1 -0.0008245 0
C 4 1.4079475 3 120.4125144 2 -0.0007124 0
C 5 1.3848260 4 121.7650311 3 -0.0002337 0
O 2 1.2657370 1 119.4856602 5 179.9998858 0
C 1 1.4000767 2 116.6062359 3 -179.9998829 0
N 8 1.3433190 1 119.9859619 6 179.9998464 0
C 9 1.3440229 8 114.1727319 1 179.9994964 0
C 10 1.4083884 9 113.7757138 8 -0.0001631 0
S 11 1.7521329 10 110.0407232 9 -0.0001398 0
C 10 1.4052603 9 126.7363287 8 179.9999893 0
C 13 1.3723157 10 118.6024004 11 -0.0031797 0
C 14 1.4012406 13 121.5427405 10 -0.0047273 0
C 15 1.3828281 14 120.5424999 13 -0.0028969 0
H 3 1.0726443 2 117.1561683 7 0.0001888 0
H 4 1.0728354 3 120.1396246 2 179.9999055 0
H 5 1.0733369 4 119.3170237 3 -179.9999272 0
H 6 1.0732729 5 119.6081888 4 -179.9998363 0
H 7 1.2000000 2 108.0043145 3 -179.9996910 0
H 13 1.0729185 10 119.7830046 11 -179.9999497 0
H 14 1.0738758 13 119.3362869 10 -179.9998656 0
H 15 1.0730873 14 119.7792989 13 179.9996614 0
H 16 1.0733694 15 120.6276406 14 179.9998028 0
$END
â Gaussian'e â äàëüíåéøåì ìèíèìèçèðóåòñÿ íîðìàëüíî (õîòÿ êàæäàÿ èòåðàöèÿ äëèòñÿ íàìíîãî äîëüøå) è îêàçûâàåòñÿ ÷òî ýòà ãåîìåòðèÿ î÷åíü áëèçêà ê îïòèìàëüíîé.

Ïîäñêàæèòå êàêèå èñïîëüçîâàòü ïàðàìåòðû äëÿ íîðìàëüíîé ñõîäèìîñòè â GAMESS?

[lmihael]

Странные Вы люди. Ведь если видно, что точности по умолчанию не хватает, зачем ее дополнительно понижать.

Прилагаю "правильный входной файл"

Спасибо большое! Это действительно работает!