24小时热门版块排行榜    

查看: 5103  |  回复: 32

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时就出现了上述的问题 不知道是怎么回事
回复此楼

» 猜你喜欢

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

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

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

auybv

金虫 (小有名气)



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

renxingye(金币+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的回帖
引用回帖:
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的回帖
引用回帖:
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
二次收敛
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出错?

结构优化出错,改初始结构比加关键词好用!个人经验。呵呵
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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 290求调剂 +4 @将就将就看 2026-03-10 8/400 2026-03-14 14:23 by 千千运气
[考研] 308 085701 四六级已过求调剂 +7 温乔乔乔乔 2026-03-12 14/700 2026-03-14 10:49 by JourneyLucky
[考研] 一志愿北京化工大学材料与化工296分求调剂 +16 稻妻小编 2026-03-09 18/900 2026-03-14 02:00 by JourneyLucky
[考研] 333求调剂 +3 球球古力 2026-03-09 3/150 2026-03-14 01:57 by JourneyLucky
[考研] 295复试调剂 +5 简木ChuFront 2026-03-09 5/250 2026-03-14 01:29 by JourneyLucky
[考研] 一志愿湖师大化学289求调剂 +6 XMCMM3.14159 2026-03-10 6/300 2026-03-14 00:28 by JourneyLucky
[考研] 0805,333求调剂 +3 112253525 2026-03-10 3/150 2026-03-13 23:42 by JourneyLucky
[考研] 材料371求调剂 +9 鳄鱼? 2026-03-11 11/550 2026-03-13 22:53 by JourneyLucky
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
[考研] 310求调剂 +3 【上上签】 2026-03-11 3/150 2026-03-13 16:16 by JourneyLucky
[考研] 302求调剂 +6 负心者当诛 2026-03-11 6/300 2026-03-13 16:11 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考研] 274求调剂0856材料化工 +12 z2839474511 2026-03-11 13/650 2026-03-13 10:39 by peike
[考研] 工科0856专硕化学工程269能调剂吗 +10 我想读研11 2026-03-10 10/500 2026-03-13 10:14 by Yuyi.
[考博] 2026年博士申请 +3 QwQwQW10 2026-03-11 3/150 2026-03-12 17:58 by gxch43
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
信息提示
请填处理意见