我发现cp2k计算的吸附能等数值偏大

reid · 发表于 Post on 2023-9-3 21:51:03

各位老铁，不知大家发现没有，cp2k计算吸附能等时，不管是ΔE还是ΔG，数值都偏大。vasp计算的一般都是1~2个eV，而cp2k有时会10个左右。不知有无解决办法

luchen · 发表于 Post on 2023-9-4 11:56:07

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

Eudaimonia · 发表于 Post on 2023-9-4 16:40:47

建议再横向对比一下不同基组的能量区别，从DZVP-MOLOPT-GTH基组使用经验来看，是存在这种倾向的

reid · 发表于 Post on 2023-9-4 22:48:05

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

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

xexlalalan · 发表于 Post on 2023-9-5 11:42:53

差这么大应该不是基组的原因。这不是偏大，这是离谱的程度了。10个eV已经是氮气三键键能的强度了，强烈怀疑你输入文件哪里有问题

reid · 发表于 Post on 2023-9-5 19:34:40

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

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

wolfli369 · 发表于 Post on 2023-9-6 08:57:38

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

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

xexlalalan · 发表于 Post on 2023-9-7 18:57:26

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

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

sobereva · 发表于 Post on 2023-9-8 01:08:14

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

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

564059617 · 发表于 Post on 2024-5-12 22:03:12

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

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

sobereva · 发表于 Post on 2024-5-13 05:05:41

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] 我发现cp2k计算的吸附能等数值偏大