计算化学公社

标题: 计算scf不收敛，错误如下，求助。 [打印本页]

作者
Author: 3044305615 时间: 2014-12-25 14:29
标题: 计算scf不收敛，错误如下，求助。
本帖最后由 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不收敛，求助原因。

作者
Author: sobereva 时间: 2014-12-25 15:14
你用了scf=qc？
把全部关键词列出来。
这才刚开始进行迭代，还远远没到不收敛呢

作者
Author: 小范范1989 时间: 2014-12-25 15:19
我来学习，我也是出现收敛问题，给出输入文件看看。

作者
Author: 小范范1989 时间: 2014-12-25 15:20

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

sob，问个问题，是不是scf=qc强制收敛，后面在计算freq的时候，是不是很容易出现虚频？谢谢sob

作者
Author: sobereva 时间: 2014-12-25 15:45

小范范1989 发表于 2014-12-25 15:20
sob，问个问题，是不是scf=qc强制收敛，后面在计算freq的时候，是不是很容易出现虚频？谢谢sob

scf=qc并不能达到强制收敛的效果，达不到收敛限而超过步数上限的时候照样会报错停止。qc只不过是比默认的DIIS收敛所需步数少（但总耗时更长）、原理更为严格，但不能保证一定能收敛。
只要qc下收敛了，那么之后的opt freq也没问题，如果没收敛，就会报错终止。所以不会造成你说的问题。
出现你说的问题，应该是指用IOp(5/13=1)使得程序即便遇到SCF没收敛也继续做下去，这样才会造成结果不准确，导致不该出现的虚频产生。

作者
Author: 小范范1989 时间: 2014-12-25 16:10

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

偶，这样啊，谢谢sob。
圣诞节快乐

作者
Author: 3044305615 时间: 2014-12-25 22:35

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

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

作者
Author: 3044305615 时间: 2014-12-25 22:36

小范范1989 发表于 2014-12-25 15:19
我来学习，我也是出现收敛问题，给出输入文件看看。

关键词已经给出

作者
Author: 3044305615 时间: 2014-12-25 22:39

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

同样学习了谢谢sob

作者
Author: sobereva 时间: 2014-12-25 22:58
就是由于你用了xqc，到了64步时普通方式迭代没收敛，自动切换到了qc方式收敛。qc方式每一步很慢，得耐心等着，从你的输出信息看才刚做了一步qc步，离达到SCF步数上限还早着。
通常xqc/qc不是推荐的解决SCF不收敛的做法，详见此贴
解决SCF不收敛问题的方法 http://sobereva.com/61

作者
Author: 3044305615 时间: 2014-12-26 09:04

sobereva 发表于 2014-12-25 22:58
就是由于你用了xqc，到了64步时普通方式迭代没收敛，自动切换到了qc方式收敛。qc方式每一步很慢，得耐心等 ...

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

作者
Author: sobereva 时间: 2014-12-26 09:07

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

不用管它，也没必要手动加它，对帮助收敛没直接用处

作者
Author: 3044305615 时间: 2014-12-26 09:57

sobereva 发表于 2014-12-26 09:07
不用管它，也没必要手动加它，对帮助收敛没直接用处

谢谢sob

作者
Author: nkallwar 时间: 2015-1-2 11:58

sobereva 发表于 2014-12-25 22:58
就是由于你用了xqc，到了64步时普通方式迭代没收敛，自动切换到了qc方式收敛。qc方式每一步很慢，得耐心等 ...

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

作者
Author: sobereva 时间: 2015-1-2 12:50

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

这样得到的几何结构可能不准确（虽然一般差得不会太多），为了准确，最好还是去掉之后继续优化一下

作者
Author: lao7 时间: 2015-1-2 21:29
用这句话试试，对我来说，如果结构没有画错，就是神器

作者
Author: lao7 时间: 2015-1-2 21:31
maxcycle=200,xqc
试试这个，好用有时候

作者
Author: lao7 时间: 2015-1-3 08:55
本帖最后由 lao7 于 2015-1-3 08:59 编辑

