计算scf不收敛，错误如下，求助。

3044305615 · 发表于 Post on 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不收敛，求助原因。

sobereva · 发表于 Post on 2014-12-25 15:14:57

你用了scf=qc？
把全部关键词列出来。
这才刚开始进行迭代，还远远没到不收敛呢

小范范1989 · 发表于 Post on 2014-12-25 15:19:46

我来学习，我也是出现收敛问题，给出输入文件看看。

小范范1989 · 发表于 Post on 2014-12-25 15:20:49

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

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

sobereva · 发表于 Post on 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没收敛也继续做下去，这样才会造成结果不准确，导致不该出现的虚频产生。

小范范1989 · 发表于 Post on 2014-12-25 16:10:12

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

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

3044305615 · 发表于 Post on 2014-12-25 22:35:47

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

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

3044305615 · 发表于 Post on 2014-12-25 22:36:13

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

关键词已经给出

3044305615 · 发表于 Post on 2014-12-25 22:39:17

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

同样学习了谢谢sob

sobereva · 发表于 Post on 2014-12-25 22:58:13

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

3044305615 · 发表于 Post on 2014-12-26 09:04:54

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

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

sobereva · 发表于 Post on 2014-12-26 09:07:35

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

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

3044305615 · 发表于 Post on 2014-12-26 09:57:46

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

谢谢sob

nkallwar · 发表于 Post on 2015-1-2 11:58:16

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

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

sobereva · 发表于 Post on 2015-1-2 12:50:17

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

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

浏览过的版块