24小时热门版块排行榜    

Znn3bq.jpeg
查看: 746  |  回复: 4
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

su571

[交流] [求助]还是L9999

又出现了L9999错误,
昨天用scan扫描了势能曲线,得到AlCd的-1价,多重度为1时的大概能量最低点,
晚上又opt(maxcyc=200,tight),结果还是不行
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

su571

谢谢两位,
K372,你说的我找不到,
liliracial,你的这几种原来都试过,都不行
不过上午改了下初始值,竟然可以了
呵呵,这几个东西算了几个星期了,终于有结果了
4楼2008-04-02 15:39:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

K372

银虫 (小有名气)


zzgyb(金币+1,VIP+0):谢谢你的帮助,欢迎继续关注计算模拟版!
我觉得不像是L9999的问题,因为9999是打印最终结果的,LZ应该看看前面哪里还有出错信息。
我以前也遇到过类似问题,前面的某一步链接(通常在最后几步)出错了但是没有中断运行,传递了错误的数据给后面的链接(通常是L9999),才导致L9999出错的
2楼2008-04-02 09:08:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliracial

银虫 (小有名气)

★ ★ ★
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收敛)。

你自己看看试试。
修炼成一块愚钝的大石头
3楼2008-04-02 09:11:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 273求调剂 +50 麦小叮当 2026-04-06 57/2850 2026-04-10 15:35 by MOF_Catal
[考研] 297求调剂 +22 ljy20040718! 2026-04-03 24/1200 2026-04-09 20:48 by yanenwang
[考研] 284求调剂 +7 让我上岸吧阿西 2026-04-09 7/350 2026-04-09 18:59 by haironglove
[考研] 311求调剂 +6 surte 2026-04-08 13/650 2026-04-09 14:00 by surte
[考研] 材料307分求大佬组收留 +17 Hll胡 2026-04-07 17/850 2026-04-09 10:53 by liuhuiying09
[考研] 一志愿华南师范大学0702物理学305调剂 +4 念常安 2026-04-07 6/300 2026-04-08 22:53 by bljnqdcc
[考研] 285求调剂 +20 哦呦呼o 2026-04-04 20/1000 2026-04-08 22:23 by yutian743
[考研] 264求调剂 +11 麦小叮当 2026-04-07 11/550 2026-04-08 16:05 by 一只好果子?
[考研] 一志愿吉大化学327求调剂 +12 王王白石 2026-04-06 13/650 2026-04-08 16:05 by luoyongfeng
[考研] 281求调剂 +10 椰子蘑菇 2026-04-06 10/500 2026-04-08 11:43 by zzucheup
[考研] 求调剂 +28 111623 2026-04-04 33/1650 2026-04-08 09:24 by 泽润东方
[考研] 生物医药调剂|SCI中科院三区一作+多项科研成果 +8 likangxing 2026-04-07 11/550 2026-04-08 00:02 by lys0704
[考研] 288环境专硕,求调材料方向 +35 lllllos 2026-04-04 39/1950 2026-04-07 23:24 by 一只好果子?
[考研] 344求调剂 +11 魏子per 2026-04-07 11/550 2026-04-07 23:01 by JourneyLucky
[考研] 305求调剂 +4 77Qi 2026-04-06 4/200 2026-04-07 20:06 by shanqishi
[考研] 软工学硕299求调剂 +6 useryy 2026-04-07 6/300 2026-04-07 09:50 by vgtyfty
[考研] 304求调剂 +4 luoye0105 2026-04-05 4/200 2026-04-06 21:05 by 木子君1218
[考研] 071000生物学调剂 +7 拉提桃 2026-04-06 7/350 2026-04-06 18:55 by 52305043001
[考研] 调剂 +4 是可乐不是可乐 2026-04-04 4/200 2026-04-04 19:41 by 唐沐儿
[考研] 266求调剂 +8 学员97LZgn 2026-04-03 8/400 2026-04-04 09:02 by 20021109
信息提示
请填处理意见