24小时热门版块排行榜

>论坛更新日志 (1273)
>虫友互识 (34)
>论文投稿 (24)
>导师招生 (23)
>考研 (21)
>基金申请 (17)
>教师之家 (17)
>硕博家园 (17)
>找工作 (17)
>考博 (16)
>公派出国 (16)
>休闲灌水 (15)
>博后之家 (14)
>论文道贺祈福 (11)
>有机交流 (11)
>文献求助 (9)

返回列表

当前主题已经存档。

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

[交流] [求助]还是L9999

又出现了L9999错误，
昨天用scan扫描了势能曲线，得到AlCd的-1价，多重度为1时的大概能量最低点，
晚上又opt(maxcyc=200,tight)，结果还是不行

回复此楼

» 猜你喜欢

1楼 2008-04-02 09:01:11

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

liliracial

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 515.5
红花: 1
帖子: 256
在线: 49.1小时
虫号: 520204
注册: 2008-03-07
专业: 催化化学

★ ★ ★
zzgyb(金币+2,VIP+0):谢谢你的帮助，欢迎继续关注计算模拟版！
zzgyb(金币+1,VIP+0):谢谢你的帮助，欢迎继续关注计算模拟版！

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

你自己看看试试。