计算化学公社

标题: MgO（100）表面弛豫时SCF不收敛怎么解决 [打印本页]

作者
Author: 科研好难 时间: 2023-11-27 10:45
标题: MgO（100）表面弛豫时SCF不收敛怎么解决
各位老师好，我用PBEsol泛函和pob-TZVP-rev2基组做MgO（100）表面弛豫时，k点设置为3，3，1并且固定了下面的原子（输入文件为MgOexp.inp），计算初scf就出现了较大的震荡，并且在50步左右维持一定的震荡规律，于是我在.inp文件里把k点改为1，2，1，MAX SCF改成了50，MIXING的ALPHA也降到了0.3，顺便把CUTOFF值增到600（输入文件为MgOexp1.inp），这样scf每步迭代时间降低很多，不过scf还是没有收敛趋势（输出文件为MgOexp1.out），请问一下各位对于我现在的情况，可以优先考虑做其他的什么尝试？结构和服务器配置如下：
(, 下载次数 Times of downloads: 4) (, 下载次数 Times of downloads: 4)

作者
Author: raolixiang 时间: 2023-11-27 11:02
slab的厚度好像要进行收敛性测试，并且可能需要固定底层的一定厚度的原子

作者
Author: sobereva 时间: 2023-11-28 03:40
这么大模型，再加上MgO是绝缘体，根本没必要设k点
SCF不收敛问题下文该说的都说了
解决CP2K中SCF不收敛的方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html）

对当前问题根本没必要用这么大的模型，用一个单胞（z方向加真空区）考虑k点做几何优化，比当前模型结构收敛容易得多而且耗时低得多得多，优化完了再扩成你想要大小的slab就完了

作者
Author: 科研好难 时间: 2023-11-28 09:59

sobereva 发表于 2023-11-28 03:40
这么大模型，再加上MgO是绝缘体，根本没必要设k点
SCF不收敛问题下文该说的都说了
解决CP2K中SCF不收敛的 ...

好的，谢谢sob老师，选择这么大的模型是因为后期要吸附稍大一点的分子结构，按着老师您的做法，对单胞（有真空区）进行优化后扩胞就可以直接吸附分子结构了是吧？

作者
Author: sobereva 时间: 2023-11-29 00:04

科研好难发表于 2023-11-28 09:59
好的，谢谢sob老师，选择这么大的模型是因为后期要吸附稍大一点的分子结构，按着老师您的做法，对单胞（ ...

是

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