计算化学公社

标题: 求助 cp2k SnO2（1 1 0）超晶胞SCF不收敛 [打印本页]

作者
Author: qczcb2 时间: 2023-9-21 09:36
标题: 求助 cp2k SnO2（1 1 0）超晶胞SCF不收敛
本帖最后由 qczcb2 于 2023-9-21 09:36 编辑

目的是导出molden所以用了较大的超晶胞。输入文件通过Multiwfn创建，计算级别为默认的PBE/DZVP-MOLOPT-SR-GTH，开启DFT-D3(BJ),SnO2块体CIF由MS导出，经cp2k优化后切（1 1 0）三层晶面（6层原子），真空层为20A，居中，最下方两层原子固定。在cp2k中优化（1 1 0）晶面单胞后，扩胞至4 2 1超胞，开启OT，开启outer SCF不收敛。将cut off加到1200以后，SCF有收敛趋势，但是收敛至10-4以后，计算的电荷密度误差增大，精度变差，最后提示 Cholesky decomposition failed. Matrix ill conditioned 。换成CG以后，最后SCF收敛也是越来越差，，电荷密度越来越离谱。
我的问题是:1,我的模型是否存在不合理的地方？如何修改？2.目前的截断能还有必要再增加吗？3.有没有必要更换基组？
附件SnO2.inp与SnO2-1.restart为块体输入输出，SnO2-3.inp与SnO2-3-1.restart为（1 1 0）晶面单胞输入输出，最后SnO2-3-3.inp与SnO2-3-3.out为4 2 1超胞OT CG算到一半停下来的结果

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