24小时热门版块排行榜    

查看: 5811  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 建议基金发布提前给出明确的时间点 +11 kulium 2026-08-21 14/700 2026-08-24 02:12 by 流流伤
[基金申请] 2026国自然函评费到账 +16 羊腰板 2026-08-21 17/850 2026-08-23 23:02 by tianxiaochun
[基金申请] 让我中一个面上吧! +12 大萍1987 2026-08-20 15/750 2026-08-23 20:44 by 大萍1987
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见