计算化学公社

 找回密码 Forget password
 注册 Register
Views: 1003|回复 Reply: 0
打印 Print 上一主题 Last thread 下一主题 Next thread

[CP2K] 求助 cp2k SnO2(1 1 0)超晶胞SCF不收敛

[复制链接 Copy URL]

39

帖子

0

威望

1016

eV
积分
1055

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 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算到一半停下来的结果

SnO2.inp

4.59 KB, 下载次数 Times of downloads: 2

块体输入

SnO2-1.restart

4.8 KB, 下载次数 Times of downloads: 0

块体输出

SnO2-3.inp

4.78 KB, 下载次数 Times of downloads: 0

晶面优化输入

SnO2-3-1.restart

5.55 KB, 下载次数 Times of downloads: 0

晶面优化输出

SnO2-3-3.inp

3.79 KB, 下载次数 Times of downloads: 2

超胞SCF输入

SnO2-3-3.out

64.28 KB, 下载次数 Times of downloads: 0

OT CG输出

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-16 20:30 , Processed in 0.188917 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list