请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 6305|回复: 32

[Gaussian/gview] 计算scf不收敛,错误如下,求助。

[复制链接]

48

帖子

0

威望

105

eV
积分
153

Level 3 能力者

发表于 2014-12-25 14:29:20 | 显示全部楼层 |阅读模式
本帖最后由 3044305615 于 2014-12-25 22:38 编辑

#p opt pbepbe/gen pseudo=read scf=xqc

Cycle  64  Pass 1  IDiag  1:
Density matrix breaks symmetry, PCut= 1.00D-04
Density has only Abelian symmetry.
RMSU=  1.75D-04    CP:  9.98D-01  6.17D-01  1.05D+00
E= -2763.13442798851     Delta-E=       -0.040685409853 Rises=F Damp=F
DIIS: error= 1.31D-02 at cycle  64 NSaved=  20.
NSaved=20 IEnMin=10 EnMin= -2763.17628199912     IErMin= 5 ErrMin= 1.07D-02
ErrMax= 1.31D-02  0.00D+00 EMaxC= 1.00D-01 BMatC= 1.01D-01 BMatP= 4.46D-02
IDIUse=2 WtCom= 0.00D+00 WtEn= 1.00D+00
EnCoef did   100 forward-backward iterations
Coeff-En:   0.000D+00 0.204D-02 0.500D-02 0.771D-01 0.284D+00 0.127D-02
Coeff-En:   0.353D-03 0.285D+00 0.134D+00 0.000D+00 0.100D-03 0.145D+00
Coeff-En:   0.567D-04 0.157D-03 0.105D-02 0.000D+00 0.109D-02 0.385D-02
Coeff-En:   0.000D+00 0.595D-01
Coeff:      0.000D+00 0.204D-02 0.500D-02 0.771D-01 0.284D+00 0.127D-02
Coeff:      0.353D-03 0.285D+00 0.134D+00 0.000D+00 0.100D-03 0.145D+00
Coeff:      0.567D-04 0.157D-03 0.105D-02 0.000D+00 0.109D-02 0.385D-02
Coeff:      0.000D+00 0.595D-01
Gap=     0.000 Goal=   None    Shift=    0.000
RMSDP=8.62D-04 MaxDP=7.88D-02 DE=-4.07D-02 OVMax= 8.76D-01
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(RPBE-PBE) =  -2763.13442799     A.U. after   65 cycles
            NFock= 64  Conv=0.86D-03     -V/T= 2.1381
KE= 2.427854839894D+03 PE=-2.662153696979D+04 EE= 1.125126781111D+04
Leave Link  502 at Thu Dec 25 12:24:57 2014, MaxMem=  3145728000 cpu:     38628.9
(Enter /data3/home/wangzg/g09/l508.exe)
Direct SCF.
Tight linear equation convergence will be used.
RPBE-PBE wavefunction.
Quadratic Convergence SCF Method.
Line search only if initial step raises energy.
Integral symmetry usage will be decided dynamically.
ReadGW:  IGet=1 IStart=         402 Next=         402 LGW=           0.
Fock matrices will be formed incrementally for  20 cycles.
Density matrix breaks symmetry, PCut= 1.00D-04
Density has only Abelian symmetry.
FoFJK:  IHMeth= 1 ICntrl=     500 DoSepK=F KAlg= 0 I1Cent=   0 FoldK=F
IRaf= 760000000 NMat=   1 IRICut=       1 DoRegI=T DoRafI=F ISym2E= 1.
FoFCou: FMM=F IPFlag=           0 FMFlag=      100000 FMFlg1=        2000
         NFxFlg=           0 DoJE=T BraDBF=F KetDBF=F FulRan=T
         wScrn=  0.000000 ICntrl=     500 IOpCl=  0 I1Cent=           0 NGrid=           0
         NMat0=    1 NMatS0=      1 NMatT0=    0 NMatD0=    1 NMtDS0=    0 NMtDT0=    0
