24小时热门版块排行榜    

查看: 621  |  回复: 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的回帖

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的回帖

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的回帖

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的回帖
相关版块跳转 我要订阅楼主 ccc13 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见