24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 5150  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 284求调剂 +6 徐同学_001 2026-04-04 12/600 2026-04-05 13:14 by 徐同学_001
[考研] 求调剂,一志愿郑州大学材料与化工专硕,英二数二342分,求老师收留 +18 v12abo 2026-04-02 20/1000 2026-04-05 11:37 by a8144223
[考研] 346分的生物与医药08600求调剂 +4 常雨阳上岸 2026-04-05 5/250 2026-04-05 10:26 by 唐沐儿
[考研] 材料与化工371求调剂 +8 陪琳看海 2026-04-04 8/400 2026-04-05 09:10 by guoweigw
[考研] 材料调剂 +9 革微桂 2026-04-04 9/450 2026-04-05 08:27 by 544594351
[考研] 085602调剂 初试总分335 +12 19123253302 2026-04-04 12/600 2026-04-05 08:08 by 544594351
[考研] 283分求调剂 +7 小聂爱学习 2026-04-03 7/350 2026-04-04 21:51 by hemengdong
[考研] 315求调剂 +13 小羊小羊_ 2026-04-02 14/700 2026-04-04 20:30 by 蓝云思雨
[考研] 一志愿C9的化学工程(085602) 340分,感觉校内调剂无望,求调剂 +9 万事宜臻 2026-04-04 9/450 2026-04-04 11:49 by 啵啵啵0119
[考研] 一志愿北京科技大学材料工程085601,求调剂 +17 cdyw 2026-04-02 18/900 2026-04-04 11:14 by w_xuqing
[考研] 297求调剂 +11 ljy20040718! 2026-04-03 13/650 2026-04-04 09:23 by 来看流星雨10
[考研] 一志愿双非085502,267分,过四级求调剂 +3 再忙也要吃饭啊 2026-04-03 3/150 2026-04-04 05:03 by gswylq
[考研] 一志愿南昌大学324求调剂 +13 hanamiko 2026-04-01 13/650 2026-04-03 18:30 by ls刘帅
[考研] 081200-11408-276学硕求调剂 +6 崔wj 2026-04-02 6/300 2026-04-03 10:19 by 蓝云思雨
[考研] 309求调剂 +14 呆菇不是戴夫 2026-04-02 14/700 2026-04-03 09:42 by 蓝云思雨
[考研] 一志愿山东大学,085600,344 +7 魏子per 2026-04-02 8/400 2026-04-02 21:12 by 百灵童888
[考研] 一志愿郑大材料工程290求调剂 +20 Youth_ 2026-03-30 20/1000 2026-04-02 14:48 by 5896
[考研] 一志愿北京科技大学材料学硕328分求调剂 +6 1段时间 2026-03-31 7/350 2026-04-02 13:57 by 3041
[考研] 261求B区调剂 +5 明仔· 2026-04-01 7/350 2026-04-02 11:17 by 邹尉尉
[考研] 一志愿中海洋材料357 +4 麦恩莉. 2026-03-30 4/200 2026-03-31 14:35 by 记事本2026
信息提示
请填处理意见