24小时热门版块排行榜    

Znn3bq.jpeg
查看: 655  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 山东高校教师考核超级无底线,员工过不下去啦 +3 qut2026 2026-04-09 4/200 2026-04-10 17:38 by doublecle
[考研] 286求调剂 +7 草木不言 2026-04-04 7/350 2026-04-10 16:15 by 柴小白
[考研] 085410-273求调剂 +4 X1999 2026-04-10 4/200 2026-04-10 15:29 by hemengdong
[考研] 296求调剂 +12 汪!?! 2026-04-08 13/650 2026-04-10 12:09 by 汪!?!
[考研] 一志愿中南大学物理学,英一66,求调剂 +4 长烟旖旎 2026-04-08 5/250 2026-04-10 10:31 by 颖果儿
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +16 努力奋斗112 2026-04-06 18/900 2026-04-10 10:06 by may_新宇
[考研] 083200 305分 求二轮调剂 不接受跨专业 +8 Claireyyyy 2026-04-09 8/400 2026-04-10 10:05 by luosha500
[考研] 一志愿双非085400电子信息344 求调剂,对材料和化学方向也感兴趣 +8 无情的小羊 2026-04-09 9/450 2026-04-10 09:30 by 松花缸1201
[考研] 297求调剂 +8 Kwgyz 2026-04-09 8/400 2026-04-09 23:22 by may_新宇
[考研] 086003调剂求助 +19 苏弋万 2026-04-09 20/1000 2026-04-09 20:48 by bljnqdcc
[考研] 一志愿华工085600 331分 +6 天下ww 2026-04-09 6/300 2026-04-09 18:59 by l_paradox
[考研] 一志愿厦大生物学332求调剂 +10 池池池池池池 2026-04-08 10/500 2026-04-09 17:10 by 独醉梦孤城
[考研] 085501机械专硕 302分 不挑专业求调剂 +5 汪某. 2026-04-09 5/250 2026-04-09 15:38 by 蒋皓禹
[考研] 085801 总分275 本科新能源 求调剂 +8 bradoner 2026-04-08 9/450 2026-04-09 13:43 by only周
[考研] 生物学328分求调剂 +9 闪电kkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
[考研] 302分求调剂 一志愿安徽大学085601 +12 zyx上岸! 2026-04-04 12/600 2026-04-07 02:09 by BruceLiu320
[考研] 复试调剂 +5 asdasdassda 2026-04-05 5/250 2026-04-06 09:32 by dongzh2009
[考研] 一志愿北交大材料工程总分358求调剂 +6 cs0106 2026-04-05 6/300 2026-04-05 16:34 by imissbao
[考研] 085600调剂 +9 东照照照 2026-04-04 9/450 2026-04-05 13:44 by ujn_zhuj
[考研] 一志愿南农090401,268,求调剂 +5 一木鸟然 2026-04-04 5/250 2026-04-04 17:07 by babysonlkd
信息提示
请填处理意见