Petite list used in FoFCou.
Defaulting to unpruned grid for atomic number  92.
Defaulting to unpruned grid for atomic number  92.
Defaulting to unpruned grid for atomic number  92.
Defaulting to unpruned grid for atomic number  92.
Defaulting to unpruned grid for atomic number  92.
Defaulting to unpruned grid for atomic number  92.
Defaulting to unpruned grid for atomic number  92.
Defaulting to unpruned grid for atomic number  92.
Defaulting to unpruned grid for atomic number  92.
Defaulting to unpruned grid for atomic number  92.
Defaulting to unpruned grid for atomic number  92.
Defaulting to unpruned grid for atomic number  92.
Defaulting to unpruned grid for atomic number  92.
Defaulting to unpruned grid for atomic number  92.
Defaulting to unpruned grid for atomic number  92.
Defaulting to unpruned grid for atomic number  92.
Iteration   1 EE= -2763.18685395436                                    Grad=2.378D-01 Max rot=    1.9217 deg.
Gradient too large for Newton-Raphson -- use scaled steepest descent instead.
Alpha Gap=    -0.027 HL deriv= 1.58D-02
     ILin= 1 X=0.000D+00 Y=-2.763186853954D+03 DE= 0.00D+00 F= -1.76D-01
Density matrix breaks symmetry, PCut= 1.00D-04
Rerun with SCF=IntRep.


没见过这种类型的scf不收敛,求助原因。

1万

帖子

25

威望

1万

eV
积分
34891

管理员

公社社长

发表于 2014-12-25 15:14:57 | 显示全部楼层
你用了scf=qc?
把全部关键词列出来。
这才刚开始进行迭代,还远远没到不收敛呢
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

1111

帖子

0

威望

3120

eV
积分
4231

Level 6 (一方通行)

发表于 2014-12-25 15:19:46 | 显示全部楼层
我来学习,我也是出现收敛问题,给出输入文件看看。
有梦想,才有实现。奋斗吧,骄傲的少年。

1111

帖子

0

威望

3120

eV
积分
4231

Level 6 (一方通行)

发表于 2014-12-25 15:20:49 | 显示全部楼层
sobereva 发表于 2014-12-25 15:14
你用了scf=qc?
把全部关键词列出来。
这才刚开始进行迭代,还远远没到不收敛呢

sob,问个问题,是不是scf=qc强制收敛,后面在计算freq的时候,是不是很容易出现虚频?谢谢sob
有梦想,才有实现。奋斗吧,骄傲的少年。

1万

帖子

25

威望

1万

eV
积分
34891

管理员

公社社长

发表于 2014-12-25 15:45:16 | 显示全部楼层
小范范1989 发表于 2014-12-25 15:20
sob,问个问题,是不是scf=qc强制收敛,后面在计算freq的时候,是不是很容易出现虚频?谢谢sob

scf=qc并不能达到强制收敛的效果,达不到收敛限而超过步数上限的时候照样会报错停止。qc只不过是比默认的DIIS收敛所需步数少(但总耗时更长)、原理更为严格,但不能保证一定能收敛。
只要qc下收敛了,那么之后的opt freq也没问题,如果没收敛,就会报错终止。所以不会造成你说的问题。
出现你说的问题,应该是指用IOp(5/13=1)使得程序即便遇到SCF没收敛也继续做下去,这样才会造成结果不准确,导致不该出现的虚频产生。
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

1111

帖子

0

威望

3120

eV
积分
4231

Level 6 (一方通行)

发表于 2014-12-25 16:10:12 | 显示全部楼层
sobereva 发表于 2014-12-25 15:45
scf=qc并不能达到强制收敛的效果,达不到收敛限而超过步数上限的时候照样会报错停止。qc只不过是比默认的 ...

偶,这样啊,谢谢sob。
圣诞节快乐
有梦想,才有实现。奋斗吧,骄傲的少年。

48

帖子

0

威望

105

eV
积分
153

