CAM-B3LYP下TiO2修饰染料分子能量计算

ca0yan9

请问一下使用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。

ca0yan9

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

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

sobereva

解决SCF不收敛问题的方法
http://sobereva.com/61

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

ca0yan9

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

谢谢Sob老师指点。