计算化学公社

标题: 请问为什么算出来的COHP和吸附能对不上呢？ [打印本页]

作者
Author: 赵寅生 时间: 2024-1-5 12:53
标题: 请问为什么算出来的COHP和吸附能对不上呢？
本帖最后由赵寅生于 2024-1-5 13:00 编辑

各位老师好，我目前有四个体系，FeN4，FeN3S， FeN4/Fe3C, FeN3S/Fe3C,我算了各中间体在Fe位点上的吸附，后续绘制了台阶图。按理说吸附能越负，吸附性越强, -ICOHP越正。但是算出来发现吸附能和ICOHP是对应不上的。在这里COHP我是用的（上旋+下旋）/2。原子对是Fe原子和最临近的O原子。
按理说，结果应该是自洽的，因为台阶图可以认为就是自由能绘制出来的，而COHP能定量反应它们的吸附强度，所以我想用COHP去对应台阶图。但是结果上来看基本是不匹配。
所以这种情况是合理的吗？还是说用COHP去理解台阶图就是不对的。
另外不同的体系之间，相同吸附物种下的COHP有比较的价值吗？
最后就是FeN4/Fe3C , FeN4/Fe3C的sigma=0.2,K点是111，（体系非常大，有100多个Fe原子，算不动）。这种精度配实验能接受吗？

作者
Author: 雾雨魔理沙 时间: 2024-3-25 08:16
我也遇到了两个体系对比，吸附能和COHP相反的结果，和楼主有同样的疑问，不知楼主解决这个问题了吗

作者
Author: xmuqyj 时间: 2024-7-24 09:29
同问，楼主解决这个问题了吗

作者
Author: masaki77 时间: 2024-10-7 15:22
遇到了同样的问题，请问楼主解决了吗？

作者
Author: PXZYYDS 时间: 2025-4-17 21:25
我也想问，楼主解决了吗？

作者
Author: sobereva 时间: 2025-4-18 11:55
COHP只是极其粗糙的理论，别太把这个当回事。如果发现说明不了问题就别盯着这个很粗糙的方法

有的是更好的考察吸附强度的指标，比如键级
使用Multiwfn对周期性体系做键级分析和NAdO分析考察成键特征
http://sobereva.com/719（http://bbs.keinsci.com/thread-47176-1-1.html）

再比如AIM拓扑分析考察临界点属性
使用Multiwfn结合CP2K做周期性体系的atom-in-molecules (AIM)拓扑分析
http://sobereva.com/717（http://bbs.keinsci.com/thread-46927-1-1.html）

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