应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】IRC出错 Convergence failure -- run terminated.
331/1返回列表
查看: 6241  |  回复: 32
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

renxingye

铁杆木虫 (小有名气)

[交流] 【求助】IRC出错 Convergence failure -- run terminated.

#p irc=(maxpoints=25,calcfc) ub3lyp/6-31g(d,p)




>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UB+HF-LYP) =  -1314.40865938     A.U. after  129 cycles
             Convg  =    0.2624D-01             -V/T =  1.8430
             S**2   =   0.0000
KE= 1.559186642165D+03 PE=-4.414967071091D+03 EE= 1.212077067922D+03
Annihilation of the first spin contaminant:
S**2 before annihilation     0.0000,   after     0.0000
Convergence failure -- run terminated.
Error termination via Lnk1e in /opt//g03/l502.exe

从过渡态的震动模式来看是比较好的！ 而且找过渡态时优化部分和频率部分都有4个Yes
但是算IRC时就出现了上述的问题 不知道是怎么回事
回复此楼

» 猜你喜欢

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

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2010-12-29 10:51:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

auybv

金虫 (小有名气)
★
renxingye(金币+1):谢谢参与
我觉得改变步长是不会有效果的，你的这个计算应该是自洽场迭代的时候出问题了。不论是opt还是irc，每走到一个新的点，都要计算能量，而量化中计算能量就是通过自洽场迭代（实际上相当于求解Schrodinger方程），这个过程默认的最大次数是128。而你的结果中显示 A.U. after  129 cycles，很明显是128圈迭代还没有满足收敛标准。解决这一问题有很多途径，比如使用scfcyc语句增大迭代次数，scfcon降低自洽场收敛标准，另外还可以使用scf=qc强制收敛。因为我不怎么遇到这种情况，所以也不清楚这些关键词是否有效，仅供参考。
赞一下(4人) 回复此楼
9楼2010-12-29 20:51:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

manson1998

金虫 (著名写手)
★
renxingye(金币+1):谢谢参与
步长改变一下看看。没准是你的这一步走到了一个死结构里
赞一下(1人) 回复此楼
2楼2010-12-29 11:44:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

06233053

荣誉版主 (文坛精英)
★
renxingye(金币+1):谢谢参与
顶楼上！
回复此楼
3楼2010-12-29 12:02:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

renxingye

铁杆木虫 (小有名气)
引用回帖:
Originally posted by manson1998 at 2010-12-29 11:44:00:
步长改变一下看看。没准是你的这一步走到了一个死结构里

试过 iop(1/8=10) 没有什么用 还是死掉了
赞一下 回复此楼
4楼2010-12-29 12:56:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

manson1998

金虫 (著名写手)
引用回帖:
Originally posted by renxingye at 2010-12-29 12:56:26:



试过 iop(1/8=10) 没有什么用 还是死掉了

第几个点出的问题？
赞一下 回复此楼
6楼2010-12-29 14:39:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

boylc789

木虫 (著名写手)
★
renxingye(金币+1):谢谢参与
顶楼上
:dnd:
回复此楼
7楼2010-12-29 15:51:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

renxingye

铁杆木虫 (小有名气)
引用回帖:
Originally posted by manson1998 at 2010-12-29 14:39:10:


第几个点出的问题？

一个是第二个点的死了

另外是第二十四个死了
赞一下 回复此楼
8楼2010-12-29 19:52:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

manson1998

金虫 (著名写手)
引用回帖:
Originally posted by renxingye at 2010-12-29 19:52:10:



一个是第二个点的死了

另外是第二十四个死了

forget it

提取第23个点no freeze opt
赞一下 回复此楼
10楼2010-12-30 08:44:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

renxingye

铁杆木虫 (小有名气)
引用回帖:
Originally posted by auybv at 2010-12-29 20:51:11:
我觉得改变步长是不会有效果的，你的这个计算应该是自洽场迭代的时候出问题了。不论是opt还是irc，每走到一个新的点，都要计算能量，而量化中计算能量就是通过自洽场迭代（实际上相当于求解Schrodinger方程），这 ...

我试过scf=(maxcycle=200,xpc) 但还是死掉了
赞一下 回复此楼
11楼2010-12-30 09:10:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xj544

木虫 (小有名气)
★
renxingye(金币+1):谢谢参与
试试加上SCF=QC
二次收敛
赞一下(1人) 回复此楼
12楼2011-01-01 00:18:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

renxingye

铁杆木虫 (小有名气)
引用回帖:
Originally posted by xj544 at 2011-01-01 00:18:36:
试试加上SCF=QC
二次收敛