Level 3 能力者

 楼主| 发表于 2014-12-25 22:35:47 | 显示全部楼层
sobereva 发表于 2014-12-25 15:14
你用了scf=qc?
把全部关键词列出来。
这才刚开始进行迭代,还远远没到不收敛呢

关键词 :#p opt pbepbe/gen pseudo=read scf=xqc

48

帖子

0

威望

105

eV
积分
153

Level 3 能力者

 楼主| 发表于 2014-12-25 22:36:13 | 显示全部楼层
小范范1989 发表于 2014-12-25 15:19
我来学习,我也是出现收敛问题,给出输入文件看看。

关键词已经给出

48

帖子

0

威望

105

eV
积分
153

Level 3 能力者

 楼主| 发表于 2014-12-25 22:39:17 | 显示全部楼层
sobereva 发表于 2014-12-25 15:45
scf=qc并不能达到强制收敛的效果,达不到收敛限而超过步数上限的时候照样会报错停止。qc只不过是比默认的 ...

同样学习了 谢谢sob

1万

帖子

25

威望

1万

eV
积分
34891

管理员

公社社长

发表于 2014-12-25 22:58:13 | 显示全部楼层
就是由于你用了xqc,到了64步时普通方式迭代没收敛,自动切换到了qc方式收敛。qc方式每一步很慢,得耐心等着,从你的输出信息看才刚做了一步qc步,离达到SCF步数上限还早着。
通常xqc/qc不是推荐的解决SCF不收敛的做法,详见此贴
解决SCF不收敛问题的方法 http://sobereva.com/61
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

48

帖子

0

威望

105

eV
积分
153

Level 3 能力者

 楼主| 发表于 2014-12-26 09:04:54 | 显示全部楼层
sobereva 发表于 2014-12-25 22:58
就是由于你用了xqc,到了64步时普通方式迭代没收敛,自动切换到了qc方式收敛。qc方式每一步很慢,得耐心等 ...

sob  老师  输出文件最后一句有个Rerun with SCF=IntRep       如果我加上scf=intrep 作用在于什么呢?对最后促进收敛you什么帮助么?  谢谢

1万

帖子

25

威望

1万

eV
积分
34891

管理员

公社社长

发表于 2014-12-26 09:07:35 | 显示全部楼层
3044305615 发表于 2014-12-26 09:04
sob  老师  输出文件最后一句有个Rerun with SCF=IntRep       如果我加上scf=intrep 作用在于什么呢?对 ...

不用管它,也没必要手动加它,对帮助收敛没直接用处
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

48

帖子

0

威望

105

eV
积分
153

Level 3 能力者

 楼主| 发表于 2014-12-26 09:57:46 | 显示全部楼层
sobereva 发表于 2014-12-26 09:07
不用管它,也没必要手动加它,对帮助收敛没直接用处

谢谢sob

66

帖子

0

威望

786

eV
积分
852

Level 4 (黑子)

发表于 2015-1-2 11:58:16 | 显示全部楼层
sobereva 发表于 2014-12-25 22:58
就是由于你用了xqc,到了64步时普通方式迭代没收敛,自动切换到了qc方式收敛。qc方式每一步很慢,得耐心等 ...

    请教sob老师,要是在几何优化用了scf=conver=6 关键词,得到的结果还需要去掉这个关键词重新计算一下验证么?
   

1万

帖子

25

威望

1万

eV
积分
34891

管理员

公社社长

发表于 2015-1-2 12:50:17 | 显示全部楼层
nkallwar 发表于 2015-1-2 11:58
请教sob老师,要是在几何优化用了scf=conver=6 关键词,得到的结果还需要去掉这个关键词重新计算一下 ...

这样得到的几何结构可能不准确(虽然一般差得不会太多),为了准确,最好还是去掉之后继续优化一下
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2018-10-19 18:13 , Processed in 0.197412 second(s), 24 queries .

快速回复 返回顶部 返回列表