计算化学公社

标题: 几何优化收敛与SCF收敛 [打印本页]

作者
Author: 最瘦二师兄 时间: 2016-5-10 15:59
标题: 几何优化收敛与SCF收敛
大家好，有个问题想请教大家，怎么区分是几何优化收敛和SCF收敛啊？优化收敛是不是在构型优化过程中遇到的啊即优化输出文件中的那四个收敛条件？而SCF收敛是指在单点能计算过程中遇到的啊？谢谢大家

作者
Author: 冰释之川 时间: 2016-5-10 19:25
简单的来说，SCF收敛是几何优化过程中【每一步】几何优化的一次循环（以RMSDP、MaxDP和DE为收敛限）；而几何优化其实相当于大循环（以Maximum Force 、RMS Force 、Maximum Displacement 和RMS Displacement 四个几何收敛判据为收敛限）

作者
Author: 最瘦二师兄 时间: 2016-5-16 09:07

冰释之川发表于 2016-5-10 19:25
简单的来说，SCF收敛是几何优化过程中【每一步】几何优化的一次循环（以RMSDP、MaxDP和DE为收敛限）；而几 ...

好的谢谢

作者
Author: 最瘦二师兄 时间: 2016-5-16 10:44

冰释之川发表于 2016-5-10 19:25
简单的来说，SCF收敛是几何优化过程中【每一步】几何优化的一次循环（以RMSDP、MaxDP和DE为收敛限）；而几 ...

那我接着请教一下  那我该如何从我的输出文件中判断它是scf不收敛还是优化不收敛呢这是我的输出文件中的最后他只告诉我收敛失败
Cycle 128  Pass 1  IDiag  1:
RMSU=  3.79D-06 CP:  1.00D+00  3.00D+00  1.34D+00  1.53D+00  1.55D+00
E= -3413.82812695882    Delta-E=    -0.000009770827 Rises=F Damp=F
DIIS: error= 1.41D-04 at cycle 125 NSaved=  20.
NSaved=20 IEnMin=20 EnMin= -3413.82812695882    IErMin=20 ErrMin= 1.41D-04
ErrMax= 1.41D-04 EMaxC= 1.00D-01 BMatC= 4.98D-06 BMatP= 4.55D-06
IDIUse=3 WtCom= 9.99D-01 WtEn= 1.41D-03
Matrix for removal  1 Erem= -3413.82695539264    Crem= 0.000D+00
Coeff-Com:  0.152D+00-0.116D+00 0.150D+00 0.688D+00-0.113D+00-0.289D+00
Coeff-Com: -0.350D-01-0.114D+01 0.360D+00-0.291D+00 0.294D+00-0.821D-01
Coeff-Com: -0.135D+00 0.390D+00 0.632D+00 0.181D+00 0.124D+00-0.499D-02
Coeff-Com:  0.809D-01 0.156D+00
Coeff-En: 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00
Coeff-En: 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00
Coeff-En: 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00
Coeff-En: 0.000D+00 0.100D+01
Coeff:    0.152D+00-0.116D+00 0.150D+00 0.687D+00-0.113D+00-0.289D+00
Coeff:    -0.350D-01-0.114D+01 0.359D+00-0.291D+00 0.294D+00-0.820D-01
Coeff:    -0.135D+00 0.390D+00 0.631D+00 0.181D+00 0.123D+00-0.498D-02
Coeff:    0.808D-01 0.157D+00
Gap=    0.007 Goal= None Shift= 0.000
RMSDP=1.22D-05 MaxDP=9.60D-04 DE=-9.77D-06 OVMax= 8.08D-03

>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(RPBE1PBE) =  -3413.82812696    A.U. after  129 cycles
         Convg  = 0.1223D-04          -V/T =  2.0183
KE= 3.352616297113D+03 PE=-2.545790568027D+04 EE= 9.907571306294D+03
Convergence failure -- run terminated.
Error termination via Lnk1e in /share/apps/gauss/g09/l502.exe at Mon May 16 02:49:09 2016.
Job cpu time:  2 days  6 hours 19 minutes 48.7 seconds.
File lengths (MBytes):  RWF= 524 Int=    0 D2E=    0 Chk=    57 Scr=    1

作者
Author: 冰释之川 时间: 2016-5-16 10:50

最瘦二师兄发表于 2016-5-16 10:44
那我接着请教一下那我该如何从我的输出文件中判断它是scf不收敛还是优化不收敛呢这是我的输出文件中的 ...

SCF不收敛

解决办法参考Sob老师的博文：http://www.sobereva.com/61

作者
Author: 最瘦二师兄 时间: 2016-5-16 14:04
本帖最后由最瘦二师兄于 2016-5-16 14:10 编辑

冰释之川发表于 2016-5-16 10:50
SCF不收敛

解决办法参考Sob老师的博文：http://www.sobereva.com/61

请问是根据这个RMSDP=1.22D-05 MaxDP=9.60D-04 DE=-9.77D-06 OVMax= 8.08D-03

>>>>>>>>>> Convergence criterion not met.因为RMSDP，MaxDP，DE是scf的收敛限而 Convergence criterion not met，从而判断出来它是scf不收敛吗

作者
Author: sobereva 时间: 2016-5-16 16:11

最瘦二师兄发表于 2016-5-16 14:04
请问是根据这个RMSDP=1.22D-05 MaxDP=9.60D-04 DE=-9.77D-06 OVMax= 8.08D-03

>>>>>>>>>> ...

这提示一看就是SCF步骤的输出，稍有使用经验就能判断。你说的那种判断也对。

作者
Author: 最瘦二师兄 时间: 2016-5-16 19:44

sobereva 发表于 2016-5-16 16:11
这提示一看就是SCF步骤的输出，稍有使用经验就能判断。你说的那种判断也对。

明白了谢谢老师

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