24小时热门版块排行榜    

查看: 2655  |  回复: 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的回帖

独孤狼

金虫 (正式写手)

独行者

【答案】应助回帖

★ ★
亍寞(金币+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的回帖
相关版块跳转 我要订阅楼主 亍寞 的主题更新
信息提示
请填处理意见