计算化学公社

标题: 过渡金属吸附H2，吸附能算出来为正，感觉很奇怪 [打印本页]

作者
Author: ZeeWong 时间: 2023-9-25 10:34
标题: 过渡金属吸附H2，吸附能算出来为正，感觉很奇怪
请教各位老师，我计算了一个过渡金属掺杂二维材料对各种气体分子的吸附能，结果非常怪异，有的高达-7 eV，H2和CO的却为正。
第一次用CP2K算周期性体系的吸附能，想请教一下问题会出在哪。
所有inp都通过Multiwfn产生。

各气体分子吸附能如下，单位eV，H2的计算输入文件见附件。
C2H6       -3.521070261
C2H4       -1.845576199
C2H2       -7.013739229
CH4             -0.640372142
CO             3.353076136
H2                2.066662045

作者
Author: sobereva 时间: 2023-9-25 11:00
CUTOFF必须保持统一
另外，算单点用的盒子应和优化整体用的盒子保持一致

作者
Author: ZeeWong 时间: 2023-9-25 11:58
非常感谢您，卢老师，我统一了截断能和盒子尺寸后大部分气体的吸附能都正常些了，可CH4的吸附能还是很奇怪，从吸附结构上来看不会这么大。能请您再看看iCH4的inp吗。

C2H6 0.448590451
C2H4 -3.568319466
C2H2 -9.132006286
CH4 -18.08744818
CO -4.606243959
H2 -5.516851042

作者
Author: sobereva 时间: 2023-9-26 07:47

ZeeWong 发表于 2023-9-25 11:58
非常感谢您，卢老师，我统一了截断能和盒子尺寸后大部分气体的吸附能都正常些了，可CH4的吸附能还是很奇怪 ...

EPS_DEFAULT、EPS_SCF这些理应都统一
你没上传优化完的结构（诸如restart文件），我不好判断
如果希望我回复，每次点击我的帖子下方的回复按钮，否则我大概率看不到

作者
Author: ZeeWong 时间: 2023-9-26 10:03

sobereva 发表于 2023-9-26 07:47
EPS_DEFAULT、EPS_SCF这些理应都统一
你没上传优化完的结构（诸如restart文件），我不好判断
如果希望 ...

感谢卢老师的回复
我统一了EPS的两个设置后，吸附能为-3.22 eV，但从吸附结构来看更像是单纯的物理吸附，总不觉能有这么大。
附件里是总的吸附结构优化后的restart文件，输入文件只统一了EPS，其他没变。

作者
Author: sobereva 时间: 2023-9-26 13:39

ZeeWong 发表于 2023-9-26 10:03
感谢卢老师的回复
我统一了EPS的两个设置后，吸附能为-3.22 eV，但从吸附结构来看更像是单纯的物理吸附 ...

对于算物理吸附能的目的，不考虑counterpoise的话，当前基组太小，BSSE效应会非常强，会严重高估吸附能。尝试改成尽可能大的比如TZV2P档次的基组，BSSE小得多

作者
Author: ZeeWong 时间: 2023-9-28 11:03

sobereva 发表于 2023-9-26 13:39
对于算物理吸附能的目的，不考虑counterpoise的话，当前基组太小，BSSE效应会非常强，会严重高估吸附能。 ...

感谢您的指导，现在吸附能看起来正常多了

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