求助：CP2K计算钴镍MOF材料吸附SO2结合能为何为正值

18920543351 · 发表于 Post on 2022-3-21 10:34:31

本帖最后由 18920543351 于 2022-3-21 12:43 编辑

各位老师好！
在计算MOF-74-Co和Ni吸附SO2结合能时，利用DZVP-SR-MOLOPT和TZVP-MOLOPT基组，计算结果Eab-Ea-Eb为正值。
下面以MOF-74-Co，基组为DZVP-SR-MOLOPT为例，输入文件已附带，望帮助！

74co-so2.inp (11.36 KB, 下载次数 Times of downloads: 23)

74co.inp (12.61 KB, 下载次数 Times of downloads: 20)

SO2.inp (2.35 KB, 下载次数 Times of downloads: 12)

wolfli369 · 发表于 Post on 2022-3-21 11:01:30

设置一下磁性，打开UKS，计算再试一下，找到稳定结构的最低能量

18920543351 · 发表于 Post on 2022-3-21 11:17:42

本帖最后由 18920543351 于 2022-3-21 12:46 编辑

wolfli369 发表于 2022-3-21 11:01
设置一下磁性，打开UKS，计算再试一下，找到稳定结构的最低能量

已经对比不同自旋多重度下能量，发现自旋极化单重态是基态，然后算完结合能也是正的。同样计算了普通单重态结合能量也是正的。

18920543351 · 发表于 Post on 2022-3-24 09:44:43

用CP2K将gamma点计算改为K点计算，OT算法改为对角化算法，得到结果也不理想。

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

[CP2K] 求助：CP2K计算钴镍MOF材料吸附SO2结合能为何为正值