计算化学公社

标题: 求助：CP2K计算钴镍MOF材料吸附SO2结合能为何为正值 [打印本页]

作者
Author: 18920543351 时间: 2022-3-21 10:34
标题: 求助：CP2K计算钴镍MOF材料吸附SO2结合能为何为正值
本帖最后由 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为例，输入文件已附带，望帮助！
(, 下载次数 Times of downloads: 23)
(, 下载次数 Times of downloads: 20)
(, 下载次数 Times of downloads: 12)

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

作者
Author: 18920543351 时间: 2022-3-21 11:17
本帖最后由 18920543351 于 2022-3-21 12:46 编辑

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

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

作者
Author: 18920543351 时间: 2022-3-24 09:44
用CP2K将gamma点计算改为K点计算，OT算法改为对角化算法，得到结果也不理想。

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