scf=xpc不是和QC一样的吗？
赞一下 回复此楼
13楼2011-01-02 12:58:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xj544

木虫 (小有名气)
renxingye(金币+1): 2011-01-03 21:39:05
引用回帖:
Originally posted by renxingye at 2011-01-02 12:58:36:

scf=xpc不是和QC一样的吗？

差不多，那你换这个试试：
在IRC 里加上maxcyc=200,不要在SCF里加

就是：#p irc=(maxpoints=25,calcfc,maxcyc=200) ub3lyp/6-31g(d,p)
赞一下 回复此楼
14楼2011-01-03 21:14:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

renxingye

铁杆木虫 (小有名气)
引用回帖:
Originally posted by xj544 at 2011-01-03 21:14:06:


差不多，那你换这个试试：
在IRC 里加上maxcyc=200,不要在SCF里加

就是：#p irc=(maxpoints=25,calcfc,maxcyc=200) ub3lyp/6-31g(d,p)

试了 还是出现原来的错误！
赞一下 回复此楼
15楼2011-01-05 12:08:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xj544

木虫 (小有名气)
引用回帖:
Originally posted by renxingye at 2011-01-05 12:08:52:



试了 还是出现原来的错误！

对了，你是用的SCF=XPC吗？高斯手册里没有这个关键字，只有XQC和QC。
你试试SCF=QC或XQC
赞一下 回复此楼
17楼2011-01-06 15:30:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

renxingye

铁杆木虫 (小有名气)
引用回帖:
Originally posted by xj544 at 2011-01-06 15:30:25:


对了，你是用的SCF=XPC吗？高斯手册里没有这个关键字，只有XQC和QC。
你试试SCF=QC或XQC

SCF=XQC 打错了 呵呵
赞一下 回复此楼
18楼2011-01-06 20:03:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xj544

木虫 (小有名气)
引用回帖:
Originally posted by renxingye at 2011-01-06 20:03:34:



SCF=XQC 打错了 呵呵

你出错的IRC 能看出趋势就不用再做了，如果看不出来，就得重新做。
能看出趋势吗？
赞一下 回复此楼
19楼2011-01-06 21:38:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

renxingye

铁杆木虫 (小有名气)
引用回帖:
Originally posted by xj544 at 2011-01-06 21:38:43:



你出错的IRC 能看出趋势就不用再做了，如果看不出来，就得重新做。
能看出趋势吗？

跑第二个点死了的 不能看

另外是第二十四个死了的可以看去趋势 比较吻合
赞一下 回复此楼
20楼2011-01-07 09:08:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xj544

木虫 (小有名气)
renxingye(金币+2): 2011-01-07 10:23:48
引用回帖:
Originally posted by renxingye at 2011-01-07 09:08:07:


跑第二个点死了的 不能看

另外是第二十四个死了的可以看去趋势 比较吻合

跑第二个点死了的，错误也是LINK9999吗？或者还是另外的错误？
第二个点就死了，有可能是输入文件有问题，你检查一下
赞一下 回复此楼
21楼2011-01-07 09:33:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

renxingye

铁杆木虫 (小有名气)
引用回帖:
Originally posted by xj544 at 2011-01-07 09:33:15:



跑第二个点死了的，错误也是LINK9999吗？或者还是另外的错误？
第二个点就死了，有可能是输入文件有问题，你检查一下

Error termination request processed by link 9999.
Error termination via Lnk1e in /opt//g03/l9999.exe

有一个是9999错误 但是不是我这两个，

9999错误怎么处理？

输入文件 我是过渡态直接做IRC的，重新算力常数的 但是后来我一个读过渡态的chk 还是9999的错误

谢谢！！！

[ Last edited by renxingye on 2011-1-7 at 20:19 ]
赞一下 回复此楼
22楼2011-01-07 10:23:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xj544

木虫 (小有名气)
引用回帖:
Originally posted by renxingye at 2011-01-07 10:23:36:



输入文件 我是过渡态直接做IRC的，重新算力常数的 但是回来我一个读过渡态的chk 还是9999的错误

那你就不要读chk,直接计算力常数，用CalcFc,不用RcFc。
赞一下 回复此楼
23楼2011-01-07 14:56:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

renxingye

铁杆木虫 (小有名气)
引用回帖:
Originally posted by xj544 at 2011-01-07 14:56:04:


那你就不要读chk,直接计算力常数，用CalcFc,不用RcFc。

#p irc=(maxpoints=25,calcfc,maxcycle=200,notrustupdate) ub3lyp/6-31g(d
,p) iop(1/8=6)

还是这样的错误：

