24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2644  |  回复: 20

LinaInverse

木虫 (小有名气)

[求助] 不进行optimize就直接报错收敛失败,求解决

rt,小可做计算有段时间也有些心得了,体系很大,124个原子,元素为C,H,B,N和重元素Ir,方法b3lyp,基组6-31g(d)+lanl2dz,-3电荷,多重度3
试过opt=tight,没用,iop(5/13)=1没试但也不想试
自洽场可以收敛,但是一圈优化都没跑就直接报错收敛失败,求大侠指点
以下是报错语句:
Symmetry not used in FoFCou.
Initial guess <Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 1.0000 <S**2>= 2.0000 S= 1.0000
Defaulting to unpruned grid for atomic number  77.
Defaulting to unpruned grid for atomic number  77.
Defaulting to unpruned grid for atomic number  77.
Defaulting to unpruned grid for atomic number  77.
Defaulting to unpruned grid for atomic number  77.
Defaulting to unpruned grid for atomic number  77.
Defaulting to unpruned grid for atomic number  77.
Defaulting to unpruned grid for atomic number  77.
EnCoef did   100 forward-backward iterations
EnCoef did     1 forward-backward iterations
EnCoef did     5 forward-backward iterations
EnCoef did     2 forward-backward iterations
Restarting incremental Fock formation.
Matrix for removal  4 Erem= -2448.69634785012     Crem= 0.000D+00
Matrix for removal  3 Erem= -2454.71182296237     Crem= 0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UB3LYP) =  -2458.97111905     A.U. after  129 cycles
             Convg  =    0.4215D-04             -V/T =  2.0365
<Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 1.0000 <S**2>= 2.0085 S= 1.0028
<L.S>= 0.000000000000E+00
Annihilation of the first spin contaminant:
S**2 before annihilation     2.0085,   after     2.0000
Convergence failure -- run terminated.
Error termination via Lnk1e in /gpfssan1/home/hyan//g09/l502.exe at Wed Jul  3 19:29:12 2013.
回复此楼
秀逗魔导师
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

Jasminer

铁杆木虫 (著名写手)


gkf高: 金币+1 2013-07-04 16:13:40
这个就是自洽场没有收敛。
2楼2013-07-04 15:41:32
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

LinaInverse

木虫 (小有名气)

引用回帖:
2楼: Originally posted by Jasminer at 2013-07-04 15:41:32
这个就是自洽场没有收敛。

自洽场不收敛会直接写出SCF Convergence failure的吧
而且这个上面已经写了SCF Done了
秀逗魔导师
3楼2013-07-04 17:07:08
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

t13340033021

金虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
引用回帖:
3楼: Originally posted by LinaInverse at 2013-07-04 17:07:08
自洽场不收敛会直接写出SCF Convergence failure的吧
而且这个上面已经写了SCF Done了...

试一下calcfc关键词
甲烷无敌
4楼2013-07-04 17:13:27
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

LinaInverse

木虫 (小有名气)

引用回帖:
4楼: Originally posted by t13340033021 at 2013-07-04 17:13:27
试一下calcfc关键词...

一圈优化都没跑,calcfc关键词能发挥作用吗?
秀逗魔导师
5楼2013-07-04 19:17:09
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

t13340033021

金虫 (正式写手)

【答案】应助回帖

★ ★
gkf高: 金币+2 2013-07-05 11:32:01
引用回帖:
5楼: Originally posted by LinaInverse at 2013-07-04 19:17:09
一圈优化都没跑,calcfc关键词能发挥作用吗?...

先解决SCF不收敛
SCF=QC
调整下初始结构
有重金属请使用ECP基组
全电子通常不容易收敛
甲烷无敌
6楼2013-07-04 19:21:09
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

LinaInverse

木虫 (小有名气)

引用回帖:
6楼: Originally posted by t13340033021 at 2013-07-04 19:21:09
先解决SCF不收敛
SCF=QC
调整下初始结构
有重金属请使用ECP基组
全电子通常不容易收敛...

主题贴上写了lanl2dz赝势基组了,你猜是用在哪个原子上?
而且问题不是SCF不收敛,如果是scf不收敛报错语句会不一样,最典型的就是经过多少多少圈优化后SCF不收敛这个句子,可是log里没有
个人经验scf=qc对付不收敛得不偿失而且也不是每次都会奏效,iop(5/13=1)对付自洽场不收敛有用但实质是跳过自洽场检验,没解决根本问题

初始结构单线态条件下优化过,正在算频率看会不会有虚频,但我们组做过很多类似配合物,构型问题恐怕不大
秀逗魔导师
7楼2013-07-04 22:48:02
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

t13340033021

金虫 (正式写手)

【答案】应助回帖

>>>>>>>>>> Convergence criterion not met.
这句话楼主看不懂吗?
帮你出谋划策的人可能会忽视一些信息,但是你自己得上心些
甲烷无敌
8楼2013-07-05 09:01:15
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

LinaInverse

木虫 (小有名气)

引用回帖:
8楼: Originally posted by t13340033021 at 2013-07-05 09:01:15
>>>>>>>>>> Convergence criterion not met.
这句话楼主看不懂吗?
帮你出谋划策的人可能会忽视一些信息,但是你自己得上心些

收敛失败和自洽场不收敛是两回事
报错语句我都遇到过,明显不同,这个不是自洽场不收敛的报错语句
秀逗魔导师
9楼2013-07-05 10:23:15
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

yypcnbe

木虫 (小有名气)

【答案】应助回帖

★ ★
gkf高: 金币+2 2013-07-05 13:50:56
引用回帖:
3楼: Originally posted by LinaInverse at 2013-07-04 17:07:08
自洽场不收敛会直接写出SCF Convergence failure的吧
而且这个上面已经写了SCF Done了...

你确定是收敛了,这句话的意思只是表示129次后SCF迭代完成了,但没有达到收敛标准,现在才达到scfcon=4,你设置的肯定不止4这么低吧,试试scfcyc=500,或者scf=qc。
10楼2013-07-05 13:37:43
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 LinaInverse 的主题更新
信息提示
请填处理意见