24小时热门版块排行榜

>论坛更新日志 (3163)
>文献求助 (217)
>虫友互识 (217)
>导师招生 (191)
>考博 (113)
>休闲灌水 (73)
>招聘信息布告栏 (64)
>硕博家园 (53)
>博后之家 (46)
>论文投稿 (45)
>基金申请 (44)
>论文道贺祈福 (40)
>考研 (39)
>公派出国 (37)
>找工作 (36)
>教师之家 (27)

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

irishes

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 299.6
散金: 5
帖子: 66
在线: 59.4小时
虫号: 508955
注册: 2008-02-22
性别: MM
专业: 质谱分析

[交流] [求助] 优化构型总是不收敛，老出LINK9999错误，加些什么关键词可以使其收敛呢，谢谢

优化构型总是不收敛，老出LINK9999错误，加些什么关键词可以使其收敛呢，谢谢

回复此楼

» 猜你喜欢

1楼 2008-04-07 11:04:54

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lixiaona158

木虫 (著名写手)

自由人

应助: 1 (幼儿园)
金币: 2664.3
散金: 160
红花: 7
帖子: 1583
在线: 158.7小时
虫号: 473345
注册: 2007-12-05
专业: 理论和计算化学

★ ★ ★ ★ ★
lei0736(金币+5,VIP+0):多谢够详细

优化构型不收敛有很多原因，不光是优化的圈数不够，下面有一些意见，你可以尝试下：

1. 使用小基组，或低级算法计算，得到scf收敛的波函数，用guess=read读初始波函数。
2. 使用scf=qc，这个计算会慢，而且需要用stable关键字来测试结果是否波函数稳定。如果这个还不收敛，会提示L508错误。
3. 改变键长，一般是缩小一点，有时会有用。
4. 计算相同体系的其他电子态，比如相应的阴离子、阳离子体系或单重态体系，得到的收敛波函数作为初始猜测进行计算。

对于优化不收敛,即L9999错误，实际上是在规定的步数内没有完成优化，即还没有找到极小值点。（或者对于过渡态优化，还没有找到过渡态）这有几种可能性：
1. 看一下能量的收敛的情况，可能正在单调减小，眼看有收敛的趋势，这样的情况下，只要加大循环的步数（opt(maxcycle=200)），可能就可以解决问题了。
2. 加大循环步数还不能解决的（循环步数有人说超过200再不收敛，再加也不会有用了，这虽然不一定绝对正确，但200步应该也差不多了），有两种可能。一是查看能量，发现能量在振荡了，且变化已经很小了，这时可能重新算一下，或者构型稍微变一下，继续优化，就可以得到收敛的结果（当然也有麻烦的，看运气和经验了）；二是构型变化太大，和你预计的差别过大，这很可能是你的初始构型太差了，优化不知道到哪里去了，这时最好检查一下初始构型，再从头优化。
3. 对于L9999快达到收敛时,考虑减小优化步长有时对于能量振荡的情况也是有用的，opt(maxstep=1).
一个建议是，对于大体系，难收敛体系，先用小基组，低精度算法优化一下，以得到较好的初始构型，再用高精度的计算接着算。如果前面的方法保留了chk文件，重新计算时需要使用 geom=allcheck 读入构型（就不必麻烦地写构型了）， guess=read（读入初始波函数，可以加快第一步SCF收敛）。