计算化学公社

标题: CAM-B3LYP下TiO2修饰染料分子能量计算 [打印本页]

作者
Author: ca0yan9 时间: 2017-11-2 11:08
标题: CAM-B3LYP下TiO2修饰染料分子能量计算
本帖最后由 ca0yan9 于 2017-11-2 13:01 编辑

请问一下使用CAM-B3LYP计算TiO2团簇表面修饰染料分子的能量，总是不收敛的情况如何处理。

已经将所有的原子进行固定（模型事先使用VASP优化出团簇结构）。是不是收敛过程的参数需要修改和调整？

SCF中60次左右的结果RMSDP=3.25D-02 MaxDP=5.23D+00 DE=-6.46D-03 OVMax= 4.96D-02，感觉离收敛还比较远，这个情况怎么处理？
Cycle  65  Pass 1  IDiag  1:
RMSU=  2.90D-02 CP:  1.00D+00  3.00D+00  3.00D+00  3.00D+00
E= -15278.4719862288    Delta-E=    -0.006456723422 Rises=F Damp=F
DIIS: error= 7.06D-04 at cycle  65 NSaved=  20.
NSaved=20 IEnMin=20 EnMin= -15278.4719862288    IErMin=15 ErrMin= 2.26D-04
ErrMax= 7.06D-04  0.00D+00 EMaxC= 1.00D-01 BMatC= 2.80D-03 BMatP= 1.54D-03
IDIUse=3 WtCom= 9.98D-01 WtEn= 2.37D-03
Coeff-Com: -0.770D-02 0.302D-02 0.283D-01 0.249D-01-0.217D-01-0.772D-01
Coeff-Com: -0.950D-02 0.176D-01-0.148D+00 0.145D+00-0.713D-01 0.489D+00
Coeff-Com: -0.310D+00 0.308D+00-0.169D+01-0.256D+01-0.254D+01 0.556D+01
Coeff-Com:  0.129D+01 0.579D+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.768D-02 0.301D-02 0.283D-01 0.248D-01-0.217D-01-0.770D-01
Coeff:    -0.948D-02 0.176D-01-0.147D+00 0.144D+00-0.711D-01 0.488D+00
Coeff:    -0.309D+00 0.307D+00-0.169D+01-0.256D+01-0.254D+01 0.555D+01
Coeff:    0.128D+01 0.580D+00
Gap=    0.061 Goal= None Shift= 0.000
Gap=    0.056 Goal= None Shift= 0.000
RMSDP=3.25D-02 MaxDP=5.23D+00 DE=-6.46D-03 OVMax= 4.96D-02。

作者
Author: ca0yan9 时间: 2017-11-2 22:06

尝试了几次，都没法收敛。

这种原子比较多的体系收敛有什么方法吗？

作者
Author: sobereva 时间: 2017-11-3 17:27
解决SCF不收敛问题的方法
http://sobereva.com/61

依次尝试。和原子数目多少没有直接关系。

作者
Author: ca0yan9 时间: 2017-11-3 21:45

sobereva 发表于 2017-11-3 17:27
解决SCF不收敛问题的方法
http://sobereva.com/61

谢谢Sob老师指点。

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