24小时热门版块排行榜    

Znn3bq.jpeg
查看: 656  |  回复: 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的回帖

张文艳

交流


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的回帖
查看全部 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的回帖

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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料085601调剂 +24 何润采123 2026-04-10 26/1300 2026-04-10 22:10 by 猪会飞
[考研] 一志愿北理工298英一数二求调剂 +13 Reframe 2026-04-10 15/750 2026-04-10 17:25 by laoshidan
[考研] 求调剂 +9 翩翩一书生 2026-04-09 9/450 2026-04-10 11:43 by wp06
[考研] 298求调剂 +8 钉叮咚冬瓜 2026-04-07 10/500 2026-04-10 11:25 by 高维春
[考研] 22408 352分求调剂0854类 +3 努力的夏末 2026-04-09 3/150 2026-04-10 10:12 by 314126402
[考研] 材料调剂 +11 一样YWY 2026-04-05 11/550 2026-04-10 09:32 by 钟洲2011
[考研] 08600生物与医药-327 +10 18755400796 2026-04-05 10/500 2026-04-10 08:14 by kangsm
[考研] 求调剂 +15 张zic 2026-04-05 16/800 2026-04-10 08:12 by kangsm
[考研] 296求调剂 +6 汪!?! 2026-04-09 6/300 2026-04-09 20:49 by lbsjt
[考研] 085400 328分 求调剂 +8 喂你一个大橙子 2026-04-09 12/600 2026-04-09 19:20 by zl8213662
[考研] 一志愿华工085600 331分 +6 天下ww 2026-04-09 6/300 2026-04-09 18:59 by l_paradox
[考研] 0703化学调剂325分 +13 15771691647 2026-04-04 15/750 2026-04-09 16:55 by 15771691647
[考研] 求机械专硕297第二批调剂 +5 拾柒12。 2026-04-08 5/250 2026-04-09 16:43 by 允当适度
[考研] 一志愿郑州大学 22408 305分求调剂 +3 安小满zzz 2026-04-08 3/150 2026-04-09 16:16 by wp06
[考研] 生物学328分求调剂 +9 闪电kkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
[考研] 材料调剂 +17 小刘同学吖吖 2026-04-06 18/900 2026-04-07 11:41 by 诗与自由
[考研] (调剂)一志愿报考哈尔滨工业大学0857资源与环境专业378分考生 +7 狠狠加油 2026-04-05 8/400 2026-04-06 16:52 by momo皓
[考研] 求调剂 +7 张.1 2026-04-05 7/350 2026-04-05 20:40 by 啵啵啵0119
[考研] 调剂 +3 李广火 2026-04-05 3/150 2026-04-05 18:57 by 蓝云思雨
[考研] 考研生物学考A区211,初试322,科目生化和生物综合,求调剂 +6 。。。54 2026-04-03 6/300 2026-04-05 14:54 by JOKER0401
信息提示
请填处理意见