24小时热门版块排行榜    

查看: 622  |  回复: 4
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 ccc13 的 4 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ccc13

[交流] 【求助】L9999不收敛的错误

计算过程,我试过加大循环的步数(opt(maxcycle=200)),不行;也曾经换用其他方法和小基组先算,然后在换大机组,也不行;发现能量在振荡变化还不小,检查一下初始构型没什么问题,也多次尝试手工调整过,再从头优化,但都失败了。
   其中:
Maximum Force    ………………        no
RMS     Force            ………………  no
Maximum Displacement   …………   yes
RMS     Displacement        ………… yes
最好的情况后两个或一个出YES,而前两个每次都是NO

请问各位有经验的高人们,这种情况该怎么解决?亟待指教,万分感激!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yyx19840628

木虫 (著名写手)

鹰的眼睛,狼的耳朵,豹的速 ...


jghe(金币+1,VIP+0):谢谢你的参与!
对于优化不收敛,即L9999错误,实际上是在规定的步数内没有完成优化,即还没有找到极小值点。(或者对于过渡态优化,还没有找到过渡态)
这有几种可能性:
1. 看一下能量的收敛的情况,可能正在单调减小,眼看有收敛的趋势,这样的情况下,只要加大循环的步数(opt(maxcycle=200)),可能就可以解决问题了。
2. 加大循环步数还不能解决的(循环步数有人说超过200再不收敛,再加也不会有用了,这虽然不一定绝对正确,但200步应该也差不多了),有两种可能。一是查看能量,发现能量在振荡了,且变化已经很小了,这时可能重新算一下,或者构型稍微变一下,继续优化,就可以得到收敛的结果(当然也有麻烦的,看运气和经验了);二是构型变化太大,和你预计的差别过大,这很可能是你的初始构型太差了,优化不知道到哪里去了,这时最好检查一下初始构型,再从头优化。
3. 对于L9999快达到收敛时,考虑减小优化步长有时对于能量振荡的情况也是有用的,opt(maxstep=1).(flyingheart )
人生的最大遗憾莫过于错误地坚持了不该坚持的,轻易地放弃了不该放弃的……
4楼2008-11-10 17:29:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

yuan888

铜虫 (小有名气)

试试这个

★ ★
lei0736(金币+1,VIP+0):谢谢
ccc13(金币+1,VIP+0):谢谢,我试试
加上iop (1/8=1,2/16=1)试试
2楼2008-11-10 12:54:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

张文艳

交流


jghe(金币+1,VIP+0):谢谢你的参与!
当是优化次数不够时,不要急于增加次数 restart,应检查每一次优化的收敛判据,看哪次 ( 第 n 次 ) 最接近收敛,截取该构型输入,加 CalcFc,重新优化;
或:加 Geom=(check,step=n),从 chk 文件中调出第 n 次构型,OPT = CalcFc。
3楼2008-11-10 16:34:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ccc13

引用回帖:
Originally posted by 张文艳 at 2008-11-10 16:34:
当是优化次数不够时,不要急于增加次数 restart,应检查每一次优化的收敛判据,看哪次 ( 第 n 次 ) 最接近收敛,截取该构型输入,加 CalcFc,重新优化;
或:加 Geom=(check,step=n),从 chk 文件中调出第 n 次构 ...

谢谢,但这个类似方法我试过,我是把能量低的那个结构坐标提出来再优化的,可还是得不到结果。
5楼2008-11-11 08:59:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见