24小时热门版块排行榜    

查看: 2744  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

亍寞

金虫 (小有名气)

[求助] opt=TS 循环步骤问题

命令这么写的,但是走完第一圈就9999出错了,请高人指教!

%mem=15000MB
%nprocshared=8
#p opt(maxcyc=300,TS,calcfc) freq b3lyp/6-311++g(2d,2p) scf(maxcyc=300,tight)

该怎么弄呢?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

亍寞

金虫 (小有名气)

引用回帖:
Originally posted by 独孤狼 at 2011-06-23 14:50:55:
L9999的错误的引发一般是因为Optimization stopped.这个可以从输出文件中看到,这个产生的原因有很多,其中一个是
-- Number of steps exceeded,  NStep= N (N=当前优化的步数)
另一个是-- Wrong numb ...

谢谢!
8楼2011-06-23 14:56:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

gmy1990

荣誉版主 (著名写手)

优秀版主优秀版主

【答案】应助回帖

亍寞(金币+2): 2011-06-23 14:51:46
看到好几个帖子说L999问题了,我把量化版的FAQ上面的贴上来
L9999错误是属于构型优化没有找到最后结果,即在规定的优化步骤内没有找到最终的合适的结构。(L9999错误)。对于优化不收敛,即L9999错误,实际上是在规定的步数内没有完成优化,即还没有找到极小值点。(或者对于过渡态优化,还没有找到过渡态),可以增加优化步骤,或是拿最后的结构 接着进行优化。此外这有几种可能性:1. 看一下能量的收敛的情况,可能正在单调减小,眼看有收敛的趋势,这样的情况下,只要加大循环的步数opt(maxcycle=200),可能就可以解决问题了。2. 加大循环步数还不能解决的(循环步数有人说超过200再不收敛,再加也不会有用了,这虽然不一定绝对正确,但200步应该也差不多了),有两种可能。一是查看能量,发现能量在振荡了,且变化已经很小了,这时可能重新算一下,或者构型稍微变一下,继续优化,就可以得到收敛的结果(当然也有麻烦的,看运气和经验了);二是构型变化太大,和你预计的差别过大,这很可能是你的初始构型太差了,优化不知道到哪里去了,这时最好检查一下初始构型,再从头优化。3. 对于L9999快达到收敛时,考虑减小优化步长有时对于能量振荡的情况也是有用的,opt(maxstep=1).一个建议是,对于大体系,难收敛体系,先用小基组,低精度算法优化一下,以得到较好的初始构型,再用高精度的计算接着算。如果前面的方法保留了chk文件,重新计算时需要使用 geom=allcheck 读入构型(就不必麻烦地写构型了), guess=read(读入初始波函数,可以加快第一步SCF收敛)。
2楼2011-06-23 12:03:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fatpig8832

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★
亍寞(金币+4): 2011-06-23 14:51:51
gmy1990(金币+3): 非常感谢你对论坛的贡献,奖励一下! 2011-06-23 15:34:08
2楼的斑竹貌似没看清楚问题...

楼主的TS只能做一步,估计是少了个 noeigen

把opt那部分写成 opt(maxcyc=300,TS,calcfc,noeigen) 再试试
3楼2011-06-23 14:17:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

亍寞

金虫 (小有名气)

引用回帖:
Originally posted by gmy1990 at 2011-06-23 12:03:30:
看到好几个帖子说L999问题了,我把量化版的FAQ上面的贴上来
L9999错误是属于构型优化没有找到最后结果,即在规定的优化步骤内没有找到最终的合适的结构。(L9999错误)。对于优化不收敛,即L9999错误,实际上是在 ...

嗯,不管怎么说还是感谢~~
4楼2011-06-23 14:48:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085602化学工程350,调剂,有没有211的 +3 利好利好. 2026-03-02 5/250 2026-03-02 21:05 by 59上学吧
[考研] 接收调剂 +6 津萌津萌 2026-03-02 14/700 2026-03-02 20:59 by 汪!?!
[考研] 化工335求调剂 +3 摸摸猫猫头 2026-03-02 3/150 2026-03-02 20:55 by 平凡的七月
[考研] 理学,工学,农学调剂,少走弯路,这里欢迎您! +5 likeihood 2026-03-02 8/400 2026-03-02 20:39 by ji493940045
[考研] 材料学硕318求调剂 +10 February_Feb 2026-03-01 10/500 2026-03-02 20:18 by hypershenger
[考研] 085600材料工程一志愿中科大总分312求调剂 +9 吃宵夜1 2026-02-28 11/550 2026-03-02 20:14 by hypershenger
[考研] 材料复试调剂 +5 学材料的点 2026-03-01 6/300 2026-03-02 20:01 by hypershenger
[考研] 调剂 +3 13853210211 2026-03-02 3/150 2026-03-02 19:59 by hypershenger
[考研] 306分材料调剂 +5 chuanzhu川烛 2026-03-01 6/300 2026-03-02 19:51 by 张晓芳0105
[考研] 261求调剂 +3 陆lh 2026-03-01 3/150 2026-03-02 19:32 by zhukairuo
[考研] 085600 英一数二272求调剂 5+6 vida_a 2026-03-01 16/800 2026-03-02 19:13 by zhukairuo
[考研] 284求调剂 +5 天下熯 2026-03-02 5/250 2026-03-02 18:38 by caszguilin
[考研] 一志愿东北大学材料专硕328,求调剂 +3 shs1083 2026-03-02 3/150 2026-03-02 17:27 by houyaoxu
[考研] 材料与化工328求调剂 +3 。,。,。,。i 2026-03-02 3/150 2026-03-02 13:09 by houyaoxu
[基金申请] 此成果不能导入原因:元数据必填信息不完整,可 进行补充。 +4 Kittylucky 2026-03-02 5/250 2026-03-02 11:07 by jurkat.1640
[考研] 化工专硕342,一志愿大连理工大学,求调剂 +6 kyf化工 2026-02-28 7/350 2026-03-02 10:56 by 无际的草原
[考研] 322求调剂 +3 熊境喆 2026-03-01 3/150 2026-03-02 08:44 by houyaoxu
[考研] 一志愿中南大学理学化学 +4 15779376950 2026-03-01 5/250 2026-03-01 19:00 by Fff-1
[考研] 307求调剂 +5 wyyyqx 2026-03-01 5/250 2026-03-01 15:21 by Fff-1
[论文投稿] Optics letters投稿被拒求助 30+3 luckyry 2026-02-26 4/200 2026-03-01 09:06 by babero
信息提示
请填处理意见