计算化学公社

标题: 求助：MgO(100)单点计算不收敛 [打印本页]

作者
Author: 风飞 时间: 2021-5-17 21:09
标题: 求助：MgO(100)单点计算不收敛
本帖最后由风飞于 2021-5-17 23:32 编辑

老师，您好，

         我首先采用QE对MgO（100）表面吸附甲苯的体系进行了结构优化，将得到的结构保存成cif文件，使用Multiwfn创建cp2K的输入文件进行单点计算，

         然后进行IGMH分析和电荷分析，但是在进行单点计算的过程中不收敛，请问可能是什么原因呢？我对应的输入文件如下：输出结果文件如下，：

主要参数设置如下，：
方法：PBE
基组：pob-TZVP
色散矫正：DFT-D3
截断能： CUTOFF 350
         REL_CUTOFF 50
SCF设置： EPS_SCF 5.0E-06
            MAX_SCF 128

作者
Author: sobereva 时间: 2021-5-17 21:54
实测开OT可以收敛

另外，这个体系应当在Multiwfn创建输入文件的时候把周期性设为XY

作者
Author: 风飞 时间: 2021-5-18 00:31
本帖最后由风飞于 2021-5-18 00:32 编辑

sobereva 发表于 2021-5-17 21:54
实测开OT可以收敛

另外，这个体系应当在Multiwfn创建输入文件的时候把周期性设为XY

老师，您好，按照您方法，对于MgO(100)体系，可以收敛了，但是对于其他体系如Al2O3（110）-T （T是指甲苯））也出现了不收敛的情况，请问对于不收敛，我们一般可以从哪些方面去调整呢？
输入文件如下：
(, 下载次数 Times of downloads: 4) ---Al2O（110）-T

作者
Author: sobereva 时间: 2021-5-18 00:42

风飞发表于 2021-5-18 00:31
老师，您好，按照您方法，对于MgO(100)体系，可以收敛了，但是对于其他体系如Al2O3（110）-T （T是指甲苯 ...

OT的情况尝试利用outerSCF
对角化的情况尝试不同收敛算法以及参数（如alpha、NBROYDEN）

作者
Author: 风飞 时间: 2021-5-19 09:18

风飞发表于 2021-5-18 00:31
老师，您好，按照您方法，对于MgO(100)体系，可以收敛了，但是对于其他体系如Al2O3（110）-T （T是指甲苯 ...

最终利用outerSCF 并且增加了单点计算SCF和OUTER SC的循环数（我的体系设置为1000，然后在700-800多的时候），可以收敛

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