计算化学公社

标题: 求助如何把MAXIT的最大值200调大 [打印本页]

作者
Author:
498746012    时间: 2021-4-8 15:53
标题: 求助如何把MAXIT的最大值200调大
各位老师好,gamess设定MAXIT最大为200,如何把它调大呀?

作者
Author:
paramecium86    时间: 2021-4-8 16:07
200圈已经是最大了。。。GAMESS的SCF就是挺坑的。除了调整各种SCF收敛的策略 也就是不停的读取之前不收敛计算的dat文件里的$VEC 部分的MO信息。多算几次试试。
作者
Author:
zjxitcc    时间: 2021-4-8 16:34
本帖最后由 zjxitcc 于 2021-4-8 16:38 编辑

如果它有收敛趋势,能量一看稳定了,就差那么几圈,尝试调大MAXIT还有希望;如果相差几十 Hartree这种严重振荡,就别想了,调大MAXIT是没用的,会一直震荡到最大圈数为止。普通GAMESS任务可以去读初猜,因此收敛性是可解决的问题。
作者
Author:
498746012    时间: 2021-4-8 16:39
zjxitcc 发表于 2021-4-8 16:34
如果它有收敛趋势,能量一看稳定了,就差那么几圈,尝试调大MAXIT还有希望;如果相差几十 Hartree这种严重 ...

老师您好,这里的能量一直下降,再有100圈估计就收敛了,但是MAXIT最大就是200
作者
Author:
498746012    时间: 2021-4-8 16:39
paramecium86 发表于 2021-4-8 16:07
200圈已经是最大了。。。GAMESS的SCF就是挺坑的。除了调整各种SCF收敛的策略 也就是不停的读取之前不收敛计 ...

谢谢老师
作者
Author:
zjxitcc    时间: 2021-4-8 16:53
498746012 发表于 2021-4-8 16:39
老师您好,这里的能量一直下降,再有100圈估计就收敛了,但是MAXIT最大就是200

你怎么可能看得出“再有100圈”你说5圈还有人信
作者
Author:
喵星大佬    时间: 2021-4-8 18:09
zjxitcc 发表于 2021-4-8 16:53
你怎么可能看得出“再有100圈”你说5圈还有人信

再有100圈确实令人吃鲸

gamess的scf日常瞬间能量变成正的
作者
Author:
498746012    时间: 2021-4-8 20:00
本帖最后由 498746012 于 2021-4-8 20:03 编辑
zjxitcc 发表于 2021-4-8 16:53
你怎么可能看得出“再有100圈”你说5圈还有人信

谢谢老师,我说错话了,我试了MAXIT=205,报错,超过200阈值
作者
Author:
498746012    时间: 2021-4-8 20:01
喵星大佬 发表于 2021-4-8 18:09
再有100圈确实令人吃鲸

gamess的scf日常瞬间能量变成正的

谢谢老师
作者
Author:
zjxitcc    时间: 2021-4-8 20:45
498746012 发表于 2021-4-8 20:00
谢谢老师,我说错话了,我试了MAXIT=205,报错,超过200阈值

那你到底是什么任务里SCF不收敛呢,普通单点计算,还是LMO-EDA?GKS-EDA?具体哪种
作者
Author:
498746012    时间: 2021-4-8 23:16
zjxitcc 发表于 2021-4-8 20:45
那你到底是什么任务里SCF不收敛呢,普通单点计算,还是LMO-EDA?GKS-EDA?具体哪种

MP2/LMO-EDA
作者
Author:
sobereva    时间: 2021-4-10 06:25
498746012 发表于 2021-4-8 20:00
谢谢老师,我说错话了,我试了MAXIT=205,报错,超过200阈值

自己改一下源码就行了,靠grep很容易搜出代码中相应的行
作者
Author:
498746012    时间: 2021-4-12 18:18
sobereva 发表于 2021-4-10 06:25
自己改一下源码就行了,靠grep很容易搜出代码中相应的行

谢谢sob老师,非常给力
作者
Author:
zjxitcc    时间: 2021-4-12 19:18
498746012 发表于 2021-4-12 18:18
谢谢sob老师,非常给力

建议换成GKS-EDA,它允许读取轨道,就不会有SCF不收敛的问题了。
作者
Author:
498746012    时间: 2021-4-13 08:42
zjxitcc 发表于 2021-4-12 19:18
建议换成GKS-EDA,它允许读取轨道,就不会有SCF不收敛的问题了。

谢谢老师,非常给力。好像只有苏培峰老师那里有GKS-EDA程序,我联系一下苏培老师要一下GKS-EDA的程序包




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3