上面，手机不太好回复。发了两次

作者
Author: nkallwar 时间: 2015-1-5 09:00

sobereva 发表于 2015-1-2 12:50
这样得到的几何结构可能不准确（虽然一般差得不会太多），为了准确，最好还是去掉之后继续优化一下

谢谢sob，不过我试了一个150个原子的双核Cr体系，用的方法和基组是M06/6-31G(d)/Lanl08(f) 做优化，发现用 scf=conver=6 完成了优化，在得到结构基础上，去掉scf=conver=6 的关键词，结果又不收敛了。（第一个link 502 过程SCF收敛了，但两个Displacement判据不收敛，优化几圈后一个link502 SCF收敛失败）。
于是我把结构做一个频率分析，准备看看存在的虚频是什么样的，好根据振动方向来微调结构，结果发现频率分析中没有虚频，这似乎说明得到的结构已经是一个稳定结构了，那这个几何结构可以使用么？（freq 作业中没有使用放宽scf收敛标准）

作者
Author: sobereva 时间: 2015-1-5 09:30
没有虚频不代表是极小点结构。收敛精度不高的情况依然可以没有虚频。只要结构处在极小点的二次元区域内就都是没有虚频的。

是否能用取决于你对结构要求的精度。不要求很准的话，勉强能用。

不收敛有很多办法解决，几何优化不收敛见《量子化学计算中帮助几何优化收敛的常用方法》（http://hi.baidu.com/sobereva/item/596c242087e820cedcf69a9d）
SCF不收敛见《解决SCF不收敛问题的方法》（http://hi.baidu.com/sobereva/item/68b8d2866c0223ebe596e0dc）

作者
Author: 3044305615 时间: 2015-1-5 20:42

lao7 发表于 2015-1-2 21:31
maxcycle=200,xqc
试试这个，好用有时候

好的谢谢我其实没太明白xqc的作用

作者
Author: nkallwar 时间: 2015-1-6 10:34

sobereva 发表于 2015-1-5 09:30
没有虚频不代表是极小点结构。收敛精度不高的情况依然可以没有虚频。只要结构处在极小点的二次元区域内就都 ...

多谢sob回复，我有一个结构做几何优化的时候，位移判据始终不收敛，力判据则是早早就收敛了，看起来是进入了一个势能平坦区，结构整体刚性还挺强的，没有大的柔链结构。那这种情况是缩小步长在周围找找还是怎么想法跳出这个势能平坦区呢？

作者
Author: sobereva 时间: 2015-1-6 10:41

nkallwar 发表于 2015-1-6 10:34
多谢sob回复，我有一个结构做几何优化的时候，位移判据始终不收敛，力判据则是早早就收敛了，看起来是进 ...

平坦区里面也有极小点，不过这个极小点可能是不显著的。
先试试gdiis或缩小步长上限

作者
Author: 3044305615 时间: 2015-1-15 10:41
本帖最后由 3044305615 于 2015-1-15 10:53 编辑

sobereva 发表于 2014-12-25 22:58
就是由于你用了xqc，到了64步时普通方式迭代没收敛，自动切换到了qc方式收敛。qc方式每一步很慢，得耐心等 ...

sob 老师  如果两个原子距离过近（限域）而导致排斥能高计算出的能量已经超出高斯给的范围了而导致的不收敛该怎么办？
输出文件太长了不知道哪部分该粘贴

命令行用了iop（5/13=1） empiricaldispersion=gd3
Variable    Old X -DE/DX Delta X Delta X Delta X    New X
                              (Linear) (Quad) (Total)
