24小时热门版块排行榜    

查看: 2213  |  回复: 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)

该怎么弄呢?
回复此楼

» 收录本帖的淘贴专辑推荐

科研资料 computational chemistry (gaussian) Gaussian

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

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的回帖

亍寞

金虫 (小有名气)

引用回帖:
Originally posted by fatpig8832 at 2011-06-23 14:17:39:
2楼的斑竹貌似没看清楚问题...

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

把opt那部分写成 opt(maxcyc=300,TS,calcfc,noeigen) 再试试

了解了,谢谢!
5楼2011-06-23 14:49:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤狼

金虫 (正式写手)

独行者

【答案】应助回帖

★ ★
亍寞(金币+2): 2011-06-23 14:56:44
gmy1990(金币+2): 非常感谢你对论坛的贡献,奖励一下! 2011-06-23 15:34:43
L9999的错误的引发一般是因为Optimization stopped.这个可以从输出文件中看到,这个产生的原因有很多,其中一个是
-- Number of steps exceeded,  NStep= N (N=当前优化的步数)
另一个是-- Wrong number of Negative eigenvalues: Desired=  1 Actual=  N (N=当前Hessian矩阵的本征值数,找过渡态需要该值为1,通过添加noeigentest参数可以取消该检查,虽然有时也能找到过渡态,但一般情况下这说明过渡态的初始猜测并不是很合理)还有一个是
Optimization aborted.
   -- No acceptable step.
这个很明显,构象优化出错了,解决方法调构象或者更改计算方法或条件……
还有很多其他的原因,我这就不一一列举了。所以解决l9999的错误最好的方法是在你输出文件中仔细找一下到底是什么原因导致计算终止,然后对症下药盲目的加maxcyc=N opt=loose opt=gdiis等等是没什么意义的
轮回的过往,注定的未来;独孤的行者,孤独的路途……
6楼2011-06-23 14:50:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤狼

金虫 (正式写手)

独行者

【答案】应助回帖

亍寞(金币+2): 2011-06-24 05:37:41
对于你这个体系,由于第一步都没算完,那很有可能是三楼说的那样初始猜测不是很好导致Hessian矩阵有2个以上的负本征值,所以优化结束,可以通过添加opt=noeigentest解决……
轮回的过往,注定的未来;独孤的行者,孤独的路途……
7楼2011-06-23 14:53:34
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 亍寞 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 对象没有,还非常想读博,难以抉择 +28 pvrw0224 2024-06-23 52/2600 2024-06-26 20:22 by pvrw0224
[基金申请] 自然基金变化?? +15 今晚推荐22 2024-06-21 25/1250 2024-06-26 19:40 by 朋友别哭010
[有机交流] 有机物的酸酐如何准确表征 10+3 方酱 2024-06-24 7/350 2024-06-26 17:53 by 宁静远行
[硕博家园] 联培博士文章第一单位署名问题交流~ +9 橙成成c 2024-06-23 21/1050 2024-06-26 16:45 by 橙成成c
[找工作] 老博士选高校 5+5 woodhead8386 2024-06-24 8/400 2024-06-26 16:35 by 大大熊呀
[硕博家园] 怎么带研究生? +24 豆豆小小2 2024-06-20 42/2100 2024-06-26 16:13 by wangyu1986w
[基金申请] 评审专家会不会很在意申请人的单位啊 +11 lancet0903 2024-06-24 13/650 2024-06-26 11:42 by 漠上藜梭
[公派出国] CSC博士联培对将来就业有用吗 +4 也就这样 2024-06-22 4/200 2024-06-26 10:46 by lxs7405
[有机交流] 酚羟基甲基化 10+4 A好运来啦啦啦 2024-06-25 6/300 2024-06-26 10:08 by 88817753
[有机交流] 过柱子,产品,杂质在是 pe:ea=100:1 也一起出来? +5 w256 2024-06-25 5/250 2024-06-26 09:31 by 小木木cc
[基金申请] 专科院校申请国自然 +14 孤单12站 2024-06-24 20/1000 2024-06-26 09:09 by chong2024
[基金申请] 博后面上到底什么时候出结果??? +7 爱学lsy 2024-06-24 7/350 2024-06-25 19:24 by 暴走的蒜泥
[第一性原理] Vasp 版权问题 10+4 竹叶青9 2024-06-22 5/250 2024-06-25 14:58 by 无所谓109
[基金申请] 演员杨幂在核心期刊发表论文,中国广播电视学期刊 +9 babu2015 2024-06-22 12/600 2024-06-25 11:05 by 我是勇敢牛牛
[教师之家] 有没有今年的影响因子? +6 jurkat.1640 2024-06-22 10/500 2024-06-24 22:48 by 潇湘之迷
[金属] 寻找钛合金热压缩代做 +4 liuyang358 2024-06-23 4/200 2024-06-24 13:50 by 搬砖狗不放弃
[基金申请] 说博后基金7月出的真打电话了吗? +12 antonysole 2024-06-24 14/700 2024-06-24 13:39 by sizhouyi
[有机交流] 生成亚胺的反应怎么能进行完全 +3 1369836 2024-06-23 3/150 2024-06-23 18:44 by hwqMSE
[基金申请] 请问三个字的中文名字,在文中到底怎么表示? 5+3 zhongyuan764 2024-06-22 6/300 2024-06-23 17:42 by 鱼翔浅底1
[论文投稿] ACS 编辑的意见 10+3 哈哈妞1993 2024-06-20 3/150 2024-06-21 17:06 by 投个论文
信息提示
请填处理意见