计算化学公社

标题: 求助用CP2K计算吸附能（结合能）的问题 [打印本页]

作者
Author: liuxp 时间: 2025-6-7 22:07
标题: 求助用CP2K计算吸附能（结合能）的问题
现在对计算结合能的方法存疑？Eads=E(AB) -EA-EB,请问一下EA和EB是各自结构优化后的能量吗？目前计算能量用TZV2P-MOLOPT-GTH级别,但是计算过EAB 和（EA+EB)能量差2.5 eV ,是不是算出来的吸附能偏大，文献的吸附步骤大部分在0.5 eV 左右，目前保持的同一级别，都开了smearing和D3？请问一下这个2.3eV 是不是不合理，不合理的话该如何改进

作者
Author: sobereva 时间: 2025-6-8 01:37
是

不说体系就光说吸附能大小根本毫无意义！认真看下文
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚、完整
http://sobereva.com/620（http://bbs.keinsci.com/thread-25787-1-1.html）
而且连泛函也不说，是什么类型吸附也不说，是否需要D3根本无从判断。而且也根本无从判断是否有开smearing的必要。

当前情况下，“这个2.3eV 是不是不合理” 你只能问算命的。

作者
Author: liuxp 时间: 2025-6-9 00:05
本帖最后由 liuxp 于 2025-6-9 00:08 编辑

sobereva 发表于 2025-6-8 01:37
是

不说体系就光说吸附能大小根本毫无意义！认真看下文

不好意思老师泛函用的PBE,元素为LaNiO3,吸附的是两个糠醛分子。请问老师这个体系需要考虑D3吗？因为看谈谈“计算时是否需要加DFT-D3色散校正？” - 思想家公社的门口：量子化学·分子模拟·二次元这个博文所以吸附一直加的D3

作者
Author: sobereva 时间: 2025-6-9 17:33

liuxp 发表于 2025-6-9 00:05
不好意思老师泛函用的PBE,元素为LaNiO3,吸附的是两个糠醛分子。请问老师这个体系需要考虑D3吗？因为看 ...

需要

作者
Author: liuxp 时间: 2025-6-9 22:50

sobereva 发表于 2025-6-9 17:33
需要

那我现在已经用的PBE＋D3＋ TZV2P-MOLOPT-SR-GTH计算的能量？算出来的吸附能在1.8 eV? 请问一下是不是有一些偏大？要考虑BSSE吗

作者
Author: sobereva 时间: 2025-6-9 23:43

liuxp 发表于 2025-6-9 22:50
那我现在已经用的PBE＋D3＋ TZV2P-MOLOPT-SR-GTH计算的能量？算出来的吸附能在1.8 eV? 请问一下是不是有 ...

如果是化学吸附，没有必然性偏大
物理吸附用counterpoise减小BSSE有益，化学吸附不需要

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