Convergence failure -- run terminated.
Error termination via Lnk1e in /opt//g03/l502.exe

怎么回事？？

谢谢
赞一下 回复此楼
24楼2011-01-12 09:27:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tangbohejin

木虫之王 (文学泰斗)
★
renxingye(金币+1):谢谢参与
引用回帖:
Originally posted by auybv at 2010-12-29 20:51:11:
我觉得改变步长是不会有效果的，你的这个计算应该是自洽场迭代的时候出问题了。不论是opt还是irc，每走到一个新的点，都要计算能量，而量化中计算能量就是通过自洽场迭代（实际上相当于求解Schrodinger方程），这 ...

顶---
回复此楼
25楼2011-01-12 11:28:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xj544

木虫 (小有名气)
引用回帖:
Originally posted by renxingye at 2011-01-12 09:27:25:




#p irc=(maxpoints=25,calcfc,maxcycle=200,notrustupdate) ub3lyp/6-31g(d
,p) iop(1/8=6)

还是这样的错误：

Convergence failure -- run terminated.
Error termination via Lnk1e in /opt ...

是l502还是l9999错误？
把link错误上面的信息也贴出来
赞一下 回复此楼
26楼2011-01-12 14:33:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

renxingye

铁杆木虫 (小有名气)
引用回帖:
Originally posted by xj544 at 2011-01-12 14:33:44:



是l502还是l9999错误？
把link错误上面的信息也贴出来

>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UB+HF-LYP) =  -1342.71493389     A.U. after  129 cycles
             Convg  =    0.1621D-01             -V/T =  1.8702
             S**2   =   0.0000
KE= 1.542956175444D+03 PE=-4.208117663644D+03 EE= 1.002083822568D+03
Annihilation of the first spin contaminant:
S**2 before annihilation     0.0000,   after     0.0000
Convergence failure -- run terminated.
Error termination via Lnk1e in /opt//g03/l502.exe

谢谢
赞一下 回复此楼
27楼2011-01-13 09:25:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xj544

木虫 (小有名气)
renxingye(金币+1): 2011-01-13 13:39:49
引用回帖:
Originally posted by renxingye at 2011-01-13 09:25:43:



>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UB+HF-LYP) =  -1342.71493389     A.U. after  129 cycles
             Convg  =    0.1621D-01            ...

像这种问题，一般有两个解决方法：
1、加上SCF=(maxcyc=N)，如N=500.
2、加SCF=QC
如果这两个方法都解决不了，那就是初始构型不太好，调整初始构型试试吧。
赞一下 回复此楼
28楼2011-01-13 09:49:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

renxingye

铁杆木虫 (小有名气)
引用回帖:
Originally posted by xj544 at 2011-01-13 09:49:07:


像这种问题，一般有两个解决方法：
1、加上SCF=(maxcyc=N)，如N=500.
2、加SCF=QC
如果这两个方法都解决不了，那就是初始构型不太好，调整初始构型试试吧。

初始构型不太好?? 这是做IRC的。可以调整初始结构？？
赞一下 回复此楼
29楼2011-01-13 13:40:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xj544

木虫 (小有名气)
引用回帖:
Originally posted by renxingye at 2011-01-13 13:40:56:



初始构型不太好?? 这是做IRC的。可以调整初始结构？？

你试试第一种方法
赞一下 回复此楼
30楼2011-01-17 08:38:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chinasd.sl

新虫 (初入文坛)
★
renxingye(金币+1): 谢谢参与
我做opt时也遇到了相同的问题，你最后怎么解决的啊，告诉一下，chinasd.sl@163.com
赞一下 回复此楼
31楼2012-03-20 10:00:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

renxingye

铁杆木虫 (小有名气)
引用回帖:
1958830楼: Originally posted by chinasd.sl at 2012-03-20 10:00:31:
我做opt时也遇到了相同的问题，你最后怎么解决的啊，告诉一下，chinasd.sl@163.com

你结构优化出错？还是IRC出错？

结构优化出错，改初始结构比加关键词好用！个人经验。呵呵
赞一下(1人) 回复此楼
32楼2012-03-20 10:59:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nedved33

金虫 (小有名气)
★
renxingye(金币+1): 谢谢参与
高斯小白想问下，怎么看第几个点出问题?
第几个点是指第几次迭代还是循环？
赞一下 回复此楼
33楼2012-03-22 19:58:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
zyxme5楼
2010-12-29 13:23   回复  
renxingye(金币+1):谢谢参与
lldddd16楼
2011-01-05 12:11   回复  
renxingye(金币+1):谢谢参与
相关版块跳转 我要订阅楼主 renxingye 的主题更新
331/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