计算化学公社

标题: 我发现cp2k计算的吸附能等数值偏大 [打印本页]

作者
Author: reid 时间: 2023-9-3 21:51
标题: 我发现cp2k计算的吸附能等数值偏大
各位老铁，不知大家发现没有，cp2k计算吸附能等时，不管是ΔE还是ΔG，数值都偏大。vasp计算的一般都是1~2个eV，而cp2k有时会10个左右。不知有无解决办法

作者
Author: luchen 时间: 2023-9-4 11:56
这没法比啊，不能说vasp一般算1-2eV，cp2k有时10个左右。应该用完全相同的结构，完全相同的k点等相关设置进行比较

作者
Author: Eudaimonia 时间: 2023-9-4 16:40
建议再横向对比一下不同基组的能量区别，从DZVP-MOLOPT-GTH基组使用经验来看，是存在这种倾向的

作者
Author: reid 时间: 2023-9-4 22:48

Eudaimonia 发表于 2023-9-4 16:40
建议再横向对比一下不同基组的能量区别，从DZVP-MOLOPT-GTH基组使用经验来看，是存在这种倾向的

您觉得哪个基组可以避免这种情况？假设可以避免，是否可以用DZVP-MOLOPT-SR-GTH结构优化，然后再用好基组算单点？

作者
Author: xexlalalan 时间: 2023-9-5 11:42
差这么大应该不是基组的原因。这不是偏大，这是离谱的程度了。10个eV已经是氮气三键键能的强度了，强烈怀疑你输入文件哪里有问题

作者
Author: reid 时间: 2023-9-5 19:34

xexlalalan 发表于 2023-9-5 11:42
差这么大应该不是基组的原因。这不是偏大，这是离谱的程度了。10个eV已经是氮气三键键能的强度了，强烈怀疑 ...

输入文件应该没问题，都是sob老师的Multiwfn产生的，并且也用了大半年了

作者
Author: wolfli369 时间: 2023-9-6 08:57

reid 发表于 2023-9-4 22:48
您觉得哪个基组可以避免这种情况？假设可以避免，是否可以用DZVP-MOLOPT-SR-GTH结构优化，然后再用好基组 ...

这是常用方法，更好的基组会有更小的误差

作者
Author: xexlalalan 时间: 2023-9-7 18:57

reid 发表于 2023-9-5 19:34
输入文件应该没问题，都是sob老师的Multiwfn产生的，并且也用了大半年了

把输入文件和结构传上来看看？

作者
Author: sobereva 时间: 2023-9-8 01:08
DZVP-MOLOPT-SR-GTH普遍是用来做结构优化的，显然算单点时候应当用尽可能好的基组，诸如TZV2P-MOLOPT-GTH。就如同做量化计算如今谁也不会就拿6-31G**算能量问题一样
没有确切前提没法说诸如CP2K算的能量偏高什么的。需要有充分的计算细节，才能再说结果是否可比。

另外，算物理吸附问题，当基组不是很大时（比如只是DZVP-MOLOPT-SR-GTH这种程度），必须考虑counterpoise校正以避免BSSE问题高估吸附强度，这CP2K直接支持的。

作者
Author: 564059617 时间: 2024-5-12 22:03

sobereva 发表于 2023-9-8 01:08
DZVP-MOLOPT-SR-GTH普遍是用来做结构优化的，显然算单点时候应当用尽可能好的基组，诸如TZV2P-MOLOPT-GTH。 ...

请问计算物理吸附时，用TZV2P-MOLOPT-GTH基组，还需要BSSE校正吗

作者
Author: sobereva 时间: 2024-5-13 05:05

564059617 发表于 2024-5-12 22:03
请问计算物理吸附时，用TZV2P-MOLOPT-GTH基组，还需要BSSE校正吗

如果有余力建议用counterpoise方式考虑BSSE校正

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