求助：MgO(100)单点计算不收敛

风飞 · 发表于 Post on 2021-5-17 21:09:44

本帖最后由风飞于 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

sobereva · 发表于 Post on 2021-5-17 21:54:03

实测开OT可以收敛

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

风飞 · 发表于 Post on 2021-5-18 00:31:14

本帖最后由风飞于 2021-5-18 00:32 编辑

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

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

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

AlT1.inp (9.51 KB, 下载次数 Times of downloads: 3) ---Al2O（110）-T

sobereva · 发表于 Post on 2021-5-18 00:42:18

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

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

风飞 · 发表于 Post on 2021-5-19 09:18:08

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

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

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

[CP2K] 求助：MgO(100)单点计算不收敛