R1       2.75816    NaN 0.00000    NaN 0.00000 2.75816
R2       4.96577    NaN 0.00000    NaN 0.00000 4.96577
R3       4.10904    NaN 0.00000    NaN 0.00000 4.10904
R4       4.96577    NaN 0.00000    NaN 0.00000 4.96577
R5       4.96577    NaN 0.00000    NaN 0.00000 4.96577
R6       4.31456    NaN 0.00000    NaN 0.00000 4.31456
R7       4.86419    NaN 0.00000    NaN 0.00000 4.86419
R8       4.31456    NaN 0.00000    NaN 0.00000 4.31456
R9       4.86418    NaN 0.00000    NaN 0.00000 4.86418
R10       4.31455    NaN 0.00000    NaN 0.00000 4.31455
R11       4.31454    NaN 0.00000    NaN 0.00000 4.31454
R12       4.31454    NaN 0.00000    NaN 0.00000 4.31454
R13       4.86418    NaN 0.00000    NaN 0.00000 4.86418
R14       4.31455    NaN 0.00000    NaN 0.00000 4.31455
R15       4.77969    NaN 0.00000    NaN 0.00000 4.77969
R16       4.22340    NaN 0.00000    NaN 0.00000 4.22340
R17       4.77969    NaN 0.00000    NaN 0.00000 4.77969
R18       4.22340    NaN 0.00000    NaN 0.00000 4.22340
还有这个表示
RMSDP=7.07D-01 MaxDP=8.51D+01 DE=    NaN OVMax= 1.00D+00

以及文件末尾
RFO could not converge Lambda in  999 iterations.
Linear search skipped for unknown reason.

作者
Author: sobereva 时间: 2015-1-15 15:19

3044305615 发表于 2015-1-15 10:41
sob 老师如果两个原子距离过近（限域）而导致排斥能高计算出的能量已经超出高斯给的范围了而导致 ...

绝对不要用5/13=1，SCF都没正常收敛，之后计算的受力可能会很邪乎
先用gview观察初始结构，把严重不合理的诸如距离过近的地方调整一下（比如用扫把型图标）

作者
Author: 3044305615 时间: 2015-1-16 16:20

sobereva 发表于 2015-1-15 15:19
绝对不要用5/13=1，SCF都没正常收敛，之后计算的受力可能会很邪乎
先用gview观察初始结构，把严重不合理 ...

sob 老师我发现我的错误是log刚开始的grimme 色散的值就是超出高斯给的范围了，为什么呢？是不是这个值没有，才导致后来的不收敛啊，您能帮忙解答一下么？但是d3不是针对1-94号元素么？

作者
Author: sobereva 时间: 2015-1-16 16:23

3044305615 发表于 2015-1-16 16:20
sob 老师我发现我的错误是log刚开始的grimme 色散的值就是超出高斯给的范围了，为什么呢？是不是这个值 ...

去掉D3再试自然就明白了。如果能正常优化完毕（此时肯定已经消除了原子间过近的接触），再重新用D3继续优化。

作者
Author: huangtf 时间: 2016-10-11 15:54
sob老师，您好！使用scf=xqc进行几何构型优化，非得把关键词stable写到输入文件(gjf)中吗？

谢谢老师！

作者
Author: sobereva 时间: 2016-10-11 16:55

huangtf 发表于 2016-10-11 15:54
sob老师，您好！使用scf=xqc进行几何构型优化，非得把关键词stable写到输入文件(gjf)中吗？

谢谢老师！

没必要

作者
Author: 彭索萍 时间: 2017-7-21 09:51
SCF=IntRep请问sob老师，这个关键词的作用？

作者
Author: sobereva 时间: 2017-7-21 22:13

彭索萍发表于 2017-7-21 09:51
SCF=IntRep请问sob老师，这个关键词的作用？

一般没什么用，不需要记住

作者
Author: 刘小英雄steven 时间: 2018-10-12 08:15
本帖最后由刘小英雄steven 于 2018-10-12 08:57 编辑

老师，柔性扫描键长，随着键长的增大，高斯会报错吗

作者
Author: sobereva 时间: 2018-10-12 09:46

刘小英雄steven 发表于 2018-10-12 08:15
老师，柔性扫描键长，随着键长的增大，高斯会报错吗

没有报错的必然性

欢迎光临计算化学公社 (http://bbs.keinsci.com/)