请教：CP2K做过渡态、电子结构等静态计算时的可靠性

卡开发发 · 发表于 Post on 2021-2-21 04:21:10

上面的deltatest我补充解释一下：使用相同的赝势情况下，abinit基于平面波可以通过调整波函数的动能截断来逼近基组极限，也就是说abinit可以认为是GTH赝势下的基组极限结果，随着zeta增加其实也并未趋于基组极限。S22测试中只包含了一部分主族元素，但过渡金属可能会出现问题，比如（DOI: 10.1021/acs.jctc.8b00957的SI）

CP2K做过渡金属体系的文章其实也并不是很多，具体原因不明，可能还有其他奇怪的例子吧。其实用什么程序无所谓，但建议各个程序使用前先对精度大致做好调研，然后对自己体系最好做好充分的测试。

邱美佳 · 发表于 Post on 2021-3-6 14:48:58

sobereva 发表于 2021-2-19 07:22
可以

另外，对于大多数波函数分析的目的，比如RDG、IRI、ELF、AIM分析等，其实TZVP-MOLOPT-GTH就够了 ...

请问sob老师cp2k如果是算DOS那要选什么级别的基组结果比较可靠呢？

sobereva · 发表于 Post on 2021-3-7 05:34:05

邱美佳发表于 2021-3-6 14:48
请问sob老师cp2k如果是算DOS那要选什么级别的基组结果比较可靠呢？

通常DZVP-MOLOPT-SR-GTH就可以，想更好点就用TZVP的

卡开发发 · 发表于 Post on 2021-3-9 09:59:05

邱美佳发表于 2021-3-6 14:48
请问sob老师cp2k如果是算DOS那要选什么级别的基组结果比较可靠呢？

建议上TZ以上的（或者小心测试），DZ可能会高估gap，参考DOI: 10.1063/1.4736998。

sobereva · 发表于 Post on 2021-3-16 08:16:10

CP2K的MOLOPT系列比同等zeta档次的其它基组精度往往高一个档，如根据原文测试，DZVP-MOLOPT-GTH算能量问题上就已经略胜于TZVP-MOLOPT-GTH，因此有较高概率DZVP-MOLOPT-SR-GTH就已经够用。不放心的话可以和pob-TZVP或TZV2P-GTH对比一下。

丁越 · 发表于 Post on 2021-10-19 14:14:43

卡开发发发表于 2021-2-18 21:50
1、CP2K的GPW是结合赝势的，所以赝波函数只有价层的且是smooth的，有可能不能正确识别NCP，也许得考虑缀加 ...

请问一下卡开发发老师，我用CP2K优化结构，用VASP计算能量，这样的方式审稿人能接受么？主要是CP2K算的很快，想利用这个优势。
其次，我算了一个含Cu、Fe的体系吸附NO、CO等，但是计算下来的结合能和VASP利用CP2K优化完的结构算结合能相差太离谱了。。肯定是VASP的结果更靠谱一点

Aridea · 发表于 Post on 2021-10-19 18:45:31

丁越发表于 2021-10-19 14:14
请问一下卡开发发老师，我用CP2K优化结构，用VASP计算能量，这样的方式审稿人能接受么？主要是CP2K算的很 ...

我最近也算了一套cp2k的弱相互作用，结合能就比较离谱，1-6eV的都有

卡开发发 · 发表于 Post on 2021-10-23 18:29:20

丁越发表于 2021-10-19 14:14
请问一下卡开发发老师，我用CP2K优化结构，用VASP计算能量，这样的方式审稿人能接受么？主要是CP2K算的很 ...

建议你可以找个稍微小一些的体系测试一下吸附能之类的数据。目前我也说不好是赝势、基组还是算法的问题，影响的因素很多，例如前面看到赝势大小芯也会显著影响精度。

卡开发发 · 发表于 Post on 2021-10-23 18:42:37

Aridea 发表于 2021-10-19 18:45
我最近也算了一套cp2k的弱相互作用，结合能就比较离谱，1-6eV的都有

因素比较多，我暂时不好说。或者可以把数据输入拿出来讨论。

Aridea · 发表于 Post on 2021-10-23 21:00:47

卡开发发发表于 2021-10-23 18:42
因素比较多，我暂时不好说。或者可以把数据输入拿出来讨论。

大大，我感觉有点像反应了，必经是带电的体系，电子转移啥的也是有可能的嘛

卡开发发 · 发表于 Post on 2021-10-23 21:23:10

Aridea 发表于 2021-10-23 21:00
大大，我感觉有点像反应了，必经是带电的体系，电子转移啥的也是有可能的嘛

带电怎么个带电法？你直接体系加1电荷？周期性边界的计算这个要非常小心。

Aridea · 发表于 Post on 2021-10-23 21:59:55

卡开发发发表于 2021-10-23 21:23
带电怎么个带电法？你直接体系加1电荷？周期性边界的计算这个要非常小心。

嗯嗯，大大我是用cp2k计算的2Dxy slab加的MT求解器，应该可以吧

卡开发发 · 发表于 Post on 2021-10-24 13:21:27

Aridea 发表于 2021-10-23 21:59
嗯嗯，大大我是用cp2k计算的2Dxy slab加的MT求解器，应该可以吧

可以。

Aridea · 发表于 Post on 2021-10-24 16:49:29

卡开发发发表于 2021-10-24 13:21
可以。

嘿嘿，当初我发帖子大大还有社长指点过了，感谢感谢~

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

[CP2K] 请教：CP2K做过渡态、电子结构等静态计算时的可靠性