CP2K计算“吸附能”过大的问题咨询

HanC · 发表于 Post on 2024-7-25 21:01:50

尊敬的各位老师好，我遇到了利用CP2K计算“对氯苯酚在金红石（110）晶面吸附”，但得到的吸附能过大的问题。我仔细查阅了相关资料，也尽可能保证了三个文件中理论方法、基组、k点以及截断能均一致，但是计算得到的吸附能过大，明显错误，如下图。请各位老师可以帮我检查下输入文件是否符合计算吸附能所对应的条件? 在输入文件中哪里出现了错误，需要如何改正？麻烦各位老师了。

HanC · 发表于 Post on 2024-8-3 11:17:28

reid 发表于 2024-8-2 16:58
都用PBE+D3+DZVP-MOLOPT-SR-GTH做结构优化，然后单点能都用PBE+D3+TZV2P-MOLOPT-GTH

好的感谢老师

reid · 发表于 Post on 2024-8-2 16:58:19

HanC 发表于 2024-7-31 20:33
感谢感谢老师回答，请问关于结构优化我有一些这样的问题：单纯的晶面和吸附有分子的晶面，我为了计算更快 ...

都用PBE+D3+DZVP-MOLOPT-SR-GTH做结构优化，然后单点能都用PBE+D3+TZV2P-MOLOPT-GTH

HanC · 发表于 Post on 2024-7-31 20:33:21

reid 发表于 2024-7-27 22:11
结构优化可以用DZVP-MOLOPT-SR-GTH，算吸附能时必须再算一遍单点能，单点能用社长说的TZV2P-MOLOPT-GTH

感谢感谢老师回答，请问关于结构优化我有一些这样的问题：单纯的晶面和吸附有分子的晶面，我为了计算更快一些，是否可以都用DZVP-MOLOPT-SR-GTH来计算，这样是否就能得到相对可靠的结构，后面我再选择更大的基组并加上D3矫正来计算吸附能，这样是否能起到降低计算时间的目的呢；还是说仅能单纯的晶面可以这样，吸附分子后的晶面就不能这样做了呢？

reid · 发表于 Post on 2024-7-27 22:11:23

结构优化可以用DZVP-MOLOPT-SR-GTH，算吸附能时必须再算一遍单点能，单点能用社长说的TZV2P-MOLOPT-GTH

HanC · 发表于 Post on 2024-7-26 18:35:50

sobereva 发表于 2024-7-26 04:00
对于色散作用不可忽略的情况下，PBE不加D3算物理吸附是开玩笑，PBE描述色散作用那么烂。
该说的都说了
...

感谢sob老师回复

sobereva · 发表于 Post on 2024-7-26 04:00:02

HanC 发表于 2024-7-25 22:45
sob老师您好，在重新编写文件的过程中，考虑到您提到的问题，我是利用Multiwfn辅助编写的CP2K的输入文件 ...

对于色散作用不可忽略的情况下，PBE不加D3算物理吸附是开玩笑，PBE描述色散作用那么烂。
该说的都说了
谈谈“计算时是否需要加DFT-D3色散校正？”
http://sobereva.com/413

HanC · 发表于 Post on 2024-7-25 23:09:06

Eudaimonia 发表于 2024-7-25 23:02
大这么多，大概率是哪里设置不同导致的，仔细检查所有输入文件会和输出文件再看一下

比较容易出问题的： ...

感谢老师解惑，我再仔细检查下您提到的问题

Eudaimonia · 发表于 Post on 2024-7-25 23:02:48

大这么多，大概率是哪里设置不同导致的，仔细检查所有输入文件会和输出文件再看一下

比较容易出问题的：自旋和UKS是否合理设置，计算二维结构晶胞大小是否一致

HanC · 发表于 Post on 2024-7-25 22:45:21

本帖最后由 HanC 于 2024-7-25 22:46 编辑

sobereva 发表于 2024-7-25 21:14
用DZVP-MOLOPT-SR-GTH那么小的基组还不考虑counterpoise校正，吸附能不高估才怪
当前是纯泛函，写AUTO_BAS ...

sob老师您好，在重新编写文件的过程中，考虑到您提到的问题，我是利用Multiwfn辅助编写的CP2K的输入文件，可能在设置过程中没有仔细注意让输入文件加入了D3校正。
在实际计算吸附能的过程中，您提到考虑counterpoise校正能够改善计算结果，那么D3是否需要考虑呢，是否需要整体和片段中均加入D3让其处于一个计算级别下呢。
您之前的讲授和帖子能够学习到，D3能够较好的来描述分子间弱相互作用，在计算吸附能的过程中，加入D3是否能更好地描述吸附能，或者说更好地描述小分子与晶体表面间相互作用以及小分子与晶体表面间的电子传递情况呢不知道这是不是能联系在一起考虑，希望老师们指导一二。

HanC · 发表于 Post on 2024-7-25 21:50:27

sobereva 发表于 2024-7-25 21:45
不考虑counterpoise校正时起码用TZV2P-MOLOPT-GTH。用TZV2P-MOLOPT-GTH再考虑counterpoise校正还有改进余地 ...

感谢sob老师

sobereva · 发表于 Post on 2024-7-25 21:45:46

不考虑counterpoise校正时起码用TZV2P-MOLOPT-GTH。用TZV2P-MOLOPT-GTH再考虑counterpoise校正还有改进余地
看MOLOPT原文里BSSE的对比

HanC · 发表于 Post on 2024-7-25 21:25:31

本帖最后由 HanC 于 2024-7-25 21:26 编辑

sobereva 发表于 2024-7-25 21:14
用DZVP-MOLOPT-SR-GTH那么小的基组还不考虑counterpoise校正，吸附能不高估才怪
当前是纯泛函，写AUTO_BAS ...

万分感谢sob老师的回复，我修改后重新计算。想追问下，计算该吸附能的情况下，最低要求什么样的基组能够满足吸附能的计算；或者说什么样的基组可以不再另加counterpoise校正，小基组加校正和较大基组哪个能更好的描述吸附能呢。

sobereva · 发表于 Post on 2024-7-25 21:14:15

用DZVP-MOLOPT-SR-GTH那么小的基组还不考虑counterpoise校正，吸附能不高估才怪
当前是纯泛函，写AUTO_BASIS RI_HFX MEDIUM作甚

哪能计算整体加D3而计算片段不加D3，根本不是同一个计算级别

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

[CP2K] CP2K计算“吸附能”过大的问题咨询